diff options
| author | simon | 2008-08-07 15:00:45 +0000 |
|---|---|---|
| committer | simon | 2008-08-07 15:00:45 +0000 |
| commit | 9633cfcb36ba18c46282c8e8a7e5dd7352212c40 (patch) | |
| tree | e766cc4a2d03d979eeba92c0053afa60aa925943 /src/includes | |
| parent | c80b732994ec39abd96da72e03fb95dfb6997422 (diff) | |
| download | scilab2c-9633cfcb36ba18c46282c8e8a7e5dd7352212c40.tar.gz scilab2c-9633cfcb36ba18c46282c8e8a7e5dd7352212c40.tar.bz2 scilab2c-9633cfcb36ba18c46282c8e8a7e5dd7352212c40.zip | |
added matrixExponential and all related sub functions
matrixExponential just work with real Double for the moment
Diffstat (limited to 'src/includes')
| -rw-r--r-- | src/includes/blas.h | 16 | ||||
| -rw-r--r-- | src/includes/lapack.h | 45 |
2 files changed, 48 insertions, 13 deletions
diff --git a/src/includes/blas.h b/src/includes/blas.h index fc12ac17..86ab62f6 100644 --- a/src/includes/blas.h +++ b/src/includes/blas.h @@ -13,6 +13,11 @@ #ifndef __BLAS_H__ #define __BLAS_H__ + +#ifndef _MACRO_C2F_ +#define _MACRO_C2F_ +#define C2F(name) name##_ +#endif /* SUBROUTINE DGEMM ( TRANSA, TRANSB, M, N, K, ALPHA, A, LDA, B, LDB, $ BETA, C, LDC ) @@ -139,10 +144,17 @@ C (NaN*0 should be NaN, not 0) * * Level 3 Blas routine. */ - +/* void dgemm_(char *TRANSA, char* TRANSB, int *M, int *N, int *K, double *ALPHA, double *A, int *LDA, double *B, int *LDB, double *BETA, - double *C, int *LDC); + double *C, int *LDC);*/ + +extern int C2F(dgemm)(); +extern int C2F(idamax)() ;/* could be transcribe easaly in c */ +extern int C2F(daxpy) () ;/* could be transcribe easaly in c */ +extern int C2F(dscal) () ;/* could be transcribe easaly in c */ +extern int C2F(dasum) () ;/* could be transcribe easaly in c */ + #endif /* !__BLAS_H__ */ diff --git a/src/includes/lapack.h b/src/includes/lapack.h index 52faccf5..1bc658ee 100644 --- a/src/includes/lapack.h +++ b/src/includes/lapack.h @@ -13,7 +13,11 @@ #ifndef __LAPACK_H__ #define __LAPACK_H__ +#ifndef _MACRO_C2F_ +#define _MACRO_C2F_ #define C2F(name) name##_ +#endif + #define getRelativeMachinePrecision() dlamch_("e", 1L) #define getOverflowThreshold() dlamch_("o", 1L) @@ -89,23 +93,42 @@ extern double dlacpy_ (char* NORM, int* M, int* N, double* A, int* LDA, /****** doubleComplex fortran function ************/ +extern double zgelsy_ () ; + + +extern double zlamch_ (); + +extern double zlange_ (); + +extern double zgetrf_ (); -extern int zlamch_ (); +extern double zgecon_ ( ) ; -extern int zlange_ (); +extern double zgetrs_ ( ) ; -extern int zgetrf_ (); +extern double zlacpy_ (); +/*extern int zgelsy_ ();*/ -extern int zgecon_ ( ) ; +/*certainly have some blas functions in */ +extern int C2F(split)(); +extern int C2F(exch)(); -extern int zgetrs_ ( ) ; +extern int C2F(balbak)(); +extern double C2F(ddot)(); +extern int C2F(pade)(); +extern int C2F(dcopy)(); +extern int C2F(dscal)(); -extern int zgelsy_ (); -/*extern doubleComplex zgelsy_ (int* M, int* N, int* NRHS, doubleComplex* A, int* LDA, - doubleComplex* B, int* LDB, int* JPVT, double* RCOND, int* RANK, - doubleComplex* WORK, int* LWORK, double* RWORK, int* INFO) ; -*/ +extern int C2F(dgeco)(); +extern int C2F(dgesl)(); +extern int C2F(coef)(); +extern int C2F(cerr)(); +extern int C2F(dclmat)(); +extern int C2F(dexpm1)(); +extern int C2F(wexpm1)(); +extern int C2F(drot)(); +extern int C2F(intexpm) (); -extern int zlacpy_ (); +extern int C2F(zcopy)(); #endif /* !__LAPACK_H__ */ |