Merge pull request #614 from PrincetonUniversity/schwarzschild_fix

felker · web-flow · commit ad811238e274 · 2024-08-13T22:38:18.000-04:00
Fixed Schwarzschild cell-width calculation
diff --git a/src/coordinates/schwarzschild.cpp b/src/coordinates/schwarzschild.cpp
@@ -354,7 +354,7 @@ void Coordinates::CenterWidth2(const int k, const int j, const int il, const int
 #pragma omp simd
   for (int i=il; i<=iu; ++i) {
     // \Delta W = r \Delta\theta
-    dx2(i) = x1v(i) * dx1f(j);
+    dx2(i) = x1v(i) * dx2f(j);
   }
   return;
 }

Original file line number	Diff line number	Diff line change
`@@ -354,7 +354,7 @@ void Coordinates::CenterWidth2(const int k, const int j, const int il, const int`
`354`	`354`	`#pragma omp simd`
`355`	`355`	`for (int i=il; i<=iu; ++i) {`
`356`	`356`	`// \Delta W = r \Delta\theta`
`357`		`- dx2(i) = x1v(i) * dx1f(j);`
	`357`	`+ dx2(i) = x1v(i) * dx2f(j);`
`358`	`358`	`}`
`359`	`359`	`return;`
`360`	`360`	`}`