[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