#956: Allow skipping setting the RHS to zero in MoL ------------------------------------+--------------------------------------- Reporter: eschnett | Owner: Type: enhancement | Status: closed Priority: minor | Milestone: Component: EinsteinToolkit thorn | Version: Resolution: fixed | Keywords: ------------------------------------+---------------------------------------
Comment (by rhaas):
Thinking about this a bit more, what happens in the following situation: * Carpet poisons new memory incl. the RHS when it is first allocated * most thorns put nothing into the ghost zones when computing the RHS, so they stay poisoned * in the buffer zones, until the full RK step is done, the poison accumulates possibly creating an overflow and/or triggering a NaN check inside MoL_Add (if it was turned on)
Is this what would happen?