aboutsummaryrefslogtreecommitdiff
path: root/src/driver/BH_diagnostics.cc
diff options
context:
space:
mode:
authorjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2003-07-04 13:10:22 +0000
committerjthorn <jthorn@f88db872-0e4f-0410-b76b-b9085cfa78c5>2003-07-04 13:10:22 +0000
commit8e5402c746dfa30f616cd38aff00c2d5a83e6e1e (patch)
treef9936ad424a930dac4f5d88460d32529bb56bd9d /src/driver/BH_diagnostics.cc
parent2a4324cff58b30936a9b2760512438423bd0e0cf (diff)
at Ed's request, add 2 new columns to the BH_diagnostics file,
giving the polar/equatorial circumference ratios git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/AHFinderDirect/trunk@1115 f88db872-0e4f-0410-b76b-b9085cfa78c5
Diffstat (limited to 'src/driver/BH_diagnostics.cc')
-rw-r--r--src/driver/BH_diagnostics.cc17
1 files changed, 11 insertions, 6 deletions
diff --git a/src/driver/BH_diagnostics.cc b/src/driver/BH_diagnostics.cc
index cbe4b59..6d262a6 100644
--- a/src/driver/BH_diagnostics.cc
+++ b/src/driver/BH_diagnostics.cc
@@ -368,8 +368,10 @@ fprintf(fileptr, "# column 14 = max z\n");
fprintf(fileptr, "# column 15 = xy-plane circumference\n");
fprintf(fileptr, "# column 16 = xz-plane circumference\n");
fprintf(fileptr, "# column 17 = yz-plane circumference\n");
-fprintf(fileptr, "# column 18 = area\n");
-fprintf(fileptr, "# column 19 = m_irreducible\n");
+fprintf(fileptr, "# column 18 = ratio of xz/xy-plane circumferences\n");
+fprintf(fileptr, "# column 19 = ratio of yz/xy-plane circumferences\n");
+fprintf(fileptr, "# column 20 = area\n");
+fprintf(fileptr, "# column 21 = m_irreducible\n");
fflush(fileptr);
return fileptr;
@@ -410,13 +412,16 @@ fprintf(fileptr,
double(min_z), double(max_z));
fprintf(fileptr,
- // {xy,xz,yz}-plane area m_irreducible
- // circumferences ====== ======
- // ====================== ====== ======
- "%#.10g\t%#.10g\t%#.10g\t%#.10g\t%#.10g\n",
+ // {xy,xz,yz}-plane xz/xy yz/xy area m_irreducible
+ // circumferences circumferende ====== ======
+ // ratios
+ // ====================== ============== ====== ======
+ "%#.10g\t%#.10g\t%#.10g\t%#.10g\t%#.10g\t%#.10g\t%#.10g\n",
double(circumference_xy),
double(circumference_xz),
double(circumference_yz),
+ double(circumference_xz / circumference_xy),
+ double(circumference_yz / circumference_xy),
double(area), double(m_irreducible));
fflush(fileptr);