The @-operator PEP (gh-4351) describes two operations that would be a valuable addition to numpy independently of the main purpose of the PEP.

Most of the functionality is already available in various places (np.dot, np.linalg.matrix_power, etc.) but no single function conforms to the PEP.