lmi0_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 Lmi0Oracle
    : public LmiOracleBase<Arr036, Mat> {
    using Base = LmiOracleBase<Arr036, Mat>;
    using Cut = std::pair<Arr036, double>;
 
  public:
    LDLTMgr _mq;  
 
  private:
    const std::vector<Mat>& m_F;  
 
  public:
    Lmi0Oracle(size_t ndim, const std::vector<Mat>& F) : _mq(ndim), m_F{F} {}
 
    auto assess_feas(const Arr036& x) -> Cut* {
        const auto n = x.size();
 
        auto getA = [&n, &x, this](size_t i, size_t j) -> double {
            auto a = 0.0;
            for (auto k = 0U; k != n; ++k) {
                a += this->m_F[k](i, j) * x[k];
            }
            return a;
        };
 
        return this->assess_impl(this->_mq, this->m_F, -1, x, getA);
    }
 
    auto operator()(const Arr036& x) -> Cut* { return assess_feas(x); }
};