diff --git a/src/linalg/decompositions.rs b/src/linalg/decompositions.rs
index e56e8cc4..a67ef389 100644
--- a/src/linalg/decompositions.rs
+++ b/src/linalg/decompositions.rs
@@ -145,9 +145,6 @@ pub fn cholesky<N, V, VS, M>(m: &M) -> Result<M, &'static str>
             if i > j {
                 out[(i, j)] = sum / out[(j, j)];
             }
-            else if i < j {
-                out[(i,j)] = N::zero();
-            }
             else if sum > N::zero() {
                 out[(i,i)] = sum.sqrt();
             }