scale_mesh_rectdom2d Module Reference

module FElib / Mesh / Rectangle 2D domain More...

Data Types
type	meshrectdom2d
	Derived type to manage a rectangular 2D computational domain. More...

Functions/Subroutines
subroutine	meshrectdom2d_init (this, negx, negy, dom_xmin, dom_xmax, dom_ymin, dom_ymax, isperiodicx, isperiodicy, refelem, nlocalmeshperprc, nprcx, nprcy, nproc, myrank)
	Initialize an object to manage a 2D rectangular computational mesh.
subroutine, public	meshrectdom2d_setuplocaldom (lcmesh, tileid, panelid, i, j, nprcx, nprcy, dom_xmin, dom_xmax, dom_ymin, dom_ymax, nex, ney)
subroutine, public	meshrectdom2d_coord_conv (x, y, xr, xs, yr, ys, vx, vy, elem)

Detailed Description

module FElib / Mesh / Rectangle 2D domain

Description

Manage mesh data of rectangle 2D domain for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshrectdom2d_init()

subroutine scale_mesh_rectdom2d::meshrectdom2d_init	(	class(meshrectdom2d), intent(inout)	this,
		integer, intent(in)	negx,
		integer, intent(in)	negy,
		real(rp), intent(in)	dom_xmin,
		real(rp), intent(in)	dom_xmax,
		real(rp), intent(in)	dom_ymin,
		real(rp), intent(in)	dom_ymax,
		logical, intent(in)	isperiodicx,
		logical, intent(in)	isperiodicy,
		type(quadrilateralelement), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in)	nprcx,
		integer, intent(in)	nprcy,
		integer, intent(in), optional	nproc,
		integer, intent(in), optional	myrank )

Initialize an object to manage a 2D rectangular computational mesh.

Parameters

[in]	negx	Number of elements in X direction (global)
[in]	negy	Number of elements in Y direction (global)
[in]	dom_xmin	Minimum X coordinate of the global domain
[in]	dom_xmax	Maximum X coordinate of the global domain
[in]	dom_ymin	Minimum Y coordinate of the global domain
[in]	dom_ymax	Maximum Y coordinate of the global domain
[in]	isperiodicx	Flag whether the domain is periodic in X direction
[in]	isperiodicy	Flag whether the domain is periodic in Y direction
[in]	refelem	Reference element for the 2D mesh
[in]	nlocalmeshperprc	Number of local meshes managed by each process
[in]	nprcx	Number of processes in X direction for domain decomposition
[in]	nprcy	Number of processes in Y direction for domain decomposition
[in]	nproc	Total number of processes (if not provided, it will be determined from the parallel environment)
[in]	myrank	Rank of the current process (if not provided, it will be determined from the parallel environment)

Definition at line 81 of file scale_mesh_rectdom2d.F90.

    
    implicit none
 
    class(MeshRectDom2D), intent(inout) :: this
    integer, intent(in) :: NeGX          !< Number of elements in X direction (global)
    integer, intent(in) :: NeGY          !< Number of elements in Y direction (global)
    real(RP), intent(in) :: dom_xmin     !< Minimum X coordinate of the global domain
    real(RP), intent(in) :: dom_xmax     !< Maximum X coordinate of the global domain
    real(RP), intent(in) :: dom_ymin     !< Minimum Y coordinate of the global domain
    real(RP), intent(in) :: dom_ymax     !< Maximum Y coordinate of the global domain
    logical, intent(in) :: isPeriodicX   !< Flag whether the domain is periodic in X direction
    logical, intent(in) :: isPeriodicY   !< Flag whether the domain is periodic in Y direction
    type(QuadrilateralElement), intent(in), target :: refElem !< Reference element for the 2D mesh
    integer, intent(in) :: NLocalMeshPerPrc !< Number of local meshes managed by each process
    integer, intent(in) :: NprcX            !< Number of processes in X direction for domain decomposition
    integer, intent(in) :: NprcY            !< Number of processes in Y direction for domain decomposition
    integer, intent(in), optional :: nproc  !< Total number of processes (if not provided, it will be determined from the parallel environment)
    integer, intent(in), optional :: myrank !< Rank of the current process (if not provided, it will be determined from the parallel environment)
 
    !-----------------------------------------------------------------------------
    
    this%NeGX = negx
    this%NeGY = negy
 
    this%xmin_gl       = dom_xmin
    this%xmax_gl       = dom_xmax
    this%ymin_gl       = dom_ymin
    this%ymax_gl       = dom_ymax
    this%dom_vol = (this%xmax_gl - this%xmin_gl) * (this%ymax_gl - this%ymin_gl)
 
    this%isPeriodicX = isperiodicx
    this%isPeriodicY = isperiodicy
 
    this%NprcX = nprcx
    this%NprcY = nprcy
 
    call meshbase2d_init( this, refelem, nlocalmeshperprc, &
      nproc, myrank )
 
    return

References scale_mesh_base2d::meshbase2d_final(), scale_mesh_base2d::meshbase2d_init(), and meshrectdom2d_setuplocaldom().

meshrectdom2d_setuplocaldom()

subroutine, public scale_mesh_rectdom2d::meshrectdom2d_setuplocaldom	(	type(localmesh2d), intent(inout)	lcmesh,
		integer, intent(in)	tileid,
		integer, intent(in)	panelid,
		integer, intent(in)	i,
		integer, intent(in)	j,
		integer, intent(in)	nprcx,
		integer, intent(in)	nprcy,
		real(rp), intent(in)	dom_xmin,
		real(rp), intent(in)	dom_xmax,
		real(rp), intent(in)	dom_ymin,
		real(rp), intent(in)	dom_ymax,
		integer, intent(in)	nex,
		integer, intent(in)	ney )

Definition at line 212 of file scale_mesh_rectdom2d.F90.

 
    use scale_meshutil_2d, only: &
      meshutil2d_genconnectivity,   &
      meshutil2d_genrectdomain,     &
      meshutil2d_buildinteriormap,  &
      meshutil2d_genpatchboundarymap
 
    use scale_localmesh_base, only: bctype_interior
    implicit none
      
    type(LocalMesh2D), intent(inout) :: lcmesh
    integer, intent(in) :: tileID
    integer, intent(in) :: panelID
    integer, intent(in) :: i, j
    integer, intent(in) :: NprcX, NprcY
    real(RP), intent(in) :: dom_xmin, dom_xmax
    real(RP), intent(in) :: dom_ymin, dom_ymax
    integer, intent(in) ::NeX, NeY
    
    class(ElementBase2D), pointer :: elem
    real(RP) :: delx, dely
    
    !-----------------------------------------------------------------------------
 
    elem => lcmesh%refElem2D
    lcmesh%tileID = tileid
    lcmesh%panelID = panelid
    !$acc update device(lcmesh%tileID, lcmesh%panelID)
    
    !--
 
    lcmesh%Ne  = nex * ney
    lcmesh%Nv  = (nex + 1)*(ney + 1)
    lcmesh%NeS = 1
    lcmesh%NeE = lcmesh%Ne
    lcmesh%NeA = lcmesh%Ne + 2*(nex + ney)
    !$acc update device(lcmesh%Ne, lcmesh%Nv, lcmesh%NeS, lcmesh%NeE, lcmesh%NeA)
 
    lcmesh%NeX = nex
    lcmesh%NeY = ney
    !$acc update device(lcmesh%NeX, lcmesh%NeY)
 
    delx = (dom_xmax - dom_xmin)/dble(nprcx)
    dely = (dom_ymax - dom_ymin)/dble(nprcy)
 
    lcmesh%xmin = dom_xmin + (i-1)*delx
    lcmesh%xmax = dom_xmin +  i   *delx
    lcmesh%ymin = dom_ymin + (j-1)*dely
    lcmesh%ymax = dom_ymin +  j   *dely
    !$acc update device(lcmesh%xmin, lcmesh%xmax, lcmesh%ymin, lcmesh%ymax)
 
    allocate( lcmesh%pos_ev(lcmesh%Nv,2) )
    allocate( lcmesh%EToV(lcmesh%Ne,elem%Nv) )
    allocate( lcmesh%EToE(lcmesh%Ne,elem%Nfaces) )
    allocate( lcmesh%EToF(lcmesh%Ne,elem%Nfaces) )
    allocate( lcmesh%BCType(lcmesh%refElem%Nfaces,lcmesh%Ne) )
    allocate( lcmesh%VMapM(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%VMapP(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%MapM(elem%NfpTot, lcmesh%Ne) )
    allocate( lcmesh%MapP(elem%NfpTot, lcmesh%Ne) )
    !$acc enter data create(lcmesh%pos_ev, lcmesh%EToV, lcmesh%EToE, lcmesh%EToF, lcmesh%BCType, &
    !$acc   lcmesh%VMapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP)
 
    lcmesh%BCType(:,:) = bctype_interior
    !$acc update device(lcmesh%BCType)
 
    !----
 
    call meshutil2d_genrectdomain( lcmesh%pos_ev, lcmesh%EToV,   & ! (out)
      lcmesh%NeX, lcmesh%xmin, lcmesh%xmax,                      & ! (in)
      lcmesh%NeY, lcmesh%ymin, lcmesh%ymax )                       ! (in)
    !$acc update device(lcmesh%pos_ev, lcmesh%EToV)
    
    !---
    call meshbase2d_setgeometricinfo( lcmesh, meshrectdom2d_coord_conv, meshrectdom2d_calc_normal )
 
    !---
 
    call meshutil2d_genconnectivity( lcmesh%EToE, lcmesh%EToF, & ! (out)
      lcmesh%EToV, lcmesh%Ne, elem%Nfaces )                      ! (in)
    !$acc update device(lcmesh%EToE, lcmesh%EToF)
 
    !---
    call meshutil2d_buildinteriormap( lcmesh%VmapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP, &
      lcmesh%pos_en, lcmesh%pos_ev, lcmesh%EToE, lcmesh%EtoF, lcmesh%EtoV,                  &
      elem%Fmask, lcmesh%Ne, elem%Np, elem%Nfp, elem%Nfaces, lcmesh%Nv )
 
    call meshutil2d_genpatchboundarymap( lcmesh%VMapB, lcmesh%MapB, lcmesh%VMapP, &
      lcmesh%pos_en, lcmesh%xmin, lcmesh%xmax, lcmesh%ymin, lcmesh%ymax,          &
      elem%Fmask, lcmesh%Ne, elem%Np, elem%Nfp, elem%Nfaces, lcmesh%Nv)
    !$acc update device(lcmesh%VMapM, lcmesh%VMapP, lcmesh%MapM, lcmesh%MapP)
    !$acc enter data copyin(lcmesh%VMapB, lcmesh%MapB)
    
    return

References scale_localmesh_base::bctype_interior, scale_mesh_base2d::meshbase2d_setgeometricinfo(), meshrectdom2d_coord_conv(), scale_meshutil_2d::meshutil2d_buildglobalmap(), scale_meshutil_2d::meshutil2d_buildinteriormap(), scale_meshutil_2d::meshutil2d_genconnectivity(), scale_meshutil_2d::meshutil2d_genpatchboundarymap(), and scale_meshutil_2d::meshutil2d_genrectdomain().

Referenced by scale_mesh_rectdom2d::meshrectdom2d::assigndomid(), scale_mesh_cubedom3d::meshcubedom3d_init(), and meshrectdom2d_init().

meshrectdom2d_coord_conv()

subroutine, public scale_mesh_rectdom2d::meshrectdom2d_coord_conv	(	real(rp), dimension(elem%np), intent(out)	x,
		real(rp), dimension(elem%np), intent(out)	y,
		real(rp), dimension(elem%np), intent(out)	xr,
		real(rp), dimension(elem%np), intent(out)	xs,
		real(rp), dimension(elem%np), intent(out)	yr,
		real(rp), dimension(elem%np), intent(out)	ys,
		real(rp), dimension(elem%nv), intent(in)	vx,
		real(rp), dimension(elem%nv), intent(in)	vy,
		type(elementbase2d), intent(in)	elem )

Definition at line 380 of file scale_mesh_rectdom2d.F90.

    implicit none
 
    type(ElementBase2D), intent(in) :: elem
    real(RP), intent(out) :: x(elem%Np), y(elem%Np)
    real(RP), intent(out) :: xr(elem%Np), xs(elem%Np), yr(elem%Np), ys(elem%Np)
    real(RP), intent(in) :: vx(elem%Nv), vy(elem%Nv)
 
    !-------------------------------------------------
    
    x(:) = vx(1) + 0.5_rp*(elem%x1(:) + 1.0_rp)*(vx(2) - vx(1))
    y(:) = vy(1) + 0.5_rp*(elem%x2(:) + 1.0_rp)*(vy(3) - vy(1))
 
    xr(:) = 0.5_rp*(vx(2) - vx(1)) !matmul(refElem%Dx1,mesh%x1(:,n))
    xs(:) = 0.0_rp                 !matmul(refElem%Dx2,mesh%x1(:,n))
    yr(:) = 0.0_rp                 !matmul(refElem%Dx1,mesh%x2(:,n))
    ys(:) = 0.5_rp*(vy(3) - vy(1)) !matmul(refElem%Dx2,mesh%x2(:,n))
 
    return

Referenced by scale_mesh_rectdom2d::meshrectdom2d::assigndomid(), and meshrectdom2d_setuplocaldom().