398     RDNode                    dEm_t, drad;  in RDIJacobian_FD()  local
401     RDRadiation(rd, &x[i], &drad);  in RDIJacobian_FD()
463       K[0][1 * 2 + 0] = hx * (a - ddiff[1].E - drad.E);  in RDIJacobian_FD()
464       K[0][1 * 2 + 1] = hx * (-ddiff[1].T - drad.T);  in RDIJacobian_FD()
469     K[1][1 * 2 + 0] = hx * (rho * dEm_t.E + drad.E);  in RDIJacobian_FD()
470     K[1][1 * 2 + 1] = hx * (rho * dEm_t.T + drad.T);  in RDIJacobian_FD()
738       RDNode                   n, nx, nt, ntx, drad, dD_R, dxD_R, dEm;  in RDIJacobian_FE()  local
741       rad = RDRadiation(rd, &n, &drad);  in RDIJacobian_FE()
746 …K[j * 2 + 0][k * 2 + 0] += (+interp[q][j] * weight[q] * (a - drad.E) * interp[q][k] + deriv[q][j] …  in RDIJacobian_FE()
747 …K[j * 2 + 0][k * 2 + 1] += (+interp[q][j] * weight[q] * (-drad.T * interp[q][k]) + deriv[q][j] * w…  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()
956     RDNode      drad, fdrad;  in RDTestDifferentiation()  local
959     rad0    = RDRadiation(rd, &n, &drad);  in RDTestDifferentiation()
968 …alPart(drad.E), (double)PetscRealPart(drad.T), (double)PetscRealPart(fdrad.E), (double)PetscRealPa…  in RDTestDifferentiation()
969                           (double)PetscRealPart(drad.T - fdrad.T)));  in RDTestDifferentiation()