@@ -40,13 +40,21 @@ void Output_HContainer<T>::write(int rx_in, int ry_in, int rz_in)
4040 int find_R = 0 ;
4141 for (int iR = 0 ; iR < size_for_loop_R; iR++)
4242 {
43+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
44+
4345 this ->_hcontainer ->loop_R (iR, rx, ry, rz);
4446 if (rx == rx_in && ry == ry_in && rz == rz_in)
4547 {
48+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
49+
4650 find_R += 1 ;
4751 this ->write_single_R (rx, ry, rz);
52+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
53+
4854 break ;
4955 }
56+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
57+
5058 }
5159 if (find_R == 0 )
5260 {
@@ -57,10 +65,18 @@ void Output_HContainer<T>::write(int rx_in, int ry_in, int rz_in)
5765template <typename T>
5866void Output_HContainer<T>::write_single_R(int rx, int ry, int rz)
5967{
68+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
69+
6070 this ->_hcontainer ->fix_R (rx, ry, rz);
71+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
72+
6173 ModuleIO::SparseMatrix<T> sparse_matrix
6274 = ModuleIO::SparseMatrix<T>(_hcontainer->get_nbasis (), _hcontainer->get_nbasis ());
75+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
76+
6377 sparse_matrix.setSparseThreshold (this ->_sparse_threshold );
78+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
79+
6480 for (int iap = 0 ; iap < this ->_hcontainer ->size_atom_pairs (); ++iap)
6581 {
6682 auto atom_pair = this ->_hcontainer ->get_atom_pair (iap);
@@ -77,6 +93,8 @@ void Output_HContainer<T>::write_single_R(int rx, int ry, int rz)
7793 }
7894 }
7995 }
96+ std::cout << __FILE__ << " " << __LINE__ << std::endl;
97+
8098 if (sparse_matrix.getNNZ () != 0 )
8199 {
82100 _ofs << rx << " " << ry << " " << rz << " " << sparse_matrix.getNNZ () << std::endl;
0 commit comments