Struct NewtonRaphsonHPolyoracle

template<typename Polytope, class bfunc>
struct NewtonRaphsonHPolyoracle

Compute intersection of H-polytope P := Ax <= b with curve p(t) = sum a_j phi_j(t) where phi_j are basis functions (e.g. polynomials) Uses Newton-Raphson to solve the transcendental equation

Template Parameters:
  • Polytope – Polytope Type

  • bfunc – feasibility constraint type

Public Types

typedef Polytope::MT MT
typedef Polytope::VT VT
typedef Polytope::NT NT
typedef Polytope::PointType Point

Public Functions

inline std::tuple<NT, Point, int> apply(NT t_prev, NT t0, NT eta, MT &A, VT &b, Polytope &P, std::vector<Point> &coeffs, bfunc phi, bfunc grad_phi, int ignore_facet = -1)