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