boundary_helper.hpp

#ifndef BOOST_REAL_BOUNDARY_HELPER_HPP
#define BOOST_REAL_BOUNDARY_HELPER_HPP

#include <algorithm>
#include "interval.hpp"

namespace boost {
    namespace real {
        namespace helper {

            bool aligned_vectors_is_lower(const std::vector<int> &lhs, const std::vector<int> &rhs) {

                // Check if lhs is lower than rhs
                auto lhs_it = lhs.cbegin();
                auto rhs_it = rhs.cbegin();
                while (rhs_it != rhs.end() && lhs_it != lhs.end() && *lhs_it == *rhs_it) {
                    ++lhs_it;
                    ++rhs_it;
                }

                if (rhs_it != rhs.cend() && lhs_it != lhs.cend()) {
                    return *lhs_it < *rhs_it;
                }

                bool lhs_all_zero = std::all_of(lhs_it, lhs.cend(), [](int i){ return i == 0; });
                bool rhs_all_zero = std::all_of(rhs_it, rhs.cend(), [](int i){ return i == 0; });

                return lhs_all_zero && !rhs_all_zero;
            }

        }
    }
}

#endif //BOOST_REAL_BOUNDARY_HELPER_HPP