Y direction for 5 scalars projection of the physical flux onto the testfunctions, with unrolled loops
do itest1=2,maxpolydegree+1 prescal(itest1-1) = ply_scalprodduallegdiff(itest1-1, itest1) & & * jacobidetstiffproj end do
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | maxpolydegree |
the maximal polynomial degree in each spatial direction. |
||
| real(kind=rk), | intent(in) | :: | length |
the length of the cubes. |
||
| real(kind=rk), | intent(inout) | :: | state(:,:,:) |
the state to alter. |
||
| real(kind=rk), | intent(in) | :: | dl_prod(2,maxpolydegree+1) |
precomputed dual legendre products: |
||
| integer, | intent(in) | :: | ielem | |||
| real(kind=rk), | intent(in) | :: | state_der(:,:) |
the state to be used to project the physical fluxes |