diff --git a/usadelndsoc/solver.py b/usadelndsoc/solver.py
index 3bafa3a311f96ef84dc062c619df7a69b7d6aa70..cb32231f051bc3c2a76023c9faefb54e17b91020 100644
--- a/usadelndsoc/solver.py
+++ b/usadelndsoc/solver.py
@@ -1692,7 +1692,7 @@ def cpr(
         state.Deltas = [solver.Delta.copy()]
         state.phis = [0]
         state.dss = [0]
-        state.Js = [solver.J_tot.A.copy()]
+        state.Js = [solver.J_tot.copy()]
 
     phi = state.phis[-1] + ds / 4
 
@@ -1705,7 +1705,7 @@ def cpr(
         )
         state.Deltas.append(solver.Delta.copy())
         state.phis.append(float(phi))
-        state.Js.append(solver.J_tot.A.copy())
+        state.Js.append(solver.J_tot.copy())
 
     for j in range(len(state.phis), max_points):
         if auto_stop and len(state.phis) > 2: