def dgen1(**kwargs): from spyctral.jacobi.quad import gq from spyctral.jacobi.eval import jacobi_polynomial [r,w] = gq(**kwargs) ps = jacobi_polynomial(r,range(N),alpha=alpha,beta=beta,scale=scale,shift=shift) return [w,ps]
def dgen1(**kwargs): from spyctral.jacobi.eval import jacobi_polynomial from spyctral.jacobi.quad import glq [x,w] = glq(**kwargs) ps = jacobi_polynomial(x,range(2*N-2),alpha=alpha,beta=beta,scale=scale,shift=shift) ps[:,0] = 0 return [w,ps]