GetFacetElements

This routine returns the face elements of a reference element.

Interface 1

Interface

INTERFACE GetFacetElements
  MODULE SUBROUTINE refelem_GetFacetElements1(refelem, ans)
    CLASS(ReferenceElement_), INTENT(IN) :: refelem
    TYPE(ReferenceElement_), INTENT(INOUT) :: ans(:)
  END SUBROUTINE refelem_GetFacetElements1
END INTERFACE GetFacetElements

example

In this example we generate face elements of a triangle by using the method called FacetElements.

PROGRAM main
  USE easifemBase
  TYPE( ReferenceTriangle_ ) :: obj
  TYPE( ReferenceElement_ ) :: facetElems( 3 )
  INTEGER( I4B ) :: ii

  CALL Initiate( obj, nsd = 2 )


  CALL GetFacetElements(obj, facetElems)

  DO ii = 1, SIZE( facetElems )
    CALL Display( facetElems( ii ),  &
      & "facetElements( " // tostring(ii) // " ) = " )
    CALL Blanklines( NOL = 2 )
  END DO
END PROGRAM main

↢ close

Interface 2

Interface

 INTERFACE GetFacetElements
  MODULE SUBROUTINE refelem_GetFacetElements2(elemType, nsd, ans)
    INTEGER(I4B), INTENT(IN) :: elemType
    INTEGER(I4B), INTENT(IN) :: nsd
    TYPE(ReferenceElement_), INTENT(INOUT) :: ans(:)
  END SUBROUTINE refelem_GetFacetElements2
END INTERFACE GetFacetElements

example

This example tests the FacetElements method.

PROGRAM main
  USE easifemBase


  block
    TYPE( ReferenceElement_ ) :: facetElems(3)
    INTEGER( I4B ) :: ii

    CALL GetFacetElements(elemType=Triangle3, nsd=2, ans=facetElems)

    DO ii = 1, SIZE( facetElems )
      CALL Display( facetElems( ii ),  &
        & "facetElements( " // tostring(ii) // " ) = " )
      CALL Blanklines( NOL = 2 )
    END DO

  end block
END PROGRAM main

↢ close