Skip to main content

RecursiveNode2D

Returns the barycentric coordinates of recursive nodes on the triangle.

Interface​

INTERFACE
  MODULE FUNCTION RecursiveNode2D(order, ipType, domain) &
    & RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: order
      !! order >= 0
    INTEGER(I4B), INTENT(IN) :: ipType
      !! interpolation point type
      !! Equidistance
      !! LobattoGaussJacobi
      !! LobattoGaussChebyshev
      !! LobattoGaussGegenbauer
      !! GaussJacobi
      !! GaussChebyshev
      !! GaussGegenbauer
    REAL(DFP), ALLOCATABLE :: ans(:, :)
      !! barycentric coordinates, in xiJ format
      !! size(ans,1) = 3 corresponding to b0, b1, b2
      !! size(ans,2) total number of points
    CHARACTER(*), OPTIONAL, INTENT(IN) :: domain
      !! unit
      !! Biunit
      !! Equilateral
  END FUNCTION RecursiveNode2D
END INTERFACE
order

Order of element.

ipType

Interpolation point type. Following values are allowed.

  • Equidistance
  • GaussJacobi
  • GaussJacobiLobatto
  • GaussChebyshev
  • GaussChebyshevLobatto
  • GaussLegendre
  • GaussLegendreLobatto
  • GaussUltraspherical
  • GaussUltrasphericalLobatto
domain

It specifies the domain of the element. It is an optional argument. It can take following values:

  • UNIT, unit segment [0,1][0,1], in this case SIZE(ans,1) is 2.
  • BIUNIT, biunit segment [−1,1][-1,1], in this case SIZE(ans, 1) is 2.
  • BARYCENTRIC, in this case SIZE(ans,1) is 3. This is also the default value.

Examples​