numpy.r_ is not a function. What is it?

by

I would argue that for all purposes r_ is a function, but one implemented by a clever hack using different syntax. Mike already explained how r_ is in reality not a function, but a class instance of RClass, which has __getitem__ implemented, so that you can use it as r_[1]. The cosmetic difference is that you use square brackets instead of curved ones, so you are not doing a function call, but you are actually indexing the object. Although this is technically true, for all purposes, it works just like a function call, but one that allows some extra syntax not allowed by a normal function.

The motivation for creating r_ probably comes from Matlab’s syntax, which allows to construct arrays in a very compact way, like x = [1:10, 15, 20:10:100]. To achieve the same in numpy, you would have to do x = np.hstack((np.arange(1,11), 15, np.arange(20,110,10))). Using colons to create ranges is not allowed in python, but they do exist in the form of the slice notation to index into a list, like L[3:5], and even A[2:10, 20:30] for multi-dimensional arrays. Under the hood, these index notation gets transformed to a call to the __getitem__ method of the object, where the colon notation gets transformed into a slice object:

In [13]: class C(object):
    ...:     def __getitem__(self, x):
    ...:         print x

In [14]: c = C()

In [15]: c[1:11, 15, 20:110:10]
(slice(1, 11, None), 15, slice(20, 110, 10))

The r_ object ‘abuses’ this fact to create a ‘function’ that accepts slice notation, which also does some additional things like concatenating everything together and returning the result, so that you can write x = np.r_[1:11, 15, 20:110:10]. The “Not a function, so takes no parameters” in the documentation is slightly misleading …

More Related Contents:

  1. Efficient evaluation of a function at every cell of a NumPy array
  2. Structure of inputs to scipy minimize function
  3. Simple Digit Recognition OCR in OpenCV-Python
  4. surface plots in matplotlib
  5. Apply function to each element of a list
  6. Best way to join / merge by range in pandas
  7. How to install python modules without root access?
  8. Why does Python code use len() function instead of a length method?
  9. Fastest way to grow a numpy numeric array
  10. How to make scipy.interpolate give an extrapolated result beyond the input range?
  11. Convert NumPy array to Python list
  12. Saving to hdf5 is very slow (Python freezing)
  13. Removing duplicate columns and rows from a NumPy 2D array
  14. How to conditionally combine two numpy arrays of the same shape
  15. scipy curve_fit doesn’t like math module
  16. How to fix ‘numpy.ndarray’ object has no attribute ‘get_figure’ when plotting subplots
  17. Make distutils look for numpy header files in the correct place
  18. LogisticRegression: Unknown label type: ‘continuous’ using sklearn in python
  19. Load just part of an image in python
  20. What does “ValueError: object too deep for desired array” mean and how to fix it?
  21. Store numpy.array in cells of a Pandas.DataFrame
  22. What’s the simplest way to extend a numpy array in 2 dimensions?
  23. Techniques for working with large Numpy arrays? [duplicate]
  24. Cython: (Why / When) Is it preferable to use Py_ssize_t for indexing?
  25. ModuleNotFoundError: No module named ‘numpy.testing.nosetester’
  26. Convert a column of timestamps into periods in pandas
  27. Consecutive, Overlapping Subsets of Array (NumPy, Python)
  28. A numpy array unexpectedly changes when changing another one despite being separate
  29. Generate random array of floats between a range
  30. Understanding scipy deconvolve

Leave a Comment