Measuring Gauss-Seidel loop-carried dependency and fixing it via loop unrolling(loiseaujc.github.io) |
Measuring Gauss-Seidel loop-carried dependency and fixing it via loop unrolling(loiseaujc.github.io) |
Written by Claude AI. Which is a problem because Claude text is riddled with engagement bait and filler, making it exhausting to read.
Maybe pedagogy was the point? I didn’t really get that from the article, maybe some intention was lost by filtering it through AI.
(On the other hand, I skimmed the article and this might be a banded system, which LAPACK can handle).
Even when one isn’t defining named constants for various kinds, it’s still good practice to use them on literal floating constants, since those default to the default kind of REAL in Fortran, not to double precision as in most later languages. A very pernicious pitfall of the language is to (say) assign or initialize X=0.1 instead of X=0.1_8 when X is double precision.