Решение задачи нахождения собственных векторов в общем виде

Все бобра! Есть ли в python функция, которая решает задачу eigenvalue decomposition в общем виде, т е для двух матриц?

Numpy.linalg.eig(a) принимает на вход только один аргумент. Я конечно понимаю что можно в эту функцию засунуть B^(-1)*A, но тогда нужно будет заморачиваться по поводу обратимости матрицы B.
Есть ли готовая функция, как в matlab/octave например?