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/zexpma.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/zexpma.c')
| -rw-r--r-- | src/matrixOperations/expm/zexpma.c | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/src/matrixOperations/expm/zexpma.c b/src/matrixOperations/expm/zexpma.c index 7ee85f17..9b378690 100644 --- a/src/matrixOperations/expm/zexpma.c +++ b/src/matrixOperations/expm/zexpma.c @@ -73,9 +73,12 @@ void zexpma(doubleComplex * in, doubleComplex * out, int _iLeadDim) zeyesa(pdblMatrixEye, _iLeadDim, _iLeadDim); - zmulma ( & zdblCst , 1 ,1, + /* zmulma ( & zdblCst , 1 ,1, pdblMatrixA , _iLeadDim, _iLeadDim, - pdblMatrixcA); + pdblMatrixcA);*/ + + for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) + pdblMatrixcA[iIndex1] = ztimess ( pdblMatrixA[iIndex1] , zdblCst ) ; /* cA = A * c */ @@ -107,9 +110,12 @@ void zexpma(doubleComplex * in, doubleComplex * out, int _iLeadDim) pdblMatrixX ); /* cX = c * X */ - zmulma ( & zdblCst , 1 ,1, + /* zmulma ( & zdblCst , 1 ,1, pdblMatrixX , _iLeadDim, _iLeadDim, - pdblMatrixcX); + pdblMatrixcX);*/ + + for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) + pdblMatrixcX[iIndex1] = ztimess ( pdblMatrixX[iIndex1] , zdblCst ) ; /* E = E + cX */ |