diff --git a/ndarray-linalg/Cargo.toml b/ndarray-linalg/Cargo.toml
index 8c5c83a8..70cb1141 100644
--- a/ndarray-linalg/Cargo.toml
+++ b/ndarray-linalg/Cargo.toml
@@ -38,7 +38,7 @@ rand = "0.8.3"
 thiserror = "2.0.0"
 
 [dependencies.ndarray]
-version = "0.16.0"
+version = "0.17.1"
 features = ["approx", "std"]
 default-features = false
 
diff --git a/ndarray-linalg/src/lobpcg/eig.rs b/ndarray-linalg/src/lobpcg/eig.rs
index e60adb04..7349c36d 100644
--- a/ndarray-linalg/src/lobpcg/eig.rs
+++ b/ndarray-linalg/src/lobpcg/eig.rs
@@ -69,7 +69,10 @@ impl<A: Float + Scalar + ScalarOperand + Lapack + PartialOrd + Default> Truncate
             lobpcg(
                 |y| self.problem.dot(&y),
                 x,
-                |mut y| y.assign(&preconditioner.dot(&y)),
+                |mut y| {
+                    let p = preconditioner.dot(&y);
+                    y.assign(&p);
+                },
                 self.constraints.clone(),
                 self.precision,
                 self.maxiter,