search

DedalusProject / dedalus

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

Construct LHS system matrix #245

Closed kharazity closed 1 year ago

kharazity commented 1 year ago

Hey everyone,

I wanted to do some numerical testing on the SEM. In particular, I want to study the condition numbers and bandwidth of the produced SE matrices. I wanted to start with a simple example like that in /examples/ldvp_2d_poisson/poisson.py.

First, how do I extract the final LHS system matrix that is used in the linear solve? Does this linear system include the boundary conditions and the Dirichlet recombination steps? What is the best way to efficiently extract this system and estimate it's conditioning?

jsoishi commented 1 year ago

Hi,

The issue tracker is not the best place for questions like these. Please ask this at our mailing list, https://groups.google.com/g/dedalus-users and we would be happy to help!

Jeff