Initiate

This subroutine initiate an instance of ReferenceHexahedron.

Interface

܀ Interface

INTERFACE Initiate
  MODULE PURE SUBROUTINE initiate_ref_Hexahedron(obj, NSD, xij)
    CLASS(ReferenceHexahedron_), INTENT(INOUT) :: obj
    INTEGER(I4B), INTENT(IN) :: NSD
    REAL(DFP), INTENT(IN), OPTIONAL :: xij(:, :)
  END SUBROUTINE initiate_ref_Hexahedron
END INTERFACE Initiate

NSD

number of spatial dimension, it should be 3.

xij

xij is the nodal coordinate of hexahedron.

• total number of rows in xij should be 3 as each row denotes a spatial dimension.
• total number of columns in xij should be 8 as each column denotes the node number.

️܀ See example

PROGRAM main
USE easifembase
USE easifemclasses
IMPLICIT NONE
TYPE(ReferenceHexahedron_) :: obj
REAL(dfp) :: avar
INTEGER(i4b), PARAMETER :: nsd = 3_I4B

CALL Initiate(obj=obj, nsd=nsd)
CALL Display(obj, 'Obj : ')

END PROGRAM main

See results

Obj : 
# ElemType : Hexahedron8
# XiDimension :: 3
# NSD : 3
# Order : 1
# EntityCounts(0) : 8
# EntityCounts(1) : 12
# EntityCounts(2) : 6
# EntityCounts(3) : 1
# Node( 1 ) : 
--------------
   -1.00000   
   -1.00000   
   -1.00000   
# Node( 2 ) : 
--------------
    1.00000   
   -1.00000   
   -1.00000   
# Node( 3 ) : 
--------------
    1.00000   
    1.00000   
   -1.00000   
# Node( 4 ) : 
--------------
   -1.00000   
    1.00000   
   -1.00000   
# Node( 5 ) : 
--------------
   -1.00000   
   -1.00000   
    1.00000   
# Node( 6 ) : 
--------------
    1.00000   
   -1.00000   
    1.00000   
# Node( 7 ) : 
--------------
   1.00000    
   1.00000    
   1.00000    
# Node( 8 ) : 
--------------
   -1.00000   
    1.00000   
    1.00000   
# Topology( 1 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    1     
# Topology( 2 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    2     
# Topology( 3 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    3     
# Topology( 4 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    4     
# Topology( 5 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    5     
# Topology( 6 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    6     
# Topology( 7 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    7     
# Topology( 8 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    8     
# Topology( 9 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    2     
# Topology( 10 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    4     
# Topology( 11 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    5     
# Topology( 12 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    3     
# Topology( 13 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    6     
# Topology( 14 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    4     
# Topology( 15 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    7     
# Topology( 16 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    4     
    8     
# Topology( 17 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    6     
# Topology( 18 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    8     
# Topology( 19 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    6     
    7     
# Topology( 20 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    7     
    8     
# Topology( 21 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    4     
    3     
    2     
# Topology( 22 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    5     
    6     
    7     
    8     
# Topology( 23 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    5     
    8     
    4     
# Topology( 24 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    2     
    3     
    7     
    6     
# Topology( 25 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    3     
    4     
    8     
    7     
# Topology( 26 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    2     
    6     
    5     
# Topology( 27 ) : 
# ElemType : Hexahedron8
# XiDim : 3
# Nptrs : 
----------
    1     
    2     
    3     
    4     
    5     
    6     
    7     
    8     

↢

Interface 2 (ReferenceHexahedron)

We can use ReferenceHexahedron function to construct an instance of ReferenceHexahedron.

܀ Interface

INTERFACE ReferenceHexahedron
  MODULE PURE FUNCTION reference_Hexahedron(NSD, xij) RESULT(obj)
    INTEGER(I4B), INTENT(IN) :: NSD
    REAL(DFP), INTENT(IN), OPTIONAL :: xij(:, :)
    TYPE(ReferenceHexahedron_) :: obj
  END FUNCTION reference_Hexahedron
END INTERFACE ReferenceHexahedron

️܀ See example

PROGRAM main
USE easifembase
USE easifemclasses
IMPLICIT NONE
TYPE(ReferenceHexahedron_) :: obj
REAL(dfp) :: avar
INTEGER(i4b), PARAMETER :: nsd = 3_I4B

obj = ReferenceHexahedron(nsd=nsd)
CALL Display(obj, 'Obj : ')

END PROGRAM main

See results

Obj : 
# ElemType : Hexahedron8
# XiDimension :: 3
# NSD : 3
# Order : 1
# EntityCounts(0) : 8
# EntityCounts(1) : 12
# EntityCounts(2) : 6
# EntityCounts(3) : 1
# Node( 1 ) : 
--------------
   -1.00000   
   -1.00000   
   -1.00000   
# Node( 2 ) : 
--------------
    1.00000   
   -1.00000   
   -1.00000   
# Node( 3 ) : 
--------------
    1.00000   
    1.00000   
   -1.00000   
# Node( 4 ) : 
--------------
   -1.00000   
    1.00000   
   -1.00000   
# Node( 5 ) : 
--------------
   -1.00000   
   -1.00000   
    1.00000   
# Node( 6 ) : 
--------------
    1.00000   
   -1.00000   
    1.00000   
# Node( 7 ) : 
--------------
   1.00000    
   1.00000    
   1.00000    
# Node( 8 ) : 
--------------
   -1.00000   
    1.00000   
    1.00000   
# Topology( 1 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    1     
# Topology( 2 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    2     
# Topology( 3 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    3     
# Topology( 4 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    4     
# Topology( 5 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    5     
# Topology( 6 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    6     
# Topology( 7 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    7     
# Topology( 8 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    8     
# Topology( 9 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    2     
# Topology( 10 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    4     
# Topology( 11 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    5     
# Topology( 12 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    3     
# Topology( 13 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    6     
# Topology( 14 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    4     
# Topology( 15 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    7     
# Topology( 16 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    4     
    8     
# Topology( 17 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    6     
# Topology( 18 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    8     
# Topology( 19 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    6     
    7     
# Topology( 20 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    7     
    8     
# Topology( 21 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    4     
    3     
    2     
# Topology( 22 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    5     
    6     
    7     
    8     
# Topology( 23 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    5     
    8     
    4     
# Topology( 24 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    2     
    3     
    7     
    6     
# Topology( 25 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    3     
    4     
    8     
    7     
# Topology( 26 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    2     
    6     
    5     
# Topology( 27 ) : 
# ElemType : Hexahedron8
# XiDim : 3
# Nptrs : 
----------
    1     
    2     
    3     
    4     
    5     
    6     
    7     
    8     

↢

Interface 3 (ReferenceHexahedron_Pointer)

The function ReferenceHexahedron_Pointer returns the pointer to newly created instance of ReferenceHexahedron.

܀ Interface

INTERFACE ReferenceHexahedron_Pointer
  MODULE PURE FUNCTION reference_Hexahedron_Pointer(NSD, xij) RESULT(obj)
    INTEGER(I4B), INTENT(IN) :: NSD
    REAL(DFP), INTENT(IN), OPTIONAL :: xij(:, :)
    CLASS(ReferenceHexahedron_), POINTER :: obj
  END FUNCTION reference_Hexahedron_Pointer
END INTERFACE ReferenceHexahedron_Pointer

️܀ See example

PROGRAM main
USE easifembase
USE easifemclasses
IMPLICIT NONE
CLASS(ReferenceElement_), pointer :: obj
REAL(dfp) :: avar
INTEGER(i4b), PARAMETER :: nsd = 3_I4B

obj => ReferenceHexahedron_Pointer(nsd=nsd)
CALL Display(obj, 'Obj : ')

END PROGRAM main

See results

Obj : 
# ElemType : Hexahedron8
# XiDimension :: 3
# NSD : 3
# Order : 1
# EntityCounts(0) : 8
# EntityCounts(1) : 12
# EntityCounts(2) : 6
# EntityCounts(3) : 1
# Node( 1 ) : 
--------------
   -1.00000   
   -1.00000   
   -1.00000   
# Node( 2 ) : 
--------------
    1.00000   
   -1.00000   
   -1.00000   
# Node( 3 ) : 
--------------
    1.00000   
    1.00000   
   -1.00000   
# Node( 4 ) : 
--------------
   -1.00000   
    1.00000   
   -1.00000   
# Node( 5 ) : 
--------------
   -1.00000   
   -1.00000   
    1.00000   
# Node( 6 ) : 
--------------
    1.00000   
   -1.00000   
    1.00000   
# Node( 7 ) : 
--------------
   1.00000    
   1.00000    
   1.00000    
# Node( 8 ) : 
--------------
   -1.00000   
    1.00000   
    1.00000   
# Topology( 1 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    1     
# Topology( 2 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    2     
# Topology( 3 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    3     
# Topology( 4 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    4     
# Topology( 5 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    5     
# Topology( 6 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    6     
# Topology( 7 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    7     
# Topology( 8 ) : 
# ElemType : Point1
# XiDim : 0
# Nptrs : 
----------
    8     
# Topology( 9 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    2     
# Topology( 10 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    4     
# Topology( 11 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    1     
    5     
# Topology( 12 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    3     
# Topology( 13 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    2     
    6     
# Topology( 14 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    4     
# Topology( 15 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    3     
    7     
# Topology( 16 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    4     
    8     
# Topology( 17 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    6     
# Topology( 18 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    5     
    8     
# Topology( 19 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    6     
    7     
# Topology( 20 ) : 
# ElemType : Line2
# XiDim : 1
# Nptrs : 
----------
    7     
    8     
# Topology( 21 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    4     
    3     
    2     
# Topology( 22 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    5     
    6     
    7     
    8     
# Topology( 23 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    5     
    8     
    4     
# Topology( 24 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    2     
    3     
    7     
    6     
# Topology( 25 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    3     
    4     
    8     
    7     
# Topology( 26 ) : 
# ElemType : Quadrangle4
# XiDim : 2
# Nptrs : 
----------
    1     
    2     
    6     
    5     
# Topology( 27 ) : 
# ElemType : Hexahedron8
# XiDim : 3
# Nptrs : 
----------
    1     
    2     
    3     
    4     
    5     
    6     
    7     
    8     

↢