146 …cient(RD rd, PetscBool limit, RDNode *n, RDNode *nx, PetscScalar *D_R, RDNode *dD_R, RDNode *dxD_R)  in RDDiffusionCoefficient()  argument
158   if (dD_R) {  in RDDiffusionCoefficient()
159     dD_R->E = -rd->c / PetscSqr(denom) * ddenom.E;  in RDDiffusionCoefficient()
160     dD_R->T = -rd->c / PetscSqr(denom) * ddenom.T;  in RDDiffusionCoefficient()
205   RDNode      n_L, nx_L, n_R, nx_R, dD_L, dxD_L, dD_R, dxD_R, dfluxL[2], dfluxR[2];  in RDDiffusion()  local
223   RDDiffusionCoefficient(rd, PETSC_TRUE, &n_R, &nx_R, &D_R, &dD_R, &dxD_R);  in RDDiffusion()
225   dfluxR[0].E = -ihx * D_R + (0.5 * dD_R.E - ihx * dxD_R.E) * nx_R.E;  in RDDiffusion()
226   dfluxR[1].E = +ihx * D_R + (0.5 * dD_R.E + ihx * dxD_R.E) * nx_R.E;  in RDDiffusion()
227   dfluxR[0].T = (0.5 * dD_R.T - ihx * dxD_R.T) * nx_R.E;  in RDDiffusion()
228   dfluxR[1].T = (0.5 * dD_R.T + ihx * dxD_R.T) * nx_R.E;  in RDDiffusion()
738       RDNode                   n, nx, nt, ntx, drad, dD_R, dxD_R, dEm;  in RDIJacobian_FE()  local
742       RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, &dD_R, &dxD_R);  in RDIJacobian_FE()
746 …q][k] + deriv[q][j] * weight[q] * ((D_R + dxD_R.E * nx.E) * deriv[q][k] + dD_R.E * nx.E * interp[q…  in RDIJacobian_FE()
747 …ad.T * interp[q][k]) + deriv[q][j] * weight[q] * (dxD_R.T * deriv[q][k] + dD_R.T * interp[q][k]) *…  in RDIJacobian_FE()