If I use a series of operations like

```
array = pm.Normal( ... )
array_reshaped = array[:, np.newaxis]
```

and then go on to use `array_reshaped`

later in the model, is that expected to cause problems for PyMC/PyTensor?

In other words, when the documentation for PyTensor says that it has tight integration with numpy, does that extend to these kinds of operations? Or do I have to use the PyTensor equivalent? And if so, what would it be?

Apologies if this isn’t clear, and thanks in advance for any guidance!