alt Hacker NewsLAPACK doesn't do sparse systems.
(On the other hand, I skimmed the article and this might be a banded system, which LAPACK can handle).
LAPACK doesn't do sparse systems.
(On the other hand, I skimmed the article and this might be a banded system, which LAPACK can handle).
True, 2D Poisson is only approximately banded. Still, I don’t think I’d bother rolling my own sparse solver. It’s a well-trod problem.