polyQuadratic_3D Function

private pure function polyQuadratic_3D(cxDir) result(phi)

This function returns matrix entries for quadratic polynomial for 3D stencil

Arguments

Type IntentOptional Attributes Name
real(kind=rk), intent(in) :: cxDir(3)

Return Value real(kind=rk), (10)


Called by

proc~~polyquadratic_3d~~CalledByGraph proc~polyquadratic_3d polyQuadratic_3D proc~build_matrixlsf_quadintp build_matrixLSF_quadIntp proc~build_matrixlsf_quadintp->proc~polyquadratic_3d proc~append_intpmatrixlsf append_intpMatrixLSF proc~append_intpmatrixlsf->proc~build_matrixlsf_quadintp interface~append~9 append interface~append~9->proc~append_intpmatrixlsf

Source Code

  pure function polyQuadratic_3D(cxDir) result (phi)
    ! --------------------------------------------------------------------------
    real(kind=rk), intent(in) :: cxDir(3)
    real(kind=rk) :: phi(10)
    ! --------------------------------------------------------------------------
    phi( 1) = 1.0_rk
    phi( 2) = cxDir(c_x)
    phi( 3) = cxDir(c_y)
    phi( 4) = cxDir(c_z)
    phi( 5) = cxDir(c_x)**2
    phi( 6) = cxDir(c_y)**2
    phi( 7) = cxDir(c_z)**2
    phi( 8) = cxDir(c_x)*cxDir(c_y)
    phi( 9) = cxDir(c_y)*cxDir(c_z)
    phi(10) = cxDir(c_z)*cxDir(c_x)

  end function polyQuadratic_3D