BTW on my box (OSX 10.8 on macbook air) this code is ~ twice slower than the numba / parakeet / pythran / cython group and ~ 3x slower than blas. Maybe it could be optimized further before merge.