Cool - I had just made the same change and was testing it - glad to see that's it!
As an aside, in general profiling, it seems that pretty much all of the programs that take a lot more
time than I'd like (film_gls - I'm looking at you!) seem to spend a huge amount of time (over 85%) in
NEWMAT, especially in the matrix multiplication routine. This isn't surprising, given that fMRI
processing is basically just a lot of linear algebra. Have you given any thought to recoding some of
the routines in NEWMAT (especially matrix multiplication and SVD) to use LAPACK and BLAS? That
would get you vector performance for free, since pretty much every OS has hand-optimized versions
of these libraries available.
Blaise
|