summaryrefslogtreecommitdiff
path: root/src/matrixOperations/chol/cchola.c
diff options
context:
space:
mode:
Diffstat (limited to 'src/matrixOperations/chol/cchola.c')
-rw-r--r--src/matrixOperations/chol/cchola.c11
1 files changed, 6 insertions, 5 deletions
diff --git a/src/matrixOperations/chol/cchola.c b/src/matrixOperations/chol/cchola.c
index 0de211be..ab576515 100644
--- a/src/matrixOperations/chol/cchola.c
+++ b/src/matrixOperations/chol/cchola.c
@@ -10,6 +10,9 @@
*
*/
+#include <malloc.h>
+#include <stdio.h>
+
#ifndef WITHOUT_LAPACK
#include "lapack.h"
#else
@@ -18,10 +21,8 @@
#include "subtraction.h"
#include "sqrt.h"
#endif
-
#include "chol.h"
-#include <malloc.h>
-#include <stdio.h>
+
void cchola(floatComplex * in, int size, floatComplex *out){
/* param in : input matrix (square matrix)
@@ -41,10 +42,10 @@ void cchola(floatComplex * in, int size, floatComplex *out){
doubleComplex* tmp;
int i=0, j=0, info=0;
- tmp=malloc((unsigned int)(size*size)*sizeof(doubleComplex));
+ tmp=(doubleComplex*)malloc((unsigned int)(size*size)*sizeof(doubleComplex));
for(i=0;i<size*size;i++) tmp[i]=DoubleComplex((double)creals(in[i]),(double)cimags(in[i]));
- C2F(zpotrf)("U", &size, tmp, &size, &info,1L);
+ zpotrf_("U", &size, tmp, &size, &info);
for(i=0;i<size*size;i++) out[i]=FloatComplex((float)zreals(tmp[i]),(float)zimags(tmp[i]));