Structure de données pour représenter le tenseur clairsemé?

Question

Qu'est-ce qu'une structure de données appropriée pour représenter un tesnor clairsemé en C++? La première option qui vient à l'esprit est un boost::unordered_map car il permet des opérations telles que la configuration rapide et la récupération d'un élément comme ci-dessous:

A(i,j,k,l) = 5 

Cependant, je voudrais aussi pouvoir faire des contractions sur un seul indice , ce qui impliquerait la sommation sur l'un des indices

C(i,j,k,m) = A(i,j,k,l)*B(l,m) 

Comment serait-il facile à mettre en œuvre cet opérateur avec un boost::unordered_map? Y a-t-il une structure de données plus appropriée?

Answer 1

Il existe des bibliothèques tenseurs disponibles, comme:

http://www.codeproject.com/KB/recipes/tensor.aspx 

et

http://cadadr.org/fm/package/ftensor.html 

Toute question avec ceux? Vous obtiendriez plus d'opérations tensorielles en utilisant une carte.