Esempio n. 1
0
    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]
Esempio n. 2
0
 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]