{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,4,22]],"date-time":"2026-04-22T08:12:05Z","timestamp":1776845525236,"version":"3.51.2"},"reference-count":30,"publisher":"American Mathematical Society (AMS)","issue":"300","license":[{"start":{"date-parts":[[2016,10,15]],"date-time":"2016-10-15T00:00:00Z","timestamp":1476489600000},"content-version":"am","delay-in-days":366,"URL":"https:\/\/www.ams.org\/publications\/copyright-and-permissions"}],"funder":[{"DOI":"10.13039\/501100003407","name":"Ministero dell\u00e2\u0080\u0099Istruzione, dell\u00e2\u0080\u0099Universit\u00c3 e della Ricerca","doi-asserted-by":"publisher","award":["PRIN 2012 N. 2012MTE38N"],"award-info":[{"award-number":["PRIN 2012 N. 2012MTE38N"]}],"id":[{"id":"10.13039\/501100003407","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100003407","name":"Ministero dell\u00e2\u0080\u0099Istruzione, dell\u00e2\u0080\u0099Universit\u00c3 e della Ricerca","doi-asserted-by":"publisher","award":["PRIN 2012 N. 2012MTE38N"],"award-info":[{"award-number":["PRIN 2012 N. 2012MTE38N"]}],"id":[{"id":"10.13039\/501100003407","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/501100003130","name":"Fonds Wetenschappelijk Onderzoek","doi-asserted-by":"publisher","award":["PRIN 2012 N. 2012MTE38N"],"award-info":[{"award-number":["PRIN 2012 N. 2012MTE38N"]}],"id":[{"id":"10.13039\/501100003130","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Math. Comp."],"abstract":"

\n We consider a linear full elliptic second order partial differential equation in a\n \n \n \n d<\/mml:mi>\n d<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n -dimensional domain,\n \n \n \n \n d<\/mml:mi>\n \n \u2265\n \n <\/mml:mo>\n 1<\/mml:mn>\n <\/mml:mrow>\n d\\ge 1<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n , approximated by isogeometric collocation methods based on uniform (tensor-product) B-splines of degrees\n \n \n \n \n p<\/mml:mi>\n :=<\/mml:mo>\n (<\/mml:mo>\n \n p<\/mml:mi>\n 1<\/mml:mn>\n <\/mml:msub>\n ,<\/mml:mo>\n \n \u2026\n \n <\/mml:mo>\n ,<\/mml:mo>\n \n p<\/mml:mi>\n d<\/mml:mi>\n <\/mml:msub>\n )<\/mml:mo>\n <\/mml:mrow>\n \\boldsymbol {p}:=(p_1,\\ldots ,p_d)<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n ,\n \n \n \n \n \n p<\/mml:mi>\n j<\/mml:mi>\n <\/mml:msub>\n \n \u2265\n \n <\/mml:mo>\n 2<\/mml:mn>\n <\/mml:mrow>\n p_j\\ge 2<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n ,\n \n \n \n \n j<\/mml:mi>\n =<\/mml:mo>\n 1<\/mml:mn>\n ,<\/mml:mo>\n \n \u2026\n \n <\/mml:mo>\n ,<\/mml:mo>\n d<\/mml:mi>\n <\/mml:mrow>\n j=1,\\ldots ,d<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n . We give a construction of the inherently non-symmetric matrices arising from this approximation technique and we perform an analysis of their spectral properties. In particular, we find and study the associated (spectral) symbol, that is, the function describing their asymptotic spectral distribution (in the Weyl sense) when the matrix-size tends to infinity or, equivalently, the fineness parameters tend to zero. The symbol is a non-negative function with a unique zero of order two at\n \n \n \n \n \n \u03b8\n \n <\/mml:mi>\n =<\/mml:mo>\n \n 0<\/mml:mn>\n <\/mml:mrow>\n <\/mml:mrow>\n \\boldsymbol {\\theta }=\\mathbf {0}<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n (where\n \n \n \n \n \n \u03b8\n \n <\/mml:mi>\n :=<\/mml:mo>\n (<\/mml:mo>\n \n \n \u03b8\n \n <\/mml:mi>\n 1<\/mml:mn>\n <\/mml:msub>\n ,<\/mml:mo>\n \n \u2026\n \n <\/mml:mo>\n ,<\/mml:mo>\n \n \n \u03b8\n \n <\/mml:mi>\n d<\/mml:mi>\n <\/mml:msub>\n )<\/mml:mo>\n <\/mml:mrow>\n \\boldsymbol {\\theta }:=(\\theta _1,\\ldots ,\\theta _d)<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n are the Fourier variables), but with infinitely many \u2018numerical zeros\u2019 for large\n \n \n \n \n \n \u2016\n \n <\/mml:mo>\n p<\/mml:mi>\n \n \n \u2016\n \n <\/mml:mo>\n \n \u221e\n \n <\/mml:mi>\n <\/mml:msub>\n <\/mml:mrow>\n \\|\\boldsymbol {p}\\|_\\infty<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n . Indeed, the symbol converges exponentially to zero with respect to\n \n \n \n \n p<\/mml:mi>\n j<\/mml:mi>\n <\/mml:msub>\n p_j<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n at all the points\n \n \n \n \n \u03b8\n \n <\/mml:mi>\n \\boldsymbol {\\theta }<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n such that\n \n \n \n \n \n \n \u03b8\n \n <\/mml:mi>\n j<\/mml:mi>\n <\/mml:msub>\n =<\/mml:mo>\n \n \u03c0\n \n <\/mml:mi>\n <\/mml:mrow>\n \\theta _j=\\pi<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n . In other words, if\n \n \n \n \n p<\/mml:mi>\n j<\/mml:mi>\n <\/mml:msub>\n p_j<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n is large, all the points\n \n \n \n \n \u03b8\n \n <\/mml:mi>\n \\boldsymbol {\\theta }<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n with\n \n \n \n \n \n \n \u03b8\n \n <\/mml:mi>\n j<\/mml:mi>\n <\/mml:msub>\n =<\/mml:mo>\n \n \u03c0\n \n <\/mml:mi>\n <\/mml:mrow>\n \\theta _j=\\pi<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n behave numerically like a zero of the symbol. The presence of the zero of order two at\n \n \n \n \n \n \u03b8\n \n <\/mml:mi>\n =<\/mml:mo>\n \n 0<\/mml:mn>\n <\/mml:mrow>\n <\/mml:mrow>\n \\boldsymbol {\\theta }=\\mathbf {0}<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n is expected because it is intrinsic in any local approximation method, such as finite differences and finite elements, of second order differential operators. However, the \u2018numerical zeros\u2019 lead to the surprising fact that, for large\n \n \n \n \n \n \u2016\n \n <\/mml:mo>\n p<\/mml:mi>\n \n \n \u2016\n \n <\/mml:mo>\n \n \u221e\n \n <\/mml:mi>\n <\/mml:msub>\n <\/mml:mrow>\n \\|\\boldsymbol {p}\\|_\\infty<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n , there is a subspace of high frequencies where the collocation matrices are ill-conditioned. This non-canonical feature is responsible for the slowdown, with respect to\n \n \n \n p<\/mml:mi>\n \\boldsymbol {p}<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n , of standard iterative methods. On the other hand, this knowledge and the knowledge of other properties of the symbol can be exploited to construct iterative solvers with convergence properties independent of the fineness parameters and of the degrees\n \n \n \n p<\/mml:mi>\n \\boldsymbol {p}<\/mml:annotation>\n <\/mml:semantics>\n <\/mml:math>\n <\/inline-formula>\n .\n <\/p>","DOI":"10.1090\/mcom\/3027","type":"journal-article","created":{"date-parts":[[2015,10,15]],"date-time":"2015-10-15T13:32:59Z","timestamp":1444915979000},"page":"1639-1680","source":"Crossref","is-referenced-by-count":29,"title":["Spectral analysis and spectral symbol of matrices in isogeometric collocation methods"],"prefix":"10.1090","volume":"85","author":[{"given":"Marco","family":"Donatelli","sequence":"first","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Carlo","family":"Garoni","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Carla","family":"Manni","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Stefano","family":"Serra-Capizzano","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Hendrik","family":"Speleers","sequence":"additional","affiliation":[],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"14","published-online":{"date-parts":[[2015,10,15]]},"reference":[{"issue":"1","key":"1","doi-asserted-by":"publisher","first-page":"186","DOI":"10.1137\/S0895479803421987","article-title":"V-cycle optimal convergence for certain (multilevel) structured linear systems","volume":"26","author":"Aric\u00f2, Antonio","year":"2004","journal-title":"SIAM J. Matrix Anal. Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0895-4798","issn-type":"print"},{"issue":"11","key":"2","doi-asserted-by":"publisher","first-page":"2075","DOI":"10.1142\/S0218202510004878","article-title":"Isogeometric collocation methods","volume":"20","author":"Auricchio, F.","year":"2010","journal-title":"Math. Models Methods Appl. Sci.","ISSN":"https:\/\/id.crossref.org\/issn\/0218-2025","issn-type":"print"},{"issue":"1","key":"3","doi-asserted-by":"publisher","first-page":"300","DOI":"10.1137\/S0036142999363188","article-title":"Superlinear convergence of conjugate gradients","volume":"39","author":"Beckermann, Bernhard","year":"2001","journal-title":"SIAM J. Numer. Anal.","ISSN":"https:\/\/id.crossref.org\/issn\/0036-1429","issn-type":"print"},{"issue":"2","key":"4","doi-asserted-by":"publisher","first-page":"746","DOI":"10.1137\/05063533X","article-title":"On the asymptotic spectrum of finite element matrix sequences","volume":"45","author":"Beckermann, Bernhard","year":"2007","journal-title":"SIAM J. Numer. Anal.","ISSN":"https:\/\/id.crossref.org\/issn\/0036-1429","issn-type":"print"},{"key":"5","series-title":"Graduate Texts in Mathematics","isbn-type":"print","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4612-0653-8","volume-title":"Matrix analysis","volume":"169","author":"Bhatia, Rajendra","year":"1997","ISBN":"https:\/\/id.crossref.org\/isbn\/0387948465"},{"key":"6","unstructured":"D. Bini, M. Capovani, O. Menchi. Metodi Numerici per l\u2019Algebra Lineare. Zanichelli (1988)."},{"key":"7","series-title":"Applied Mathematical Sciences","isbn-type":"print","volume-title":"A practical guide to splines","volume":"27","author":"de Boor, Carl","year":"2001","ISBN":"https:\/\/id.crossref.org\/isbn\/0387953663"},{"key":"8","series-title":"Universitext","isbn-type":"print","doi-asserted-by":"publisher","DOI":"10.1007\/978-1-4612-1426-7","volume-title":"Introduction to large truncated Toeplitz matrices","author":"B\u00f6ttcher, Albrecht","year":"1999","ISBN":"https:\/\/id.crossref.org\/isbn\/0387985700"},{"key":"9","series-title":"Wavelet Analysis and its Applications","isbn-type":"print","volume-title":"An introduction to wavelets","volume":"1","author":"Chui, Charles K.","year":"1992","ISBN":"https:\/\/id.crossref.org\/isbn\/0121745848"},{"key":"10","doi-asserted-by":"crossref","unstructured":"J. A. Cottrell, T. J. R. Hughes, Y. Bazilevs. Isogeometric Analysis: Toward Integration of CAD and FEA. John Wiley & Sons (2009).","DOI":"10.1002\/9780470749081"},{"key":"11","doi-asserted-by":"publisher","first-page":"230","DOI":"10.1016\/j.cma.2014.06.001","article-title":"Robust and optimal multi-iterative techniques for IgA Galerkin linear systems","volume":"284","author":"Donatelli, Marco","year":"2015","journal-title":"Comput. Methods Appl. Mech. Engrg.","ISSN":"https:\/\/id.crossref.org\/issn\/0045-7825","issn-type":"print"},{"key":"12","doi-asserted-by":"publisher","first-page":"1120","DOI":"10.1016\/j.cma.2014.11.036","article-title":"Robust and optimal multi-iterative techniques for IgA collocation linear systems","volume":"284","author":"Donatelli, Marco","year":"2015","journal-title":"Comput. Methods Appl. Mech. Engrg.","ISSN":"https:\/\/id.crossref.org\/issn\/0045-7825","issn-type":"print"},{"key":"13","unstructured":"M. Donatelli, C. Garoni, C. Manni, S. Serra-Capizzano, H. Speleers. Symbol-based multigrid methods for Galerkin B-spline isogeometric analysis. Tech. Report TW650 (2014), Dept. Computer Science, KU Leuven."},{"key":"14","unstructured":"C. Garoni. Structured matrices coming from PDE approximation theory: spectral analysis, spectral symbol and design of fast iterative solvers. Ph.D. Thesis in Mathematics of Computation, University of Insubria, Como, Italy (2014). Available online at \\url{http:\/\/hdl.handle.net\/10277\/568}."},{"issue":"4","key":"15","doi-asserted-by":"publisher","first-page":"751","DOI":"10.1007\/s00211-013-0600-2","article-title":"On the spectrum of stiffness matrices arising from isogeometric analysis","volume":"127","author":"Garoni, Carlo","year":"2014","journal-title":"Numer. Math.","ISSN":"https:\/\/id.crossref.org\/issn\/0029-599X","issn-type":"print"},{"issue":"2","key":"16","doi-asserted-by":"publisher","first-page":"213","DOI":"10.1007\/s00020-014-2157-6","article-title":"Tools for determining the asymptotic spectral distribution of non-Hermitian perturbations of Hermitian matrix-sequences and applications","volume":"81","author":"Garoni, Carlo","year":"2015","journal-title":"Integral Equations Operator Theory","ISSN":"https:\/\/id.crossref.org\/issn\/0378-620X","issn-type":"print"},{"issue":"1","key":"17","doi-asserted-by":"publisher","first-page":"84","DOI":"10.1016\/j.jat.2006.05.002","article-title":"The asymptotic properties of the spectrum of nonsymmetrically perturbed Jacobi matrix sequences","volume":"144","author":"Golinskii, Leonid","year":"2007","journal-title":"J. Approx. Theory","ISSN":"https:\/\/id.crossref.org\/issn\/0021-9045","issn-type":"print"},{"issue":"5-8","key":"18","doi-asserted-by":"publisher","first-page":"301","DOI":"10.1016\/j.cma.2008.12.004","article-title":"Efficient quadrature for NURBS-based isogeometric analysis","volume":"199","author":"Hughes, T. J. R.","year":"2010","journal-title":"Comput. Methods Appl. Mech. Engrg.","ISSN":"https:\/\/id.crossref.org\/issn\/0045-7825","issn-type":"print"},{"issue":"3","key":"19","doi-asserted-by":"publisher","first-page":"939","DOI":"10.1137\/S0036142995285034","article-title":"Preconditioning Chebyshev spectral collocation by finite-difference operators","volume":"34","author":"Kim, Sang Dong","year":"1997","journal-title":"SIAM J. Numer. Anal.","ISSN":"https:\/\/id.crossref.org\/issn\/0036-1429","issn-type":"print"},{"issue":"1","key":"20","doi-asserted-by":"publisher","first-page":"3","DOI":"10.1137\/S0036144504445376","article-title":"Convergence analysis of Krylov subspace iterations with methods from potential theory","volume":"48","author":"Kuijlaars, Arno B. J.","year":"2006","journal-title":"SIAM Rev.","ISSN":"https:\/\/id.crossref.org\/issn\/1095-7200","issn-type":"print"},{"issue":"1","key":"21","doi-asserted-by":"publisher","first-page":"142","DOI":"10.1006\/jath.2001.3617","article-title":"Asymptotic zero distribution of orthogonal polynomials with discontinuously varying recurrence coefficients","volume":"113","author":"Kuijlaars, A. B. J.","year":"2001","journal-title":"J. Approx. Theory","ISSN":"https:\/\/id.crossref.org\/issn\/0021-9045","issn-type":"print"},{"key":"22","doi-asserted-by":"publisher","first-page":"170","DOI":"10.1016\/j.cma.2013.07.017","article-title":"Isogeometric collocation: cost comparison with Galerkin methods and extension to adaptive hierarchical NURBS discretizations","volume":"267","author":"Schillinger, Dominik","year":"2013","journal-title":"Comput. Methods Appl. Mech. Engrg.","ISSN":"https:\/\/id.crossref.org\/issn\/0045-7825","issn-type":"print"},{"key":"23","doi-asserted-by":"crossref","first-page":"109","DOI":"10.1016\/S0024-3795(97)00231-0","article-title":"On the extreme eigenvalues of Hermitian (block) Toeplitz matrices","volume":"270","author":"Serra, Stefano","year":"1998","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"},{"issue":"3","key":"24","doi-asserted-by":"publisher","first-page":"461","DOI":"10.1007\/s002110050400","article-title":"The rate of convergence of Toeplitz based PCG methods for second order nonlinear boundary value problems","volume":"81","author":"Serra, Stefano","year":"1999","journal-title":"Numer. Math.","ISSN":"https:\/\/id.crossref.org\/issn\/0029-599X","issn-type":"print"},{"issue":"1-3","key":"25","doi-asserted-by":"publisher","first-page":"121","DOI":"10.1016\/S0024-3795(00)00311-6","article-title":"Distribution results on the algebra generated by Toeplitz sequences: a finite-dimensional approach","volume":"328","author":"Serra Capizzano, S.","year":"2001","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"},{"key":"26","doi-asserted-by":"publisher","first-page":"371","DOI":"10.1016\/S0024-3795(02)00504-9","article-title":"Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations","volume":"366","author":"Serra Capizzano, S.","year":"2003","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"},{"issue":"1","key":"27","doi-asserted-by":"publisher","first-page":"180","DOI":"10.1016\/j.laa.2006.04.012","article-title":"The GLT class as a generalized Fourier analysis and applications","volume":"419","author":"Serra-Capizzano, Stefano","year":"2006","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"},{"key":"28","doi-asserted-by":"publisher","first-page":"41","DOI":"10.1016\/S0024-3795(02)00719-X","article-title":"Analysis of preconditioning strategies for collocation linear systems","volume":"369","author":"Serra Capizzano, Stefano","year":"2003","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"},{"issue":"6","key":"29","doi-asserted-by":"publisher","first-page":"1962","DOI":"10.1137\/S0036142997322722","article-title":"Spectral analysis of Hermite cubic spline collocation systems","volume":"36","author":"Sun, Weiwei","year":"1999","journal-title":"SIAM J. Numer. Anal.","ISSN":"https:\/\/id.crossref.org\/issn\/0036-1429","issn-type":"print"},{"issue":"1-3","key":"30","doi-asserted-by":"publisher","first-page":"91","DOI":"10.1016\/S0024-3795(97)10079-9","article-title":"Locally Toeplitz sequences: spectral properties and applications","volume":"278","author":"Tilli, P.","year":"1998","journal-title":"Linear Algebra Appl.","ISSN":"https:\/\/id.crossref.org\/issn\/0024-3795","issn-type":"print"}],"container-title":["Mathematics of Computation"],"original-title":[],"language":"en","link":[{"URL":"http:\/\/www.ams.org\/mcom\/2016-85-300\/S0025-5718-2015-03027-0\/S0025-5718-2015-03027-0.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"},{"URL":"https:\/\/www.ams.org\/mcom\/2016-85-300\/S0025-5718-2015-03027-0\/S0025-5718-2015-03027-0.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,4,21]],"date-time":"2026-04-21T18:49:14Z","timestamp":1776797354000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.ams.org\/mcom\/2016-85-300\/S0025-5718-2015-03027-0\/"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2015,10,15]]},"references-count":30,"journal-issue":{"issue":"300","published-print":{"date-parts":[[2016,7]]}},"alternative-id":["S0025-5718-2015-03027-0"],"URL":"https:\/\/doi.org\/10.1090\/mcom\/3027","archive":["CLOCKSS","Portico"],"relation":{},"ISSN":["1088-6842","0025-5718"],"issn-type":[{"value":"1088-6842","type":"electronic"},{"value":"0025-5718","type":"print"}],"subject":[],"published":{"date-parts":[[2015,10,15]]}}}