Statevector.evolve fails for repeated calls on non-qubit statevectors

[chriseclectic]

Environment

• Qiskit Terra version: 0.22
• Python version: All
• Operating system: All

What is happening?

Repeat calls of Statevector/DensityMatrix.evolve on subsystems for non-qubit states raises numpy axis errors.

How can we reproduce the issue?

import qiskit.quantum_info as qi
import numpy as np

mat = np.eye(3)

state = qi.random_statevector((3, 3, 3), seed=123)
state.evolve(mat, [0]).evolve(mat, [1])

Some long traceback to...
ValueError: axes don't match array

This doesn't happen for qubit subsystems.

What should happen?

No error

Any suggestions?

I think I have identified the bug in the OpShape.compose method where an attribute is not being set correctly for non-qubit subsystems.