MappedScalar2D Derived Type

  1. SELF_MappedScalar_2D.f90
  2. SELF_MappedScalar_2D
  3. MappedScalar2D

type, public, extends(MappedScalar2D_t) :: MappedScalar2D


Inherits

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

Contents

Variables

Type-Bound Procedures

Source Code

MappedScalar2D

Components

TypeVisibilityAttributesNameInitial
integer, public :: M
integer, public :: N
integer(kind=c_size_t), public :: alloc_jas =0
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

High-water-mark backing store for the arrays above (AMR Stage 6b; see SELF_DataPool). The public members are pointers remapped onto the leading part of these pools at the exact logical shape, so every extent, stride and shape() stays exactly as before while Resize can reuse the storage across an adaptation epoch. Nothing should index the pools.
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(:)

Bytes currently allocated for jas_gpu, so Resize can reuse it (Stage 6b). The capacity counters for the five arrays inherited from Scalar2D are declared there.
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
real(kind=prec), public, pointer, contiguous:: pool_avgBoundary(:)=> null()
real(kind=prec), public, pointer, contiguous:: pool_boundary(:)=> null()
real(kind=prec), public, pointer, contiguous:: pool_boundarynormal(:)=> null()
real(kind=prec), public, pointer, contiguous:: pool_extBoundary(:)=> null()
real(kind=prec), public, pointer, contiguous:: pool_interior(:)=> null()

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

procedure, public :: MapArrays => MapArrays_Scalar2D_t

  • public subroutine MapArrays_Scalar2D_t(this, Np, nVar, nElem)

    Size the backing pools for (Np,nVar,nElem) and remap the public arrays onto them at the exact logical shape. Bounds-remapping a rank-1 contiguous target to a rank-4 pointer is what keeps every extent and stride identical to a plain allocate while allowing the pool underneath to be larger and reused. See SELF_DataPool.

    Arguments

    TypeIntentOptionalAttributesName
    class(Scalar2D_t), intent(inout) :: this
    integer, intent(in) :: Np
    integer, intent(in) :: nVar
    integer, intent(in) :: nElem

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 :: Resize => Resize_MappedScalar2D

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

    Rebind to a new element count, reusing host pools and device buffers where they fit (AMR Stage 6b), and without uploading zeros the way Init does.

    Read more…

    Arguments

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

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
    !! Bytes currently allocated for jas_gpu, so Resize can reuse it (Stage 6b). The capacity
    !! counters for the five arrays inherited from Scalar2D are declared there.
    integer(c_size_t) :: alloc_jas = 0
    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 :: Resize => Resize_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