subroutine sourcemethod_DGModel3D(this)
implicit none
class(DGModel3D),intent(inout) :: this
! Local
integer :: i,j,k,iel
real(prec) :: s(1:this%nvar),dsdx(1:this%nvar,1:3)
call gpuCheck(hipMemcpy(c_loc(this%solution%interior), &
this%solution%interior_gpu,sizeof(this%solution%interior), &
hipMemcpyDeviceToHost))
call gpuCheck(hipMemcpy(c_loc(this%solutiongradient%interior), &
this%solutiongradient%interior_gpu,sizeof(this%solutiongradient%interior), &
hipMemcpyDeviceToHost))
do concurrent(i=1:this%solution%N+1,j=1:this%solution%N+1, &
k=1:this%solution%N+1,iel=1:this%mesh%nElem)
s = this%solution%interior(i,j,k,iel,1:this%nvar)
dsdx = this%solutionGradient%interior(i,j,k,iel,1:this%nvar,1:3)
this%source%interior(i,j,k,iel,1:this%nvar) = this%source3d(s,dsdx)
enddo
call gpuCheck(hipMemcpy(this%source%interior_gpu, &
c_loc(this%source%interior), &
sizeof(this%source%interior), &
hipMemcpyHostToDevice))
endsubroutine sourcemethod_DGModel3D