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