Skip to main content

FEVariable example 43c

!!! note "" This example tests / operator. Matrix, Quadrature

Use association​

  • [[FEVariable_]]

Usage​

!!! note "Import modules and declare variable"

PROGRAM main
  USE easifemBase
  IMPLICIT NONE
  TYPE(FEVariable_) :: obj

!!! note "constant / constant"

  obj = QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant ) &
      / QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant )
  CALL Display(obj, "constant / constant")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant ) &
      / 1.0_DFP
  CALL Display(obj, "constant / constant")

!!! note "space / constant"

  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableSpace ) &
      / QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant )
  CALL Display(obj, "space / constant")
  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableSpace ) &
      / 1.0_DFP
  CALL Display(obj, "space / constant")

!!! note "time / constant"

  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableTime ) &
      / QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant )
  CALL Display(obj, "time / constant")
  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableTime ) &
      / 1.0_DFP
  CALL Display(obj, "time / constant")

!!! note "spacetime / constant"

  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 24.0_DFP), [3,2,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableSpaceTime ) &
      / QuadratureVariable( reshape(arange(1.0_DFP, 6.0_DFP), [3,2]), &
      & typeFEVariableMatrix,    &
      & typeFEVariableConstant )
  CALL Display(obj, "spacetime / constant")
  call display("=================================")
  obj = QuadratureVariable( reshape(arange(1.0_DFP, 24.0_DFP), [3,2,2,2]), &
      & typeFEVariableMatrix, &
      & typeFEVariableSpaceTime ) &
      / 1.0_DFP
  CALL Display(obj, "spacetime / constant")
END PROGRAM main