diff options
| author | simon | 2008-08-11 12:32:00 +0000 |
|---|---|---|
| committer | simon | 2008-08-11 12:32:00 +0000 |
| commit | b074501da9007198d67b37dd3527eeef2bb51dc8 (patch) | |
| tree | b93ba517a3fa0e2b107f28803e248b7f29d18aaa /src/matrixOperations/expm/dexpma.c | |
| parent | 652e93fcb2d27ad9daddd5db06a502010d173253 (diff) | |
| download | scilab2c-b074501da9007198d67b37dd3527eeef2bb51dc8.tar.gz scilab2c-b074501da9007198d67b37dd3527eeef2bb51dc8.tar.bz2 scilab2c-b074501da9007198d67b37dd3527eeef2bb51dc8.zip | |
all test passed now but problem of precision
Diffstat (limited to 'src/matrixOperations/expm/dexpma.c')
| -rw-r--r-- | src/matrixOperations/expm/dexpma.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/matrixOperations/expm/dexpma.c b/src/matrixOperations/expm/dexpma.c index 9bfc4ea2..3bb3156b 100644 --- a/src/matrixOperations/expm/dexpma.c +++ b/src/matrixOperations/expm/dexpma.c @@ -59,7 +59,7 @@ void dexpma (double* in, double* out, int _iLeadDim){ for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) pdblMatrixA[iIndex1] = in[iIndex1] / dblS2 ; - + /* Pade approximation for exp(A)*/ /*X = A */ /*C2F(dcopy)(&iSquare, pdblMatrixA, &iOne, pdblMatrixX, &iOne );*/ @@ -89,7 +89,7 @@ void dexpma (double* in, double* out, int _iLeadDim){ for(iLoop1 = 2 ; iLoop1 <= iMax ; iLoop1++) { dblCst = dblCst * (iMax - iLoop1 + 1 ) / (iLoop1 * (2 * iMax - iLoop1 + 1)); - + /*Temp = X*/ /*C2F(dcopy)(&iSquare, pdblMatrixX, &iOne, pdblMatrixTemp, &iOne);*/ for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) |