diff options
Diffstat (limited to 'src/c/matrixOperations/expm')
-rw-r--r-- | src/c/matrixOperations/expm/cexpma.c | 2 | ||||
-rw-r--r-- | src/c/matrixOperations/expm/dexpma.c | 4 | ||||
-rw-r--r-- | src/c/matrixOperations/expm/sexpma.c | 4 | ||||
-rw-r--r-- | src/c/matrixOperations/expm/zexpma.c | 2 |
4 files changed, 6 insertions, 6 deletions
diff --git a/src/c/matrixOperations/expm/cexpma.c b/src/c/matrixOperations/expm/cexpma.c index 7da0fc70..e7c9f0b2 100644 --- a/src/c/matrixOperations/expm/cexpma.c +++ b/src/c/matrixOperations/expm/cexpma.c @@ -69,7 +69,7 @@ void cexpma(floatComplex * in, floatComplex * out, int _iLeadDim) pfltMatrixX[iIndex1] = pfltMatrixA[iIndex1] ; - ceyea(pfltMatrixEye, _iLeadDim, _iLeadDim); + ceyea(pfltMatrixEye, _iLeadDim, _iLeadDim); /*cmulma ( & cfltCst , 1 ,1, pfltMatrixA , _iLeadDim, _iLeadDim, diff --git a/src/c/matrixOperations/expm/dexpma.c b/src/c/matrixOperations/expm/dexpma.c index d9b42327..8939813d 100644 --- a/src/c/matrixOperations/expm/dexpma.c +++ b/src/c/matrixOperations/expm/dexpma.c @@ -67,7 +67,7 @@ void dexpma (double* in, double* out, int _iLeadDim){ for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) pdblMatrixX[iIndex1] = pdblMatrixA[iIndex1] ; - deyea(pdblMatrixEye, _iLeadDim, _iLeadDim); + deyea(pdblMatrixEye, _iLeadDim, _iLeadDim); /*cA = A * c*/ @@ -106,7 +106,7 @@ void dexpma (double* in, double* out, int _iLeadDim){ pdblMatrixcX[iIndex1] = pdblMatrixX[iIndex1] * dblCst ; /*E = E + cX*/ - dadda ( out, iSquare , pdblMatrixcX , iSquare , out ) ; + dadda ( out, iSquare , pdblMatrixcX , iSquare , out ) ; if(iFlag == 1) /*D = D + cX*/ { diff --git a/src/c/matrixOperations/expm/sexpma.c b/src/c/matrixOperations/expm/sexpma.c index cd969981..b265a212 100644 --- a/src/c/matrixOperations/expm/sexpma.c +++ b/src/c/matrixOperations/expm/sexpma.c @@ -68,7 +68,7 @@ void sexpma (float* in, float* out, int _iLeadDim){ for ( iIndex1 = 0 ; iIndex1 < iSquare ; iIndex1++ ) pfltMatrixX[iIndex1] = pfltMatrixA[iIndex1] ; - seyea(pfltMatrixEye, _iLeadDim, _iLeadDim); + seyea(pfltMatrixEye, _iLeadDim, _iLeadDim); /*cA = A * c*/ @@ -110,7 +110,7 @@ void sexpma (float* in, float* out, int _iLeadDim){ pfltMatrixcX[iIndex1] = pfltMatrixX[iIndex1] * (float) dblCst ; /*E = E + cX*/ - sadda ( out, iSquare , pfltMatrixcX , iSquare , out ) ; + sadda ( out, iSquare , pfltMatrixcX , iSquare , out ) ; if(iFlag == 1) /*D = D + cX*/ { diff --git a/src/c/matrixOperations/expm/zexpma.c b/src/c/matrixOperations/expm/zexpma.c index 3de35e02..680f68ab 100644 --- a/src/c/matrixOperations/expm/zexpma.c +++ b/src/c/matrixOperations/expm/zexpma.c @@ -71,7 +71,7 @@ void zexpma(doubleComplex * in, doubleComplex * out, int _iLeadDim) pdblMatrixX[iIndex1] = pdblMatrixA[iIndex1] ; - zeyea(pdblMatrixEye, _iLeadDim, _iLeadDim); + zeyea(pdblMatrixEye, _iLeadDim, _iLeadDim); /* zmulma ( & zdblCst , 1 ,1, pdblMatrixA , _iLeadDim, _iLeadDim, |