Added GEAM implementation.
More robust debugging, fixed problem with GEAM, need to apply the same concept across BLAS implementations.
Lots of updates. Appears to be working (worked out the algorithm by hand and matched results).

Unfortunately the performance on it still varies *wildly*. As a result I
am looking at implementing the non-optimized algorithm as well.