scale_mesh_hierarchy_2d Module Reference

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

Data Types
type	meshhierarchy2d
	Derived type for mesh hierarchy in 2D domain. More...
type	meshhierarchylevel2d
	Derived type to represent mesh hierarchy level in 2D domain. More...
type	meshhierarchylocalmgdata2d
	Derived type to manage 2D local mesh data for multigrid. More...

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

Detailed Description

module FElib / Mesh / 2D domain

Description

Manage mesh hierarchy of 2D domain for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshhierarchy2d_init()

subroutine scale_mesh_hierarchy_2d::meshhierarchy2d_init	(	class(meshhierarchy2d), intent(inout), target	this,
		class(meshbase2d), 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, intent(in)	h_level_num )

Initialize an object for mesh hierarchy in 2D 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

Definition at line 105 of file scale_mesh_hierarchy_2d.F90.

    implicit none
    class(MeshHierarchy2D), intent(inout), target :: this
    class(MeshBase2D), 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 :: poly_lev
    integer :: h_lev
 
    integer :: ldom_id
 
    type(MeshHierarchyLevel2D), pointer :: level_ptr
    class(MeshBase2D), pointer :: mesh2D_ptr
 
    class(MeshPtr2D), pointer :: fine_mesh_ptr
    class(MeshPtr2D), pointer :: coarse_mesh_ptr
    !-------------------------------------------------------------
 
    call meshhierarchybase_init( this, p_level_num, h_level_num )
 
    !- Setup p-mesh hierarchy
 
    if ( porder_list(p_level_num) /= 1 .and. h_level_num > 0 ) then
      log_info('MeshHierarchy2D_Init',*) 'Currently, only p=1 is supported for the coarsest p-mesh level when h_LEVEL_NUM > 0. Check!'
      call prc_abort
    end if
 
    !- Setup p-mesh hierarchy
 
    allocate( this%elem2D_list(p_level_num) )
    do poly_lev=1, p_level_num
      call this%elem2D_list(poly_lev)%Init( 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 (meshrectdom2d)
        call construct_rectdom2d_mesh( this, this%p_mesh_list(poly_lev)%ptr, &
          parent_mesh%NeGX, parent_mesh%NeGY,     &
          parent_mesh, this%elem2D_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
 
    allocate( this%h_mesh_list(this%NUM_hMG_LEVEL) )
 
    this%h_mesh_list(mesh_hierarchy_hmg_finest_level)%ptr => this%p_mesh_list(this%NUM_pMG_LEVEL)%ptr
    do h_lev=mesh_hierarchy_hmg_finest_level + 1, this%NUM_hMG_LEVEL
      select type(parent_mesh)
      type is (meshrectdom2d)
        call construct_rectdom2d_mesh( this, this%h_mesh_list(h_lev)%ptr, &
          negx_list(h_lev), negy_list(h_lev),                             &
          parent_mesh, this%elem2D_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
 
    return

References scale_mesh_hierarchy_base::mesh_hierarchy_hmg_finest_level, scale_mesh_hierarchy_base::mesh_hierarchy_pmg_finest_level, scale_mesh_hierarchy_base::mesh_hierarchy_type_hmg, scale_mesh_hierarchy_base::mesh_hierarchy_type_pmg, scale_mesh_hierarchy_base::meshhierarchy_construct_pmg_mat1d(), scale_mesh_hierarchy_base::meshhierarchybase_final(), scale_mesh_hierarchy_base::meshhierarchybase_init(), scale_mesh_hierarchy_base::meshhierarchylocalmgdatabase_final(), and scale_mesh_hierarchy_base::meshhierarchylocalmgdatabase_init().