Abstract
It has been observed that the task of matrix assembly in Isogeometric Analysis (IGA) is more challenging than in the case of traditional finite element methods. The additional difficulties associated with IGA are caused by the increased degree and the larger supports of the functions that occur in the integrals defining the matrix elements. Recently we introduced an interpolation-based approach that approximately transforms the integrands into piecewise polynomials and uses look-up tables to evaluate their integrals [17, 18]. The present paper relies on this earlier work and proposes to use tensor methods to accelerate the assembly process further. More precisely, we show how to represent the matrices that occur in IGA as sums of a small number of Kronecker products of auxiliary matrices that are defined by univariate integrals. This representation, which is based on a low-rank tensor approximation of certain parts of the integrands, makes it possible to achieve a significant speedup of the assembly process without compromising the overall accuracy of the simulation.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Similar content being viewed by others
Notes
- 1.
Geometry + Simulation Modules, see www.gs.jku.at/gismo, also [10].
References
Antolin, P., Buffa, A., Calabrò, F., Martinelli, M., Sangalli, G.: Efficient matrix computation for tensor-product isogeometric analysis: the use of sum factorization. Comp. Meth. Appl. Mech. Engrg. 285, 817–828 (2015)
Auricchio, F., Calabrò, F., Hughes, T., Reali, A., Sangalli, G.: A simple algorithm for obtaining nearly optimal quadrature rules for NURBS-based isogeometric analysis. Comp. Meth. Appl. Mech. Engrg. 249–252, 15–27 (2012)
Bazilevs, Y., Beirão da Veiga, L., Cottrell, J.A., Hughes, T.J.R., Sangalli, G.: Isogeometric analysis: approximation, stability and error estimates for \(h\)-refined meshes. Math. Models Methods Appl. Sci. 16(7), 1031–1090 (2006)
Calabrò, F., Manni, C., Pitolli, F.: Computation of quadrature rules for integration with respect to refinable functions on assigned nodes. Appl. Numer. Math. 90, 168–189 (2015)
Cottrell, J.A., Hughes, T.J.R., Bazilevs, Y.: Isogeometric Analysis: Toward Integration of CAD and FEA. Wiley, Chichester (2009)
De Boor, C.: A Practical Guide to Splines. Applied Mathematical Sciences. Springer, Berlin (2001)
Hackbusch, W.: Tensor Spaces and Numerical Tensor Calculus. Springer, Berlin (2012)
Hughes, T., Cottrell, J., Bazilevs, Y.: Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Comp. Meth. Appl. Mech. Engrg. 194(39–41), 4135–4195 (2005)
Hughes, T., Reali, A., Sangalli, G.: Efficient quadrature for NURBS-based isogeometric analysis. Comp. Meth. Appl. Mech. Engrg. 199(5–8), 301–313 (2010)
Jüttler, B., Langer, U., Mantzaflaris, A., Moore, S.E., Zulehner, W.: Geometry + simulation modules: implementing isogeometric analysis. PAMM 14(1), 961–962 (2014)
Karatarakis, A., Karakitsios, P., Papadrakakis, M.: Computation of the isogeometric analysis stiffness matrix on GPU. In: Papadrakakis, M., Kojic, M., Tuncer, I., (eds.) Proceedings of the 3rd South-East European Conference on Computational Mechanics (SEECCM) (2013). www.eccomasproceedings.org/cs2013
Karatarakis, A., Karakitsios, P., Papadrakakis, M.: GPU accelerated computation of the isogeometric analysis stiffness matrix. Comp. Meth. Appl. Mech. Engrg. 269, 334–355 (2014)
Khoromskij, B.N.: \({\cal O}(d \log n)\)-quantics approximation of \(N\)-\(d\) tensors in high-dimensional numerical modeling. Constr. Appr. 34(2), 257–280 (2011)
Khoromskij, B.N.: Tensor-structured numerical methods in scientific computing: survey on recent advances. Chemometr. Intell. Lab. Syst. 110(1), 1–19 (2012)
Khoromskij, B.N.: Tensor numerical methods for multidimensional PDEs: theoretical analysis and initial applications. In: Proceedings of ESAIM, pp. 1–28 (2014)
Kleiss, S., Pechstein, C., Jüttler, B., Tomar, S.: IETI - isogeometric tearing and interconnecting. Comp. Meth. Appl. Mech. Engrg. 247–248, 201–215 (2012)
Mantzaflaris, A., Jüttler, B.: Exploring matrix generation strategies in isogeometric analysis. In: Floater, M., Lyche, T., Mazure, M.-L., Mørken, K., Schumaker, L.L. (eds.) MMCS 2012. LNCS, vol. 8177, pp. 364–382. Springer, Heidelberg (2014)
Mantzaflaris, A., Jüttler, B.: Integration by interpolation and look-up for Galerkin-based isogeometric analysis. Comp. Methods Appl. Mech. Engrg. 284, 373–400 (2015). Isogeometric Analysis Special Issue
Patzák, B., Rypl, D.: Study of computational efficiency of numerical quadrature schemes in the isogeometric analysis. In: Proceedings of the 18 th International Conference on Engineering Mechanics, EM 2012, pp. 1135–1143 (2012)
Schillinger, D., Evans, J., Reali, A., Scott, M., Hughes, T.: Isogeometric collocation: cost comparison with Galerkin methods and extension to adaptive hierarchical NURBS discretizations. Comp. Meth. Appl. Mech. Engrg. 267, 170–232 (2013)
Schillinger, D., Hossain, S., Hughes, T.: Reduced Bézier element quadrature rules for quadratic and cubic splines in isogeometric analysis. Comp. Meth. Appl. Mech. Engrg. 277, 1–45 (2014)
Acknowledgement
This research was supported by the National Research Network “Geometry + Simulation” (NFN S117), funded by the Austrian Science Fund (FWF).
Author information
Authors and Affiliations
Corresponding author
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2015 Springer International Publishing Switzerland
About this paper
Cite this paper
Mantzaflaris, A., Jüttler, B., Khoromskij, B.N., Langer, U. (2015). Matrix Generation in Isogeometric Analysis by Low Rank Tensor Approximation. In: Boissonnat, JD., et al. Curves and Surfaces. Curves and Surfaces 2014. Lecture Notes in Computer Science(), vol 9213. Springer, Cham. https://doi.org/10.1007/978-3-319-22804-4_24
Download citation
DOI: https://doi.org/10.1007/978-3-319-22804-4_24
Published:
Publisher Name: Springer, Cham
Print ISBN: 978-3-319-22803-7
Online ISBN: 978-3-319-22804-4
eBook Packages: Computer ScienceComputer Science (R0)