Implementing the Bartels–Stewart algorithm in Eigen3?
问题 In the past when I've needed to solve the Sylvester equation, AX + XB = C , I've used scipy 's function, solve_sylvester [1], which apparently works by using the Bartels-Stewart algorithm to get things into upper triangular form, and then solving the equation using lapack . I now need to solve the equation using eigen . eigen provides an function, matrix_function_solve_triangular_sylvester [2], which seems by the documentation to be similar to the lapack function which scipy calls. I'm