diff --git a/q3d/firedrakeplus/solver.py b/q3d/firedrakeplus/solver.py
index a2a10d4..2e30fb5 100644
--- a/q3d/firedrakeplus/solver.py
+++ b/q3d/firedrakeplus/solver.py
@@ -224,9 +224,6 @@ def _non_graddesc_solve(times_list, q_soln, bcs, solver_parameters, newton_param
         except ConvergenceError:
             pr.fail('Convergence error')
             return False
-        
-        # add the energy of q_soln to the energies
-        energies.append(compute_energy(q_soln))
 
         # print energy
         pr.iter_info_verbose(f'NON GD SOLVE COMPLETED', f'energy = {energies[-1]}', i=len(energies))