search

danfortunato / ultraSEM

The ultraspherical spectral element method
MIT License
29 stars 3 forks source link

Implement O(p^4) leaf solver #1

Closed danfortunato closed 5 years ago

danfortunato commented 5 years ago

Building the solution operator in ultraSEMLeaf/initialize requires solving a p^2 x p^2 almost-banded-block-banded linear system (of bandwidth p) with p righthand sides. The current implementation uses backslash, which is O(p^5). We can achieve O(p^4) by using Woodbury.

danfortunato commented 5 years ago

Fixed with b405166304b795b9d32d7f12e6f77637b1683be1.