Skip to main content
Log in

Spectral analysis of matrices in Galerkin methods based on generalized B-splines with high smoothness

  • Published:
Numerische Mathematik Aims and scope Submit manuscript

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.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8

Similar content being viewed by others

Notes

  1. 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.

  2. 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.

  3. A different ordering has been used in [14] in the bivariate setting.

References

  1. 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)

    Article  MathSciNet  MATH  Google Scholar 

  2. Beckermann, B., Kuijlaars, A.B.J.: Superlinear convergence of conjugate gradients. SIAM J. Numer. Anal. 39, 300–329 (2001)

    Article  MathSciNet  MATH  Google Scholar 

  3. 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)

    Article  MathSciNet  MATH  Google Scholar 

  4. Carnicer, J.M., Mainar, E., Peña, J.M.: Critical length for design purposes and extended Chebyshev spaces. Constr. Approx. 20, 55–71 (2004)

    Article  MathSciNet  MATH  Google Scholar 

  5. Chui, C.K.: An Introduction to Wavelets. Academic Press, New York (1992)

    MATH  Google Scholar 

  6. Costantini, P., Lyche, T., Manni, C.: On a class of weak Tchebycheff systems. Numer. Math. 101, 333–354 (2005)

    Article  MathSciNet  MATH  Google Scholar 

  7. Cottrell, J.A., Hughes, T.J.R., Bazilevs, Y.: Isogeometric Analysis: Toward Integration of CAD and FEA. Wiley, New York (2009)

    Book  Google Scholar 

  8. 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)

    Article  MathSciNet  MATH  Google Scholar 

  9. 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)

    Article  MathSciNet  Google Scholar 

  10. 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)

    Article  MathSciNet  Google Scholar 

  11. 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)

  12. 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)

  13. 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)

  14. 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)

    Article  MathSciNet  MATH  Google Scholar 

  15. 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)

  16. 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)

  17. Golinskii, L., Serra-Capizzano, S.: The asymptotic properties of the spectrum of nonsymmetrically perturbed Jacobi matrix sequences. J. Approx. Theory 144, 84–102 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  18. Grenander, U., Szegö, G.: Toeplitz Forms and Their Applications, 2nd edn. Chelsea, New York (1984)

    MATH  Google Scholar 

  19. 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)

    Article  MathSciNet  MATH  Google Scholar 

  20. 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)

    Article  MathSciNet  MATH  Google Scholar 

  21. Kuijlaars, A.B.J.: Convergence analysis of Krylov subspace iterations with methods from potential theory. SIAM Rev. 48, 3–40 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  22. Kvasov, B.I., Sattayatham, P.: GB-splines of arbitrary order. J. Comput. Appl. Math. 104, 63–88 (1999)

    Article  MathSciNet  MATH  Google Scholar 

  23. 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)

    Article  MathSciNet  MATH  Google Scholar 

  24. Manni, C., Pelosi, F., Sampoli, M.L.: Isogeometric analysis in advection–diffusion problems: tension splines approximation. J. Comput. Appl. Math. 236, 511–528 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  25. Manni, C., Reali, A., Speleers, H.: Isogeometric collocation methods with generalized B-splines. Comput. Math. Appl. 70, 1659–1675 (2015)

    Article  MathSciNet  Google Scholar 

  26. Serra, S.: Multi-iterative methods. Comput. Math. Appl. 26, 65–87 (1993)

    Article  MathSciNet  MATH  Google Scholar 

  27. 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)

    Article  MathSciNet  MATH  Google Scholar 

  28. Serra-Capizzano, S.: Generalized locally Toeplitz sequences: spectral analysis and applications to discretized partial differential equations. Linear Algebra Appl. 366, 371–402 (2003)

    Article  MathSciNet  MATH  Google Scholar 

  29. Serra-Capizzano, S.: The GLT class as a generalized Fourier analysis and applications. Linear Algebra Appl. 419, 180–233 (2006)

    Article  MathSciNet  MATH  Google Scholar 

  30. Serra-Capizzano, S., Tablino Possio, C.: Analysis of preconditioning strategies for collocation linear systems. Linear Algebra Appl. 369, 41–75 (2003)

    Article  MathSciNet  MATH  Google Scholar 

  31. 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)

    Article  MathSciNet  MATH  Google Scholar 

  32. Tilli, P.: A note on the spectral distribution of Toeplitz matrices. Linear Multilinear Algebra 45, 147–159 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  33. Tyrtyshnikov, E.E.: A unifying approach to some old and new theorems on distribution and clustering. Linear Algebra Appl. 232, 1–43 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  34. Unser, M., Blu, T.: Cardinal exponential splines: Part I—Theory and filtering algorithms. IEEE Trans. Signal Proc. 53, 1425–1438 (2005)

    Article  MathSciNet  Google Scholar 

  35. Wang, G., Fang, M.: Unified and extended form of three types of splines. J. Comput. Appl. Math. 216, 498–508 (2008)

    Article  MathSciNet  MATH  Google Scholar 

Download references

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

Authors

Corresponding author

Correspondence to Hendrik Speleers.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

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

Download citation

  • Received:

  • Revised:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00211-016-0796-z

Mathematics Subject Classification

Navigation