subroutine sourcemethod_DGModel1D(this)
implicit none
class(DGModel1D),intent(inout) :: this
! Local
integer :: iel
integer :: i
real(prec) :: f(1:this%solution%nvar),dfdx(1:this%solution%nvar)
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 iel = 1,this%mesh%nelem
do i = 1,this%solution%interp%N+1
f = this%solution%interior(i,iel,1:this%solution%nvar)
dfdx = this%solutionGradient%interior(i,iel,1:this%solution%nvar)
this%source%interior(i,iel,1:this%solution%nvar) = &
this%source1d(f,dfdx)
enddo
enddo
call gpuCheck(hipMemcpy(this%source%interior_gpu, &
c_loc(this%source%interior), &
sizeof(this%source%interior), &
hipMemcpyHostToDevice))
endsubroutine sourcemethod_DGModel1D