A low-rank solver for conforming multipatch Isogeometric Analysis

Abstract

In this paper, we propose an innovative isogeometric low-rank solver for the linear elasticity model problem, specifically designed to allow multipatch domains. Our approach splits the domain into subdomains, each formed by the union of neighboring patches. Within each subdomain, we employ Tucker low-rank matrices and vectors to approximate the system matrices and right-hand side vectors, respectively. This enables the construction of local approximate fast solvers. These local solvers are then combined into an overlapping Schwarz preconditioner, which is utilized in a truncated preconditioned conjugate gradient method. Numerical experiments demonstrate the significant memory storage benefits and a uniformly bounded number of iterations with respect to both mesh size and spline degree.