Abstract. We construct and analyze an overlapping Schwarz preconditioner for elliptic problems discretized with isogeometric analysis. The preconditioner is based on partitioning the domain of the problem into overlapping subdomains, solving local isogeometric problems on these subdomains, and solving an additional coarse isogeometric problem associated with the subdomain mesh. We develop an h-analysis of the preconditioner, showing in particular that the resulting algorithm is scalable and its convergence rate depends linearly on the ratio between subdomain and "overlap sizes" for fixed polynomial degree p and regularity k of the basis functions. Numerical results in two-and three-dimensional tests show the good convergence properties of the preconditioner with respect to the isogeometric discretization parameters h, p, k, number of subdomains N , overlap size, and also jumps in the coefficients of the elliptic operator.Key words. domain decomposition methods, overlapping Schwarz, scalable preconditioners, isogeometric analysis, finite elements, NURBS AMS subject classifications. 65N55, 65N30, 65F10 DOI. 10.1137/110833476 1. Introduction. Isogeometric analysis (IGA) based on NURBS (nonuniform rational B-splines) was introduced by Hughes, Cottrell, and Bazilevs in [27] as an innovative numerical methodology for the analysis of PDE problems, allowing for an exact description of CAD-type geometries. NURBS are a standard in the computeraided design (CAD) community mainly because such spline functions allow excellent representations of free-form surfaces, and there are very efficient algorithms for their evaluation, refinement, and derefinement. In the isogeometric framework, the NURBS basis functions representing the CAD geometry are also used as the basis for the discrete solution space of PDEs, following an isoparametric paradigm. In addition to exact representation of CAD geometries, another advantage of using NURBS basis functions is the higher control on the regularity of the discrete space. For instance, spaces of global C k regularity are easily built, thus allowing for fewer degrees of freedom, better performance in case of vibrations, easier approximation of higher order problems, and other advantages. IGA methodologies have been summarized in the recent book [18] and studied in, e.g., [2,4,9,10,19,23,28,29,11,5,8]. IGA methods are having a growing impact on fields as diverse as fluid dynamics [6,7,40,15,26], structural mechanics [3,1,12,20,30,39], and electromagnetics [17,16].The discrete systems produced by isogeometric methods are better conditioned than the systems produced by standard finite elements or finite differences, but their conditioning can still degenerates rapidly for decreasing mesh size h, increasing poly-