Linear solver now works with all boundaries, under all conditions (flat, not flat, internal, surface)