Skip to main content

FEVariable example 41d

!!! note "" This example tests / operator. Scalar Quadrature

Use association​

  • [[FEVariable_]]

Usage​

!!! note "Import modules and declare variable"

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

!!! note "constant / space"

  call display("=================================")
  obj = QuadratureVariable( 1.0_DFP, &
    & typeFEVariableScalar,    &
    & typeFEVariableConstant ) &
    / QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
    & typeFEVariableScalar, &
    & typeFEVariableSpace )
  CALL Display(obj, "constnat / space")

!!! note "space / space"

  call display("=================================")
  obj = QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
      & typeFEVariableScalar, &
      & typeFEVariableSpace ) &
      /  QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
      & typeFEVariableScalar, &
      & typeFEVariableSpace )
  CALL Display(obj, "space / space")

!!! note "constant / time"

  call display("=================================")
  obj = QuadratureVariable( 1.0_DFP, &
    & typeFEVariableScalar,    &
    & typeFEVariableConstant ) &
    /  QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
    & typeFEVariableScalar, &
    & typeFEVariableTime )
  CALL Display(obj, "constant / time")

!!! note "time / time"

  call display("=================================")
  obj = QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
      & typeFEVariableScalar, &
      & typeFEVariableTime ) &
      / QuadratureVariable( arange(1.0_DFP, 3.0_DFP), &
      & typeFEVariableScalar, &
      & typeFEVariableTime )
  CALL Display(obj, "time / time")

!!! note "constant / spacetime"

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

!!! note "spacetime / spacetime"

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