241 static PetscErrorCode RDGetLocalArrays(RD rd, TS ts, Vec X, Vec Xdot, PetscReal *Theta, PetscReal *…  in RDGetLocalArrays()  argument
261   if (istheta && rd->endpoint) PetscCall(TSThetaGetTheta(ts, Theta));  in RDGetLocalArrays()
262   else *Theta = 1.;  in RDGetLocalArrays()
265 …PetscCall(VecWAXPY(*X0loc, -(*Theta) * (*dt), *Xloc_t, *Xloc));     /* back out the value at the s…  in RDGetLocalArrays()
316   PetscReal     hx, Theta, dt;  in RDIFunction_FD()  local
321   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FD()
332     rad = (1. - Theta) * RDRadiation(rd, &x0[i], 0) + Theta * RDRadiation(rd, &x[i], 0);  in RDIFunction_FD()
348       PetscScalar D_R, bcTheta = rd->bcmidpoint ? Theta : 1.;  in RDIFunction_FD()
369 …PetscScalar diff = (1. - Theta) * RDDiffusion(rd, hx, x0, i, 0) + Theta * RDDiffusion(rd, hx, x, i…  in RDIFunction_FD()
384   PetscReal     hx, Theta, dt;  in RDIJacobian_FD()  local
389   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FD()
409       dEm_t.E = dEm1.E / (Theta * dt);  in RDIJacobian_FD()
410       dEm_t.T = dEm1.T / (Theta * dt);  in RDIJacobian_FD()
432       PetscScalar D, bcTheta = rd->bcmidpoint ? Theta : 1.;  in RDIJacobian_FD()
443 …K[0][1 * 2 + 0] = (bcTheta / Theta) * hx * (1. - 2. * D * (-1. / hx) - 2. * nx.E * dD.E + 2. * nx.…  in RDIJacobian_FD()
444         K[0][1 * 2 + 1] = (bcTheta / Theta) * hx * (-2. * nx.E * dD.T);  in RDIJacobian_FD()
445 …K[0][2 * 2 + 0] = (bcTheta / Theta) * hx * (-2. * D * (1. / hx) - 2. * nx.E * dD.E - 2. * nx.E * d…  in RDIJacobian_FD()
448         K[0][1 * 2 + 0] = -1. / Theta;  in RDIJacobian_FD()
449         K[0][2 * 2 + 0] = 1. / Theta;  in RDIJacobian_FD()
455       K[0][0 * 2 + 0] = -1. / Theta;  in RDIJacobian_FD()
456       K[0][1 * 2 + 0] = 1. / Theta;  in RDIJacobian_FD()
632   PetscReal     hx, Theta, dt, weight[5], interp[5][2], deriv[5][2];  in RDIFunction_FE()  local
637   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIFunction_FE()
657       rad = (1. - Theta) * RDRadiation(rd, &n0, 0) + Theta * RDRadiation(rd, &n, 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()
680       PetscReal   ratio, bcTheta = rd->bcmidpoint ? Theta : 1.;  in RDIFunction_FE()
718   PetscReal     hx, Theta, dt, weight[5], interp[5][2], deriv[5][2];  in RDIJacobian_FE()  local
724   PetscCall(RDGetLocalArrays(rd, ts, X, Xdot, &Theta, &dt, &X0loc, &x0, &Xloc, &x, &Xloc_t, &xdot));  in RDIJacobian_FE()
762       n.E  = (1 - Theta) * x0[0].E + Theta * x[0].E;  in RDIJacobian_FE()
763       n.T  = (1 - Theta) * x0[0].T + Theta * x[0].T;  in RDIJacobian_FE()