Avoid PPCG dependency on hegvd wrapper change

Author: cheerly-pku

source/source_base/module_container/base/third_party/lapack.h

@@ -228,7 +228,7 @@ void hegvd(const int itype, const char jobz, const char uplo, const int n,
             float* a, const int lda,
             float* b, const int ldb, float* w,
             float* work, int lwork, float* rwork, int lrwork,
-            int* iwork, int liwork, int& info)
+            int* iwork, int liwork, int info)
 {
     // call the fortran routine
     ssygvd_(&itype, &jobz, &uplo, &n,
@@ -242,7 +242,7 @@ void hegvd(const int itype, const char jobz, const char uplo, const int n,
             double* a, const int lda,
             double* b, const int ldb, double* w,
             double* work, int lwork, double* rwork, int lrwork,
-            int* iwork, int liwork, int& info)
+            int* iwork, int liwork, int info)
 {
     // call the fortran routine
     dsygvd_(&itype, &jobz, &uplo, &n,
@@ -255,7 +255,7 @@ void hegvd(const int itype, const char jobz, const char uplo, const int n,
             std::complex<float>* a, const int lda,
             std::complex<float>* b, const int ldb, float* w,
             std::complex<float>* work, int lwork, float* rwork, int lrwork,
-            int* iwork, int liwork, int& info)
+            int* iwork, int liwork, int info)
 {
     // call the fortran routine
     chegvd_(&itype, &jobz, &uplo, &n,
@@ -269,7 +269,7 @@ void hegvd(const int itype, const char jobz, const char uplo, const int n,
             std::complex<double>* a, const int lda,
             std::complex<double>* b, const int ldb, double* w,
             std::complex<double>* work, int lwork, double* rwork, int lrwork,
-            int* iwork, int liwork, int& info)
+            int* iwork, int liwork, int info)
 {
     // call the fortran routine
     zhegvd_(&itype, &jobz, &uplo, &n,

source/source_hsolver/diago_ppcg.cpp

@@ -40,7 +40,7 @@ Real max_generalized_residual(
 template <typename T>
 struct PpcgLapack
 {
-    using Real = typename ct::kernels::lapack_hegvd<T, base_device::DEVICE_CPU>::Real;
+    using Real = typename ct::kernels::lapack_hegvx<T, base_device::DEVICE_CPU>::Real;
 
     static void heevd(int n, T* a, Real* w)
     {
@@ -50,8 +50,8 @@ struct PpcgLapack
     static void hegvd(int n, T* a, T* b, Real* w)
     {
         std::vector<T> eigen_vec(n * n, T(0));
-        ct::kernels::lapack_hegvd<T, base_device::DEVICE_CPU>()(
-            n, n, a, b, w, eigen_vec.data());
+        ct::kernels::lapack_hegvx<T, base_device::DEVICE_CPU>()(
+            n, n, a, b, n, w, eigen_vec.data());
         std::copy(eigen_vec.begin(), eigen_vec.end(), a);
     }