Attempt to add streams and improve BLAS performance.

Streams aren't giving anything, likely because BLAS
calls are kernel-launch bound.  Will investigate more
later.

Dt = D^T and dgemm('T','N' is best, so dgemms changed
accordingly.  #ifdef'ed sections are attempts to coalece
the gemm loop; don't work yet.