MappedScalar2D Derived Type

  1. SELF_MappedScalar_2D.f90
  2. SELF_MappedScalar_2D
  3. MappedScalar2D

type, public, extends(MappedScalar2D_t) :: MappedScalar2D


Inherits

type~~mappedscalar2d~2~~InheritsGraph type~mappedscalar2d~2 MappedScalar2D c_ptr c_ptr type~mappedscalar2d~2->c_ptr jas_gpu, mortarBuff_gpu, halo_sendbuf_gpu, halo_recvbuf_gpu type~mappedscalar2d_t MappedScalar2D_t type~mappedscalar2d~2->type~mappedscalar2d_t type~semquad SEMQuad type~mappedscalar2d_t->type~semquad geometry type~scalar2d Scalar2D type~mappedscalar2d_t->type~scalar2d type~semquad->type~scalar2d nScale, J type~vector2d Vector2D type~semquad->type~vector2d x, nHat type~tensor2d Tensor2D type~semquad->type~tensor2d dxds, dsdx type~scalar2d_t Scalar2D_t type~scalar2d->type~scalar2d_t type~vector2d_t Vector2D_t type~vector2d->type~vector2d_t type~self_dataobj SELF_DataObj type~scalar2d_t->type~self_dataobj type~tensor2d_t Tensor2D_t type~tensor2d->type~tensor2d_t type~vector2d_t->type~self_dataobj type~lagrange Lagrange type~self_dataobj->type~lagrange interp EquationParser EquationParser type~self_dataobj->EquationParser eqn type~metadata Metadata type~self_dataobj->type~metadata meta type~tensor2d_t->type~self_dataobj type~lagrange_t Lagrange_t type~lagrange->type~lagrange_t

Contents

Variables

Type-Bound Procedures

Source Code

MappedScalar2D

Components

TypeVisibilityAttributesNameInitial
integer, public :: M
integer, public :: N
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:):: avgBoundary
character(len=3), public :: backend ="cpu"
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:):: boundary
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:):: boundarynormal
type(EquationParser), public, allocatable:: eqn(:)
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:):: extBoundary
type(SEMQuad), public, pointer:: geometry=> null()
logical, public :: geometry_associated =.false.
integer, public :: halo_inflight =0
integer, public :: halo_nactive =0
type(c_ptr), public :: halo_recvbuf_gpu =c_null_ptr
integer, public, allocatable:: halo_reqs(:)
type(c_ptr), public :: halo_sendbuf_gpu =c_null_ptr
logical, public :: halo_static_done =.false.
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:):: interior
type(Lagrange), public, pointer:: interp
type(c_ptr), public :: jas_gpu
type(Metadata), public, allocatable:: meta(:)
real(kind=prec), public, allocatable, dimension(:,:,:,:):: mortarBuff
type(c_ptr), public :: mortarBuff_gpu =c_null_ptr
integer, public :: nElem
integer, public :: nVar

Type-Bound Procedures

procedure, public :: AssociateGeometry => AssociateGeometry_MappedScalar2D_t

procedure, public :: AverageSides => AverageSides_Scalar2D_t

procedure, public :: BoundaryInterp => BoundaryInterp_Scalar2D_t

procedure, public :: DissociateGeometry => DissociateGeometry_MappedScalar2D_t

procedure, public :: Free => Free_MappedScalar2D

generic, public :: Gradient => Gradient_Scalar2D_t

  • public subroutine Gradient_Scalar2D_t(this, df)

    Arguments

    TypeIntentOptionalAttributesName
    class(Scalar2D_t), intent(in) :: this
    real(kind=prec), intent(inout) :: df(1:this%N+1,1:this%N+1,1:this%nelem,1:this%nvar,1:2)

generic, public :: GridInterp => GridInterp_Scalar2D_t

  • public subroutine GridInterp_Scalar2D_t(this, f)

    Arguments

    TypeIntentOptionalAttributesName
    class(Scalar2D_t), intent(in) :: this
    real(kind=prec), intent(inout) :: f(1:this%M+1,1:this%M+1,1:this%nelem,1:this%nvar)

procedure, public :: Init => Init_MappedScalar2D

  • public subroutine Init_MappedScalar2D(this, interp, nVar, nElem)

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(out) :: this
    type(Lagrange), intent(in), target:: interp
    integer, intent(in) :: nVar
    integer, intent(in) :: nElem

procedure, private :: MPIMortarExchangeAsync => MPIMortarExchangeAsync_MappedScalar2D

  • public subroutine MPIMortarExchangeAsync_MappedScalar2D(this, mesh)

    GPU-resident analogue of the base-class mortar message posting; messages are posted on device memory (GPU-aware MPI), following MPIExchangeAsync.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(Mesh2D), intent(inout) :: mesh

generic, public :: MappedDGGradient => MappedDGGradient_MappedScalar2D_t, MappedDGGradient_MappedScalar2D

procedure, private :: MappedDGGradient_MappedScalar2D

generic, public :: MappedGradient => MappedGradient_MappedScalar2D_t, MappedGradient_MappedScalar2D

  • public subroutine MappedGradient_MappedScalar2D_t(this, df)

    Calculates the gradient of a function using the strong form of the gradient in mapped coordinates.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D_t), intent(in) :: this
    real(kind=prec), intent(out) :: df(1:this%N+1,1:this%N+1,1:this%nelem,1:this%nvar,1:2)

  • public subroutine MappedGradient_MappedScalar2D(this, df)

    Calculates the gradient of a function using the strong form of the gradient in mapped coordinates.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(c_ptr), intent(out) :: df

procedure, private :: MappedGradient_MappedScalar2D

  • public subroutine MappedGradient_MappedScalar2D(this, df)

    Calculates the gradient of a function using the strong form of the gradient in mapped coordinates.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(c_ptr), intent(out) :: df

procedure, public :: MortarExchange => MortarExchange_MappedScalar2D

  • public subroutine MortarExchange_MappedScalar2D(this, mesh)

    GPU implementation of the mortar exchange (see the base class for the algorithm) : traces are staged, reoriented, restricted, and projected in device memory with the SELF_Mortar kernels.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(Mesh2D), intent(inout) :: mesh

procedure, public :: SetDescription => SetDescription_DataObj

  • public subroutine SetDescription_DataObj(this, ivar, description)

    Set the description of the ivar-th variable

    Arguments

    TypeIntentOptionalAttributesName
    class(SELF_DataObj), intent(inout) :: this
    integer, intent(in) :: ivar
    character, intent(in) :: description

generic, public :: SetEquation => SetEquation_DataObj

  • public subroutine SetEquation_DataObj(this, ivar, eqnChar)

    Sets the equation parser for the ivar-th variable

    Arguments

    TypeIntentOptionalAttributesName
    class(SELF_DataObj), intent(inout) :: this
    integer, intent(in) :: ivar
    character, intent(in) :: eqnChar

procedure, public :: SetInteriorFromEquation => SetInteriorFromEquation_MappedScalar2D

  • public subroutine SetInteriorFromEquation_MappedScalar2D(this, geometry, time)

    Sets the this % interior attribute using the eqn attribute, geometry (for physical positions), and provided simulation time.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(SEMQuad), intent(in) :: geometry
    real(kind=prec), intent(in) :: time

procedure, public :: SetName => SetName_DataObj

  • public subroutine SetName_DataObj(this, ivar, name)

    Set the name of the ivar-th variable

    Arguments

    TypeIntentOptionalAttributesName
    class(SELF_DataObj), intent(inout) :: this
    integer, intent(in) :: ivar
    character, intent(in) :: name

procedure, public :: SetUnits => SetUnits_DataObj

  • public subroutine SetUnits_DataObj(this, ivar, units)

    Set the units of the ivar-th variable

    Arguments

    TypeIntentOptionalAttributesName
    class(SELF_DataObj), intent(inout) :: this
    integer, intent(in) :: ivar
    character, intent(in) :: units

procedure, public :: SideExchange => SideExchange_MappedScalar2D

procedure, public :: SideExchangeFinish => SideExchangeFinish_MappedScalar2D

  • public subroutine SideExchangeFinish_MappedScalar2D(this, mesh)

    Complete the aggregated halo exchange started by SideExchangeStart: wait on the persistent requests and unpack the received traces into extBoundary, applying side flips. The host blocks in MPI_Waitall (which also drives MPI progress) while previously launched device kernels continue to execute. No-op if no exchange is in flight.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(Mesh2D), intent(inout) :: mesh

procedure, public :: SideExchangeStart => SideExchangeStart_MappedScalar2D

  • public subroutine SideExchangeStart_MappedScalar2D(this, mesh, nactive)

    Begin the aggregated halo exchange and launch the local (same-rank) side exchange kernel.

    Read more…

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedScalar2D), intent(inout) :: this
    type(Mesh2D), intent(inout) :: mesh
    integer, intent(in), optional :: nactive

procedure, public :: UpdateDevice => UpdateDevice_Scalar2D_t

procedure, public :: UpdateHost => UpdateHost_Scalar2D_t

generic, public :: WriteHDF5 => WriteHDF5_MPI_Scalar2D_t, WriteHDF5_Scalar2D_t

  • public subroutine WriteHDF5_MPI_Scalar2D_t(this, fileId, group, elemoffset, nglobalelem)

    Arguments

    TypeIntentOptionalAttributesName
    class(Scalar2D_t), intent(in) :: this
    integer(kind=HID_T), intent(in) :: fileId
    character, intent(in) :: group
    integer, intent(in) :: elemoffset
    integer, intent(in) :: nglobalelem

  • public subroutine WriteHDF5_Scalar2D_t(this, fileId, group)

    Arguments

    TypeIntentOptionalAttributesName
    class(Scalar2D_t), intent(in) :: this
    integer(kind=HID_T), intent(in) :: fileId
    character, intent(in) :: group

Source Code

  type,extends(MappedScalar2D_t),public :: MappedScalar2D

    type(c_ptr) :: jas_gpu ! jacobian weighted scalar for gradient calculation
    type(c_ptr) :: mortarBuff_gpu = c_null_ptr ! mortar trace staging (lazy allocation)

    ! Packed device buffers and persistent MPI requests for the aggregated
    ! halo exchange. The side tables are shared across fields and live on
    ! mesh%decomp; the buffers and requests are per-field (sized by this
    ! field's variable count) and created lazily on the first exchange.
    type(c_ptr) :: halo_sendbuf_gpu = c_null_ptr ! packed device send buffer
    type(c_ptr) :: halo_recvbuf_gpu = c_null_ptr ! packed device receive buffer
    integer,allocatable :: halo_reqs(:) ! persistent requests; receives in 1:nnbr, sends in nnbr+1:2*nnbr
    integer :: halo_nactive = 0 ! variable count baked into halo_reqs
    integer :: halo_inflight = 0 ! variable count of the exchange in flight (0 = none)
    logical :: halo_static_done = .false. ! all variables have been exchanged at least once

  contains
    procedure,public :: Init => Init_MappedScalar2D
    procedure,public :: Free => Free_MappedScalar2D

    procedure,public :: SetInteriorFromEquation => SetInteriorFromEquation_MappedScalar2D

    procedure,public :: SideExchange => SideExchange_MappedScalar2D
    procedure,public :: SideExchangeStart => SideExchangeStart_MappedScalar2D
    procedure,public :: SideExchangeFinish => SideExchangeFinish_MappedScalar2D

    procedure,public :: MortarExchange => MortarExchange_MappedScalar2D
    procedure,private :: MPIMortarExchangeAsync => MPIMortarExchangeAsync_MappedScalar2D

    generic,public :: MappedGradient => MappedGradient_MappedScalar2D
    procedure,private :: MappedGradient_MappedScalar2D

    generic,public :: MappedDGGradient => MappedDGGradient_MappedScalar2D
    procedure,private :: MappedDGGradient_MappedScalar2D

  endtype MappedScalar2D