Skip to main content

FEVariable example 33c

!!! note "" This example tests * operator. Matrix, Nodal

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