Open dbischof90 opened 1 year ago
For multidimensional FDM calculations in C#, we can have more or less banded diagonal matrices. Initializing those is very handy in Python's scipy; for example see https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.diags.html
scipy
Building sparse matrices with diagonal values (and vectors) already seems to be implemented. Is it a lot of trouble to implement a more general diagonal storage?
For multidimensional FDM calculations in C#, we can have more or less banded diagonal matrices. Initializing those is very handy in Python's
scipy; for example see https://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.diags.htmlBuilding sparse matrices with diagonal values (and vectors) already seems to be implemented. Is it a lot of trouble to implement a more general diagonal storage?