I am writing a scientific code in c which requires numerous tensor products, and higher-order tensor operations as well, e.g. contraction, kronecker product, etc. I am planning on implementing this myself using efficient algorithms as in the literature and employing a good BLAS library. I did some looking about and haven’t really found anything for c (plenty for c++, python etc.). Am I missing something or has this really not been well implemented in c?
Thanks.
Perhaps you could investigate tensor extensions to the GSL — GNU Scientific Library.