HydroBase should initialize excision mask before MoL_PostStep since it is used in MoL_PostStep

Create issue
Issue #1310 closed
Christian Reisswig created an issue

Some routines in GRHydro query whether a point is excised using hydrobase's excision mask. HydroBase initializes the excision mask to zero in certain groups. However, this initialization should happen before any other code reads the excision mask. Since GRHydro reads the excision mask in Con2Prim, which is scheduled in MoL_PostStep, the mask should be set before MoL_PostStep!

Most notably, this is required in post_recover_variables (since the mask is not checkpointed and thus we critically rely on a correct mask initialzation), and also in post_regrid, where new points must be correctly initialized.

I attach a simple patch that fixes this issue.

Keyword: HydroBase

Comments (7)

  1. Frank Löffler
    • changed status to open
    • removed comment

    I wonder why this patch is necessary. It is correct, and should be applied, but doesn't Cactus schedule things sorted alphabetically anymore? HydroBase should have been scheduled before MoL.

  2. Log in to comment