The cost of continuity: Performance of iterative solvers on isogeometric finite elements
In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional...
| Main Authors: | , , , |
|---|---|
| Format: | Journal Article |
| Published: |
2013
|
| Online Access: | http://hdl.handle.net/20.500.11937/51388 |
| _version_ | 1848758685580394496 |
|---|---|
| author | Collier, N. Dalcin, L. Pardo, D. Calo, Victor |
| author_facet | Collier, N. Dalcin, L. Pardo, D. Calo, Victor |
| author_sort | Collier, N. |
| building | Curtin Institutional Repository |
| collection | Online Access |
| description | In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional finite element spaces, and Cp-1 B-splines, which represent maximum continuity We provide theoretical estimates for the increase in cost of the matrix-vector product as well as for the construction and application of black-box preconditioners. We accompany these estimates with numerical results and study their sensitivity to various grid parameters such as element size h and polynomial order of approximation p in addition to the aforementioned continuity of the basis. Finally, we present timing results for a range of preconditioning options for the Laplace problem. We conclude that the matrix-vector product operation is at most 33p2/8 times more expensive for the more continuous space, although for moderately low p, this number is significantly reduced. Moreover, if static condensation is not employed, this number further reduces to at most a value of 8, even for high p. Preconditioning options can be up to p3 times more expensive to set up, although this difference significantly decreases for some popular preconditioners such as incomplete LU factorization. © 2013 Society for Industrial and Applied Mathematics. |
| first_indexed | 2025-11-14T09:47:55Z |
| format | Journal Article |
| id | curtin-20.500.11937-51388 |
| institution | Curtin University Malaysia |
| institution_category | Local University |
| last_indexed | 2025-11-14T09:47:55Z |
| publishDate | 2013 |
| recordtype | eprints |
| repository_type | Digital Repository |
| spelling | curtin-20.500.11937-513882018-03-29T09:08:25Z The cost of continuity: Performance of iterative solvers on isogeometric finite elements Collier, N. Dalcin, L. Pardo, D. Calo, Victor In this paper we study how the use of a more continuous set of basis functions affects the cost of solving systems of linear equations resulting from a discretized Galerkin weak form. Specifically, we compare performance of linear solvers when discretizing using Co B-splines, which span traditional finite element spaces, and Cp-1 B-splines, which represent maximum continuity We provide theoretical estimates for the increase in cost of the matrix-vector product as well as for the construction and application of black-box preconditioners. We accompany these estimates with numerical results and study their sensitivity to various grid parameters such as element size h and polynomial order of approximation p in addition to the aforementioned continuity of the basis. Finally, we present timing results for a range of preconditioning options for the Laplace problem. We conclude that the matrix-vector product operation is at most 33p2/8 times more expensive for the more continuous space, although for moderately low p, this number is significantly reduced. Moreover, if static condensation is not employed, this number further reduces to at most a value of 8, even for high p. Preconditioning options can be up to p3 times more expensive to set up, although this difference significantly decreases for some popular preconditioners such as incomplete LU factorization. © 2013 Society for Industrial and Applied Mathematics. 2013 Journal Article http://hdl.handle.net/20.500.11937/51388 10.1137/120881038 restricted |
| spellingShingle | Collier, N. Dalcin, L. Pardo, D. Calo, Victor The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title | The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title_full | The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title_fullStr | The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title_full_unstemmed | The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title_short | The cost of continuity: Performance of iterative solvers on isogeometric finite elements |
| title_sort | cost of continuity: performance of iterative solvers on isogeometric finite elements |
| url | http://hdl.handle.net/20.500.11937/51388 |