Fixed error in the Non-Lapack part

1 files changed, 7 insertions, 0 deletions

diff --git a/src/matrixOperations/chol/cchola.c b/src/matrixOperations/chol/cchola.c
index 4ebf473e..0de211be 100644
--- a/src/matrixOperations/chol/cchola.c
+++ b/src/matrixOperations/chol/cchola.c
@@ -77,6 +77,13 @@ void cchola(floatComplex * in, int size, floatComplex *out){
 		tmp = FloatComplex(creals(in[i*size+i])-accu,0);
 		out[i*size+i]=csqrts(tmp);	
 	}
+	
+	/*Zeros in the lower triangular part*/
+	for (i=0;i<size;i++){
+		for (j=i+1;j<size;j++){
+			out[j+i*size]=FloatComplex(0,0);
+		}
+	}
 
 #endif	
 }