Stocking large numbers into numpy array

by

As Nils Werner already mentioned, numpy’s native ctypes cannot save numbers that large, but python itself can since the int objects use an arbitrary length implementation.
So what you can do is tell numpy not to convert the numbers to ctypes but use the python objects instead. This will be slower, but it will work.

In [14]: x = np.array([18,30,31,31,15], dtype=object)

In [15]: 150**x
Out[15]: 
array([1477891880035400390625000000000000000000L,
       191751059232884086668491363525390625000000000000000000000000000000L,
       28762658884932613000273704528808593750000000000000000000000000000000L,
       28762658884932613000273704528808593750000000000000000000000000000000L,
       437893890380859375000000000000000L], dtype=object)

In this case the numpy array will not store the numbers themselves but references to the corresponding int objects. When you perform arithmetic operations they won’t be performed on the numpy array but on the objects behind the references.
I think you’re still able to use most of the numpy functions with this workaround but they will definitely be a lot slower than usual.

But that’s what you get when you’re dealing with numbers that large 😀
Maybe somewhere out there is a library that can deal with this issue a little better.

Just for completeness, if precision is not an issue, you can also use floats:

In [19]: x = np.array([18,30,31,31,15], dtype=np.float64)

In [20]: 150**x
Out[20]: 
array([  1.47789188e+39,   1.91751059e+65,   2.87626589e+67,
         2.87626589e+67,   4.37893890e+32])

More Related Contents:

  1. What are the advantages of NumPy over regular Python lists?
  2. ValueError: setting an array element with a sequence
  3. Numpy `logical_or` for more than two arguments
  4. What is the difference between contiguous and non-contiguous arrays?
  5. Calculate mean across dimension in a 2D array
  6. Quick way to upsample numpy array by nearest neighbor tiling [duplicate]
  7. How does NumPy’s transpose() method permute the axes of an array?
  8. Count consecutive occurences of values varying in length in a numpy array
  9. Concatenate two NumPy arrays vertically
  10. Getting the indices of several elements in a NumPy array at once
  11. How to do n-D distance and nearest neighbor calculations on numpy arrays
  12. Multiplying across in a numpy array
  13. initialize a numpy array
  14. How to conditionally combine two numpy arrays of the same shape
  15. What does dtype=object mean while creating a numpy array?
  16. Numpy: efficient way to generate combinations from given ranges
  17. Selecting close matches from one array based on another reference array
  18. Repeat each values of an array different times
  19. Input and output numpy arrays to h5py
  20. Why does numpy have a corresponding function for many ndarray methods?
  21. Numpy sub-array assignment with advanced, mixed indexing
  22. Efficiently sorting a numpy array in descending order?
  23. Numpy: Divide each row by a vector element
  24. Setting the fmt option in numpy.savetxt
  25. Access n-th dimension in python [duplicate]
  26. Numpy shuffle multidimensional array by row only, keep column order unchanged
  27. Randomly select from numpy array
  28. How to “scale” a numpy array?
  29. Generate random array of floats between a range
  30. Numpy `ValueError: operands could not be broadcast together with shape …`

Leave a Comment