42 lambda, rhocp, bd, dt, n)
43 class(
ax_t),
intent(in) :: Ax
44 type(
mesh_t),
intent(inout) :: msh
45 type(
space_t),
intent(inout) :: Xh
46 type(
field_t),
intent(inout) :: s
47 type(
field_t),
intent(inout) :: s_res
48 type(
field_t),
intent(inout) :: f_Xh
49 type(
coef_t),
intent(inout) :: c_Xh
50 type(
field_t),
intent(in) :: lambda
51 real(kind=
rp),
intent(in) :: rhocp
52 real(kind=
rp),
intent(in) :: bd
53 real(kind=
rp),
intent(in) :: dt
54 integer,
intent(in) :: n
57 call copy(c_xh%h1, lambda%x, n)
58 call cfill(c_xh%h2, rhocp * bd / dt, n)
62 call ax%compute(s_res%x, s%x, c_xh, msh, xh)
65 s_res%x(i,1,1,1) = (-s_res%x(i,1,1,1)) + f_xh%x(i,1,1,1)
Defines a Matrix-vector product.
subroutine, public cfill(a, c, n)
Set all elements to a constant c .
subroutine, public copy(a, b, n)
Copy a vector .
integer, parameter, public rp
Global precision used in computations.
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.
Defines a function space.
Base type for a matrix-vector product providing .
Coefficients defined on a given (mesh, ) tuple. Arrays use indices (i,j,k,e): element e,...
Abstract type to compute scalar residual.
Wrapper type for the routine to compute the scalar residual on the CPU.
The function space for the SEM solution fields.