Discussion :

[bougnbie]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[wiztricks]

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