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