scale_multigrid_fieldset_base Module Reference

module FElib / Multigrid / Field set base More...

Data Types
type	mgfieldset2d
	Derived type for 2D multigrid field set. More...
type	mgfieldset3d
	Derived type for 3D multigrid field set. More...

Functions/Subroutines
subroutine	mgfieldset2d_init (this, mesh2d, aux_scalar_num, aux_vec_num, level_id)
	Initialize an object to manage multigrid field set in 2D.

Detailed Description

module FElib / Multigrid / Field set base

Description

Manage field sets for multigrid method

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

mgfieldset2d_init()

subroutine scale_multigrid_fieldset_base::mgfieldset2d_init	(	class(mgfieldset2d), intent(inout), target	this,
		class(meshbase2d), intent(in)	mesh2d,
		integer, intent(in)	aux_scalar_num,
		integer, intent(in)	aux_vec_num,
		integer, intent(in)	level_id )

Initialize an object to manage multigrid field set in 2D.

Definition at line 98 of file scale_multigrid_fieldset_base.F90.

    use scale_mesh_rectdom2d, only: meshrectdom2d
    implicit none
    class(MGFieldSet2D), intent(inout), target :: this
    class(MeshBase2D), intent(in) :: mesh2D
    integer, intent(in) :: aux_scalar_num
    integer, intent(in) :: aux_vec_num
    integer, intent(in) :: level_id
 
    integer :: iv
 
    type(MeshFieldCommRectDom2D), pointer :: var_comm2D
    type(MeshFieldCommRectDom2D), pointer :: aux_comm2D
 
    type(MeshFieldCommCubedSphereDom2D), pointer :: var_comm2D_cs
    type(MeshFieldCommCubedSphereDom2D), pointer :: aux_comm2D_cs
    !---------------------------------
 
    this%level_id = level_id
    
    call this%dq%Init( "dq", "1", mesh2d )    
    call this%f%Init( "f", "1", mesh2d )
    call this%res%Init( "res", "1", mesh2d )
 
    this%aux_var_num_tot = aux_scalar_num + 2*aux_vec_num
    if ( this%aux_var_num_tot > 0 ) then
      allocate( this%aux_var(this%aux_var_num_tot) )
      do iv=1, this%aux_var_num_tot
        call this%aux_var(iv)%Init( "aux_var", "1", mesh2d )
      end do
    end if
    
    select type(mesh2d)
    type is (meshrectdom2d)
      allocate(var_comm2d)
      call var_comm2d%Init( 1, 0, 0, mesh2d )
      this%var_comm_ptr => var_comm2d
    type is (meshcubedspheredom2d)
      allocate(var_comm2d_cs)
      call var_comm2d_cs%Init( 1, 0, 0, mesh2d )
      this%var_comm_ptr => var_comm2d_cs
    end select
 
    if ( this%aux_var_num_tot > 0 ) then
      select type(mesh2d)
      type is (meshrectdom2d)
        allocate( aux_comm2d )
        call aux_comm2d%Init( aux_scalar_num, aux_vec_num, 0, mesh2d )
        this%aux_comm_ptr => aux_comm2d
      type is (meshcubedspheredom2d)
        allocate( aux_comm2d_cs )
        call aux_comm2d_cs%Init( aux_scalar_num, aux_vec_num, 0, mesh2d )
        this%aux_comm_ptr => aux_comm2d_cs
      end select
    end if
  
    call this%Dx%Init( mesh2d%refElem2D%Dx1, storage_format='ELL')
    call this%Dy%Init( mesh2d%refElem2D%Dx2, storage_format='ELL')
    call this%Lift%Init( mesh2d%refElem2D%Lift, storage_format='ELL')
    return