CalculateEntropy_DGModel3D Subroutine

  1. SELF_DGModel3D.f90
  2. SELF_DGModel3D
  3. CalculateEntropy_DGModel3D

public subroutine CalculateEntropy_DGModel3D(this)

Arguments

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

Calls

proc~~calculateentropy_dgmodel3d~~CallsGraph proc~calculateentropy_dgmodel3d CalculateEntropy_DGModel3D interface~hipmemcpy hipMemcpy proc~calculateentropy_dgmodel3d->interface~hipmemcpy proc~gpucheck gpuCheck proc~calculateentropy_dgmodel3d->proc~gpucheck mpi_allreduce mpi_allreduce proc~calculateentropy_dgmodel3d->mpi_allreduce

Contents

Source Code

CalculateEntropy_DGModel3D

Source Code

  subroutine CalculateEntropy_DGModel3D(this)
    implicit none
    class(DGModel3D),intent(inout) :: this
    ! Local
    integer :: iel,i,j,k,ierror
    real(prec) :: e,jac
    real(prec) :: s(1:this%nvar)

    call gpuCheck(hipMemcpy(c_loc(this%solution%interior), &
                            this%solution%interior_gpu,sizeof(this%solution%interior), &
                            hipMemcpyDeviceToHost))

    e = 0.0_prec
    do iel = 1,this%geometry%nelem
      do k = 1,this%solution%interp%N+1
        do j = 1,this%solution%interp%N+1
          do i = 1,this%solution%interp%N+1
            jac = abs(this%geometry%J%interior(i,j,k,iel,1))
            s = this%solution%interior(i,j,k,iel,1:this%nvar)
            e = e+this%entropy_func(s)*jac
          enddo
        enddo
      enddo
    enddo

    if(this%mesh%decomp%mpiEnabled) then
      call mpi_allreduce(e, &
                         this%entropy, &
                         1, &
                         this%mesh%decomp%mpiPrec, &
                         MPI_SUM, &
                         this%mesh%decomp%mpiComm, &
                         iError)
    else
      this%entropy = e
    endif

  endsubroutine CalculateEntropy_DGModel3D