[Commits] [svn:einsteintoolkit] IllinoisGRMHD/trunk/src/ (Rev. 8)
zachetie at gmail.com
zachetie at gmail.com
Mon Oct 7 11:25:07 CDT 2013
User: zetienne
Date: 2013/10/07 11:25 AM
Modified:
/trunk/src/
driver_conserv_to_prims.C
Log:
Add code description & license info in header to Conserv-to-Prims solver
File Changes:
Directory: /trunk/src/
======================
File [modified]: driver_conserv_to_prims.C
Delta lines: +27 -3
===================================================================
--- trunk/src/driver_conserv_to_prims.C 2013-10-07 15:24:04 UTC (rev 7)
+++ trunk/src/driver_conserv_to_prims.C 2013-10-07 16:25:07 UTC (rev 8)
@@ -1,6 +1,30 @@
-// Here is a wrapper for the 2d solver described in Noble et al.
-// It is meant to be an interface with their code.
-// Note that it assumes a simple gamma law for the moment. No hybrid.
+/* We evolve forward in time a set of functions called the
+ * "conservative variables" forward in time, and any time
+ * the conserv's are updated, we must solve for the
+ * primitive variables (rho, pressure, velocities) using a
+ * Newton-Raphson technique, before reconstructing &
+ * evaluating the RHSs of the MHD equations again.
+ *
+ * This file contains the driver routine for this Newton-
+ * Raphson solver. Truncation errors in conservative
+ * variables can lead to no solutions in primitive
+ * variables. We correct for these errors in the code
+ * through a number of tricks described in the appendices
+ * of http://arxiv.org/pdf/1112.0568.pdf.
+ *
+ * This is a wrapper for the 2d solver of Noble et al. See
+ * harm_utoprim_2d.c for references and copyright notice
+ * for that solver. This wrapper was primarily written by
+ * Zachariah Etienne & Yuk Tung Liu, in 2011-2013.
+ *
+ * For optimal compatibility, this wrapper is licensed under
+ * the GPL v2 or any later version.
+ *
+ * Note that this code assumes a simple gamma law for the
+ * moment, though it would be easy to extend to a piecewise
+ * polytrope. */
+
+
#include <iostream>
#include <iomanip>
#include <fstream>
More information about the Commits
mailing list