diff options
author | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-14 15:37:35 +0000 |
---|---|---|
committer | rhaas <rhaas@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-14 15:37:35 +0000 |
commit | a8d0fcb877d375f4544460cec1bf8be6bccb30fa (patch) | |
tree | 9dab0c5aa34ac5e631f7f4e35669014afa0c0ee1 /param.ccl | |
parent | 9f3b48b486078faf82a04b2a39af415be7cdec87 (diff) |
GRHydro: correct atmosphere handling after primitve reconstruction
I have fixed a bug in GRHydro's atmosphere handling after reconstruction.
After reconstruction, the reconstructed plus and minus face values are tested
for whether they drop below atmosphere level and are then reset. In
particular, this means that plus and minus face values for cell i don't
generally coincide with the minus/plus face values of the corresponding cells
i-1 and i+1. If the reconstruction previously found that the Riemann problem
was trivial between, say rhominus(i+1) and rhoplus(i), when changing
rhoplus(i), it is not anymore! However, the mask "trivial_rp" testing for a
trivial Riemann problem is not changed! This leads to inconsistent behavior at
the boundary of the atmosphere.
The symptom was a drift in the center of mass of a static and perfectly
symmetric TOV star.
Patch by Christian Reisswig.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@332 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'param.ccl')
0 files changed, 0 insertions, 0 deletions