146 static void RDDiffusionCoefficient(RD rd, PetscBool limit, RDNode *n, RDNode *nx, PetscScalar *D_R,…  in RDDiffusionCoefficient()  argument
157   *D_R      = rd->c / denom;  in RDDiffusionCoefficient()
206   PetscScalar D_L, D_R, fluxL, fluxR;  in RDDiffusion()  local
223   RDDiffusionCoefficient(rd, PETSC_TRUE, &n_R, &nx_R, &D_R, &dD_R, &dxD_R);  in RDDiffusion()
224   fluxR       = D_R * nx_R.E;  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()
348       PetscScalar D_R, bcTheta = rd->bcmidpoint ? Theta : 1.;  in RDIFunction_FD()  local
357         RDDiffusionCoefficient(rd, rd->bclimit, &n, &nx, &D_R, 0, 0);  in RDIFunction_FD()
358         f[0].E = hx * (n.E - 2. * D_R * nx.E - rd->Eapplied);  in RDIFunction_FD()
651       PetscScalar Em_t, rad, D_R, D0_R;  in RDIFunction_FE()  local
669       RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0);  in RDIFunction_FE()
671 …f[i + j].E += (deriv[q][j] * weight[q] * ((1. - Theta) * D0_R * n0x.E + Theta * D_R * nx.E) + inte…  in RDIFunction_FE()
679       PetscScalar D_R, D_R_bc;  in RDIFunction_FE()  local
687       RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0);  in RDIFunction_FE()
689       ratio = PetscRealPart(D_R / D_R_bc);  in RDIFunction_FE()
736       PetscScalar              D_R;  in RDIJacobian_FE()  local
742       RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, &dD_R, &dxD_R);  in RDIJacobian_FE()
746 …j] * weight[q] * (a - drad.E) * interp[q][k] + deriv[q][j] * weight[q] * ((D_R + dxD_R.E * nx.E) *…  in RDIJacobian_FE()
758       PetscScalar D_R, D_R_bc;  in RDIJacobian_FE()  local
766       RDDiffusionCoefficient(rd, PETSC_TRUE, &n, &nx, &D_R, 0, 0);  in RDIJacobian_FE()
768       ratio = PetscRealPart(D_R / D_R_bc);  in RDIJacobian_FE()