Matrice creuse, usage?

Bonjour à tous,
j'utilise numpy pour faire du calcul, et dans le cadre de mes travaux je dois construire des matrices à plusieurs dimensions de taille non négligeable. Il s'agit de matrices extrêmement creuses.
Par exemple, j'ai fait

Code:

---

1 2	a= np.zeros(((1000,1000,1000)))

---

Et ça passe, par contre avec une dimension supplémentaire:

Code:

---

1 2	a= np.zeros((((1000,1000,1000,1000))))

---

j'ai une memory error (je suis sous python 64-bit je précise) ce qui est bien normal vu la taille du truc.

Aussi, j'aimerais savoir s'il existe un moyen de gérer les matrices creuses sous numpy, de manière pas trop compliquée ?
En effet, je dois être capable d'effectuer de nombreuses manipulations sur les matrices creuses ainsi obtenues (multiplication, modification d'une valeur particulière, etc).

Merci

Citation:

---

Envoyé par bougnbie

Aussi, j'aimerais savoir s'il existe un moyen de gérer les matrices creuses sous numpy, de manière pas trop compliquée ?

---

Pas à ma connaissance.
Il existe un module sparse dans SciPy.
Et des bouts de solutions +/- partiels que vous trouverez sans peine avec Google ou autre.

- W