# Hypermatrices?

Hi all,

I want to store and manipulate arrays of $n\times n\times n$ numbers $$x_{ijk}\quad\text{where}\quad 0\leq i, j, k\leq n-1$$ with Sage.

Here is what I have been doing so far:

`import numpy as np`

`n = 7`

`X = np.zeros((n, n, n))`

`X[2][4][1] = 17`

etc...

Remark: So $X$ is a so called 3-hypermatrix.

**Question:** Are NumPy arrays my best option when handling hypermatrices with Sage?

**Precision**: I will need to manipulate lists of such hypermatrices, and even hypermatrices whose elements will be hypermatrices, e.g.
$$X_{ijk}^{(abcd)}$$
and loop over such lists, e.g. for a ..., for b ..., for c ..., do ... with the 3-hypermatrix $X^{(abcd)}$.