lmi_old_oracle.hpp

Go to the documentation of this file.

 
// -*- coding: utf-8 -*-
#pragma once
 
#include <utility>  // for move
#include <vector>
 
#include "lmi_oracle_base.hpp"
 
template <typename Arr036, typename Mat = Arr036> class LmiOldOracle
    : public LmiOracleBase<Arr036, Mat> {
    using Base = LmiOracleBase<Arr036, Mat>;
    using Cut = std::pair<Arr036, double>;
 
    LDLTMgr _mgr;
    const std::vector<Mat>& m_F;
    Mat m_F0;
 
  public:
    LmiOldOracle(size_t ndim, const std::vector<Mat>& F, Mat B)
        : _mgr{ndim}, m_F{F}, m_F0{std::move(B)} {}
    auto assess_feas(const Arr036& x) -> Cut* {
        const auto n = x.size();
 
        Mat A{this->m_F0};
        for (auto k = 0U; k != n; ++k) {
            for (auto i = 0U; i != this->_mgr._n; ++i) {
                for (auto j = 0U; j != this->_mgr._n; ++j) {
                    A(i, j) -= this->m_F[k](i, j) * x[k];
                }
            }
        }
 
        auto getA = [&A](size_t i, size_t j) { return A(i, j); };
        return this->assess_impl(this->_mgr, this->m_F, +1, x, getA);
    }
 
    auto operator()(const Arr036& x) -> Cut* { return assess_feas(x); }
};