computes coefficients of interpolation such that quantity at position x,y can be computed as
u(x,y) = sum_{i=1,intorder+1,j=1,intorder+1} coeffs(i,j) u(i,j)
i,j is Cartesian ginterpolation grid with x_i=0...intorder, y_j=0...intorder x,y is coordinate of to be interpolated point with respect to southwest point of interpolation stamp maesured in grid distance
INPUT: intorder: order of interpolation x: x-coordinate w.r.t southwest node point y: y-coordinate w.r.t southwest node point
OUTPUT: coeffs: coefficients of interpolating polynomial
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | intorder | |||
| real(kind=FP), | intent(in) | :: | x | |||
| real(kind=FP), | intent(in) | :: | y | |||
| real(kind=FP), | intent(out), | dimension(intorder+1,intorder+1) | :: | coeffs |