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