Work in using the Matrix package which directly wraps LAPACK routines into common %*%, chol(), t() etc methods .

Drastic slowdown with this version!! Something is very, very, very wrong.

Runtime for rNIW.snappy2.5 :
   user  system elapsed
206.853   0.683 208.451
TODO: what's up with Inverse Wisharts??
Slow as mud because S4 is slow as mud:

                             self.time self.pct total.time total.pct
"standardGeneric"                13.40     6.58     151.30     74.30
".Call"                          12.60     6.19      16.58      8.14
"match"                          11.70     5.75      27.40     13.46
"possibleExtends"                 7.94     3.90      18.30      8.99
".identC"                         6.60     3.24       9.82      4.82
"exists"                          6.04     2.97       6.40      3.14
"validObject"                     5.84     2.87      78.88     38.74
".getClassFromCache"              5.82     2.86      11.44      5.62
"assign"                          5.30     2.60       5.30      2.60
"is"                              5.24     2.57      32.96     16.19
"getClassDef"                     5.16     2.53      15.16      7.44
"ls"                              5.04     2.47      15.12      7.42
"doTryCatch"                      4.44     2.18     203.64    100.00
"paste"                           4.28     2.10       4.50      2.21
"tryCatch"                        3.82     1.88     203.64    100.00
"get"                             3.80     1.87       3.80      1.87
".cacheGenericTable"              3.58     1.76       5.70      2.80
"initialize"                      3.26     1.60      96.54     47.41
"as.vector"                       2.98     1.46      20.70     10.16
".getGeneric"                     2.90     1.42      14.24      6.99
"as"                              2.88     1.41     114.62     56.29
"loadMethod"                      2.74     1.35      15.08      7.41
"structure"                       2.16     1.06       4.30      2.11
"tryCatchList"                    2.06     1.01     203.64    100.00
"%*%"                             2.00     0.98       9.50      4.67
...

doTryCatch and paste and tryCatchList might be my calls, but "standardGeneric" and "possibleExtends" are defs not.
There's a way to precache the method: use selectMethod
Hilarious fail: `Matrix` made everything slower, not faster.
Ah, I discovered that crossprod() was built in already--Matrix's version was just an implementation, not a necessity.

Update all algorithms to use it.
Kill snappy2.5, which is a pile of fail.
Instead of snappy3.5, plans for all the paths forward.
How to sample Inverse Wisharts.