# The N-dimensional array (ndarray)¶

An ndarray is a (usually fixed-size) multidimensional container of items of the same type and size. The number of dimensions and items in an array is defined by its shape, which is a tuple of N non-negative integers that specify the sizes of each dimension. The type of items in the array is specified by a separate data-type object (dtype), one of which is associated with each ndarray.

As with other container objects in Python, the contents of an ndarray can be accessed and modified by indexing or slicing the array (using, for example, N integers), and via the methods and attributes of the ndarray.

Different ndarrays can share the same data, so that changes made in one ndarray may be visible in another. That is, an ndarray can be a “view” to another ndarray, and the data it is referring to is taken care of by the “base” ndarray.

Example

A 2-dimensional array of size 2 x 3, composed of 4-byte integer elements:

>>> x = np.array([[1, 2, 3], [4, 5, 6]], np.int32)
>>> type(x)
<class 'mxnet.numpy.ndarray'>
>>> x.shape
(2, 3)
>>> x.dtype
dtype('int32')


The array can be indexed using Python container-like syntax:

>>> # The element of x in the *second* row, *third* column, namely, 6.
>>> x[1, 2]
array(6, dtype=int32)  # this is different than the official NumPy which returns a np.int32 object


For example slicing can produce views of the array if the elements to be sliced is continguous in memory:

>>> y = x[1,:]
>>> y
array([9, 5, 6], dtype=int32)  # this also changes the corresponding element in x
>>> x
array([[1, 2, 3],
[9, 5, 6]], dtype=int32)


## Constructing arrays¶

New arrays can be constructed using the routines detailed in Array creation routines, and also by using the low-level ndarray constructor:

 ndarray ndarray(handle, writable=True):

## Indexing arrays¶

Arrays can be indexed using an extended Python slicing syntax, array[selection].

## Internal memory layout of an ndarray¶

An instance of class ndarray consists of a contiguous one-dimensional segment of computer memory (owned by the array, or by some other object), combined with an indexing scheme that maps N integers into the location of an item in the block. The ranges in which the indices can vary is specified by the shape of the array. How many bytes each item takes and how the bytes are interpreted is defined by the data-type object associated with the array.

Note

mxnet.numpy.ndarray currently only supports storing elements in C-order/row-major and contiguous memory space. The following content on explaining a variety of memory layouts of an ndarray are copied from the official NumPy documentation as a comprehensive reference.

A segment of memory is inherently 1-dimensional, and there are many different schemes for arranging the items of an N-dimensional array in a 1-dimensional block. NumPy is flexible, and ndarray objects can accommodate any strided indexing scheme. In a strided scheme, the N-dimensional index $$(n_0, n_1, ..., n_{N-1})$$ corresponds to the offset (in bytes):

$n_{\mathrm{offset}} = \sum_{k=0}^{N-1} s_k n_k$

from the beginning of the memory block associated with the array. Here, $$s_k$$ are integers which specify the strides of the array. The column-major order (used, for example, in the Fortran language and in Matlab) and row-major order (used in C) schemes are just specific kinds of strided scheme, and correspond to memory that can be addressed by the strides:

$s_k^{\mathrm{column}} = \mathrm{itemsize} \prod_{j=0}^{k-1} d_j , \quad s_k^{\mathrm{row}} = \mathrm{itemsize} \prod_{j=k+1}^{N-1} d_j .$

where $$d_j$$ = self.shape[j].

Both the C and Fortran orders are contiguous, i.e., single-segment, memory layouts, in which every part of the memory block can be accessed by some combination of the indices.

While a C-style and Fortran-style contiguous array, which has the corresponding flags set, can be addressed with the above strides, the actual strides may be different. This can happen in two cases:

1. If self.shape[k] == 1 then for any legal index index[k] == 0. This means that in the formula for the offset $$n_k = 0$$ and thus $$s_k n_k = 0$$ and the value of $$s_k$$ = self.strides[k] is arbitrary.

2. If an array has no elements (self.size == 0) there is no legal index and the strides are never used. Any array with no elements may be considered C-style and Fortran-style contiguous.

Point 1. means that self and self.squeeze() always have the same contiguity and aligned flags value. This also means that even a high dimensional array could be C-style and Fortran-style contiguous at the same time.

An array is considered aligned if the memory offsets for all elements and the base offset itself is a multiple of self.itemsize. Understanding memory-alignment leads to better performance on most hardware.

Note

Points (1) and (2) are not yet applied by default. Beginning with NumPy 1.8.0, they are applied consistently only if the environment variable NPY_RELAXED_STRIDES_CHECKING=1 was defined when NumPy was built. Eventually this will become the default.

You can check whether this option was enabled when your NumPy was built by looking at the value of np.ones((10,1), order='C').flags.f_contiguous. If this is True, then your NumPy has relaxed strides checking enabled.

Warning

It does not generally hold that self.strides[-1] == self.itemsize for C-style contiguous arrays or self.strides[0] == self.itemsize for Fortran-style contiguous arrays is true.

Data in new ndarrays is in the row-major (C) order, unless otherwise specified, but, for example, basic array slicing often produces views in a different scheme.

Note

Several algorithms in NumPy work on arbitrarily strided arrays. However, some algorithms require single-segment arrays. When an irregularly strided array is passed in to such algorithms, a copy is automatically made.

## Array attributes¶

Array attributes reflect information that is intrinsic to the array itself. Generally, accessing an array through its attributes allows you to get and sometimes set intrinsic properties of the array without creating a new array. The exposed attributes are the core parts of an array and only some of them can be reset meaningfully without creating a new array. Information on each attribute is given below.

### Memory layout¶

The following attributes contain information about the memory layout of the array:

 ndarray.shape Tuple of array dimensions. ndarray.ndim Number of array dimensions. ndarray.size Number of elements in the array.

### Data type¶

The data type object associated with the array can be found in the dtype attribute:

 ndarray.dtype Data-type of the array’s elements.

## Array methods¶

An ndarray object has many methods which operate on or with the array in some fashion, typically returning an array result. These methods are briefly explained below. (Each method’s docstring has a more complete description.)

For the following methods there are also corresponding functions in numpy: all(), any(), argmax(), argmin(), argpartition(), argsort(), choose(), clip(), compress(), copy(), cumprod(), cumsum(), diagonal(), imag(), max, mean(), min, nonzero(), partition(), prod(), ptp(), put(), ravel(), real(), repeat(), reshape(), round, searchsorted(), sort(), squeeze(), std(), sum(), swapaxes(), take(), trace(), transpose(), var().

### Array conversion¶

 ndarray.item(*args) Copy an element of an array to a standard Python scalar and return it. ndarray.copy([order]) Return a coyp of the array, keeping the same device. ndarray.astype(dtype[, order, casting, …]) Copy of the array, cast to a specified type.

### Shape manipulation¶

For reshape, resize, and transpose, the single tuple argument may be replaced with n integers which will be interpreted as an n-tuple.

 ndarray.reshape(*args, **kwargs) Returns a copy of the array with a new shape. ndarray.transpose(*axes) Permute the dimensions of an array. ndarray.swapaxes(axis1, axis2) Return a copy of the array with axis1 and axis2 interchanged. ndarray.flatten([order]) Return a copy of the array collapsed into one dimension. ndarray.squeeze([axis]) Remove single-dimensional entries from the shape of a.

### Item selection and manipulation¶

For array methods that take an axis keyword, it defaults to None. If axis is None, then the array is treated as a 1-D array. Any other value for axis represents the dimension along which the operation should proceed.

 ndarray.nonzero() Return the indices of the elements that are non-zero. ndarray.take(indices[, axis, mode]) Convenience fluent method for take(). ndarray.repeat(repeats[, axis]) Repeat elements of an array. ndarray.argsort([axis, descending, stable]) Convenience fluent method for argsort(). ndarray.sort([axis, descending, stable]) Convenience fluent method for sort().

### Calculation¶

Many of these methods take an argument named axis. In such cases,

• If axis is None (the default), the array is treated as a 1-D array and the operation is performed over the entire array. This behavior is also the default if self is a 0-dimensional array or array scalar. (An array scalar is an instance of the types/classes float32, float64, etc., whereas a 0-dimensional array is an ndarray instance containing precisely one array scalar.)

• If axis is an integer, then the operation is done over the given axis (for each 1-D subarray that can be created along the given axis).

Example of the axis argument

A 3-dimensional array of size 3 x 3 x 3, summed over each of its three axes

>>> x
array([[[ 0,  1,  2],
[ 3,  4,  5],
[ 6,  7,  8]],
[[ 9, 10, 11],
[12, 13, 14],
[15, 16, 17]],
[[18, 19, 20],
[21, 22, 23],
[24, 25, 26]]])
>>> x.sum(axis=0)
array([[27, 30, 33],
[36, 39, 42],
[45, 48, 51]])
>>> # for sum, axis is the first keyword, so we may omit it,
>>> # specifying only its value
>>> x.sum(0), x.sum(1), x.sum(2)
(array([[27, 30, 33],
[36, 39, 42],
[45, 48, 51]]),
array([[ 9, 12, 15],
[36, 39, 42],
[63, 66, 69]]),
array([[ 3, 12, 21],
[30, 39, 48],
[57, 66, 75]]))


The parameter dtype specifies the data type over which a reduction operation (like summing) should take place. The default reduce data type is the same as the data type of self. To avoid overflow, it can be useful to perform the reduction using a larger data type.

For several methods, an optional out argument can also be provided and the result will be placed into the output array given. The out argument must be an ndarray and have the same number of elements. It can have a different data type in which case casting will be performed.

 ndarray.max([axis, out, keepdims]) Return the maximum along a given axis. ndarray.argmax([axis, out, keepdims]) Return indices of the maximum values along the given axis. ndarray.min([axis, out, keepdims]) Convenience fluent method for min(). ndarray.argmin([axis, out, keepdims]) Return indices of the minium values along the given axis. ndarray.clip([min, max, out]) Return an array whose values are limited to [min, max]. ndarray.sum([axis, dtype, out, keepdims]) Return the sum of the array elements over the given axis. ndarray.mean([axis, dtype, out, keepdims]) Returns the average of the array elements along given axis. ndarray.prod([axis, dtype, out, keepdims]) Return the product of the array elements over the given axis. ndarray.cumsum([axis, dtype, out]) Return the cumulative sum of the elements along the given axis. ndarray.var([axis, dtype, out, correction, …]) Returns the variance of the array elements, along given axis. ndarray.std([axis, dtype, out, correction, …]) Returns the standard deviation of the array elements along given axis. ndarray.round([decimals, out]) Convenience fluent method for round(). ndarray.all([axis, out, keepdims]) ndarray.any([axis, out, keepdims])

## Arithmetic, matrix multiplication, and comparison operations¶

Arithmetic and comparison operations on ndarrays are defined as element-wise operations, and generally yield ndarray objects as results.

Each of the arithmetic operations (+, -, *, /, //, %, divmod(), ** or pow(), <<, >>, &, ^, |, ~) and the comparisons (==, <, >, <=, >=, !=) is equivalent to the corresponding universal function (or ufunc for short) in NumPy.

Comparison operators:

 ndarray.__lt__(other) x.__lt__(y) <=> x < y ndarray.__le__(other) x.__le__(y) <=> x <= y ndarray.__gt__(other) x.__gt__(y) <=> x > y ndarray.__ge__(other) x.__ge__(y) <=> x >= y ndarray.__eq__(other) x.__eq__(y) <=> x == y ndarray.__ne__(other) x.__ne__(y) <=> x != y

Truth value of an array (bool()):

 ndarray.__bool__()

Note

Truth-value testing of an array invokes ndarray.__bool__(), which raises an error if the number of elements in the array is larger than 1, because the truth value of such arrays is ambiguous.

Unary operations:

 ndarray.__neg__() x.__neg__() <=> -x ndarray.__abs__() ndarray.__invert__() x.__invert__() <=> ~x

Arithmetic:

 ndarray.__add__(other) x.__add__(y) <=> x + y ndarray.__sub__(other) x.__sub__(y) <=> x - y ndarray.__mul__(other) x.__mul__(y) <=> x * y ndarray.__truediv__(other) x.__truediv__(y) <=> x / y ndarray.__mod__(other) x.__mod__(y) <=> x % y ndarray.__pow__(other) x.__pow__(y) <=> x ** y ndarray.__and__(other) x.__and__(y) <=> x & y ndarray.__or__(other) x.__or__(y) <=> x | y ndarray.__xor__(other) x.__xor__(y) <=> x ^ y

Note

• Any third argument to pow() is silently ignored, as the underlying ufunc takes only two arguments.

• The three division operators are all defined; div is active by default, truediv is active when __future__ division is in effect.

• Because ndarray is a built-in type (written in C), the __r{op}__ special methods are not directly defined.

• The functions called to implement many arithmetic special methods for arrays can be modified using __array_ufunc__.

Arithmetic, in-place:

 ndarray.__iadd__(other) x.__iadd__(y) <=> x += y ndarray.__isub__(other) x.__isub__(y) <=> x -= y ndarray.__imul__(other) x.__imul__(y) <=> x *= y ndarray.__itruediv__(other) x.__itruediv__(y) <=> x /= y ndarray.__imod__(other) x.__imod__(y) <=> x %= y ndarray.__iand__(other) x.__iand__(y) <=> x &= y ndarray.__ior__(other) x.__ior__(y) <=> x |= y ndarray.__ixor__(other) x.__ixor__(y) <=> x ^= y

Warning

In place operations will perform the calculation using the precision decided by the data type of the two operands, but will silently downcast the result (if necessary) so it can fit back into the array. Therefore, for mixed precision calculations, A {op}= B can be different than A = A {op} B. For example, suppose a = ones((3,3)). Then, a += 3j is different than a = a + 3j: while they both perform the same computation, a += 3 casts the result to fit back in a, whereas a = a + 3j re-binds the name a to the result.

Matrix Multiplication:

 ndarray.__matmul__(other) x.__matmul__(y) <=> x @ y

## Special methods¶

For standard library functions:

 ndarray.__reduce__() Helper for pickle. ndarray.__setstate__(state)

Basic customization:

 ndarray.__new__ Create and return a new object.

Container customization: (see Indexing)

 ndarray.__len__() Number of elements along the first axis. ndarray.__getitem__(key) Return self[key]. ndarray.__setitem__(key, value) Sets self[key] to value.

Conversion; the operations index(), int() and float(). They work only on arrays that have one element in them and return the appropriate scalar.

 ndarray.__index__() ndarray.__int__() ndarray.__float__()

String representations:

 ndarray.__str__() Returns a string representation of the array. ndarray.__repr__() Returns a string representation of the array.