Merge branch 'fix_exact_riemann_pstar_test' into finish_riemann_cleaning

Util/exact_riemann/README.md

@@ -12,3 +12,14 @@ https://ui.adsabs.harvard.edu/abs/2015ApJS..216...31Z/abstract
 and more details are given there.
 
 To build the solver, simply type 'make' in this directory.
+
+
+Notes:
+
+* ``test2`` is a vacuum
+
+  The velocity in the star region should be 0, but it seems that
+  slight differences in the rarefaction integration from the left and
+  from the right give different roundoff, resulting in a small, finite
+  velocity.  This can be controlled a bit by tightening the tolerance
+  in the rarefaction integration.

Util/exact_riemann/riemann_star_state.H

@@ -142,13 +142,11 @@ riemann_star_state(const amrex::Real rho_l, const amrex::Real u_l, const amrex::
 
         amrex::Real pstar_new = pstar - Z_l * Z_r * (ustar_r - ustar_l) / (Z_l + Z_r);
 
-        // estimate the error in the current star solution
-        amrex::Real err1 = std::abs(ustar_r - ustar_l);
-        amrex::Real err2 = pstar_new - pstar;
+        // estimate the error in the current pstar solution
 
-        std::cout << "iter = " << iter << " " << err1 / std::max(std::abs(ustar_l), std::abs(ustar_r))  << " " << err2 / pstar << std::endl;
+        amrex::Real error = std::abs(pstar_new - pstar);
 
-        if (err2 < tol * pstar) {
+        if (error < tol * pstar) {
             converged = true;
         }