def dgen2(**kwargs): from spyctral.laguerre.quad import pgq from spyctral.laguerre.eval import laguerre_function [r,w] = pgq(**kwargs) ps = laguerre_function(r,range(N),alpha=alpha,scale=scale,shift=shift) return [w,ps]
def dgen1(**kwargs): from spyctral.laguerre.eval import laguerre_function from spyctral.laguerre.weights import sqrt_weight [x,w] = pgq(**kwargs) ps = laguerre_function(x,range(2*N),alpha=alpha,scale=scale,shift=shift) ps[:,0] = 0 ps = (ps.T*sqrt_weight(x,alpha=alpha,scale=scale,shift=shift)).T return [w,ps]