| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| class(ECDGModel3D), | intent(inout) | :: | this |
subroutine CalculateTendency_ECDGModel3D(this)
implicit none
class(ECDGModel3D),intent(inout) :: this
! Local
integer :: ndof
call this%solution%BoundaryInterp()
call this%solution%SideExchange(this%mesh)
call this%PreTendency()
call this%SetBoundaryCondition()
if(this%gradient_enabled) then
call this%CalculateSolutionGradient()
call this%SetGradientBoundaryCondition()
call this%solutionGradient%AverageSides()
endif
call this%SourceMethod()
call this%BoundaryFlux()
call this%TwoPointFluxMethod()
call this%twoPointFlux%MappedDivergence(this%fluxDivergence%interior_gpu)
call ECDGSurfaceContribution_3D_gpu( &
this%flux%boundarynormal_gpu, &
this%geometry%J%interior_gpu, &
this%solution%interp%bMatrix_gpu, &
this%solution%interp%qWeights_gpu, &
this%fluxDivergence%interior_gpu, &
this%solution%interp%N,this%solution%nVar,this%mesh%nElem)
ndof = this%solution%nVar* &
this%mesh%nElem* &
(this%solution%interp%N+1)* &
(this%solution%interp%N+1)* &
(this%solution%interp%N+1)
call CalculateDSDt_gpu(this%fluxDivergence%interior_gpu, &
this%source%interior_gpu, &
this%dSdt%interior_gpu,ndof)
endsubroutine CalculateTendency_ECDGModel3D