488 static void RDEvaluate(PetscReal interp[][2], PetscReal deriv[][2], PetscInt q, const RDNode x[], P…  in RDEvaluate()
496     n->E += interp[q][j] * x[i + j].E;  in RDEvaluate()
497     n->T += interp[q][j] * x[i + j].T;  in RDEvaluate()
506 …drature(RD rd, PetscReal hx, PetscInt *nq, PetscReal weight[], PetscReal interp[][2], PetscReal de…  in RDGetQuadrature()
617       interp[q][j] = refinterp[q][j];  in RDGetQuadrature()
632   PetscReal     hx, Theta, dt, weight[5], interp[5][2], deriv[5][2];  in RDIFunction_FE()  local
646   PetscCall(RDGetQuadrature(rd, hx, &nq, weight, interp, deriv));  in RDIFunction_FE()
653       RDEvaluate(interp, deriv, q, x, i, &n, &nx);  in RDIFunction_FE()
654       RDEvaluate(interp, deriv, q, x0, i, &n0, &n0x);  in RDIFunction_FE()
655       RDEvaluate(interp, deriv, q, xdot, i, &nt, &ntx);  in RDIFunction_FE()
671 …q][j] * weight[q] * ((1. - Theta) * D0_R * n0x.E + Theta * D_R * nx.E) + interp[q][j] * weight[q] …  in RDIFunction_FE()
672         f[i + j].T += interp[q][j] * weight[q] * (rho * Em_t + rad);  in RDIFunction_FE()
718   PetscReal     hx, Theta, dt, weight[5], interp[5][2], deriv[5][2];  in RDIJacobian_FE()  local
727   PetscCall(RDGetQuadrature(rd, hx, &nq, weight, interp, deriv));  in RDIJacobian_FE()
739       RDEvaluate(interp, deriv, q, x, i, &n, &nx);  in RDIJacobian_FE()
740       RDEvaluate(interp, deriv, q, xdot, i, &nt, &ntx);  in RDIJacobian_FE()
746 …= (+interp[q][j] * weight[q] * (a - drad.E) * interp[q][k] + deriv[q][j] * weight[q] * ((D_R + dxD…  in RDIJacobian_FE()
747 …k * 2 + 1] += (+interp[q][j] * weight[q] * (-drad.T * interp[q][k]) + deriv[q][j] * weight[q] * (d…  in RDIJacobian_FE()
748           K[j * 2 + 1][k * 2 + 0] += interp[q][j] * weight[q] * drad.E * interp[q][k];  in RDIJacobian_FE()
749 …K[j * 2 + 1][k * 2 + 1] += interp[q][j] * weight[q] * (a * rd->rho * dEm.T + drad.T) * interp[q][k…  in RDIJacobian_FE()