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 pgq(N,s=1.,t=1.,scale=1.,shift=0.):
from spyctral.mapjpoly.maps import weight
[x,w] = gq(N,s=s,t=t,scale=scale,shift=shift)
# Weight accordingly
w /= weight(x,s=s,t=t,scale=scale,shift=shift)
return [x,w]
def gq(N,s=1.,t=1.,scale=1.,shift=0.):
from spyctral.common.maps import physical_scaleshift as pss
from spyctral.mapjpoly.maps import st_to_ab, r_to_x
from spyctral.jacobi.quad import gq
[alpha,beta]=st_to_ab(s,t)
[r,w] = gq(N,alpha=alpha,beta=beta)
x = r_to_x(r)
#x *= scale
#x += shift
pss(x,scale=scale,shift=shift)
return [x,w]
def dgen1(**kwargs):
from spyctral.jacobi.eval import jacobi_polynomial
[x,w] = gq(**kwargs)
ps = jacobi_polynomial(x,range(2*N),alpha=alpha,beta=beta,scale=scale,shift=shift)
ps[:,0] = 0
return [w,ps]
