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