interpolate_replace Module Subroutine

public impure module subroutine interpolate_replace(xold, xnew, fx, idxx)

Uses

Interpolate f(xold) -> f(xnew). Only consider xnew values that are contained within xold for now. Xnew returns untouched, but fx is overwritten with interpolated values on the grid xnew(idxx)

Arguments

Type IntentOptional Attributes Name
real(kind=dp), intent(in) :: xold(:)

The grid of values on which our function has been evaluated

real(kind=dp), intent(in) :: xnew(:)

The new grid of values that we want our function to be evaluated on.

real(kind=dp), intent(inout), allocatable :: fx(:)

On input: evaluated function f(xold) On output: the evaluated function f(xnew)

integer, intent(out), optional, allocatable :: idxx(:)

The indices of values of xnew that are used to evaluate fx