Thanks! Looks like the solver hits a bump while accelerating on a high gear. That should not be a problem, though. I'll take a look at the code, but it will take time.Meanwhile, if you have time, you can add "#define DEBUG" in the beginning of msadams.c and recompile GSL. This should make msadams print out _lots_ of debug info and data, that would help decipher this issue. Interesting stuff happens between t=5116.13 and t=5116.27.
Yeah, wasn't quite expecting a 1 GB output file. Attached is a truncated section that has the steps in question. I've prefixed "PROGRAM" to the output lines that are not coming out of msadams.c
M
stdout.dat.trunc.gz
Description: GNU Zip compressed data
