CellBasis

This function returns the cell basis for hexahedron element.

Interface

܀ Interface

 INTERFACE CellBasis_Hexahedron
  MODULE PURE FUNCTION CellBasis_Hexahedron1(  &
    & n1, &
    & n2, &
    & n3, &
    & x,   &
    & y,   &
    & z)   &
    & RESULT(ans)
    INTEGER(I4B), INTENT(IN) :: n1
    !! order along axis 1
    !! it should be greater than 1
    INTEGER(I4B), INTENT(IN) :: n2
    !! order along axis 2
    !! it should be greater than 1
    INTEGER(I4B), INTENT(IN) :: n3
    !! order along axis 3
    !! it should be greater than 1
    REAL(DFP), INTENT(IN) :: x(:), y(:), z(:)
    !! point of evaluation
    !! these points should be between [-1, 1].
    REAL(DFP) :: ans( &
      & SIZE(x), &
      & (n1 - 1_I4B) * (n2 - 1_I4B) * (n3 - 1_I4B))
  END FUNCTION CellBasis_Hexahedron1
END INTERFACE CellBasis_Hexahedron

️܀ See example

PROGRAM main
USE easifemBase
USE easifemClasses
IMPLICIT NONE

REAL(DFP), ALLOCATABLE :: xxx(:, :, :), yyy(:, :, :), zzz(:, :, :), &
& x(:), y(:), z(:), basisValue(:, :), val(:, :, :)
TYPE(VTKPlot_) :: avtk
TYPE(string) :: label(1)
INTEGER(I4B) :: ii

x = linspace(-1.0_DFP, 1.0_DFP, 11)
y = linspace(-1.0_DFP, 1.0_DFP, 11)
z = linspace(-1.0_DFP, 1.0_DFP, 11)
CALL MeshGrid(xxx, yyy, zzz, x, y, z)
x = RESHAPE(xxx, [SIZE(xxx)])
y = RESHAPE(yyy, [SIZE(yyy)])
z = RESHAPE(zzz, [SIZE(zzz)])

basisValue = CellBasis_Hexahedron(2, 2, 2, x, y, z)

DO ii = 1, SIZE(label)
  label(ii) = tostring(ii)
END DO

CALL avtk%Plot(xxx, yyy, zzz, basisValue, label, "cellBasis.vts")

END PROGRAM main

↢