DedalusProject / dedalus

A flexible framework for solving PDEs with modern spectral methods.
http://dedalus-project.org/
GNU General Public License v3.0
489 stars 115 forks source link

Matsolver updates #251

Closed kburns closed 1 year ago

kburns commented 1 year ago

This PR adds some matsolver updates, updating the interface with umfpack and adding a flexible base class for SuperLU factorized solves that can easily be subclassed to customize all options. Included are NATURAL and COLAMD regular and transposed solvers, since those seem to be the most commonly used.

Relatedly, we should compare NATURAL and COLAMD transposed for all the examples, and see if there's any reason not to switch the default to COLAMD transposed.

@bpbrown @lecoanet @kyle-augustson

kburns commented 1 year ago

Merging this in because it's useful. But we should keep in mind the idea of changing the default.