loiseaujc

3 ms·

1.
by loiseaujc 17d ago
Author here. It is Fortran indeed, but there is nothing specific to Fortran in what I discuss. Had you written the kernel in C or C++, you'd run in exactly the same problems.
2.
by loiseaujc 17d ago
Author here. It actually is a block tridiagonal matrix with tridiagonal blocks. The test problem only has a quarter millions of unknowns so, for a production run, I wouldn't bother write my own sparse solver either. Here, it is done ma
3.
by loiseaujc 17d ago
Author here. Glad you liked it. Claude has been used to help me structure the post, but not the writing itself (except maybe proof-reading it). All the techniques I use are pretty standard tricks to retain someone's attention (in parti
4.
by loiseaujc 17d ago
Author here. The second-order accurate finite difference approximation of the 2D Poisson operator is a block tridiagonal matrix with tridiagonal blocks. There are ways to handle it with LAPACK, but you'll still need special care. And 5
5.
by loiseaujc 17d ago
Author here. Performance-wise, it will never ever be a serious solver. No matter how hard you try. And for this particular matrix, (2D Poisson, block tridiagonal with tridiagonal blocks, quarter millions unknows), LAPACK isn't gonna do
6.
by loiseaujc 17d ago
Author here. Sure it does, and we'll get there! I just need to be done with my teaching duties a bit.
7.
by loiseaujc 17d ago
Author here. It is indeed specific to this 2D Poisson matrix, although loop-unrolling + substitution is a very general technique in performance engineering. Eventually, whenever I'm done teaching classes, we'll get to specialized
8.
by loiseaujc 17d ago
Author here. I did use Claude to help me structure the architecture of the post a bit, but not the text. These are idioms I use all the time as an educator to keep easily-bored students engaged as much as possible. Really sorry that this tr
9.
Measuring Gauss-Seidel loop-carried dependency and fixing it via loop unrolling(loiseaujc.github.io)
2 points by loiseaujc 19d ago|0 comments
10.
The HydroGym reinforcement learning platform for fluid dynamics(nature.com)
1 points by loiseaujc 26d ago|0 comments
11.
by loiseaujc 1y ago
OP here. > This is a big red flag. If they have little background in programming, teach good standards of Python. That is a big issue, I'll admit it. And actually, in my Uni, I am a very strong proponent of making an intro to progr
12.
by loiseaujc 1y ago
OP here. I actually love Julia and acknowledge its ecosystem can be fantastic for learning purposes. Just see the Computational Thinking course at MIT. I'm not dismissing Julia. Actually, the first lines of my conclusions are > In t
13.
by loiseaujc 1y ago
OP here. 0-based vs 1-based indexing is one these hot debates I don't want get in as I've mentioned in the post. Yet, point is, if you take pretty much any math textbook, indexing starts at 1, whether we like it or not. And for ma