Modify without lapack part

1 files changed, 7 insertions, 0 deletions

diff --git a/src/matrixOperations/chol/dchola.c b/src/matrixOperations/chol/dchola.c
index a35f90b2..a627b236 100644
--- a/src/matrixOperations/chol/dchola.c
+++ b/src/matrixOperations/chol/dchola.c
@@ -58,5 +58,12 @@ void dchola(double * in, int size, double *out){
 		out[i*size+i]=sqrt(in[i*size+i]-accu);	
 	}
 
+	/*Zeros in the lower triangular part*/
+	for (i=0;i<size;i++){
+		for (j=i+1;j<size;j++){
+			out[j+i*size]=0;
+		}
+	}
+
 #endif
 }