I have a sparse, square, symmetric matrix with the following structure: (Let's say the size of the matrix is $N \times N$)
Here, the area under the blue stripes is the non-zero elements. Could someone tell me if there is an algorithm to invert this kind of matrix that is simple yet more efficient than Gaussian elimination and LU decomposition? Thank you in advance.