955 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])  in g3_vu()  argument
963       g3[((c * Nc + c) * dim + d) * dim + d] += nu;  in g3_vu()
964       g3[((c * Nc + d) * dim + d) * dim + c] += nu;  in g3_vu()
1023 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])  in g3_conduct_vu()  argument
1033       g3[((c * Nc + c) * dim + d) * dim + d] += mu / Re; // gradU  in g3_conduct_vu()
1035       g3[((c * Nc + d) * dim + d) * dim + c] += mu / Re; // gradU transpose  in g3_conduct_vu()
1037       g3[((c * Nc + d) * dim + c) * dim + d] -= 2.0 / 3.0 * mu / Re;  in g3_conduct_vu()
1065 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])  in g3_wT()  argument
1070   for (d = 0; d < dim; ++d) g3[d * dim + d] = alpha;  in g3_wT()
1108 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[])  in g3_conduct_wT()  argument
1115   for (d = 0; d < dim; ++d) g3[d * dim + d] = k / Pe;  in g3_conduct_wT()