@@ -46,8 +46,9 @@ void lapackEigen(int &npw, std::vector<std::complex<double>> &hm, double *e, boo
4646 char tmp_c1 = ' V' , tmp_c2 = ' U' ;
4747 zheev_ (&tmp_c1, &tmp_c2, &npw, hm.data (), &npw, e, work2, &lwork, rwork, &info);
4848 end = clock ();
49- if (outtime)
49+ if (outtime) {
5050 std::cout << " Lapack Run time: " << (double )(end - start) / CLOCKS_PER_SEC << " S" << std::endl;
51+ }
5152 delete[] rwork;
5253 delete[] work2;
5354}
@@ -78,7 +79,8 @@ class DiagoBPCGPrepare
7879 // calculate eigenvalues by LAPACK;
7980 double *e_lapack = new double [npw];
8081 auto ev = DIAGOTEST::hmatrix;
81- if (mypnum == 0 ) lapackEigen (npw, ev, e_lapack, false );
82+ if (mypnum == 0 ) { lapackEigen (npw, ev, e_lapack, false );
83+ }
8284 // initial guess of psi by perturbing lapack psi
8385 ModuleBase::ComplexMatrix psiguess (nband, npw);
8486 std::default_random_engine p (1 );
@@ -248,7 +250,14 @@ TEST(DiagoBPCGTest, readH)
248250 // read Hamilt matrix from file data-H
249251 std::vector<std::complex <double >> hm;
250252 std::ifstream ifs;
251- ifs.open (" H-KPoints-Si64.dat" );
253+ std::string filename = " H-KPoints-Si64.dat" ;
254+ ifs.open (filename);
255+ // open file and check status
256+ if (!ifs.is_open ())
257+ {
258+ std::cout << " Error opening file " << filename << std::endl;
259+ exit (1 );
260+ }
252261 DIAGOTEST::readh (ifs, hm);
253262 ifs.close ();
254263 int dim = DIAGOTEST::npw;
@@ -280,7 +289,8 @@ int main(int argc, char **argv)
280289
281290 testing::InitGoogleTest (&argc, argv);
282291 ::testing::TestEventListeners &listeners = ::testing::UnitTest::GetInstance ()->listeners ();
283- if (myrank != 0 ) delete listeners.Release (listeners.default_result_printer ());
292+ if (myrank != 0 ) { delete listeners.Release (listeners.default_result_printer ());
293+ }
284294
285295 int result = RUN_ALL_TESTS ();
286296 if (myrank == 0 && result != 0 )
0 commit comments