module FElib / Mesh / 3D domain More...

Data Types
type	meshhierarchy3d
	Derived type for mesh hierarchy in 3D domain. More...
type	meshhierarchylevel3d
	Derived type to represent mesh hierarchy level in 3D domain. More...
type	meshhierarchylocalmgdata3d
	Derived type to manage 3D local mesh data for multigrid. More...

Functions/Subroutines
subroutine	meshhierarchy3d_init (this, parent_mesh, porder_list, p_level_num, negx_list, negy_list, negz_list, h_level_num)
	Initialize an object for mesh hierarchy in 3D domain.

Detailed Description

module FElib / Mesh / 3D domain

Description: Manage mesh hierarchy of 3D domain for element-based methods

Author: Yuta Kawai, Team SCALE

Function/Subroutine Documentation

◆ meshhierarchy3d_init()

subroutine scale_mesh_hierarchy_3d::meshhierarchy3d_init	(	class(meshhierarchy3d), intent(inout), target	this,
		class(meshbase3d), intent(in), target	parent_mesh,
		integer, dimension(p_level_num), intent(in)	porder_list,
		integer, intent(in)	p_level_num,
		integer, dimension(h_level_num), intent(in)	negx_list,
		integer, dimension(h_level_num), intent(in)	negy_list,
		integer, dimension(h_level_num), intent(in)	negz_list,
		integer, intent(in)	h_level_num )

Initialize an object for mesh hierarchy in 3D domain.

Parameters

[in]	parent_mesh	Pointer to the finest mesh
[in]	p_level_num	Number of p-mesh levels
[in]	porder_list	Polynomial order list for p-mesh levels
[in]	h_level_num	Number of h-mesh levels
[in]	negx_list	Number of elements in x-direction for h-mesh levels
[in]	negy_list	Number of elements in y-direction for h-mesh levels
[in]	negz_list	Number of elements in z-direction for h-mesh levels

Definition at line 105 of file scale_mesh_hierarchy_3d.F90.

    implicit none
    class(MeshHierarchy3D), intent(inout), target :: this
    class(MeshBase3D), intent(in), target :: parent_mesh  !< Pointer to the finest mesh
    integer, intent(in) :: p_LEVEL_NUM                    !< Number of p-mesh levels
    integer, intent(in) :: porder_list(p_LEVEL_NUM)       !< Polynomial order list for p-mesh levels
    integer, intent(in) :: h_LEVEL_NUM                    !< Number of h-mesh levels
    integer, intent(in) :: NeGX_list(h_LEVEL_NUM)         !< Number of elements in x-direction for h-mesh levels
    integer, intent(in) :: NeGY_list(h_LEVEL_NUM)         !< Number of elements in y-direction for h-mesh levels
    integer, intent(in) :: NeGZ_list(h_LEVEL_NUM)         !< Number of elements in z-direction for h-mesh levels
 
    integer :: poly_lev
    integer :: h_lev
 
    integer :: ldom_id
 
    type(MeshHierarchyLevel3D), pointer :: level_ptr
    class(MeshBase3D), pointer :: mesh3D_ptr
 
    class(MeshPtr3D), pointer :: fine_mesh_ptr
    class(MeshPtr3D), pointer :: coarse_mesh_ptr
 
    integer :: NeGZ
    !-------------------------------------------------------------
 
    call meshhierarchybase_init( this, p_level_num, h_level_num )
 
    !--
    if ( parent_mesh%refElem3D%PolyOrder_h /= parent_mesh%refElem3D%PolyOrder_v ) then
      log_info('MeshHierarchy3D_Init',*) 'Currently, PolyOrder_h should equal PolyOrder_v. Check!'
      call prc_abort
    end if
    if ( porder_list(p_level_num) /= 1 .and. h_level_num > 0 ) then
      log_info('MeshHierarchy3D_Init',*) 'Currently, only p=1 is supported for the coarsest p-mesh level. Check!'
      call prc_abort
    end if
    if ( h_level_num > 0 ) then
      negz = negz_list(1)
      do h_lev=2, h_level_num
        if ( negz_list(h_lev) /= negz ) then
          log_info('MeshHierarchy3D_Init',*) 'Currently, the number of elements should be same in the vertical direction. Check!'
          call prc_abort
        end if
      end do
    end if
 
    !- Setup p-mesh hierarchy
 
    allocate( this%elem3D_list(p_level_num) )
    do poly_lev=1, p_level_num
      call this%elem3D_list(poly_lev)%Init( porder_list(poly_lev), porder_list(poly_lev), .false. )
    end do
 
    allocate( this%p_mesh_list(this%NUM_pMG_LEVEL) )
 
    this%p_mesh_list(mesh_hierarchy_pmg_finest_level)%ptr => parent_mesh
    do poly_lev=mesh_hierarchy_pmg_finest_level + 1, this%NUM_pMG_LEVEL
      select type(parent_mesh)
      type is (meshcubedom3d)
        call construct_cubedom3d_mesh( this, this%p_mesh_list(poly_lev)%ptr,    &
          parent_mesh%NeGX, parent_mesh%NeGY, parent_mesh%NeGZ, parent_mesh%FZ, &
          parent_mesh, this%elem3D_list(poly_lev) )
      end select
    end do
 
    allocate( this%p_level(this%NUM_pMG_LEVEL) )
 
    do poly_lev=mesh_hierarchy_pmg_finest_level, this%NUM_pMG_LEVEL
      call this%p_level(poly_lev)%Init( poly_lev, this%p_mesh_list(poly_lev)%ptr, &
        this%p_mesh_list, this%NUM_pMG_LEVEL, mesh_hierarchy_type_pmg )
    end do
 
    !- Setup h-mesh hierarchy
 
    if ( this%NUM_hMG_LEVEL > 0 ) then
      allocate( this%h_mesh_list(this%NUM_hMG_LEVEL) )
 
      this%h_mesh_list(mesh_hierarchy_hmg_finest_level)%ptr => parent_mesh
      do h_lev=mesh_hierarchy_hmg_finest_level + 1, this%NUM_hMG_LEVEL
        select type(parent_mesh)
        type is (meshcubedom3d)
          call construct_cubedom3d_mesh( this, this%h_mesh_list(h_lev)%ptr,       &
            negx_list(h_lev), negy_list(h_lev), negz_list(h_lev), parent_mesh%FZ, &
            parent_mesh, this%elem3D_list(p_level_num) )
        end select
      end do
 
      allocate( this%h_level(this%NUM_hMG_LEVEL) )
 
      do h_lev=mesh_hierarchy_hmg_finest_level, this%NUM_hMG_LEVEL
        call this%h_level(h_lev)%Init( h_lev, this%h_mesh_list(h_lev)%ptr, &
          this%h_mesh_list, this%NUM_hMG_LEVEL, mesh_hierarchy_type_hmg )
      end do
    end if
 
    return

Data Types

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ meshhierarchy3d_init()