@@ -14,7 +14,7 @@ namespace LR_Util
1414 auto dr_real = DMR_real.get_DMR_vector ()[is];
1515 assert (dr != nullptr );
1616 assert (dr_real != nullptr );
17- for (int ia = 0 ;ia < nat;ia++)
17+ for (int ia = 0 ;ia < nat;ia++) {
1818 for (int ja = 0 ;ja < nat;ja++)
1919 {
2020 auto ap = dr->find_pair (ia, ja);
@@ -25,9 +25,11 @@ namespace LR_Util
2525 auto dR = ap->get_R_index (iR);
2626 auto ptr = ap->get_HR_values (iR).get_pointer ();
2727 auto ptr_real = ap_real->get_HR_values (dR.x , dR.y , dR.z ).get_pointer ();
28- for (int i = 0 ;i < ap->get_size ();++i) ptr_real[i] = (get_imag ? ptr[i].imag () : ptr[i].real ());
28+ for (int i = 0 ;i < ap->get_size ();++i) { ptr_real[i] = (get_imag ? ptr[i].imag () : ptr[i].real ());
29+ }
2930 }
3031 }
32+ }
3133 }
3234 }
3335
@@ -37,7 +39,7 @@ namespace LR_Util
3739 const char & type)
3840 {
3941 bool get_imag = (type == ' I' || type == ' i' );
40- for (int ia = 0 ;ia < nat;ia++)
42+ for (int ia = 0 ;ia < nat;ia++) {
4143 for (int ja = 0 ;ja < nat;ja++)
4244 {
4345 auto ap = HR.find_pair (ia, ja);
@@ -48,8 +50,10 @@ namespace LR_Util
4850 auto dR = ap->get_R_index (iR);
4951 auto ptr = ap->get_HR_values (iR).get_pointer ();
5052 auto ptr_real = ap_real->get_HR_values (dR.x , dR.y , dR.z ).get_pointer ();
51- for (int i = 0 ;i < ap->get_size ();++i) get_imag ? ptr[i].imag (ptr_real[i]) : ptr[i].real (ptr_real[i]);
53+ for (int i = 0 ;i < ap->get_size ();++i) { get_imag ? ptr[i].imag (ptr_real[i]) : ptr[i].real (ptr_real[i]);
54+ }
5255 }
5356 }
57+ }
5458 }
5559}
0 commit comments