334       PetscScalar Em0, Em1;  in RDIFunction_FD()  local
336       RDMaterialEnergy(rd, &x[i], &Em1, NULL);  in RDIFunction_FD()
337       Em_t = (Em1 - Em0) / dt;  in RDIFunction_FD()
404       PetscScalar Em0, Em1;  in RDIJacobian_FD()  local
407       RDMaterialEnergy(rd, &x[i], &Em1, &dEm1);  in RDIJacobian_FD()
659         PetscScalar Em0, Em1;  in RDIFunction_FE()  local
661         RDMaterialEnergy(rd, &n, &Em1, 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
884     rd->MaterialEnergy(rd, &n, &Em1, 0);  in RDTestDifferentiation()
885     fdEm.E = (Em1 - Em0) / epsilon;  in RDTestDifferentiation()
888     rd->MaterialEnergy(rd, &n, &Em1, 0);  in RDTestDifferentiation()
889     fdEm.T = (Em1 - Em0) / (T0 * epsilon);  in RDTestDifferentiation()