875 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_qu()  argument
878   for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0;  in g1_qu()
893 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_vu()  argument
902         if (c == d) g1[(c * NcJ + d) * dim + e] += u[e];  in g1_vu()
917 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_conduct_qu()  argument
923   for (d = 0; d < dim; ++d) g1[d * dim + d] = p_th / u[uOff[TEMP]];  in g1_conduct_qu()
940 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_conduct_qT()  argument
946   for (d = 0; d < dim; ++d) g1[d] = -p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d];  in g1_conduct_qT()
1006 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_conduct_vu()  argument
1017         if (c == d) g1[(c * NcJ + d) * dim + e] += p_th / u[uOff[TEMP]] * u[uOff[VEL] + e];  in g1_conduct_vu()
1059 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_wT()  argument
1062   for (d = 0; d < dim; ++d) g1[d] = u[uOff[0] + d];  in g1_wT()
1098 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[])  in g1_conduct_wT()  argument
1105   for (d = 0; d < dim; ++d) g1[d] += c_p * p_th / u[uOff[TEMP]] * u[uOff[VEL] + d];  in g1_conduct_wT()