488 static void RDEvaluate(PetscReal interp[][2], PetscReal deriv[][2], PetscInt q, const RDNode x[], P…  in RDEvaluate()
498     nx->E += deriv[q][j] * x[i + j].E;  in RDEvaluate()
499     nx->T += deriv[q][j] * x[i + j].T;  in RDEvaluate()
506 …RD rd, PetscReal hx, PetscInt *nq, PetscReal weight[], PetscReal interp[][2], PetscReal deriv[][2])  in RDGetQuadrature()
618       deriv[q][j]  = refderiv[q][j] / hx;  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 …f[i + j].E += (deriv[q][j] * weight[q] * ((1. - Theta) * D0_R * n0x.E + Theta * D_R * nx.E) + inte…  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 …[j] * weight[q] * (a - drad.E) * interp[q][k] + deriv[q][j] * weight[q] * ((D_R + dxD_R.E * nx.E) …  in RDIJacobian_FE()
747 …+interp[q][j] * weight[q] * (-drad.T * interp[q][k]) + deriv[q][j] * weight[q] * (dxD_R.T * deriv[…  in RDIJacobian_FE()