search

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

Linalg improvements #276

Closed kburns closed 8 months ago

kburns commented 8 months ago

This PR includes some improvements to linear algebra using cython routines.

kburns commented 8 months ago

For 2D RBC at (Nx, Nz) = (256, 512), this speeds things up by about 5%, and for 3D HOT at 256^3 it speeds things up by about 4%. Bigger differences on other more anisotropic simulations where inverse Jacobi added a lot.

kburns commented 8 months ago

The old behavior, if desired, is available through the OLD_CSR_MATVECS config option.