Hmmm... that sounds good, still not sure how Tamas' FFAs would work. I was preparing to suggest the *other* route that Tamas mentions on his page, basically using the GSL data structures thinly wrapped in a Lisp object (to tie it into the GC). I guess both have their pros and cons.
BTW, I have one question, purely out of curiosity. Even if you use an FFA, GSL is going to
expect data in the format the GSL developers designed, e.g.:
00 01 02 03 XX XX XX XX
10 11 12 13 XX XX XX XX
20 21 22 23 XX XX XX XX
Which means
that either the Lisp array will have extra junk to carry around, and
hence new matrix operations will need to be defined, or some
C side trickery must be performed to translate it, which is what we are
trying to avoid, right? How do you deal
with this (or plan to)?
Best,
Zach