template <typename T>
Constraint_sol<T> feasible(const Constraint_sys<T>& constraints, size_t sol_size);