scale_mesh_hierarchy_base Module Reference

module FElib / Mesh / Hierarchy base More...

Data Types
type	meshhierarchybase
	Base type for mesh hierarchy. More...
type	meshhierarchylevelbase
type	meshhierarchylocalmgdatabase
	Base type to manage local mesh data for multigrid. More...

Functions/Subroutines
subroutine, public	meshhierarchybase_init (this, p_level_num, h_level_num)
	Initialize a base object for mesh hierarchy.
subroutine, public	meshhierarchybase_final (this)
	Finalize a base object for mesh hierarchy.
subroutine, public	meshhierarchy_construct_pmg_mat1d (mat1d, np_i, np_o)
	Construct a 1D p-multigrid matrix used for transfer between different p-levels.
subroutine, public	meshhierarchylocalmgdatabase_init (this, lcmesh)
	Initialize a base object to manage local multigrid data for mesh hierarchy.
subroutine, public	meshhierarchylocalmgdatabase_final (this)
	Finalize a base object to manage local multigrid data for mesh hierarchy.

Variables
integer, parameter, public	mesh_hierarchy_pmg_finest_level = 1
	Finest level index in p-MG.
integer, parameter, public	mesh_hierarchy_hmg_finest_level = 1
	Finest level index in h-MG.
integer, parameter, public	mesh_hierarchy_type_pmg = 1
	Type ID of mesh hierarchy: p-MG.
integer, parameter, public	mesh_hierarchy_type_hmg = 2
	Type ID of mesh hierarchy: h-MG.

Detailed Description

module FElib / Mesh / Hierarchy base

Description

A base module to provide a mesh hierarchy

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshhierarchybase_init()

subroutine, public scale_mesh_hierarchy_base::meshhierarchybase_init	(	class(meshhierarchybase), intent(inout)	this,
		integer, intent(in)	p_level_num,
		integer, intent(in)	h_level_num )

Initialize a base object for mesh hierarchy.

Parameters

[in]	p_level_num	Number of p-mesh levels
[in]	h_level_num	Number of h-mesh levels

Definition at line 83 of file scale_mesh_hierarchy_base.F90.

    implicit none
    class(MeshHierarchyBase), intent(inout) :: this
    integer, intent(in) :: p_LEVEL_NUM   !< Number of p-mesh levels
    integer, intent(in) :: h_LEVEL_NUM   !< Number of h-mesh levels
    !-------------------------------------------------------------
 
    this%NUM_TOT_LEVEL = 0
    this%NUM_hMG_LEVEL = h_level_num
    this%NUM_pMG_LEVEL = p_level_num
 
    this%NUM_TOT_LEVEL = this%NUM_hMG_LEVEL + this%NUM_pMG_LEVEL
    return

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

meshhierarchybase_final()

subroutine, public scale_mesh_hierarchy_base::meshhierarchybase_final

(

class(meshhierarchybase), intent(inout)

this

)

Finalize a base object for mesh hierarchy.

Definition at line 102 of file scale_mesh_hierarchy_base.F90.

    implicit none
    class(MeshHierarchyBase), intent(inout) :: this
    !-------------------------------------------------------------
    return

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

meshhierarchy_construct_pmg_mat1d()

subroutine, public scale_mesh_hierarchy_base::meshhierarchy_construct_pmg_mat1d	(	real(rp), dimension(np_o,np_i), intent(out)	mat1d,
		integer, intent(in)	np_i,
		integer, intent(in)	np_o )

Construct a 1D p-multigrid matrix used for transfer between different p-levels.

Parameters

[in]	np_i	Number of points with input DOF
[in]	np_o	Number of points with output DOF
[out]	mat1d	Transfer matrix with p-MG

Definition at line 111 of file scale_mesh_hierarchy_base.F90.

    use scale_polynomial, only: &
      polynomial_genlagrangepoly, &
      polynomial_gengausslegendrept,        &
      polynomial_gengausslegendreptintweight
    use scale_element_line, only: lineelement
    implicit none
    integer, intent(in) :: np_i               !< Number of points with input DOF
    integer, intent(in) :: np_o               !< Number of points with output DOF
    real(RP), intent(out) :: mat1D(np_o,np_i) !< Transfer matrix with p-MG
 
    type(LineElement) :: elem1D_i
    type(LineElement) :: elem1D_o
 
    real(RP), allocatable :: lag_i(:,:)
    real(RP), allocatable :: lag_o(:,:)
 
    integer :: np_int
    real(RP), allocatable :: int_w(:)
    real(RP), allocatable :: int_pts(:)
 
    integer :: pi, po
    !---------------------------------------------------------------------
 
    call elem1d_i%Init( np_i-1, .false. )
    call elem1d_o%Init( np_o-1, .false. )
 
    !-
    np_int = max(np_i,np_o)
    
    allocate( int_pts(np_int) )
    int_pts(:) = polynomial_gengausslegendrept(np_int)
 
    allocate( int_w(np_int) )     
    int_w(:) = polynomial_gengausslegendreptintweight(np_int)
 
    !-
    allocate( lag_i(np_int,np_i) )
    lag_i(:,:) = polynomial_genlagrangepoly(elem1d_i%PolyOrder, elem1d_i%x1, int_pts)
 
    allocate( lag_o(np_int,np_o) )
    lag_o(:,:) = polynomial_genlagrangepoly(elem1d_o%PolyOrder, elem1d_o%x1, int_pts)
 
    do pi=1, np_i
    do po=1, np_o
      mat1d(po,pi) = sum( int_w(:) * lag_i(:,pi) * lag_o(:,po) )
    end do
    end do
    mat1d(:,:) = matmul( elem1d_o%invM, mat1d )
 
    !-
    call elem1d_i%Final()
    call elem1d_o%Final()
 
    return

References scale_polynomial::polynomial_gengausslegendrept(), scale_polynomial::polynomial_gengausslegendreptintweight(), and scale_polynomial::polynomial_genlagrangepoly().

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

meshhierarchylocalmgdatabase_init()

subroutine, public scale_mesh_hierarchy_base::meshhierarchylocalmgdatabase_init	(	class(meshhierarchylocalmgdatabase), intent(inout)	this,
		class(localmeshbase), intent(in)	lcmesh )

Initialize a base object to manage local multigrid data for mesh hierarchy.

Definition at line 171 of file scale_mesh_hierarchy_base.F90.

    implicit none
    class(MeshHierarchyLocalMGDataBase), intent(inout) :: this
    class(LocalMeshBase), intent(in) :: lcmesh
    !-------------------------------------------------------------
 
    allocate( this%Ic2f(lcmesh%refElem%Np,lcmesh%refElem%Nv,lcmesh%NeA) )
    allocate( this%Ic2f_emap(lcmesh%Ne) )
 
    allocate( this%CoarseLocalMesh_tileIDlist(1) )
    allocate( this%CoarseLocalMesh_lcdomIDlist(1) )
    return

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

meshhierarchylocalmgdatabase_final()

subroutine, public scale_mesh_hierarchy_base::meshhierarchylocalmgdatabase_final

(

class(meshhierarchylocalmgdatabase), intent(inout)

this

)

Finalize a base object to manage local multigrid data for mesh hierarchy.

Definition at line 188 of file scale_mesh_hierarchy_base.F90.

    implicit none
    class(MeshHierarchyLocalMGDataBase), intent(inout) :: this
    !------------------------------------------------------------- 
    deallocate( this%Ic2f )
    deallocate( this%Ic2f_emap )
    deallocate( this%CoarseLocalMesh_tileIDlist )
    deallocate( this%CoarseLocalMesh_lcdomIDlist )
    return

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

Variable Documentation

mesh_hierarchy_pmg_finest_level

integer, parameter, public scale_mesh_hierarchy_base::mesh_hierarchy_pmg_finest_level = 1

Finest level index in p-MG.

Definition at line 65 of file scale_mesh_hierarchy_base.F90.

  integer, public, parameter :: MESH_HIERARCHY_pMG_FINEST_LEVEL = 1  !< Finest level index in p-MG

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

mesh_hierarchy_hmg_finest_level

integer, parameter, public scale_mesh_hierarchy_base::mesh_hierarchy_hmg_finest_level = 1

Finest level index in h-MG.

Definition at line 66 of file scale_mesh_hierarchy_base.F90.

  integer, public, parameter :: MESH_HIERARCHY_hMG_FINEST_LEVEL = 1  !< Finest level index in h-MG

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

mesh_hierarchy_type_pmg

integer, parameter, public scale_mesh_hierarchy_base::mesh_hierarchy_type_pmg = 1

Type ID of mesh hierarchy: p-MG.

Definition at line 68 of file scale_mesh_hierarchy_base.F90.

  integer, public, parameter :: MESH_HIERARCHY_TYPE_pMG         = 1  !< Type ID of mesh hierarchy: p-MG

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().

mesh_hierarchy_type_hmg

integer, parameter, public scale_mesh_hierarchy_base::mesh_hierarchy_type_hmg = 2

Type ID of mesh hierarchy: h-MG.

Definition at line 69 of file scale_mesh_hierarchy_base.F90.

  integer, public, parameter :: MESH_HIERARCHY_TYPE_hMG         = 2  !< Type ID of mesh hierarchy: h-MG

Referenced by scale_mesh_hierarchy_2d::meshhierarchy2d_init(), and scale_mesh_hierarchy_3d::meshhierarchy3d_init().