BUG: Previous idea was wrong - boundary faces are 0 anyhow.

pmgbergen · Aug 27, 2024 · dcd3000 · dcd3000
1 parent 4ea2b34
commit dcd3000
Showing 1 changed file with 2 additions and 6 deletions.
diff --git a/src/darsia/utils/fv.py b/src/darsia/utils/fv.py
@@ -71,13 +71,9 @@ def __init__(
                 # After all, this is identical to the arithmetic avg of cell volumes
                 # connected to the face. This requires careful handling of the boundary faces.
                 # NOTE: This assumes equidistant grid spacing.
+                # NOTE: Fluxes on boundary are 0 and excluded.
                 volume_scaling = np.ones(grid.num_faces, dtype=float)
-                volume_scaling[grid.exterior_faces] = 0.5
-                mass_matrix = sps.diags(
-                    np.prod(grid.voxel_size)
-                    * volume_scaling
-                    * np.ones(grid.num_faces, dtype=float)
-                )
+                mass_matrix = sps.diags(np.prod(grid.voxel_size) * volume_scaling)
             else:
                 raise NotImplementedError