fix an illegal call of get_ngk()

maki49 · Sep 1, 2023 · da4aa52 · da4aa52
1 parent 1779e38
commit da4aa52
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/source/module_psi/psi.cpp b/source/module_psi/psi.cpp
@@ -281,6 +281,7 @@ template <typename T, typename Device> int Psi<T, Device>::get_current_nbas() co
 
 template <typename T, typename Device> const int& Psi<T, Device>::get_ngk(const int ik_in) const
 {
+    if (!this->ngk) return this->nbasis;
     return this->ngk[ik_in];
 }