Before revision:   To prevent spurious oscillations caused by irregular intersections of mesh cells, a penalization ensuring the stability of the approach and defining implicitly the extension to host domains is added.
After revision:  Thereby, an expensive remeshing and adaptation of the sparse matrix data structure are avoided and the computations are accelerated. To prevent spurious oscillations caused by irregular intersections of mesh cells, a penalization ensuring the stability of the approach and defining implicitly the extension to host domains is added.