[Numpy-discussion] sparse array data

Wolfgang Kerzendorf wkerzendorf at gmail.com
Wed May 2 12:16:28 EDT 2012


Hi all,

I'm currently writing a code that needs three dimensional data (for the physicists it's dimensions are atom, ion, level). The problem is that not all combinations do exist (a sparse array). Sparse matrices in scipy only deal with two dimensions. The operations that I need to do on those are running functions like exp(item/constant) on all of the items. I also want to sum them up in the last dimension. What's the best way to make a class that takes this kind of data and does the required operations fast. Maybe some phycisists have implemented these things already. Any thoughts?

Cheers
   Wolfgang


More information about the NumPy-Discussion mailing list