|
1 | 1 | #ifndef LAPACK_HPP |
2 | 2 | #define LAPACK_HPP |
3 | | - |
| 3 | +#include <iostream> |
4 | 4 | extern "C" |
5 | 5 | { |
6 | 6 | // ================================================================================= |
@@ -49,10 +49,10 @@ extern "C" |
49 | 49 | // ================================================================================= |
50 | 50 | // gvx |
51 | 51 | void dsygvx_(const int* itype, const char* jobz, const char* range, const char* uplo, |
52 | | - const int* n, double* A, const int* lda, double* B, const int* ldb, |
53 | | - const double* vl, const double* vu, const int* il, const int* iu, |
54 | | - const double* abstol, const int* m, double* w, double* Z, const int* ldz, |
55 | | - double* work, int* lwork, int*iwork, int* ifail, int* info); |
| 52 | + const int* n, double* A, const int* lda, double* B, const int* ldb, |
| 53 | + const double* vl, const double* vu, const int* il, const int* iu, |
| 54 | + const double* abstol, const int* m, double* w, double* Z, const int* ldz, |
| 55 | + double* work, const int* lwork, int* iwork, int* ifail, int* info); |
56 | 56 |
|
57 | 57 | void chegvx_(const int* itype,const char* jobz,const char* range,const char* uplo, |
58 | 58 | const int* n,std::complex<float> *a,const int* lda,std::complex<float> *b, |
@@ -424,8 +424,8 @@ class LapackWrapper |
424 | 424 | int* ifail, |
425 | 425 | int& info) |
426 | 426 | { |
427 | | - // dsygvx_(&itype, &jobz, &range, &uplo, &n, a, &lda, b, &ldb, &vl, |
428 | | - // &vu, &il,&iu, &abstol, &m, w, z, &ldz, work, &lwork, rwork, iwork, ifail, &info); |
| 427 | + dsygvx_(&itype, &jobz, &range, &uplo, &n, a, &lda, b, &ldb, &vl, |
| 428 | + &vu, &il, &iu, &abstol, &m, w, z, &ldz, work, &lwork, iwork, ifail, &info); |
429 | 429 | } |
430 | 430 |
|
431 | 431 | // wrap function of fortran lapack routine xhegvx ( pointer version ). |
|
0 commit comments