[ET Trac] [Einstein Toolkit] #623: use GRHydro_polytrope_handle whenever prim2conpolytrope is called

Einstein Toolkit trac-noreply at einsteintoolkit.org
Fri Oct 14 03:49:00 CDT 2011


#623: use GRHydro_polytrope_handle whenever prim2conpolytrope is called
--------------------+-------------------------------------------------------
 Reporter:  rhaas   |       Owner:     
     Type:  defect  |      Status:  new
 Priority:  minor   |   Milestone:     
Component:  Other   |     Version:     
 Keywords:          |  
--------------------+-------------------------------------------------------
 Currently GRHydro uses grhydro_eos_handle in both Primitive2Conservative
 and Primitive2ConservativePoly and relies on grhydro_eos_handle being an
 actual handle to the 2d polytype in the later.
 This means that calling Primitive2ConservativePoly from within a run with
 grhydro_eos_type == "General" uses the wrong handle (namely the general
 eos). This is not normally a problem since Primitive2ConservativePoly is
 only ever called if grhydro_eos_type == "polytype", with the exception of
 EoSChangeGammaK which always calls Primitive2ConservativePoly (and has to
 call a polytype routine).

 There are basically two ways to make this work in this case:
 * use GRHydro_polytrope_handle in Primitive2ConservativePoly
 * call prim2conPoly from within EoSChangeGammaK

 The first option has the disadvantage that this is not what
 Conservative2Primitive does. The second has the disadvantage that
 EoSChangeGammaK has to also recompute Yecons (and any other new
 conservatives we might define) on top of calling prim2conPoly.

 Any comments?

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/623>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit


More information about the Trac mailing list