diff options
| author | torset | 2009-02-16 15:31:16 +0000 |
|---|---|---|
| committer | torset | 2009-02-16 15:31:16 +0000 |
| commit | 2e41ca898e113e9fe0cb7c27c60a5f0ac2a74517 (patch) | |
| tree | 8f78c74b06936ed4c445274bdbbbc74e650282cd /src/matrixOperations/chol/zchola.c | |
| parent | 92d3899c14d576c20762d7860f5e2f137923da0c (diff) | |
| download | scilab2c-2e41ca898e113e9fe0cb7c27c60a5f0ac2a74517.tar.gz scilab2c-2e41ca898e113e9fe0cb7c27c60a5f0ac2a74517.tar.bz2 scilab2c-2e41ca898e113e9fe0cb7c27c60a5f0ac2a74517.zip | |
Fixed error in the Non-Lapack part
Diffstat (limited to 'src/matrixOperations/chol/zchola.c')
| -rw-r--r-- | src/matrixOperations/chol/zchola.c | 12 |
1 files changed, 10 insertions, 2 deletions
diff --git a/src/matrixOperations/chol/zchola.c b/src/matrixOperations/chol/zchola.c index ce6543ff..83661dc5 100644 --- a/src/matrixOperations/chol/zchola.c +++ b/src/matrixOperations/chol/zchola.c @@ -10,7 +10,7 @@ * */ -#ifdef WITHOUT_LAPACK +#ifndef WITHOUT_LAPACK #include "lapack.h" #else #include "multiplication.h" @@ -29,7 +29,7 @@ void zchola(doubleComplex * in, int size, doubleComplex *out){ param out : output upper triangular matrix */ -#ifdef WITHOUT_LAPACK +#ifndef WITHOUT_LAPACK int i=0, j=0, info=0; for(i=0;i<size*size;i++) out[i]=in[i]; @@ -65,6 +65,14 @@ void zchola(doubleComplex * in, int size, doubleComplex *out){ tmp = DoubleComplex(zreals(in[i*size+i])-accu,0); out[i*size+i]=zsqrts(tmp); } + + /*Zeros in the lower triangular part*/ + for (i=0;i<size;i++){ + for (j=i+1;j<size;j++){ + out[j+i*size]=DoubleComplex(0,0); + } + } + #endif |