Skip to main content

InterpolationPoint

This routine returns the interpolation points on triangle.

Interface​

INTERFACE
  MODULE FUNCTION InterpolationPoint_Triangle(order, ipType, &
    & layout, xij) RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: order
    !! order
    INTEGER(I4B), INTENT(IN) :: ipType
    !! interpolation point type
    REAL(DFP), OPTIONAL, INTENT(IN) :: xij(:, :)
    !! Coord of domain in xij format
    CHARACTER(*), INTENT(IN) :: layout
    !! local node numbering layout, always VEFC
    REAL(DFP), ALLOCATABLE :: ans(:, :)
    !! xij coordinates
  END FUNCTION InterpolationPoint_Triangle
END INTERFACE
xij
  • xij contains nodal coordinates of triangle in xij format.
  • SIZE(xij,1) = nsd, and SIZE(xij,2)=3
  • If xij is absent then unit triangle is assumed
ipType

  • ipType is interpolation point type, it can take following values

  • Equidistance, uniformly/evenly distributed points

  • BlythPozChebyshev

  • BlythPozLegendre

  • GaussLegendreLobatto, which is same as IsaacLegendre

  • GaussChebyshevLobatto, which is same as IsaacChebyshev

  • GaussJacobiLobatto

  • GaussUltrasphericalLobatto

  • IsaacChebyshev

  • IsaacLegendre

  • ChenBabuska TODO

  • Hesthaven TODO

  • Feket TODO

layout
  • layout specifies the arrangement of points. The nodes are always returned in “VEFC” format (vertex, edge, face, cell). 1:3 are vertex points, then edge, and then internal nodes. The internal nodes also follow the same convention. Please read Gmsh manual on this topic.