CholeskyUtilities.hpp

Go to the documentation of this file.

#ifndef ZONOOPT_CHOLESKY_UTILITIES_HPP_
#define ZONOOPT_CHOLESKY_UTILITIES_HPP_
 
#include "Eigen/Sparse"
#include "Eigen/Dense"
 
namespace ZonoOpt::detail {
 
    struct LDLT_data
    {
        Eigen::SparseMatrix<zono_float> L;
        Eigen::DiagonalMatrix<zono_float, -1> Dinv;
        Eigen::PermutationMatrix<-1, -1, int> P, Pinv;
        bool factorized = false;
    };
 
    void get_LDLT_data(const Eigen::SimplicialLDLT<Eigen::SparseMatrix<zono_float>>& solver, LDLT_data& data);
 
    Eigen::Vector<zono_float, -1> solve_LDLT(const LDLT_data& data, const Eigen::Vector<zono_float, -1>& b);
 
    void affine_set_projection(Eigen::Ref<Eigen::Vector<zono_float, -1>> z, const Eigen::SparseMatrix<zono_float>& A,
        const Eigen::SparseMatrix<zono_float>& AT, const Eigen::Vector<zono_float, -1>& b, const LDLT_data& AAT_ldlt);
 
} // end namespace detail
// end namespace ZonoOpt
 
#endif