Skip to content

Commit

Permalink
Formatting
Browse files Browse the repository at this point in the history
  • Loading branch information
@fluidnumerics-joe
fluidnumerics-joe committed Nov 26, 2024
1 parent cd31b3f commit 0c3932e
Show file tree
Hide file tree
Showing 17 changed files with 3,063 additions and 3,063 deletions.
416 changes: 208 additions & 208 deletions src/SELF_DomainDecomposition_t.f90
View file

Large diffs are not rendered by default.

314 changes: 157 additions & 157 deletions src/SELF_Mesh_1D.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -26,186 +26,186 @@

module SELF_Mesh_1D

use SELF_Constants
use SELF_Lagrange
use SELF_Data
use SELF_Scalar_1D
use SELF_SupportRoutines
use SELF_HDF5
use SELF_Mesh
use SELF_Constants
use SELF_Lagrange
use SELF_Data
use SELF_Scalar_1D
use SELF_SupportRoutines
use SELF_HDF5
use SELF_Mesh

! External Libs !
use HDF5
! External Libs !
use HDF5

use iso_c_binding
use iso_c_binding

implicit none
implicit none

type, extends(SEMMesh) :: Mesh1D
integer, pointer, dimension(:, :) :: elemInfo
real(prec), pointer, dimension(:) :: nodeCoords
integer, pointer, dimension(:) :: globalNodeIDs
integer, pointer, dimension(:, :) :: BCType
character(LEN=255), allocatable :: BCNames(:)
integer, dimension(2) :: bcid = 0 ! Boundary conditions for the left and right endpoints
type,extends(SEMMesh) :: Mesh1D
integer,pointer,dimension(:,:) :: elemInfo
real(prec),pointer,dimension(:) :: nodeCoords
integer,pointer,dimension(:) :: globalNodeIDs
integer,pointer,dimension(:,:) :: BCType
character(LEN=255),allocatable :: BCNames(:)
integer,dimension(2) :: bcid = 0 ! Boundary conditions for the left and right endpoints

contains
procedure, public :: Init => Init_Mesh1D
procedure, public :: Free => Free_Mesh1D
generic, public :: StructuredMesh => UniformBlockMesh_Mesh1D
procedure, private :: UniformBlockMesh_Mesh1D
procedure, public :: ResetBoundaryConditionType => ResetBoundaryConditionType_Mesh1D
contains
procedure,public :: Init => Init_Mesh1D
procedure,public :: Free => Free_Mesh1D
generic,public :: StructuredMesh => UniformBlockMesh_Mesh1D
procedure,private :: UniformBlockMesh_Mesh1D
procedure,public :: ResetBoundaryConditionType => ResetBoundaryConditionType_Mesh1D

procedure, public :: Write_Mesh => Write_Mesh1D
procedure,public :: Write_Mesh => Write_Mesh1D

end type Mesh1D
endtype Mesh1D

contains

subroutine Init_Mesh1D(this, nElem, nNodes, nBCs)
implicit none
class(Mesh1D), intent(out) :: this
integer, intent(in) :: nElem
integer, intent(in) :: nNodes
integer, intent(in) :: nBCs

this%nGeo = 1
this%nElem = nElem
this%nGlobalElem = nElem
this%nNodes = nNodes
this%nCornerNodes = nElem*2
this%nUniqueNodes = 0
this%nBCs = nBCs
this%bcid = 0

allocate (this%elemInfo(1:4, 1:nElem))
allocate (this%nodeCoords(1:nNodes))
allocate (this%globalNodeIDs(1:nNodes))
allocate (this%BCType(1:4, 1:nBCs))

allocate (this%BCNames(1:nBCs))
call this%decomp%Init()

end subroutine Init_Mesh1D

subroutine Free_Mesh1D(this)
implicit none
class(Mesh1D), intent(inout) :: this

this%nElem = 0
this%nNodes = 0
this%nCornerNodes = 0
this%nUniqueNodes = 0
this%nBCs = 0
deallocate (this%elemInfo)
deallocate (this%nodeCoords)
deallocate (this%globalNodeIDs)
deallocate (this%BCType)
deallocate (this%BCNames)
call this%decomp%Free()

end subroutine Free_Mesh1D

subroutine UniformBlockMesh_Mesh1D(this, nElem, x)
implicit none
class(Mesh1D), intent(out) :: this
integer, intent(in) :: nElem
real(prec), intent(in) :: x(1:2)
! Local
integer :: iel, ngeo
integer :: nid, nNodes
integer :: i
real(prec) :: xU(1:nElem + 1)
type(Lagrange), target :: linearInterp
type(Lagrange), target :: nGeoInterp
type(Scalar1D) :: xLinear
type(Scalar1D) :: xGeo

ngeo = 1

nNodes = nElem*(nGeo + 1)
call this%Init(nElem, nNodes, 2)
this%quadrature = GAUSS_LOBATTO

! Set the hopr_nodeCoords
xU = UniformPoints(x(1), x(2), 1, nElem + 1)

call linearInterp%Init(1, GAUSS_LOBATTO, &
nGeo, GAUSS_LOBATTO)

call nGeoInterp%Init(nGeo, GAUSS_LOBATTO, &
nGeo, GAUSS_LOBATTO)

! Create a linear interpolant to interpolate to nGeo grid
call xLinear%Init(linearInterp, 1, nElem)
call xGeo%Init(nGeoInterp, 1, nElem)

do iel = 1, nElem
xLinear%interior(1:2, iel, 1) = xU(iel:iel + 1)
end do

call xLinear%GridInterp(xGeo%interior)

! Set the element information
nid = 1
do iel = 1, nElem
this%elemInfo(1, iel) = selfLineLinear ! Element Type
this%elemInfo(2, iel) = 1 ! Element Zone
this%elemInfo(3, iel) = nid ! Node Index Start
do i = 1, nGeo + 1
this%nodeCoords(nid) = xGeo%interior(i, iel, 1)
nid = nid + 1
end do
this%elemInfo(4, iel) = nid - 1 ! Node Index End
end do

call xLinear%Free()
call xGeo%Free()
call linearInterp%Free()
call nGeoInterp%Free()

end subroutine UniformBlockMesh_Mesh1D

subroutine ResetBoundaryConditionType_Mesh1D(this, leftbc, rightbc)
subroutine Init_Mesh1D(this,nElem,nNodes,nBCs)
implicit none
class(Mesh1D),intent(out) :: this
integer,intent(in) :: nElem
integer,intent(in) :: nNodes
integer,intent(in) :: nBCs

this%nGeo = 1
this%nElem = nElem
this%nGlobalElem = nElem
this%nNodes = nNodes
this%nCornerNodes = nElem*2
this%nUniqueNodes = 0
this%nBCs = nBCs
this%bcid = 0

allocate(this%elemInfo(1:4,1:nElem))
allocate(this%nodeCoords(1:nNodes))
allocate(this%globalNodeIDs(1:nNodes))
allocate(this%BCType(1:4,1:nBCs))

allocate(this%BCNames(1:nBCs))
call this%decomp%Init()

endsubroutine Init_Mesh1D

subroutine Free_Mesh1D(this)
implicit none
class(Mesh1D),intent(inout) :: this

this%nElem = 0
this%nNodes = 0
this%nCornerNodes = 0
this%nUniqueNodes = 0
this%nBCs = 0
deallocate(this%elemInfo)
deallocate(this%nodeCoords)
deallocate(this%globalNodeIDs)
deallocate(this%BCType)
deallocate(this%BCNames)
call this%decomp%Free()

endsubroutine Free_Mesh1D

subroutine UniformBlockMesh_Mesh1D(this,nElem,x)
implicit none
class(Mesh1D),intent(out) :: this
integer,intent(in) :: nElem
real(prec),intent(in) :: x(1:2)
! Local
integer :: iel,ngeo
integer :: nid,nNodes
integer :: i
real(prec) :: xU(1:nElem+1)
type(Lagrange),target :: linearInterp
type(Lagrange),target :: nGeoInterp
type(Scalar1D) :: xLinear
type(Scalar1D) :: xGeo

ngeo = 1

nNodes = nElem*(nGeo+1)
call this%Init(nElem,nNodes,2)
this%quadrature = GAUSS_LOBATTO

! Set the hopr_nodeCoords
xU = UniformPoints(x(1),x(2),1,nElem+1)

call linearInterp%Init(1,GAUSS_LOBATTO, &
nGeo,GAUSS_LOBATTO)

call nGeoInterp%Init(nGeo,GAUSS_LOBATTO, &
nGeo,GAUSS_LOBATTO)

! Create a linear interpolant to interpolate to nGeo grid
call xLinear%Init(linearInterp,1,nElem)
call xGeo%Init(nGeoInterp,1,nElem)

do iel = 1,nElem
xLinear%interior(1:2,iel,1) = xU(iel:iel+1)
enddo

call xLinear%GridInterp(xGeo%interior)

! Set the element information
nid = 1
do iel = 1,nElem
this%elemInfo(1,iel) = selfLineLinear ! Element Type
this%elemInfo(2,iel) = 1 ! Element Zone
this%elemInfo(3,iel) = nid ! Node Index Start
do i = 1,nGeo+1
this%nodeCoords(nid) = xGeo%interior(i,iel,1)
nid = nid+1
enddo
this%elemInfo(4,iel) = nid-1 ! Node Index End
enddo

call xLinear%Free()
call xGeo%Free()
call linearInterp%Free()
call nGeoInterp%Free()

endsubroutine UniformBlockMesh_Mesh1D

subroutine ResetBoundaryConditionType_Mesh1D(this,leftbc,rightbc)
!! This method can be used to reset all of the boundary elements
!! boundary condition type to the desired value.
!!
!! Note that ALL physical boundaries will be set to have this boundary
!! condition
implicit none
class(Mesh1D), intent(inout) :: this
integer, intent(in) ::leftbc, rightbc
implicit none
class(Mesh1D),intent(inout) :: this
integer,intent(in) ::leftbc,rightbc

this%bcid(1) = leftbc
this%bcid(2) = rightbc
this%bcid(1) = leftbc
this%bcid(2) = rightbc

end subroutine ResetBoundaryConditionType_Mesh1D
endsubroutine ResetBoundaryConditionType_Mesh1D

subroutine Write_Mesh1D(this, meshFile)
! Writes mesh output in HOPR format (serial IO only)
implicit none
class(Mesh1D), intent(inout) :: this
character(*), intent(in) :: meshFile
! Local
integer(HID_T) :: fileId
subroutine Write_Mesh1D(this,meshFile)
! Writes mesh output in HOPR format (serial IO only)
implicit none
class(Mesh1D),intent(inout) :: this
character(*),intent(in) :: meshFile
! Local
integer(HID_T) :: fileId

call Open_HDF5(meshFile, H5F_ACC_RDWR_F, fileId)
call Open_HDF5(meshFile,H5F_ACC_RDWR_F,fileId)

call WriteAttribute_HDF5(fileId, 'nElems', this%nElem)
call WriteAttribute_HDF5(fileId, 'Ngeo', this%nGeo)
call WriteAttribute_HDF5(fileId, 'nBCs', this%nBCs)
call WriteAttribute_HDF5(fileId,'nElems',this%nElem)
call WriteAttribute_HDF5(fileId,'Ngeo',this%nGeo)
call WriteAttribute_HDF5(fileId,'nBCs',this%nBCs)

call WriteArray_HDF5(fileId, 'BCType', this%bcType)
call WriteArray_HDF5(fileId,'BCType',this%bcType)

! Read local subarray of ElemInfo
call WriteArray_HDF5(fileId, 'ElemInfo', this%elemInfo)
! Read local subarray of ElemInfo
call WriteArray_HDF5(fileId,'ElemInfo',this%elemInfo)

! Read local subarray of NodeCoords and GlobalNodeIDs
call WriteArray_HDF5(fileId, 'NodeCoords', this%nodeCoords)
call WriteArray_HDF5(fileId, 'GlobalNodeIDs', this%globalNodeIDs)
! Read local subarray of NodeCoords and GlobalNodeIDs
call WriteArray_HDF5(fileId,'NodeCoords',this%nodeCoords)
call WriteArray_HDF5(fileId,'GlobalNodeIDs',this%globalNodeIDs)

call Close_HDF5(fileID)
call Close_HDF5(fileID)

end subroutine Write_Mesh1D
endsubroutine Write_Mesh1D

end module SELF_Mesh_1D
endmodule SELF_Mesh_1D
Loading

0 comments on commit 0c3932e

Please sign in to comment.