lmi_oracle_base.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <memory>   // for unique_ptr
#include <utility>  // for pair, forward
#include <vector>
 
#include "ldlt_mgr.hpp"
 
template <typename Arr036, typename Mat = Arr036> class LmiOracleBase {
  protected:
    using Cut = std::pair<Arr036, double>;
 
    std::unique_ptr<Cut> cut = std::make_unique<Cut>();
 
    template <typename LDLT, typename Fn> auto assess_impl(LDLT& mgr, const std::vector<Mat>& F,
                                                           const int sign, const Arr036& x,
                                                           Fn&& getA) -> Cut* {
        const auto n = x.size();
        if (mgr.factor(std::forward<Fn>(getA))) {
            return nullptr;
        }
        const auto ep = mgr.witness();  // call before sym_quad() !!!
        Arr036 g{x};
        for (auto i = 0U; i != n; ++i) {
            g[i] = sign * mgr.sym_quad(F[i]);
        }
        this->cut->first = std::move(g);
        this->cut->second = std::move(ep);
        return this->cut.get();
    }
};