Efficient mass and stiffness matrix assembly via weighted Gaussian quadrature rules for Bsplines
Access Status
Authors
Date
2017Collection
Type
Metadata
Show full item recordAbstract
Calabro et al. (2017) changed the paradigm of the mass and stiffness computation from the traditional elementwise assembly to a rowwise concept, showing that the latter one offers integration that may be orders of magnitude faster. Considering a Bspline basis function as a nonnegative measure, each mass matrix row is integrated by its own quadrature rule with respect to that measure. Each rule is easy to compute as it leads to a linear system of equations, however, the quadrature rules are of the NewtonCotes type, that is, they require a number of quadrature points that is equal to the dimension of the spline space. In this work, we propose weighted quadrature rules of Gaussian type which require the minimum number of quadrature points while guaranteeing exactness of integration with respect to the weight function. The weighted Gaussian rules arise as solutions of nonlinear systems of equations. We derive rules for the mass and stiffness matrices for uniform $C^1$ quadratic and $C^2$ cubic isogeometric discretizations. Our rules further reduce the number of quadrature points by a factor of $(\frac{p+1}{2p+1})^d$ when compared to Calabro et al. (2017), $p$ being the polynomial degree and $d$ the dimension of the problem, and consequently reduce the computational cost of the mass and stiffness matrix assembly by a similar factor.
Citation
Source Title
Department
Related items
Showing items related by title, author, creator and subject.

Barton, M.; Calo, Victor (2016)We introduce optimal quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. Using the homotopy continuation concept (Barton and Calo, 2016) that ...

Barton, M.; Calo, Victor (2017)© 2016 Elsevier LtdWe introduce Gaussian quadrature rules for spline spaces that are frequently used in Galerkin discretizations to build mass and stiffness matrices. By definition, these spaces are of even degrees. The ...

Barton, M.; Calo, Victor (2016)We introduce a new concept for generating optimal quadrature rules for splines. To generate an optimal quadrature rule in a given (target) spline space, we build an associated source space with known optimal quadrature ...