I'm working on translating a discontinuous Galerkin code from MATLAB to C and I'm at the final point where I need to solve a sparse system. I've taken a course in C before but I'm very rusty and wasn't all that good to begin with.
I've used LAPACK before but was wondering if there are better options out there. My professor uses PetSc but I was looking at the documentation and have no idea how to install it let alone use it. The "tutorial" section on the ANL site is written at a pretty high level already.
Also, for sparse solvers in C, is there a "special way" that the array has to be stored? Or can I just input my arrays as is into the solver and get the solution?
Here's a plot of the matrix structure