diff options
| -rw-r--r-- | src/matrixOperations/division/zrdivma.c | 24 |
1 files changed, 11 insertions, 13 deletions
diff --git a/src/matrixOperations/division/zrdivma.c b/src/matrixOperations/division/zrdivma.c index 16fb0398..ae51dd7a 100644 --- a/src/matrixOperations/division/zrdivma.c +++ b/src/matrixOperations/division/zrdivma.c @@ -12,7 +12,6 @@ #include "matrixDivision.h" -#include <stdio.h> #include <string.h> void zrdivma( doubleComplex* in1, int lines1, int columns1 , @@ -99,17 +98,17 @@ void zrdivma( doubleComplex* in1, int lines1, int columns1 , cNorm = 'N'; C2F(zgetrs)(&cNorm, &columns2, &lines1, poAf, &columns2, pIpiv, poBt, &columns2, &iInfo); { - int i,j,ij,ji; - for(j = 0 ; j < lines2 ; j++) - { - for(i = 0 ; i < lines1 ; i++) + int i,j,ij,ji; + for(j = 0 ; j < lines2 ; j++) { - ij = i + j * lines1; - ji = j + i * max(lines2, columns2); - out[ij] = DoubleComplex ( zreals( poBt[ji]) , -zimags ( poBt[ji])); - printf ( "\n\t\t\t<debug>%e + %e\n " , zreals( poBt[ji]) , -zimags ( poBt[ji])); - } - } + for(i = 0 ; i < lines1 ; i++) + { + ij = i + j * lines1; + ji = j + i * max(lines2, columns2); + out[ij] = DoubleComplex ( zreals( poBt[ji]) , -zimags ( poBt[ji])); + + } + } } iExit = 1; @@ -144,8 +143,7 @@ void zrdivma( doubleComplex* in1, int lines1, int columns1 , ij = i + j * lines1; ji = j + i * max(lines2, columns2); out[ij] = DoubleComplex ( zreals( poBt[ji]) , -zimags ( poBt[ji])); - printf ( "\n\t\t\t<debug>%e + %e\n " , zreals( poBt[ji]) , -zimags ( poBt[ji])); - } + } } } } |