Popis: |
In this paper we develop asymptotically optimal algorithms for fast computations with the discrete harmonic Poincare–Steklov operators (Dirichlet–Neumann mapping) for interior and exterior problems in the presence of a nested mesh refinement. Our approach is based on the multilevel interface solver applied to the Schur complement reduction onto the nested refined interface associated with a nonmatching decomposition of a polygon by rectangular substructures. This paper extends methods from Khoromskij and Prossdorf (1995), where the finite element approximations of interior problems on quasi‐uniform grids have been considered. For both interior and exterior problems, the matrix–vector multiplication with the compressed Schur complement matrix on the interface is shown to have a complexity of the order O(Nr log3Nu), where Nr = O((1 + pr) Nu) is the number of degrees of freedom on the polygonal boundary under consideration, Nu is the boundary dimension of a finest quasi‐uniform level and pr ⩾ 0 defines the refinement depth. The corresponding memory needs are estimated by O(Nr logqNu), where q = 2 or q = 3 in the case of interior and exterior problems, respectively. |