def dweighted_jacobi_function(x,ns,s=1.,t=1.,scale=1.,shift=0.): from spyctral.mapjpoly.maps import sqrt_weight, dsqrt_weight w = sqrt_weight(x,s=s,t=t,scale=scale,shift=shift) dw = dsqrt_weight(x,s=s,t=t,scale=scale,shift=shift) PB = jacobi_function(x,ns,s=s,t=t,scale=scale,shift=shift) dPB = djacobi_function(x,ns,s=s,t=t,scale=scale,shift=shift) # Yay product rule return (PB.T*dw).T + (dPB.T*w).T
def weighted_jacobi_function(x,ns,s=1.,t=1.,scale=1.,shift=0.): from spyctral.mapjpoly.maps import sqrt_weight w = sqrt_weight(x,s=s,t=t,scale=scale,shift=shift) return (jacobi_function(x,ns,s=s,t=t,scale=scale,shift=shift).T*w).T