MappedVector3D Derived Type

  1. SELF_MappedVector_3D.f90
  2. SELF_MappedVector_3D
  3. MappedVector3D

type, public, extends(MappedVector3D_t) :: MappedVector3D


Inherits

type~~mappedvector3d~2~~InheritsGraph type~mappedvector3d~2 MappedVector3D c_ptr c_ptr type~mappedvector3d~2->c_ptr halo_sendbuf_gpu, halo_recvbuf_gpu, mortarBuff_gpu type~mappedvector3d_t MappedVector3D_t type~mappedvector3d~2->type~mappedvector3d_t type~semhex SEMHex type~mappedvector3d_t->type~semhex geometry type~vector3d~2 Vector3D type~mappedvector3d_t->type~vector3d~2 type~semhex->type~vector3d~2 x, nHat, xMesh type~tensor3d~2 Tensor3D type~semhex->type~tensor3d~2 dxds, dsdx type~scalar3d~2 Scalar3D type~semhex->type~scalar3d~2 nScale, J type~lagrange Lagrange type~semhex->type~lagrange meshToModel type~vector3d_t Vector3D_t type~vector3d~2->type~vector3d_t type~tensor3d_t Tensor3D_t type~tensor3d~2->type~tensor3d_t type~self_dataobj SELF_DataObj type~vector3d_t->type~self_dataobj type~scalar3d_t Scalar3D_t type~scalar3d~2->type~scalar3d_t type~lagrange_t Lagrange_t type~lagrange->type~lagrange_t type~tensor3d_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 type~scalar3d_t->type~self_dataobj

Contents

Variables

Type-Bound Procedures

Source Code

MappedVector3D

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

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() is unchanged while Resize can reuse the storage across an adaptation epoch. Nothing should index the pools directly.
type(EquationParser), public, allocatable:: eqn(:)
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:,:,:):: extBoundary
type(SEMHex), public, pointer:: geometry=> null()
logical, public :: geometry_associated =.false.
type(c_ptr), public :: halo_recvbuf_gpu =c_null_ptr
type(c_ptr), public :: halo_sendbuf_gpu =c_null_ptr
real(kind=prec), public, pointer, contiguous, dimension(:,:,:,:,:,:):: interior
type(Lagrange), public, pointer:: interp
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_MappedVector3D_t

procedure, public :: AverageSides => AverageSides_Vector3D_t

procedure, public :: BoundaryInterp => BoundaryInterp_Vector3D_t

generic, public :: Curl => Curl_Vector3D_t

  • public subroutine Curl_Vector3D_t(this, curlf)

    Arguments

    TypeIntentOptionalAttributesName
    class(Vector3D_t), intent(in) :: this
    real(kind=prec), intent(out) :: curlf(1:this%N+1,1:this%N+1,1:this%N+1,1:this%nelem,1:this%nvar,1:3)

procedure, public :: DissociateGeometry => DissociateGeometry_MappedVector3D_t

generic, public :: Divergence => Divergence_Vector3D_t

  • public subroutine Divergence_Vector3D_t(this, df)

    Arguments

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

procedure, public :: Free => Free_MappedVector3D

generic, public :: Gradient => Gradient_Vector3D_t

  • public subroutine Gradient_Vector3D_t(this, df)

    Arguments

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

generic, public :: GridInterp => GridInterp_Vector3D_t

  • public subroutine GridInterp_Vector3D_t(this, f)

    Arguments

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

    (Output) Array of function values, defined on the target grid

procedure, public :: Init => Init_Vector3D_t

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

    Arguments

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

procedure, public :: MPIExchangeAsync => MPIExchangeAsync_MappedVector3D

  • public subroutine MPIExchangeAsync_MappedVector3D(this, mesh)

    Post the aggregated halo exchange: one MPI_Irecv/MPI_Isend pair per neighboring rank, carrying every (side,variable,component) boundary trace shared with that rank in a single packed device buffer. The boundary array is laid out with the component index outermost, so the pack/unpack kernels treat the vector as 3*nvar scalar variables. Packed buffers are allocated on first use; the shared side tables are built by SideExchange before this is called.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(Mesh3D), intent(inout) :: mesh

procedure, private :: MPIMortarExchangeAsync => MPIMortarExchangeAsync_MappedVector3D

  • public subroutine MPIMortarExchangeAsync_MappedVector3D(this, mesh)

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

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(Mesh3D), intent(inout) :: mesh

procedure, private :: MPIMortarFluxAsync => MPIMortarFluxAsync_MappedVector3D

  • public subroutine MPIMortarFluxAsync_MappedVector3D(this, mesh)

    Posts the one-directional messages for MortarFluxCollect on device memory : each remote small face sends its boundaryNormal trace to the big face's rank.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(Mesh3D), intent(inout) :: mesh

procedure, public :: MapArrays => MapArrays_Vector3D_t

  • public subroutine MapArrays_Vector3D_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. See SELF_DataPool for why this is done with pointer remapping rather than by over-allocating the element dimension.

    Arguments

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

generic, public :: MappedDGDivergence => MappedDGDivergence_MappedVector3D_t, MappedDGDivergence_MappedVector3D

  • public subroutine MappedDGDivergence_MappedVector3D_t(this, df)

    Computes the divergence of a 3-D vector using the weak form On input, the attribute of the vector is assigned and the attribute is set to the physical directions of the vector. This method will project the vector onto the contravariant basis vectors.

    Arguments

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

  • public subroutine MappedDGDivergence_MappedVector3D(this, df)

    Computes the divergence of a 3-D vector using the weak form On input, the attribute of the vector is assigned and the attribute is set to the physical directions of the vector. This method will project the vector onto the contravariant basis vectors.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(in) :: this
    type(c_ptr), intent(inout) :: df

procedure, private :: MappedDGDivergence_MappedVector3D

  • public subroutine MappedDGDivergence_MappedVector3D(this, df)

    Computes the divergence of a 3-D vector using the weak form On input, the attribute of the vector is assigned and the attribute is set to the physical directions of the vector. This method will project the vector onto the contravariant basis vectors.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(in) :: this
    type(c_ptr), intent(inout) :: df

generic, public :: MappedDivergence => MappedDivergence_MappedVector3D_t, MappedDivergence_MappedVector3D

procedure, private :: MappedDivergence_MappedVector3D

procedure, public :: MortarExchange => MortarExchange_MappedVector3D

  • public subroutine MortarExchange_MappedVector3D(this, mesh)

    GPU implementation of the vector mortar exchange; the kernels treat the (variable, direction) pairs as 3*nvar independent trace lines.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(Mesh3D), intent(inout) :: mesh

procedure, public :: MortarFluxCollect => MortarFluxCollect_MappedVector3D

  • public subroutine MortarFluxCollect_MappedVector3D(this, mesh)

    GPU implementation of MortarFluxCollect (see the base class for the algorithm and conservation statement). Stages the small faces' boundaryNormal traces in the mortar buffer, then overwrites the big face's integrand on device.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(Mesh3D), intent(inout) :: mesh

procedure, public :: Resize => Resize_MappedVector3D

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

    Rebind to a new element count, reusing host pools and device buffers where they fit (AMR Stage 6b). Inherits the Vector3D resize (host pools plus the five device buffers) and adds this class's mortar buffers, which are sized by mesh%nMortars rather than nElem: nMortars changes with the mesh independently of the element count, so a stale buffer would silently under-size the next mortar exchange. They are invalidated here and lazily re-created at the correct size on next use.

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), 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, SetEquation_Vector3D_t

  • 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

  • public subroutine SetEquation_Vector3D_t(this, idir, ivar, eqnChar)

    Sets the equation parser for the idir direction and ivar-th variable

    Arguments

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

procedure, public :: SetInteriorFromEquation => SetInteriorFromEquation_MappedVector3D

  • public subroutine SetInteriorFromEquation_MappedVector3D(this, geometry, time)

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

    Arguments

    TypeIntentOptionalAttributesName
    class(MappedVector3D), intent(inout) :: this
    type(SEMHex), 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_MappedVector3D

procedure, public :: UpdateDevice => UpdateDevice_Vector3D_t

procedure, public :: UpdateHost => UpdateHost_Vector3D_t

generic, public :: WriteHDF5 => WriteHDF5_MPI_Vector3D_t, WriteHDF5_Vector3D_t

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

    Arguments

    TypeIntentOptionalAttributesName
    class(Vector3D_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_Vector3D_t(this, fileId, group)

    Arguments

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

Source Code

  type,extends(MappedVector3D_t),public :: MappedVector3D

    ! Packed device buffers for the aggregated MPI halo exchange. The side
    ! tables are shared across fields and live on mesh%decomp; the buffers
    ! are per-field (sized for 3*nvar variables, since all three vector
    ! components are exchanged) and allocated 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
    type(c_ptr) :: mortarBuff_gpu = c_null_ptr ! mortar trace staging (lazy allocation)

  contains

    procedure,public :: Resize => Resize_MappedVector3D
    procedure,public :: Free => Free_MappedVector3D
    procedure,public :: SideExchange => SideExchange_MappedVector3D
    procedure,public :: MPIExchangeAsync => MPIExchangeAsync_MappedVector3D
    procedure,public :: MortarExchange => MortarExchange_MappedVector3D
    procedure,private :: MPIMortarExchangeAsync => MPIMortarExchangeAsync_MappedVector3D
    procedure,public :: MortarFluxCollect => MortarFluxCollect_MappedVector3D
    procedure,private :: MPIMortarFluxAsync => MPIMortarFluxAsync_MappedVector3D

    generic,public :: MappedDivergence => MappedDivergence_MappedVector3D
    procedure,private :: MappedDivergence_MappedVector3D

    generic,public :: MappedDGDivergence => MappedDGDivergence_MappedVector3D
    procedure,private :: MappedDGDivergence_MappedVector3D

    procedure,public :: SetInteriorFromEquation => SetInteriorFromEquation_MappedVector3D

  endtype MappedVector3D