58 const PetscInt VEL  = 0;  variable
740   for (d = 0; d < dim; ++d) f0[0] += p_th / u[uOff[TEMP]] * u_x[uOff_x[VEL] + d * dim + d];  in f0_conduct_q()
743 …for (d = 0; d < dim; ++d) f0[0] -= p_th / (u[uOff[TEMP]] * u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x…  in f0_conduct_q()
757     f0[c] += u_t[uOff[VEL] + c];  in f0_v()
759     for (d = 0; d < dim; ++d) f0[c] += u[uOff[VEL] + d] * u_x[uOff_x[VEL] + c * dim + d];  in f0_v()
775   for (d = 0; d < dim; ++d) f0[d] = rho * S * u_t[uOff[VEL] + d];  in f0_conduct_v()
779     for (d = 0; d < dim; ++d) f0[c] += rho * u[uOff[VEL] + d] * u_x[uOff_x[VEL] + c * dim + d];  in f0_conduct_v()
814   for (c = 0; c < Nc; ++c) u_div += PetscRealPart(u_x[uOff_x[VEL] + c * dim + c]);  in f1_conduct_v()
818 … = 0; d < dim; ++d) f1[c * dim + d] += coef * (u_x[uOff_x[VEL] + c * dim + d] + u_x[uOff_x[VEL] + …  in f1_conduct_v()
835   for (d = 0; d < dim; ++d) f0[0] += u[uOff[VEL] + d] * u_x[uOff_x[TEMP] + d];  in f0_w()
850 …for (d = 0; d < dim; ++d) f0[0] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d] * u_x[uOff_x[TEMP…  in f0_conduct_w()
937 …+= p_th / PetscSqr(u[uOff[TEMP]]) * (-u_x[uOff_x[VEL] + d * dim + d] + 2.0 / u[uOff[TEMP]] * u[uOf…  in g0_conduct_qT()
946   for (d = 0; d < dim; ++d) g1[d] = -p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d];  in g1_conduct_qT()
979   for (d = 0; d < dim; ++d) g0[d] -= p_th * S / PetscSqr(u[uOff[TEMP]]) * u_t[uOff[VEL] + d];  in g0_conduct_vT()
983 …d < dim; ++d) g0[c] -= p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x[uOff_x[VEL] + c * d…  in g0_conduct_vT()
1002 …  for (d = 0; d < dim; ++d) g0[c * Nc + d] += p_th / u[uOff[TEMP]] * u_x[uOff_x[VEL] + c * Nc + d];  in g0_conduct_vu()
1017         if (c == d) g1[(c * NcJ + d) * dim + e] += p_th / u[uOff[TEMP]] * u[uOff[VEL] + e];  in g1_conduct_vu()
1095 …for (d = 0; d < dim; ++d) g0[0] -= c_p * p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d] * u_x[u…  in g0_conduct_wT()
1105   for (d = 0; d < dim; ++d) g1[d] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d];  in g1_conduct_wT()
1189 …, "top wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in UniformBoundaryConditions()
1191 …bottom wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in UniformBoundaryConditions()
1193 …"right wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in UniformBoundaryConditions()
1195 … "left wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in UniformBoundaryConditions()
1224     PetscCall(PetscDSSetResidual(ds, VEL, f0_v, f1_v));  in SetupProblem()
1228     PetscCall(PetscDSSetJacobian(ds, VEL, VEL, g0_vu, g1_vu, NULL, g3_vu));  in SetupProblem()
1229     PetscCall(PetscDSSetJacobian(ds, VEL, PRES, NULL, NULL, g2_vp, NULL));  in SetupProblem()
1230     PetscCall(PetscDSSetJacobian(ds, PRES, VEL, NULL, g1_qu, NULL, NULL));  in SetupProblem()
1231     PetscCall(PetscDSSetJacobian(ds, TEMP, VEL, g0_wu, NULL, NULL, NULL));  in SetupProblem()
1236       PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_quadratic_v, 0, NULL));  in SetupProblem()
1239       exactFuncs[VEL]    = quadratic_u;  in SetupProblem()
1242       exactFuncs_t[VEL]  = quadratic_u_t;  in SetupProblem()
1249       PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_cubic_v, 0, NULL));  in SetupProblem()
1252       exactFuncs[VEL]    = cubic_u;  in SetupProblem()
1255       exactFuncs_t[VEL]  = cubic_u_t;  in SetupProblem()
1262       PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_cubic_trig_v, 0, NULL));  in SetupProblem()
1265       exactFuncs[VEL]    = cubic_trig_u;  in SetupProblem()
1268       exactFuncs_t[VEL]  = cubic_trig_u_t;  in SetupProblem()
1277       exactFuncs[VEL]    = taylor_green_u;  in SetupProblem()
1280       exactFuncs_t[VEL]  = taylor_green_u_t;  in SetupProblem()
1291     PetscCall(PetscDSSetResidual(ds, VEL, f0_conduct_v, f1_conduct_v));  in SetupProblem()
1295     PetscCall(PetscDSSetJacobian(ds, VEL, VEL, g0_conduct_vu, g1_conduct_vu, NULL, g3_conduct_vu));  in SetupProblem()
1296     PetscCall(PetscDSSetJacobian(ds, VEL, PRES, NULL, NULL, g2_conduct_vp, NULL));  in SetupProblem()
1297     PetscCall(PetscDSSetJacobian(ds, VEL, TEMP, g0_conduct_vT, NULL, NULL, NULL));  in SetupProblem()
1298     PetscCall(PetscDSSetJacobian(ds, PRES, VEL, g0_conduct_qu, g1_conduct_qu, NULL, NULL));  in SetupProblem()
1300     PetscCall(PetscDSSetJacobian(ds, TEMP, VEL, g0_conduct_wu, NULL, NULL, NULL));  in SetupProblem()
1305 … PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_conduct_quadratic_v, 0, NULL));  in SetupProblem()
1309       exactFuncs[VEL]    = quadratic_u;  in SetupProblem()
1312       exactFuncs_t[VEL]  = quadratic_u_t;  in SetupProblem()
1320       PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_conduct_pipe_v, 0, NULL));  in SetupProblem()
1324       exactFuncs[VEL]    = pipe_u;  in SetupProblem()
1327       exactFuncs_t[VEL]  = pipe_u_t;  in SetupProblem()
1335 …PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, VEL, 0, 0, f0_conduct_bd_pipe_v, 0, NULL)…  in SetupProblem()
1339 …PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, VEL, 0, 0, f0_conduct_bd_pipe_v, 0, NULL)…  in SetupProblem()
1343 …, "top wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in SetupProblem()
1346 …bottom wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in SetupProblem()
1351 …PetscCall(PetscWeakFormSetIndexResidual(wf, NULL, 0, VEL, 0, 1, f0_conduct_pipe_wiggly_v, 0, NULL)…  in SetupProblem()
1355       exactFuncs[VEL]    = pipe_wiggly_u;  in SetupProblem()
1358       exactFuncs_t[VEL]  = pipe_wiggly_u_t;  in SetupProblem()
1366 …PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, VEL, 0, 0, f0_conduct_bd_pipe_wiggly_v, 0…  in SetupProblem()
1370 …PetscCall(PetscWeakFormSetIndexBdResidual(wf, label, id, VEL, 0, 0, f0_conduct_bd_pipe_wiggly_v, 0…  in SetupProblem()
1374 …, "top wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in SetupProblem()
1377 …bottom wall velocity", label, 1, &id, VEL, 0, NULL, (PetscVoidFn *)exactFuncs[VEL], (PetscVoidFn *…  in SetupProblem()
1411   PetscCall(PetscDSSetExactSolution(ds, VEL, exactFuncs[VEL], ctx));  in SetupProblem()
1414   PetscCall(PetscDSSetExactSolutionTimeDerivative(ds, VEL, exactFuncs_t[VEL], ctx));  in SetupProblem()
1433   PetscCall(DMGetField(dm, VEL, NULL, (PetscObject *)&fe));  in CreateCellDM()
1485   PetscCall(DMSetField(dm, VEL, NULL, (PetscObject)fe[VEL]));  in SetupDiscretization()
1495   PetscCall(PetscFEDestroy(&fe[VEL]));  in SetupDiscretization()