[ET Trac] [Einstein Toolkit] #1841: GRHydro tracers broken

Einstein Toolkit trac-noreply at einsteintoolkit.org
Mon Nov 23 10:20:55 CST 2015


#1841: GRHydro tracers broken
------------------------------------+---------------------------------------
  Reporter:  I.Hawke@…              |       Owner:                     
      Type:  defect                 |      Status:  new                
  Priority:  unset                  |   Milestone:                     
 Component:  EinsteinToolkit thorn  |     Version:  development version
Resolution:                         |    Keywords:  GRHydro            
------------------------------------+---------------------------------------

Comment (by I.Hawke@…):

 This is not the only problem with tracers, as I've since found.

 Currently all reconstruction of tracers is done on the primitive
 variables, but the conserved are never updated. This can be fixed by
 adding

   if (evolve_tracer .ne. 0) then
     call Prim2ConservativeTracer(CCTK_PASS_FTOF)
   end if

 at the end of GRHydro_Reconstruct.F90, but *only* if the use_cxx_code
 parameter is set to no.

 To truly fix tracers, the C++ version GRHydro_Reconstruct.cc will also
 need a P2C_Tracer routine, which needs calling at the end.

-- 
Ticket URL: <https://trac.einsteintoolkit.org/ticket/1841#comment:1>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit


More information about the Trac mailing list