
For more information about solving both elliptic and parabolic problems in irregular domains using the finite difference method in conjunction with level set functions, we refer the reader to [1]}, [2]} (and references therein). Once the finite difference stencil has been formed, the resulting system of linear equations is solved using LU decomposition.
