src/sh.c @@ -371,7 +371,7 @@ ts_int i,j; for(i=0;i<sph->l;i++){ for(j=0;j<2*i+1;j++){ sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j]; sph->sumUlm2[i][j]+=sph->ulm[i][j]* sph->ulm[i][j]; } } sph->N++;