Skip to main content

Overview

RefLine_ defines the reference element for line element. This element will be used for constructing finite elements on lines (1D). RefLine_ is a subclass of AbstractRefElem_.

  • To see the full structure of this class go to RefLine_.

Initiate

Use Initiate to construct an instance . An example is shows below:

PROGRAM main
USE easifemBase
USE easifemClasses
IMPLICIT NONE
TYPE(RefLine_) :: obj
REAL(DFP) :: avar, xij(1, 2)
CALL obj%Initiate(nsd=1, baseInterpol="Lagrange", baseContinuity="H1")
CALL display(obj%MdEncode(), 'Line2 : '//char_lf//char_lf)
END PROGRAM main
See results

Line2 :

Element typeLine2
Xidimension1
NSD1
tPoints2
tLines1
tSurfaces0
tVolumes0
BaseContinuityH1
BaseInterpolationLagrangeInterpolation

Nodal Coordinates:

x-11

PointTopology( 1 ) :

Element typePoint1
Xidimension0
Nptrs1

PointTopology( 2 ) :

Element typePoint1
Xidimension0
Nptrs2

EdgeTopology( 1 ) :

Element typeLine2
Xidimension1
Nptrs12

Getting the name

To get the element name use GetName.

Getting facet elements

To get the facet elements use GetFacetElements. An example is given below:

PROGRAM main
USE easifemBase
USE easifemClasses
IMPLICIT NONE
TYPE(RefLine_) :: obj
TYPE(AbstractRefElementPointer_), ALLOCATABLE :: faceElem(:)
INTEGER(I4B) :: ii
TYPE(String) :: astr

REAL(DFP) :: avar, xij(1, 2)
CALL obj%Initiate(nsd=1, baseInterpol="Lagrange", baseContinuity="H1")
CALL obj%GetFacetElements(faceElem)

DO ii = 1, SIZE(faceElem)
astr = "**Facet Element "//tostring(ii)//"**"//char_lf//char_lf// &
& faceElem(ii)%ptr%MdEncode()
CALL display(astr, "")
END DO
END PROGRAM main
See results

Facet Element 1

Element typePoint1
Xidimension0
NSD1
tPoints1
tLines0
tSurfaces0
tVolumes0
BaseContinuityH1
BaseInterpolationLagrangeInterpolation

Nodal Coordinates:

x-1

PointTopology( 1 ) :

Element typePoint1
Xidimension0
Nptrs1

Facet Element 2

Element typePoint1
Xidimension0
NSD1
tPoints1
tLines0
tSurfaces0
tVolumes0
BaseContinuityH1
BaseInterpolationLagrangeInterpolation

Nodal Coordinates:

x1

PointTopology( 1 ) :

Element typePoint1
Xidimension0
Nptrs1

Getting nodal coordinates

To get nodal coordinates of element use GetNodeCoord.