Abstract
We present a first step towards the spectral analysis of matrices arising from IgA Galerkin methods based on hyperbolic and trigonometric GB-splines. Second order differential problems with constant coefficients are considered and discretized by means of sequences of both nested and non-nested spline spaces. We prove that there always exists an asymptotic eigenvalue distribution which can be compactly described by a symbol, just like in the polynomial case. There is a complete similarity between the symbol expressions in the hyperbolic, trigonometric and polynomial cases. This results in similar spectral features of the corresponding matrices. We also analyze the IgA discretization based on trigonometric GB-splines for the eigenvalue problem related to the univariate Laplace operator. We prove that, for non-nested spaces, the phase parameter can be exploited to improve the spectral approximation with respect to the polynomial case. As part of the analysis, we derive several Fourier properties of cardinal GB-splines.
Similar content being viewed by others
Notes
The functions \(N_{i,1}^{U,V}\) may also depend on p because of the definition of \(\widetilde{U}_i, \widetilde{V}_i\), but we omit the parameter p in order to avoid a heavier notation. Moreover, in the most interesting cases (polynomial, hyperbolic and trigonometric GB-splines) the functions \( N_{i, 1}^{U,V}\) are independent of p.
For polynomial B-splines, a degree p in Galerkin approximation corresponds to a degree \(2p+1\) in collocation, see [11, Remark 3.2]. We have followed the same rule here.
A different ordering has been used in [14] in the bivariate setting.
References
Aricó, A., Donatelli, M., Serra-Capizzano, S.: V-cycle optimal convergence for certain (multilevel) structured linear systems. SIAM J. Matrix Anal. Appl. 26, 186–214 (2004)
Beckermann, B., Kuijlaars, A.B.J.: Superlinear convergence of conjugate gradients. SIAM J. Numer. Anal. 39, 300–329 (2001)
Beirão da Veiga, L., Buffa, A., Rivas, J., Sangalli, G.: Some estimates for \(h\)–\(p\)–\(k\)-refinement in isogeometric analysis. Numer. Math. 118, 271–305 (2011)
Carnicer, J.M., Mainar, E., Peña, J.M.: Critical length for design purposes and extended Chebyshev spaces. Constr. Approx. 20, 55–71 (2004)
Chui, C.K.: An Introduction to Wavelets. Academic Press, New York (1992)
Costantini, P., Lyche, T., Manni, C.: On a class of weak Tchebycheff systems. Numer. Math. 101, 333–354 (2005)
Cottrell, J.A., Hughes, T.J.R., Bazilevs, Y.: Isogeometric Analysis: Toward Integration of CAD and FEA. Wiley, New York (2009)
Cottrell, J.A., Reali, A., Bazilevs, Y., Hughes, T.J.R.: Isogeometric analysis of structural vibrations. Comput. Methods Appl. Mech. Eng. 195, 5257–5296 (2006)
Donatelli, M., Garoni, C., Manni, C., Serra-Capizzano, S., Speleers, H.: Robust and optimal multi-iterative techniques for IgA Galerkin linear systems. Comput. Methods Appl. Mech. Eng. 284, 230–264 (2015)
Donatelli, M., Garoni, C., Manni, C., Serra-Capizzano, S., Speleers, H.: Robust and optimal multi-iterative techniques for IgA collocation linear systems. Comput. Methods Appl. Mech. Eng. 284, 1120–1146 (2015)
Donatelli, M., Garoni, C., Manni, C., Serra-Capizzano, S., Speleers, H.: Spectral analysis and spectral symbol of matrices in isogeometric collocation methods. Math. Comput. (to appear)
Donatelli, M., Garoni, C., Manni, C., Serra-Capizzano, S., Speleers, H.: Symbol-based multigrid methods for Galerkin B-spline isogeometric analysis. In: Technical Report TW650, Department of Computer Science, KU Leuven (2014)
Garoni, C., Hughes, T.J.R., Reali, A., Serra-Capizzano, S., Speleers, H.: Smoothness versus polynomial degree: why IgA outperforms FEA in the spectral approximation (in preparation)
Garoni, C., Manni, C., Pelosi, F., Serra-Capizzano, S., Speleers, H.: On the spectrum of stiffness matrices arising from isogeometric analysis. Numer. Math. 127, 751–799 (2014)
Garoni, C., Manni, C., Serra-Capizzano, S., Sesana, D., Speleers, H.: Spectral analysis and spectral symbol of matrices in isogeometric Galerkin methods. Math. Comput. (to appear)
Garoni, C., Serra-Capizzano, S.: The theory of generalized locally Toeplitz sequences: a review, an extension, and a few representative applications. In: Technical Report 2015–023, Department of Information Technology, Uppsala University (2015)
Golinskii, L., Serra-Capizzano, S.: The asymptotic properties of the spectrum of nonsymmetrically perturbed Jacobi matrix sequences. J. Approx. Theory 144, 84–102 (2007)
Grenander, U., Szegö, G.: Toeplitz Forms and Their Applications, 2nd edn. Chelsea, New York (1984)
Hughes, T.J.R., Cottrell, J.A., Bazilevs, Y.: Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comput. Methods Appl. Mech. Eng. 194, 4135–4195 (2005)
Hughes, T.J.R., Evans, J.A., Reali, A.: Finite element and NURBS approximations of eigenvalue, boundary-value, and initial-value problems. Comput. Methods Appl. Mech. Eng. 272, 290–320 (2014)
Kuijlaars, A.B.J.: Convergence analysis of Krylov subspace iterations with methods from potential theory. SIAM Rev. 48, 3–40 (2006)
Kvasov, B.I., Sattayatham, P.: GB-splines of arbitrary order. J. Comput. Appl. Math. 104, 63–88 (1999)
Manni, C., Pelosi, F., Sampoli, M.L.: Generalized B-splines as a tool in isogeometric analysis. Comput. Methods Appl. Mech. Eng. 200, 867–881 (2011)
Manni, C., Pelosi, F., Sampoli, M.L.: Isogeometric analysis in advection–diffusion problems: tension splines approximation. J. Comput. Appl. Math. 236, 511–528 (2011)
Manni, C., Reali, A., Speleers, H.: Isogeometric collocation methods with generalized B-splines. Comput. Math. Appl. 70, 1659–1675 (2015)
Serra, S.: Multi-iterative methods. Comput. Math. Appl. 26, 65–87 (1993)
Serra-Capizzano, S.: The rate of convergence of Toeplitz based PCG methods for second order nonlinear boundary value problems. Numer. Math. 81, 461–495 (1999)
Serra-Capizzano, S.: Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations. Linear Algebra Appl. 366, 371–402 (2003)
Serra-Capizzano, S.: The GLT class as a generalized Fourier analysis and applications. Linear Algebra Appl. 419, 180–233 (2006)
Serra-Capizzano, S., Tablino Possio, C.: Analysis of preconditioning strategies for collocation linear systems. Linear Algebra Appl. 369, 41–75 (2003)
Serra-Capizzano, S., Tilli, P.: On unitarily invariant norms of matrix-valued linear positive operators. J. Taylor & Francis Group Inequal. Appl. 7, 309–330 (2002)
Tilli, P.: A note on the spectral distribution of Toeplitz matrices. Linear Multilinear Algebra 45, 147–159 (1998)
Tyrtyshnikov, E.E.: A unifying approach to some old and new theorems on distribution and clustering. Linear Algebra Appl. 232, 1–43 (1996)
Unser, M., Blu, T.: Cardinal exponential splines: Part I—Theory and filtering algorithms. IEEE Trans. Signal Proc. 53, 1425–1438 (2005)
Wang, G., Fang, M.: Unified and extended form of three types of splines. J. Comput. Appl. Math. 216, 498–508 (2008)
Acknowledgments
This work was partially supported by INdAM-GNCS Gruppo Nazionale per il Calcolo Scientifico, by the MIUR ‘Futuro in Ricerca 2013’ Programme through the project DREAMS, and by the ‘Uncovering Excellence’ Programme of the University of Rome ‘Tor Vergata’ through the project DEXTEROUS.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Roman, F., Manni, C. & Speleers, H. Spectral analysis of matrices in Galerkin methods based on generalized B-splines with high smoothness. Numer. Math. 135, 169–216 (2017). https://doi.org/10.1007/s00211-016-0796-z
Received:
Revised:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00211-016-0796-z