when we are working with an irregular domain in solution of PDEs using finite difference method, can we divide the irregular domain into several rectangular sub-domains and then solve the problem?
the boundary condition is known and when the firs sub-domain is solved, the obtained results are employed as boundary condition for the next sub-domain
This idea is similar to the numerical trapezoidal integration.