Change default matsolvers to colamd

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

Change default matsolvers to colamd #271

Closed kburns closed 9 months ago

kburns commented 10 months ago

The COLAMD ordering seems to be much better at reducing fill-in than NATURAL, particularly for multi-dim. coupled problems. I'm not aware of cases with substantial performance degradations, but haven't looked thoroughly yet...