#include "MUQ/Modeling/LinearAlgebra/LyapunovSolver.h"

#include <Eigen/Eigenvalues>

using namespace muq::Modeling;


void LyapunovSolver::compute(Eigen::MatrixXd const& A, Eigen::MatrixXd const& Q)
{

    // Compute the Schur complement of A
    Eigen::RealSchur<Eigen::MatrixXd> 
    

}

