@@ -304,7 +304,8 @@ void Wavefunc_in_pw::produce_local_basis_in_pw(const UnitCell& ucell,
304304/* for(int is_N = 0; is_N < 2; is_N++)*/ // for rotate base
305305 for (int is_N = 0 ; is_N < 1 ; is_N++)
306306 {
307- if (L==0 && is_N==1 ) continue ;
307+ if (L==0 && is_N==1 ) { continue ;
308+ }
308309 if (ucell.atoms [it].ncpp .has_so )
309310 {
310311 const double j = std::abs (double (L+is_N) - 0.5 );
@@ -332,15 +333,16 @@ void Wavefunc_in_pw::produce_local_basis_in_pw(const UnitCell& ucell,
332333 std::complex <double > fup,fdown;
333334 // int nc;
334335 // This routine creates two functions only in the case j=l+1/2 or exit in the other case
335- if (fabs (j-L+0.5 )<1e-4 ) continue ;
336+ if (fabs (j-L+0.5 )<1e-4 ) { continue ;
337+ }
336338 delete[] chiaux;
337339 chiaux = new double [npw];
338340 // Find the functions j= l- 1/2
339- if (L==0 )
341+ if (L==0 ) {
340342 for (int ig=0 ;ig<npw;ig++){
341343 chiaux[ig] = flq[ig];
342344 }
343- else
345+ } else
344346 {
345347 /* for(int ib = 0;ib < ucell.atoms[it].nchi;ib++)
346348 {
@@ -368,8 +370,9 @@ void Wavefunc_in_pw::produce_local_basis_in_pw(const UnitCell& ucell,
368370 for (int m = 0 ;m<2 *L+1 ;m++)
369371 {
370372 const int lm = L*L +m;
371- if (iwall + 2 * L + 1 > ucell.natomwfc )
373+ if (iwall + 2 * L + 1 > ucell.natomwfc ) {
372374 ModuleBase::WARNING_QUIT (" this->wf.atomic_wfc()" , " error: too many wfcs" );
375+ }
373376 for (int ig = 0 ; ig < npw; ig++)
374377 {
375378 aux[ig] = sk[ig] * ylm (lm,ig) * chiaux[ig];
0 commit comments