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