@@ -194,21 +194,33 @@ void unkOverlap_lcao::init(const UnitCell& ucell,
194194 }
195195 }
196196
197- for (auto & co1: center2_orb11)
198- for (auto & co2: co1.second )
199- for (auto & co3: co2.second )
200- for (auto & co4: co3.second )
201- for (auto & co5: co4.second )
202- for (auto & co6: co5.second )
197+ for (auto & co1: center2_orb11) {
198+ for (auto & co2: co1.second ) {
199+ for (auto & co3: co2.second ) {
200+ for (auto & co4: co3.second ) {
201+ for (auto & co5: co4.second ) {
202+ for (auto & co6: co5.second ) {
203203 co6.second .init_radial_table ();
204+ }
205+ }
206+ }
207+ }
208+ }
209+ }
204210
205- for (auto & co1: center2_orb21_r)
206- for (auto & co2: co1.second )
207- for (auto & co3: co2.second )
208- for (auto & co4: co3.second )
209- for (auto & co5: co4.second )
210- for (auto & co6: co5.second )
211+ for (auto & co1: center2_orb21_r) {
212+ for (auto & co2: co1.second ) {
213+ for (auto & co3: co2.second ) {
214+ for (auto & co4: co3.second ) {
215+ for (auto & co5: co4.second ) {
216+ for (auto & co6: co5.second ) {
211217 co6.second .init_radial_table ();
218+ }
219+ }
220+ }
221+ }
222+ }
223+ }
212224
213225 rcut_orb_.resize (orb.get_ntype ());
214226 for (int it = 0 ; it < orb.get_ntype (); ++it) {
@@ -439,8 +451,9 @@ void unkOverlap_lcao::cal_orb_overlap(const UnitCell& ucell)
439451 // if ( !pv.in_this_processor(iw1,iw2) ) continue;
440452
441453 // iw1 and iw2 never have overlap
442- if (orb1_orb2_R[iw1][iw2].empty ())
454+ if (orb1_orb2_R[iw1][iw2].empty ()) {
443455 continue ;
456+ }
444457
445458 int atomType1 = iw2it (ucell,iw1);
446459 int ia1 = iw2ia (ucell,iw1);
0 commit comments