const double* x, const int* incx,
const double* b, double* y, const int* incy);
// Compute eigenvalues and vectors
-void dgeev_ (char* jobvl, char* jobvr,
- int* n, double* A, int* lda,
+void dgeev_ (const char* jobvl, const char* jobvr,
+ const int* n, double* A, const int* lda,
double* lambda_re, double* lambda_im,
- double* vl, int* ldvl,
- double* vr, int* ldva,
- double* work, int* lwork,
+ double* vl, const int* ldvl,
+ double* vr, const int* ldva,
+ double* work, const int* lwork,
int* info);
// Compute eigenvalues and vectors (expert)
-void dgeevx_ (char* balanc, char* jobvl, char* jobvr, char* sense,
- int* n, double* A, int* lda,
+void dgeevx_ (const char* balanc, const char* jobvl, const char* jobvr,
+ const char* sense,
+ const int* n, double* A, const int* lda,
double* lambda_re, double* lambda_im,
- double* vl, int* ldvl,
- double* vr, int* ldvr,
+ double* vl, const int* ldvl,
+ double* vr, const int* ldvr,
int* ilo, int* ihi,
double* scale, double* abnrm,
double* rconde, double* rcondv,
- double* work, int* lwork,
+ double* work, const int* lwork,
int* iwork, int* info);
// Compute singular value decomposition
void dgesvd_ (int* jobu, int* jobvt,
- int* n, int* m, double* A, int* lda,
+ const int* n, const int* m, double* A, const int* lda,
double* s,
- double* u, int* ldu,
- double* vt, int* ldvt,
- double* work, int* lwork,
+ double* u, const int* ldu,
+ double* vt, const int* ldvt,
+ double* work, const int* lwork,
int* info);