22 Nov 2004 rkarl   » (Journeyer)

Took another wack at some semi-Lagrangian Navier-Stokes code thats I abaondoned some time back. Basically, it worked, but I never bothered to make the LHS symmetric. Instead, being lazy, I just solved with BiCGSTAB_l.

So, after a few really stupid errors, it now works again, solving with conjugate gradients. Unfortunatly, the bloody thing now takes forever to converge. Arg! Guess this means that I have to write preconditioning code. Oof.

It is some consolation that now, at least, each iteration is about 3x as fast

Latest blog entries     Older blog entries

New Advogato Features

New HTML Parser: The long-awaited libxml2 based HTML parser code is live. It needs further work but already handles most markup better than the original parser.

Keep up with the latest Advogato features by reading the Advogato status blog.

If you're a C programmer with some spare time, take a look at the mod_virgule project page and help us with one of the tasks on the ToDo list!