Accidentally left squared norm in places.

EXP-code · Aug 19, 2023 · 18d4acc · 18d4acc
1 parent 05360b6
commit 18d4acc
Showing 1 changed file with 9 additions and 9 deletions.
diff --git a/src/CBrockDisk.cc b/src/CBrockDisk.cc
@@ -593,13 +593,13 @@ void CBrockDisk::determine_fields_at_point_polar
 
     get_dens_coefs(l,expcoef->row(2*l - 1), pc);
     get_dens_coefs(l,expcoef->row(2*l    ), ps);
-    dens += (pc*cosm[0][l] + ps*sinm[0][l]) * M_SQRT2 / sqnorm(l,n);
+    dens += (pc*cosm[0][l] + ps*sinm[0][l]) * M_SQRT2;
 
     get_pot_coefs(l,expcoef->row(2*l - 1), pc, dpc);
     get_pot_coefs(l,expcoef->row(2*l    ), ps, dps);
-    potl += (pc *cosm[0][l] + ps *sinm[0][l]) * M_SQRT2 / sqnorm(l,n);
-    potr += (dpc*cosm[0][l] + dps*sinm[0][l]) * M_SQRT2 / sqnorm(l,n);
-    potp += (-pc*sinm[0][l] + ps *cosm[0][l]) * M_SQRT2 * l / sqnorm(l,n);
+    potl += (pc *cosm[0][l] + ps *sinm[0][l]) * M_SQRT2;
+    potr += (dpc*cosm[0][l] + dps*sinm[0][l]) * M_SQRT2;
+    potp += (-pc*sinm[0][l] + ps *cosm[0][l]) * M_SQRT2 * l;
   }
 
   *tdens0 /= scale*scale*scale;
@@ -620,8 +620,8 @@ void CBrockDisk::get_pot_coefs(int l, const Eigen::VectorXd& coef,
   pp = dpp = 0.0;
 
   for (int i=0; i<nmax; i++) {
-    pp  += potd[0](l, i) * coef[i] / sqnorm(l,n);
-    dpp += dpot[0](l, i) * coef[i] / sqnorm(l,n);
+    pp  += potd[0](l, i) * coef[i] / sqnorm(l,i);
+    dpp += dpot[0](l, i) * coef[i] / sqnorm(l,i);
   }
 
   p = -pp;
@@ -638,8 +638,8 @@ void CBrockDisk::get_pot_coefs_safe
   pp = dpp = 0.0;
 
   for (int i=0; i<nmax; i++) {
-    pp  += potd1(l, i) * coef[i] / sqnorm(l,n);
-    dpp += dpot1(l, i) * coef[i] / sqnorm(l,n);
+    pp  += potd1(l, i) * coef[i] / sqnorm(l,i);
+    dpp += dpot1(l, i) * coef[i] / sqnorm(l,i);
   }
 
   p  = -pp;
@@ -653,7 +653,7 @@ void CBrockDisk::get_dens_coefs
   double pp = 0.0;
 
   for (int i=0; i<nmax; i++)
-    pp  += dend(l, i) * coef[i] / sqnorm(l,n);
+    pp  += dend(l, i) * coef[i] / sqnorm(l,i);
 
   p = pp;
 }