bary_weights Subroutine

public subroutine bary_weights(n, w)

Calculates the barycentric weights for a given set of nodes. As of now, assumues CGL nodes. TODO: extend stably for the general case (CGL, CGR, CG)

Arguments

Type IntentOptional Attributes Name
integer, intent(in) :: n

Number of interpolation nodes = n+1.

real(kind=rp), intent(out), dimension(:) :: w

(n+1,). Barycentric weights