File richardson_extrapolation.hpp

template<typename Point, typename NT, typename Polytope, typename func> struct RichardsonExtrapolationODESolver

Public Types

typedef std::vector<Point> pts

typedef std::vector<Polytope*> bounds

typedef std::vector<NT> coeffs

typedef std::vector<coeffs> scoeffs

typedef std::vector<pts> ptsv

typedef std::vector<ptsv> ptsm

typedef Polytope::VT VT

Public Functions

inline RichardsonExtrapolationODESolver(NT initial_time, NT step, pts initial_state, func oracle, bounds boundaries)

inline void initialize_solver()

inline void step(int k, bool accepted)

inline void print_state()

inline void steps(int num_steps, bool accepted)

inline Point get_state(int index)

inline void set_state(int index, Point p)

Public Members

unsigned int dim

const unsigned int MAX_TRIES = 5

NT eta

NT eta_temp

NT t

NT t_prev

NT tol = 1e-7

NT error = NT(-1)

NT den

Point num

Point y

VT Ar

VT Av

RKODESolver<Point, NT, Polytope, func> *solver

func F

bounds Ks

pts xs

pts xs_prev

ptsm A

bool flag

Point x_prev_bound

int prev_facet = -1