I have a pentadiagonal system and I was hoping that there are fast implementations of solvers related to this.
I was trying to use BandedMatrices.jl and then do LU decomposition. But in this case it returns only 1 banded matrix that I need to pivot. If i know specify that I do not want pivoting, then it doesn’t generate banded matrices at all