fix bug

haozhihan · haozhihan · commit a3817e4983c0 · 2024-12-26T10:06:41.000Z
diff --git a/source/module_hamilt_pw/hamilt_pwdft/operator_pw/velocity_pw.cpp b/source/module_hamilt_pw/hamilt_pwdft/operator_pw/velocity_pw.cpp
@@ -48,7 +48,7 @@ void Velocity::act
 ) const
 {
     ModuleBase::timer::tick("Operator", "Velocity");
-    const int npw = psi_in->get_current_nbas();
+    const int npw = psi_in->get_ngk(this->ik);
     const int max_npw = psi_in->get_nbasis() / psi_in->npol;
     const int npol = psi_in->npol;
     const std::complex<double>* tmpsi_in = psi0;

Original file line number	Diff line number	Diff line change
`@@ -48,7 +48,7 @@ void Velocity::act`
`48`	`48`	`) const`
`49`	`49`	`{`
`50`	`50`	`ModuleBase::timer::tick("Operator", "Velocity");`
`51`		`- const int npw = psi_in->get_current_nbas();`
	`51`	`+ const int npw = psi_in->get_ngk(this->ik);`
`52`	`52`	`const int max_npw = psi_in->get_nbasis() / psi_in->npol;`
`53`	`53`	`const int npol = psi_in->npol;`
`54`	`54`	`const std::complex<double>* tmpsi_in = psi0;`