32 lambda, rhocp, bd, dt, n)
33 class(ax_t),
intent(in) :: Ax
34 type(mesh_t),
intent(inout) :: msh
35 type(space_t),
intent(inout) :: Xh
36 type(field_t),
intent(inout) :: s
37 type(field_t),
intent(inout) :: s_res
38 type(field_t),
intent(inout) :: f_Xh
39 type(coef_t),
intent(inout) :: c_Xh
40 real(kind=rp),
intent(in) :: lambda
41 real(kind=rp),
intent(in) :: rhocp
42 real(kind=rp),
intent(in) :: bd
43 real(kind=rp),
intent(in) :: dt
44 integer,
intent(in) :: n
48 c_xh%h1(i,1,1,1) = lambda
51 c_xh%h2(i,1,1,1) = rhocp * (bd / dt)
55 call ax%compute(s_res%x, s%x, c_xh, msh, xh)
58 s_res%x(i,1,1,1) = (-s_res%x(i,1,1,1)) + f_xh%x(i,1,1,1)
Residuals in the scalar equation (CPU version).
subroutine scalar_residual_cpu_compute(Ax, s, s_res, f_Xh, c_Xh, msh, Xh, lambda, rhocp, bd, dt, n)
Compute the residual.
Defines the residual for the scalar transport equation.
Abstract type to compute scalar residual.
Wrapper type for the routine to compute the scalar residual on the CPU.