diff options
Diffstat (limited to 'src/matrixOperations/division/zrdiva.c')
| -rw-r--r-- | src/matrixOperations/division/zrdiva.c | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/src/matrixOperations/division/zrdiva.c b/src/matrixOperations/division/zrdiva.c index f3b2eac3..43ed4d0f 100644 --- a/src/matrixOperations/division/zrdiva.c +++ b/src/matrixOperations/division/zrdiva.c @@ -112,7 +112,7 @@ int zrdiva( doubleComplex* in1, int lines1, int columns1 , dblRcond = sqrt(dblEps); cNorm = 'F'; iMax = Max(lines2, columns2); - printf ( "columns2 : %d\n" , columns2 ) ; + memset(pJpvt, 0x00,(unsigned int) sizeof(int) *(unsigned int) lines2); C2F(zgelsy)(&columns2, &lines2, &lines1, poAt, &columns2, poBt, &iMax, pJpvt, &dblRcond, &pRank[0], poDwork, &iWork, pRwork, &iInfo); @@ -133,7 +133,8 @@ int zrdiva( 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])); + } } } } |