85     CeedScalar Tau_d[3], stab[5][3], body_force[5] = {0.}, U_dot[5];  in MassFunction_Newtonian()  local
86     UnpackState_U(s_dot.U, U_dot);  in MassFunction_Newtonian()
88     Stabilization(context, s, Tau_d, grad_s, U_dot, body_force, stab);  in MassFunction_Newtonian()
200     CeedScalar Tau_d[3], stab[5][3], U_dot[5] = {0};  in RHSFunction_Newtonian()  local
202     Stabilization(context, s, Tau_d, grad_s, U_dot, body_force, stab);  in RHSFunction_Newtonian()
265     CeedScalar Tau_d[3], stab[5][3], U_dot[5] = {0}, qi_dot[5];  in IFunction_Newtonian()  local
268     UnpackState_U(s_dot.U, U_dot);  in IFunction_Newtonian()
270     for (CeedInt j = 0; j < 5; j++) v[j][i] = wdetJ * (U_dot[j] - body_force[j]);  in IFunction_Newtonian()
280     Stabilization(context, s, Tau_d, grad_s, U_dot, body_force, stab);  in IFunction_Newtonian()
369     CeedScalar dstab[5][3], U_dot[5] = {0};  in IJacobian_Newtonian()  local
370     for (CeedInt j = 0; j < 5; j++) U_dot[j] = context->ijacobian_time_shift * dU[j];  in IJacobian_Newtonian()
371     Stabilization(context, s, Tau_d, grad_ds, U_dot, dbody_force, dstab);  in IJacobian_Newtonian()