Beispiel #1
0
    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]
Beispiel #2
0
 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]