Вот такую библиотечку сделал.

array

Обертка над стандартными массивами Си.

Поддерживает тензорное индексирование:

c_array(x)[_i][_j] = c_array(y)[_i][_k] * c_array(z)[_k][_j] + c_array(z)[_i][_j];

Срезы

c_array(x)(_j,_i) = c_array(y)(_i,_j);

Производительностью, увы, похвастать не могу. На умножении матриц на Интеле отстаю в 2 раза

Выкидывать жалко, может кому пригодится...