scale_mesh_rectdom2d Module Reference

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

Data Types
type	meshrectdom2d

Functions/Subroutines
subroutine	meshrectdom2d_init (this, negx, negy, dom_xmin, dom_xmax, dom_ymin, dom_ymax, isperiodicx, isperiodicy, refelem, nlocalmeshperprc, nprcx, nprcy, nproc, myrank)
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 )

Definition at line 75 of file scale_mesh_rectdom2d.F90.

    
    implicit none
 
    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
 
    !-----------------------------------------------------------------------------
    
    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 205 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
    
    !--
 
    lcmesh%Ne  = nex * ney
    lcmesh%Nv  = (nex + 1)*(ney + 1)
    lcmesh%NeS = 1
    lcmesh%NeE = lcmesh%Ne
    lcmesh%NeA = lcmesh%Ne + 2*(nex + ney)
 
    lcmesh%NeX = nex
    lcmesh%NeY = 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
 
    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) )
 
    lcmesh%BCType(:,:) = bctype_interior
 
    !----
 
    call meshutil2d_genrectdomain( lcmesh%pos_ev, lcmesh%EToV,   & ! (out)
      lcmesh%NeX, lcmesh%xmin, lcmesh%xmax,                      & ! (in)
      lcmesh%NeY, lcmesh%ymin, lcmesh%ymax )                       ! (in)
    
    !---
    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)
 
    !---
    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)
    
    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 359 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().