13 integer,
intent(in) :: nl, nelv, ldim
14 real(kind=
rp),
intent(inout) :: e(nl**ldim, nelv)
15 real(kind=
rp),
intent(inout) :: r(nl**ldim, nelv)
16 real(kind=
rp),
intent(inout) :: s(nl*nl,2,ldim, nelv)
17 real(kind=
rp),
intent(inout) :: d(nl**ldim, nelv)
21 if(.not. ldim .eq. 3)
then
23 call tnsr2d_el_cpu(e(1,ie), nl, r(1,ie), nl, s(1,2,1,ie), s(1,1,2,ie))
25 r(i,ie) = d(i,ie) * e(i,ie)
27 call tnsr2d_el_cpu(e(1,ie), nl, r(1,ie), nl, s(1,1,1,ie), s(1,2,2,ie))
32 s(1,2,1,ie), s(1,1,2,ie), s(1,1,3,ie))
34 r(i,ie) = d(i,ie) * e(i,ie)
37 s(1,1,1,ie), s(1,2,2,ie), s(1,2,3,ie))
subroutine, public fdm_do_fast_cpu(e, r, s, d, nl, ldim, nelv)
integer, parameter, public rp
Global precision used in computations.
subroutine, public tnsr2d_el_cpu(v, nv, u, nu, A, Bt)
subroutine, public tnsr3d_el_cpu(v, nv, u, nu, A, Bt, Ct)