diff options
Diffstat (limited to 'packages/base/src/Internal/C/lapack-aux.c')
-rw-r--r-- | packages/base/src/Internal/C/lapack-aux.c | 12 |
1 files changed, 4 insertions, 8 deletions
diff --git a/packages/base/src/Internal/C/lapack-aux.c b/packages/base/src/Internal/C/lapack-aux.c index 177d373..ff7ad92 100644 --- a/packages/base/src/Internal/C/lapack-aux.c +++ b/packages/base/src/Internal/C/lapack-aux.c | |||
@@ -55,10 +55,6 @@ typedef float complex TCF; | |||
55 | #define NODEFPOS 2006 | 55 | #define NODEFPOS 2006 |
56 | #define NOSPRTD 2007 | 56 | #define NOSPRTD 2007 |
57 | 57 | ||
58 | inline int mod (int a, int b); | ||
59 | |||
60 | inline int64_t mod_l (int64_t a, int64_t b); | ||
61 | |||
62 | //////////////////////////////////////////////////////////////////////////////// | 58 | //////////////////////////////////////////////////////////////////////////////// |
63 | void asm_finit() { | 59 | void asm_finit() { |
64 | #ifdef i386 | 60 | #ifdef i386 |
@@ -247,7 +243,7 @@ int svd_l_C(OCMAT(a),OCMAT(u), DVEC(s),OCMAT(v)) { | |||
247 | ldvt = q; | 243 | ldvt = q; |
248 | } | 244 | } |
249 | }DEBUGMSG("svd_l_C"); | 245 | }DEBUGMSG("svd_l_C"); |
250 | 246 | ||
251 | double *rwork = (double*) malloc(5*q*sizeof(double)); | 247 | double *rwork = (double*) malloc(5*q*sizeof(double)); |
252 | CHECK(!rwork,MEM); | 248 | CHECK(!rwork,MEM); |
253 | integer lwork = -1; | 249 | integer lwork = -1; |
@@ -423,7 +419,7 @@ int eig_l_R(ODMAT(a),ODMAT(u), CVEC(s),ODMAT(v)) { | |||
423 | //////////////////// symmetric real eigensystem //////////// | 419 | //////////////////// symmetric real eigensystem //////////// |
424 | 420 | ||
425 | int dsyev_(char *jobz, char *uplo, integer *n, doublereal *a, | 421 | int dsyev_(char *jobz, char *uplo, integer *n, doublereal *a, |
426 | integer *lda, doublereal *w, doublereal *work, integer *lwork, | 422 | integer *lda, doublereal *w, doublereal *work, integer *lwork, |
427 | integer *info); | 423 | integer *info); |
428 | 424 | ||
429 | int eig_l_S(int wantV,DVEC(s),ODMAT(v)) { | 425 | int eig_l_S(int wantV,DVEC(s),ODMAT(v)) { |
@@ -665,7 +661,7 @@ int linearSolveLSC_l(OCMAT(a),OCMAT(b)) { | |||
665 | int dgelss_(integer *m, integer *n, integer *nrhs, | 661 | int dgelss_(integer *m, integer *n, integer *nrhs, |
666 | doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal * | 662 | doublereal *a, integer *lda, doublereal *b, integer *ldb, doublereal * |
667 | s, doublereal *rcond, integer *rank, doublereal *work, integer *lwork, | 663 | s, doublereal *rcond, integer *rank, doublereal *work, integer *lwork, |
668 | integer *info); | 664 | integer *info); |
669 | 665 | ||
670 | int linearSolveSVDR_l(double rcond,ODMAT(a),ODMAT(b)) { | 666 | int linearSolveSVDR_l(double rcond,ODMAT(a),ODMAT(b)) { |
671 | integer m = ar; | 667 | integer m = ar; |
@@ -955,7 +951,7 @@ int schur_l_R(ODMAT(u), ODMAT(s)) { | |||
955 | int zgees_(char *jobvs, char *sort, L_fp select, integer *n, | 951 | int zgees_(char *jobvs, char *sort, L_fp select, integer *n, |
956 | doublecomplex *a, integer *lda, integer *sdim, doublecomplex *w, | 952 | doublecomplex *a, integer *lda, integer *sdim, doublecomplex *w, |
957 | doublecomplex *vs, integer *ldvs, doublecomplex *work, integer *lwork, | 953 | doublecomplex *vs, integer *ldvs, doublecomplex *work, integer *lwork, |
958 | doublereal *rwork, logical *bwork, integer *info); | 954 | doublereal *rwork, logical *bwork, integer *info); |
959 | 955 | ||
960 | int schur_l_C(OCMAT(u), OCMAT(s)) { | 956 | int schur_l_C(OCMAT(u), OCMAT(s)) { |
961 | integer m = sr; | 957 | integer m = sr; |