Skip to main content

FEVariable example 22b

!!! note "" This example tests subtraction operator. Vector, Nodal

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 = NodalVariable( arange(1.0_DFP, 3.0_DFP), &
    & typeFEVariableVector,    &
    & typeFEVariableConstant ) &
    - NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
    & typeFEVariableVector, &
    & typeFEVariableSpace )
  CALL Display(obj, "constant-space")

!!! note "space-space"

  call display("=================================")
  obj = NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
                      & typeFEVariableVector, &
                      & typeFEVariableSpace ) &
      - NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
                      & typeFEVariableVector, &
                      & typeFEVariableSpace )
  CALL Display(obj, "space-space")

!!! note "constant-time"

  call display("=================================")
  obj = NodalVariable( arange(1.0_DFP, 3.0_DFP), &
    & typeFEVariableVector,    &
    & typeFEVariableConstant ) &
    - NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
    & typeFEVariableVector, &
    & typeFEVariableTime )
  CALL Display(obj, "constant-time")

!!! note "time-time"

  call display("=================================")
  obj = NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
      & typeFEVariableVector, &
      & typeFEVariableTime ) &
      - NodalVariable( reshape(arange(1.0_DFP, 6.0_DFP),[3,2]), &
      & typeFEVariableVector, &
      & typeFEVariableTime )
  CALL Display(obj, "time-time")

!!! note "constant-spacetime"

  call display("=================================")
  obj = NodalVariable( arange(1.0_DFP, 3.0_DFP), &
    & typeFEVariableVector,    &
    & typeFEVariableConstant ) &
    - NodalVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
    & typeFEVariableVector, &
    & typeFEVariableSpaceTime )
  CALL Display(obj, "constant-spacetime")

!!! note "spacetime-spacetime"

  call display("=================================")
  obj = NodalVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
    & typeFEVariableVector, &
    & typeFEVariableSpaceTime ) &
    - NodalVariable( reshape(arange(1.0_DFP, 12.0_DFP), [3,2,2]), &
    & typeFEVariableVector, &
    & typeFEVariableSpaceTime )
  CALL Display(obj, "spacetime-spacetime")
END PROGRAM main