Commits

Sam Skillman committed 99d895b Merge

Merged in gbryan/enzo-dev-gb (pull request #124)

Correction in the computation of the gravitational potential RHS for the MultiGrid solver

  • Participants
  • Parent commits 3ef8b77, 30751bf

Comments (0)

Files changed (1)

File src/enzo/Grid_SolveForPotential.C

   float InverseVolumeElement = 1;
   for (dim = 0; dim < GridRank; dim++) {
     size *= GravitatingMassFieldDimension[dim];
-    InverseVolumeElement *= (GravitatingMassFieldDimension[dim]);
+    InverseVolumeElement *= (GravitatingMassFieldDimension[dim]-1);
   }
   tol_dim = max(sqrt(float(size))*1e-6, tol_dim);