Free_MappedScalar2D Subroutine

  1. SELF_MappedScalar_2D.f90
  2. SELF_MappedScalar_2D
  3. Free_MappedScalar2D

public subroutine Free_MappedScalar2D(this)

Arguments

TypeIntentOptionalAttributesName
class(MappedScalar2D), intent(inout) :: this

Calls

proc~~free_mappedscalar2d~~CallsGraph proc~free_mappedscalar2d Free_MappedScalar2D proc~free_scalar2d_t Free_Scalar2D_t proc~free_mappedscalar2d->proc~free_scalar2d_t proc~gpucheck gpuCheck proc~free_mappedscalar2d->proc~gpucheck interface~hipfree hipFree proc~free_mappedscalar2d->interface~hipfree mpi_request_free mpi_request_free proc~free_mappedscalar2d->mpi_request_free mpi_finalized mpi_finalized proc~free_mappedscalar2d->mpi_finalized

Contents

Source Code

Free_MappedScalar2D

Source Code

  subroutine Free_MappedScalar2D(this)
    implicit none
    class(MappedScalar2D),intent(inout) :: this
    ! Local
    integer :: n,iError
    logical :: mpiIsFinalized

    ! Host storage is owned by the pools in the parent type (Stage 6b), so the parent Free
    ! releases it rather than deallocating these pointers.
    call Free_Scalar2D_t(this)

    call gpuCheck(hipFree(this%interior_gpu))
    call gpuCheck(hipFree(this%boundary_gpu))
    call gpuCheck(hipFree(this%extBoundary_gpu))
    call gpuCheck(hipFree(this%avgBoundary_gpu))
    ! boundarynormal_gpu is allocated in Init alongside the four above and must be released
    ! here with them. Omitting it leaks (N+1)*4*nElem*2*nvar reals on every Free/Init cycle,
    ! which the adaptive loop performs once per epoch for each of the model's five
    ! MappedScalar2D fields. Compare Free_Scalar2D in src/gpu/SELF_Scalar_2D.f90 and
    ! Free_Scalar3D in src/gpu/SELF_Scalar_3D.f90, which both free it.
    call gpuCheck(hipFree(this%boundarynormal_gpu))
    call gpuCheck(hipFree(this%jas_gpu))
    if(c_associated(this%mortarBuff_gpu)) then
      call gpuCheck(hipFree(this%mortarBuff_gpu))
      this%mortarBuff_gpu = c_null_ptr
    endif

    if(c_associated(this%halo_sendbuf_gpu)) call gpuCheck(hipFree(this%halo_sendbuf_gpu))
    if(c_associated(this%halo_recvbuf_gpu)) call gpuCheck(hipFree(this%halo_recvbuf_gpu))
    this%halo_sendbuf_gpu = c_null_ptr
    this%halo_recvbuf_gpu = c_null_ptr

    if(allocated(this%halo_reqs)) then
      ! Persistent requests can only be released while MPI is still
      ! initialized; if the mesh (and its MPI finalization) was freed first,
      ! MPI has reclaimed them already.
      call MPI_FINALIZED(mpiIsFinalized,iError)
      if(.not. mpiIsFinalized) then
        do n = 1,size(this%halo_reqs)
          call MPI_REQUEST_FREE(this%halo_reqs(n),iError)
        enddo
      endif
      deallocate(this%halo_reqs)
    endif
    this%halo_nactive = 0
    this%halo_inflight = 0
    this%halo_static_done = .false.

  endsubroutine Free_MappedScalar2D