scale_mesh_linedom1d Module Reference

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

Data Types
type	meshlinedom1d

Functions/Subroutines
subroutine, public	meshlinedom1d_init (this, neg, dom_xmin, dom_xmax, refelem, nlocalmeshperprc, nproc, myrank, fx)
	Initialize an object to manage a 1D domain.
subroutine, public	meshlinedom1d_final (this)
	Finalize an object to manage a 1D domain.
subroutine, public	meshlinedom1d_generate (this)
	Generate meshes for the 1D domain.

Detailed Description

module FElib / Mesh / 1D domain

Description

Manage mesh data of 1D domain for element-based methods

Author

Yuta Kawai, Team SCALE

Function/Subroutine Documentation

meshlinedom1d_init()

subroutine, public scale_mesh_linedom1d::meshlinedom1d_init	(	class(meshlinedom1d), intent(inout)	this,
		integer, intent(in)	neg,
		real(rp), intent(in)	dom_xmin,
		real(rp), intent(in)	dom_xmax,
		type(lineelement), intent(in), target	refelem,
		integer, intent(in)	nlocalmeshperprc,
		integer, intent(in), optional	nproc,
		integer, intent(in), optional	myrank,
		real(rp), dimension(neg+1), intent(in), optional	fx )

Initialize an object to manage a 1D domain.

Parameters

[in]	neg	Total number of elements in the 1D domain
[in]	dom_xmin	Minimum coordinate of the 1D domain
[in]	dom_xmax	Maximum coordinate of the 1D domain
[in]	refelem	Reference element for the 1D mesh
[in]	nlocalmeshperprc	Number of local meshes managed by each process
[in]	nproc	Total number of processes (if not provided, it will be determined internally)
[in]	myrank	Rank of the current process (if not provided, it will be determined internally)
[in]	fx	Coordinates of the element boundaries (if not provided, they will be generated internally)

Definition at line 62 of file scale_mesh_linedom1d.F90.

    
    implicit none
 
    class(MeshLineDom1D), intent(inout) :: this
    integer, intent(in) :: NeG                       !< Total number of elements in the 1D domain
    real(RP), intent(in) :: dom_xmin                 !< Minimum coordinate of the 1D domain
    real(RP), intent(in) :: dom_xmax                 !< Maximum coordinate of the 1D domain
    type(LineElement), intent(in), target :: refElem !< Reference element for the 1D mesh
    integer, intent(in) :: NLocalMeshPerPrc          !< Number of local meshes managed by each process
    integer, intent(in), optional :: nproc           !< Total number of processes (if not provided, it will be determined internally)
    integer, intent(in), optional :: myrank          !< Rank of the current process (if not provided, it will be determined internally)
    real(RP), intent(in), optional :: FX(NeG+1)      !< Coordinates of the element boundaries (if not provided, they will be generated internally)
    !-----------------------------------------------------------------------------
 
    call meshbase1d_init(this, &
      neg,                                    &
      dom_xmin, dom_xmax,                     &
      refelem, nlocalmeshperprc,              &
      nproc, myrank, fx )
 
      this%dom_vol = (this%xmax_gl - this%xmin_gl)
 

References scale_mesh_base1d::meshbase1d_init().

Referenced by scale_mesh_linedom1d::meshlinedom1d::generate().

meshlinedom1d_final()

subroutine, public scale_mesh_linedom1d::meshlinedom1d_final

(

class(meshlinedom1d), intent(inout)

this

)

Finalize an object to manage a 1D domain.

Definition at line 92 of file scale_mesh_linedom1d.F90.

    implicit none
    class(MeshLineDom1D), intent(inout) :: this
    !-----------------------------------------------------------------------------
  
    call meshbase1d_final(this)
 

References scale_mesh_base1d::meshbase1d_final().

Referenced by scale_mesh_linedom1d::meshlinedom1d::generate().

meshlinedom1d_generate()

subroutine, public scale_mesh_linedom1d::meshlinedom1d_generate

(

class(meshlinedom1d), intent(inout), target

this

)

Generate meshes for the 1D domain.

Definition at line 103 of file scale_mesh_linedom1d.F90.

    implicit none
    class(MeshLineDom1D), intent(inout), target :: this
            
    integer :: n
    integer :: p
    type(LocalMesh1D), pointer :: mesh
 
    integer :: tileID_table(this%LOCAL_MESH_NUM, this%PRC_NUM)
    integer :: panelID_table(this%LOCAL_MESH_NUM*this%PRC_NUM)
    integer :: pi_table(this%LOCAL_MESH_NUM*this%PRC_NUM)
 
    !integer :: TILE_NUM_PER_PANEL
    integer :: tileID
    
    !-----------------------------------------------------------------------------
 
   ! TILE_NUM_PER_PANEL = this%LOCAL_MESH_NUM_global / 1
    
    
    !--- Construct the connectivity of patches  (only master node)
 
    call meshbase1d_assigndomid( this,    & ! (in)
        tileid_table, panelid_table,      & ! (out)
        pi_table )                          ! (out)
 
    !--- Setup local meshes managed by my process
 
    do n=1, this%LOCAL_MESH_NUM
      mesh => this%lcmesh_list(n)
      tileid = tileid_table(n, mesh%PRC_myrank+1)
      call meshbase1d_setuplocaldom( mesh,              & ! (inout)
         tileid,  panelid_table(tileid),                & ! (in)
         pi_table(tileid), this%Nprc,                   & ! (in)
         this%xmin_gl, this%xmax_gl,                    & ! (in)
         this%NeG / this%Nprc, this%FX(:) )               ! (in)
 
      !---
      ! write(*,*) "** my_rank=", mesh%PRC_myrank
      ! write(*,*) " tileID:", mesh%tileID
      ! write(*,*) " pnlID:", mesh%panelID, "-- i (within a panel)=", pi_table(tileID)
      ! write(*,*) " local mesh:", n, "( total", this%LOCAL_MESH_NUM, ")"
      ! write(*,*) " panel_connect:", this%tilePanelID_globalMap(:,mesh%tileID)
      ! write(*,*) " tile_connect:", this%tileID_globalMap(:,mesh%tileID)
      ! write(*,*) " face_connect:", this%tileFaceID_globalMap(:,mesh%tileID)
      ! write(*,*) " domain size"
      ! write(*,*) "   NeX:", mesh%Ne
      ! write(*,*) "   [X]:",  mesh%xmin, mesh%xmax
    end do
 
    this%isGenerated = .true.
 
    return

References scale_mesh_base1d::meshbase1d_assigndomid(), and scale_mesh_base1d::meshbase1d_setuplocaldom().

Referenced by scale_mesh_linedom1d::meshlinedom1d::generate().