scale_mesh_linedom1d.F90

Go to the documentation of this file.

!-------------------------------------------------------------------------------
!> module FElib / Mesh / 1D domain
!!
!! @par Description
!!      Manage mesh data of 1D domain for element-based methods
!!
!! @author Yuta Kawai, Team SCALE
!<
!-------------------------------------------------------------------------------
#include "scaleFElib.h"
module scale_mesh_linedom1d
 
  !-----------------------------------------------------------------------------
  !
  !++ used modules
  !
  use scale_precision
 
  use scale_mesh_base1d, only: &
    meshbase1d, meshbase1d_init, meshbase1d_final, &
    meshbase1d_assigndomid, meshbase1d_setgeometricinfo, meshbase1d_setuplocaldom
 
  use scale_localmesh_1d, only: localmesh1d
  use scale_element_base, only: elementbase1d
  use scale_element_line, only: lineelement
 
  !-----------------------------------------------------------------------------
  implicit none
  private
 
  !-----------------------------------------------------------------------------
  !
  !++ Public type & procedure
  ! 
  type, extends(meshbase1d), public :: meshlinedom1d
  contains
    procedure :: init => meshlinedom1d_init
    procedure :: final => meshlinedom1d_final
    procedure :: generate => meshlinedom1d_generate
  end type meshlinedom1d
 
  public :: meshlinedom1d_init, meshlinedom1d_final
  public :: meshlinedom1d_generate
  
  !-----------------------------------------------------------------------------
  !
  !++ Public parameters & variables
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private procedure
  !
  
  !-----------------------------------------------------------------------------
  !
  !++ Private parameters & variables
  !
 
contains
  !> Initialize an object to manage a 1D domain
  subroutine meshlinedom1d_init(this,       & ! (inout)
    neg,                                    & ! (in)
    dom_xmin, dom_xmax,                     & ! (in)
    refelem, nlocalmeshperprc,              & ! (in)
    nproc, myrank, fx )                       ! (in)
    
    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)
 
  end subroutine meshlinedom1d_init
 
  !> Finalize an object to manage a 1D domain
  subroutine meshlinedom1d_final( this ) ! (inout)
    implicit none
    class(meshlinedom1d), intent(inout) :: this
    !-----------------------------------------------------------------------------
  
    call meshbase1d_final(this)
 
  end subroutine meshlinedom1d_final
  
  !> Generate meshes for the 1D domain
!OCL SERIAL
  subroutine meshlinedom1d_generate( this ) ! (inout)
    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
  end subroutine meshlinedom1d_generate
  
  !-------------------------------------------------
end module scale_mesh_linedom1d