This function returns matrix entries for quadratic polynomial for 1D stencil
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | cxDir(3) |
pure function polyQuadratic_1D(cxDir) result (phi)
! --------------------------------------------------------------------------
real(kind=rk), intent(in) :: cxDir(3)
real(kind=rk) :: phi(3)
! --------------------------------------------------------------------------
phi(1) = 1.0_rk
phi(2) = cxDir(c_x)
phi(3) = cxDir(c_x)**2
end function polyQuadratic_1D