diff --git a/src/compas_fd/solvers/fd_solver.py b/src/compas_fd/solvers/fd_solver.py
index 7371e89a..01229b60 100644
--- a/src/compas_fd/solvers/fd_solver.py
+++ b/src/compas_fd/solvers/fd_solver.py
@@ -17,6 +17,7 @@ def solve(self) -> None:
         nd = self.numdata
         b = nd.p[nd.free] - nd.Df.dot(nd.xyz[nd.fixed])
         nd.xyz[nd.free] = spsolve(nd.Di, b)
+        nd.reset_xyz()
 
     @property
     def is_converged(self) -> bool: