Iterative solution for layered geometries