interpol_coeffs Subroutine

public subroutine interpol_coeffs(intorder, x, y, coeffs)

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

Arguments

Type IntentOptional 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

Calls

proc~~interpol_coeffs~~CallsGraph proc~interpol_coeffs interpol_coeffs dgetrf dgetrf proc~interpol_coeffs->dgetrf dgetri dgetri proc~interpol_coeffs->dgetri proc~get_stderr get_stderr proc~interpol_coeffs->proc~get_stderr sgetrf sgetrf proc~interpol_coeffs->sgetrf sgetri sgetri proc~interpol_coeffs->sgetri

Called by

proc~~interpol_coeffs~~CalledByGraph proc~interpol_coeffs interpol_coeffs proc~create_cart_map_matrix create_cart_map_matrix proc~create_cart_map_matrix->proc~interpol_coeffs proc~create_polar_map_matrix create_polar_map_matrix proc~create_polar_map_matrix->proc~interpol_coeffs