Just a curious onlooker - don't mean to criticize. A lot of this seems to be recreating stuff that Fortran 95/2003 does natively, but I guess this is for C/C++ people?
The library also implements the algorithms in three backends (CUDA, OpenCL and native CPU) using the same API. We'll be adding support for SSE/AVX/NEON to make it more performance portable inthe future.
EDIT: The CUDA and OpenCL backends will obviously be faster than MKL. We'll be adding SSE / AVX support at some point which'll make the CPU backend faster as well.
We wrote some blogs on this:
-http://arrayfire.com/triangle-counting-in-graphs-on-the-gpu-...
-http://arrayfire.com/triangle-counting-in-graphs-on-the-gpu-...
We plan on looking at additional algorithms in the future.