334       PetscScalar Em0, Em1;  in RDIFunction_FD()  local
335       RDMaterialEnergy(rd, &x0[i], &Em0, NULL);  in RDIFunction_FD()
337       Em_t = (Em1 - Em0) / dt;  in RDIFunction_FD()
404       PetscScalar Em0, Em1;  in RDIJacobian_FD()  local
406       RDMaterialEnergy(rd, &x0[i], &Em0, NULL);  in RDIJacobian_FD()
659         PetscScalar Em0, Em1;  in RDIFunction_FE()  local
660         RDMaterialEnergy(rd, &n0, &Em0, NULL);  in RDIFunction_FE()
662         Em_t = (Em1 - Em0) / dt;  in RDIFunction_FE()
878     PetscScalar T0 = 1000., T1 = T0 * (1. + epsilon), Em0, Em1;  in RDTestDifferentiation()  local
881     rd->MaterialEnergy(rd, &n, &Em0, &dEm);  in RDTestDifferentiation()
885     fdEm.E = (Em1 - Em0) / epsilon;  in RDTestDifferentiation()
889     fdEm.T = (Em1 - Em0) / (T0 * epsilon);  in RDTestDifferentiation()