scale_mesh_base Module Reference

module FElib / Mesh / Base More...

Data Types
type	meshbase
interface	meshbase_get_localmesh
type	meshdiminfo

Functions/Subroutines
subroutine, public	meshbase_init (this, ndimtype, refelem, nlocalmeshperprc, nsidetile, nprocs)
subroutine, public	meshbase_final (this)
subroutine, public	meshbase_setgeometricinfo (mesh, ndim)

Detailed Description

module FElib / Mesh / Base

Description

Base module to mangage meshes for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshbase_init()

subroutine, public scale_mesh_base::meshbase_init	(	class(meshbase), intent(inout)	this,
		integer, intent(in)	ndimtype,
		class(elementbase), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in)	nsidetile,
		integer, intent(in), optional	nprocs )

Definition at line 90 of file scale_mesh_base.F90.

    
    use scale_prc, only: &
      prc_nprocs
    implicit none
 
    class(MeshBase), intent(inout) :: this
    integer, intent(in) :: ndimtype
    class(ElementBase), intent(in), target :: refElem
    integer, intent(in) :: NLocalMeshPerPrc
    integer, intent(in) :: NsideTile
    integer, intent(in), optional :: nprocs
 
    integer :: n
    !-----------------------------------------------------------------------------
    
    if (present(nprocs)) then
      this%PRC_NUM = nprocs
    else
      this%PRC_NUM = prc_nprocs
    end if
 
    this%LOCAL_MESH_NUM        = nlocalmeshperprc
    this%LOCAL_MESH_NUM_global = this%PRC_NUM * this%LOCAL_MESH_NUM
 
    this%refElem => refelem
        
    allocate( this%tileID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tileFaceID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tilePanelID_globalMap(nsidetile, this%LOCAL_MESH_NUM_global) )
    allocate( this%tileID_global2localMap(this%LOCAL_MESH_NUM_global) )
    allocate( this%PRCRank_globalMap(this%LOCAL_MESH_NUM_global) )
    allocate( this%dimInfo(ndimtype) )
    
    this%isGenerated = .false.
 
    return

Referenced by scale_mesh_base1d::meshbase1d_init(), scale_mesh_base2d::meshbase2d_init(), scale_mesh_base3d::meshbase3d_init(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().

meshbase_final()

subroutine, public scale_mesh_base::meshbase_final

(

class(meshbase), intent(inout)

this

)

Definition at line 132 of file scale_mesh_base.F90.

    implicit none
    class(MeshBase), intent(inout) :: this
    !-----------------------------------------------------------------------------
  
    if ( allocated(this%tileID_globalMap) ) then
      deallocate( this%tileID_globalMap )
      deallocate( this%tileFaceID_globalMap )
      deallocate( this%tilePanelID_globalMap )
      deallocate( this%tileID_global2localMap )
      deallocate( this%PRCRank_globalMap )
      deallocate( this%dimInfo )
    end if
 
    return

Referenced by scale_mesh_base1d::meshbase1d_final(), scale_mesh_base2d::meshbase2d_final(), scale_mesh_base3d::meshbase3d_final(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().

meshbase_setgeometricinfo()

subroutine, public scale_mesh_base::meshbase_setgeometricinfo	(	class(localmeshbase), intent(inout)	mesh,
		integer, intent(in)	ndim )

Definition at line 150 of file scale_mesh_base.F90.

    implicit none
    
    class(LocalMeshBase), intent(inout) :: mesh
    integer, intent(in) :: ndim
 
    class(ElementBase), pointer :: refElem
    !-----------------------------------------------------------------------------
 
    refelem => mesh%refElem
 
    allocate( mesh%pos_en(refelem%Np,mesh%Ne,ndim) )
    !allocate( mesh%fx(refElem%Nfaces*refElem%Nfp,mesh%Ne) )
    !allocate( mesh%fy(refElem%Nfaces*refElem%Nfp,mesh%Ne) )
    allocate( mesh%normal_fn(refelem%NfpTot,mesh%Ne,ndim) )
    allocate( mesh%sJ(refelem%NfpTot,mesh%Ne) )
    allocate( mesh%J(refelem%Np,mesh%Ne) )
    allocate( mesh%Fscale(refelem%NfpTot,mesh%Ne) )
    allocate( mesh%Escale(refelem%Np,mesh%Ne,ndim,ndim) )
    allocate( mesh%Gsqrt(refelem%Np,mesh%Ne) )
    allocate( mesh%G_ij(refelem%Np,mesh%Ne, ndim,ndim) )
    allocate( mesh%GIJ (refelem%Np,mesh%Ne, ndim,ndim) )
 
    return

Referenced by scale_mesh_base1d::meshbase1d_setgeometricinfo(), and scale_mesh_base::meshbase_get_localmesh::meshbase_get_localmesh().