#1311: SetMask_SphericalSurface sets mask incorrectly for multiple spheres ----------------------------------------+----------------------------------- Reporter: reisswig@… | Owner: Type: defect | Status: new Priority: blocker | Milestone: ET_2013_05 Component: EinsteinToolkit thorn | Version: development version Keywords: SetMask_SphericalSurface | ----------------------------------------+----------------------------------- SetMask_SphericalSurface loops over all spheres and sets points to "excised" within the sphere radius, and to "not excised" outside. Unfortunately, it clears the mask for each sphere individually, thus overwriting "excised" points from a previous sphere!
The attached patch fixes this by first clearing the mask everywhere before looping over spheres, and then setting points to "excised" if inside a given sphere.