162   const PetscScalar  eqss[]   = {0.5, 1., 0.5}; /* equations rescaling for a symmetric Jacobian */  in C_0()  local
164 …for (PetscInt k = 0; k < 3; k++) f0[k] = eqss[k] * (u_t[uOff[C_FIELD_ID] + k] - outerp[k] + alpha …  in C_0()
185   const PetscScalar  eqss[]   = {0.5, 1., 1., 0.5, 1., 0.5};  in C_0_3d()  local
187 …for (PetscInt k = 0; k < 6; k++) f0[k] = eqss[k] * (u_t[uOff[C_FIELD_ID] + k] - outerp[k] + alpha …  in C_0_3d()
205   const PetscScalar eqss[] = {0.5, 1., 0.5};  in JC_0_c0c0()  local
209 …for (PetscInt j = 0; j < 3; j++) J[k * 3 + j] = eqss[k] * (alpha * dfnorm * dC[j] * u[uOff[C_FIELD…  in JC_0_c0c0()
211     J[k * 3 + k] += eqss[k] * (alpha * fnorm + u_tShift);  in JC_0_c0c0()
232   const PetscScalar eqss[] = {0.5, 1., 1., 0.5, 1., 0.5};  in JC_0_c0c0_3d()  local
236 …for (PetscInt j = 0; j < 6; j++) J[k * 6 + j] = eqss[k] * (alpha * dfnorm * dC[j] * u[uOff[C_FIELD…  in JC_0_c0c0_3d()
238     J[k * 6 + k] += eqss[k] * (alpha * fnorm + u_tShift);  in JC_0_c0c0_3d()
246   const PetscScalar  eqss[] = {0.5, 1., 0.5};  in JC_0_c0p1()  local
248   J[0] = -2 * gradp[0] * eqss[0];  in JC_0_c0p1()
251   J[2] = -gradp[1] * eqss[1];  in JC_0_c0p1()
252   J[3] = -gradp[0] * eqss[1];  in JC_0_c0p1()
255   J[5] = -2 * gradp[1] * eqss[2];  in JC_0_c0p1()
261   const PetscScalar  eqss[] = {0.5, 1., 1., 0.5, 1., 0.5};  in JC_0_c0p1_3d()  local
263   J[0] = -2 * gradp[0] * eqss[0];  in JC_0_c0p1_3d()
267   J[3] = -gradp[1] * eqss[1];  in JC_0_c0p1_3d()
268   J[4] = -gradp[0] * eqss[1];  in JC_0_c0p1_3d()
271   J[6] = -gradp[2] * eqss[2];  in JC_0_c0p1_3d()
273   J[8] = -gradp[0] * eqss[2];  in JC_0_c0p1_3d()
276   J[10] = -2 * gradp[1] * eqss[3];  in JC_0_c0p1_3d()
280   J[13] = -gradp[2] * eqss[4];  in JC_0_c0p1_3d()
281   J[14] = -gradp[1] * eqss[4];  in JC_0_c0p1_3d()
285   J[17] = -2 * gradp[2] * eqss[5];  in JC_0_c0p1_3d()