lag_intrp_grd2 Subroutine

public subroutine lag_intrp_grd2(nx, ny, npx, npy, matx, maty, f, fintrp, rwrk)

Interpolates nodal function values to a grid using barycentric lagrange interpolation in two variables.

Arguments

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

Number of nodes along x = nx+1 and along y = ny+1

integer, intent(in) :: ny

Number of nodes along x = nx+1 and along y = ny+1

integer, intent(in) :: npx

Number of interpolating points along x = npx+1 and along y = npy+1

integer, intent(in) :: npy

Number of interpolating points along x = npx+1 and along y = npy+1

real(kind=rp), intent(in), dimension(:,:) :: matx

(npx+1,nx+1). Interpolation matrix along x

real(kind=rp), intent(in), dimension(:,:) :: maty

(npy+1,npy+1). Interpolation matrix along y

real(kind=rp), intent(in), dimension(:) :: f

((nx+1)*(ny+1),). Value of the function at the nodes

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

((npx+1)*(npy+1),). Value of the function at the nodes

real(kind=rp), intent(inout), dimension(:) :: rwrk

Workspace array of minimum size (npx+1)*(ny+1).