refactor polynomial array-like operations to allow numpy arrays over polynomials
action matrix solver working via eigenvalues only