Add experimental PW LOBPCG eigensolver

source/source_estate/elecstate_print.cpp

@@ -58,6 +58,7 @@ void print_scf_iterinfo(const std::string& ks_solver,
            {"scalapack_gvx", "GV"},
            {"cusolver", "CU"},
            {"bpcg", "BP"},
+           {"lobpcg", "LB"},
            {"pexsi", "PE"},
            {"cusolvermp", "CM"},
            {"sdft", "CT"}}; // CT = Chebyshev Trace, for pure SDFT (nbands=0) where no H diagonalization is performed

source/source_hsolver/CMakeLists.txt

@@ -4,6 +4,7 @@ list(APPEND objects
     diago_david.cpp
     diago_dav_subspace.cpp
     diago_bpcg.cpp
+    diago_lobpcg.cpp
     para_linear_transform.cpp
     hsolver_pw.cpp
     hsolver_lcaopw.cpp