motzkin_burger.hpp

Go to the documentation of this file.

/*****************************************************************************
    
    motzkin_burger.hpp -- Motzkin-Burger algorithm for double-description of
        a polyhedral cone.

    This file is a part of the Arageli library.

    Copyright (C) Nikolai Yu. Zolotykh, 1999--2006
    Copyright (C) Sergey S. Lyalin, 2006
    University of Nizhni Novgorod, Russia

*****************************************************************************/

#ifndef _ARAGELI_motzkin_burger_hpp_
#define _ARAGELI_motzkin_burger_hpp_

#include "config.hpp"

#include <vector>


namespace Arageli
{

namespace ctrl
{

struct skeleton_motzkin_burger_idler
{
    class abort : public ctrl::abort {};

    template <typename A, typename F, typename Q, typename E>
    void preamble (const A& a, const F& f, const Q& q, const E& e) const {}

    void begin_gauss () const {}

    void end_gauss () const {}

    template <typename I>
    void find_non_zero (const I& i) const {}

    void zero_row () const {}

    template <typename I, typename J>
    void swap_cols_rows (const I& i, const J& j) const {}

    template <typename J>
    void eliminate_col (const J& j) const {}

    template <typename A, typename F, typename Q>
    void show_matrices (const A& a, const F& f, const Q& q) const {}

    void begin_motzkin () const {}

    void end_motzkin () const {}

    template <typename I1, typename I2>
    void select_col (const I1& current_column, const I2& new_column) const {}

    void zero_solution () const {}

    void corollary_inequality () const {}

    void begin_row_balancing () const {}

    template <typename I, typename J>
    void balanced_rows (const I& j_p, const J& j_m) const {}

    void end_row_balancing () const {}

    void delete_negates () const {}

    template <typename A, typename F, typename Q, typename E>
    void conclusion (const A& a, const F& f, const Q& q, const E& e) const {}
};

} // namespace ctrl



template
<
    typename A, typename F,
    typename Q, typename E,
    typename Ctrler
>
void skeleton_motzkin_burger
(
    A& a, F& f, Q& q, E& e,
    Ctrler ctrler
);



template <typename A, typename F, typename Q, typename E>
inline void skeleton_motzkin_burger (A& a, F& f, Q& q, E& e)
{
    return skeleton_motzkin_burger
        (a, f, q, e, ctrl::skeleton_motzkin_burger_idler());
}



template
<
    typename A, typename F,
    typename Q, typename E,
    typename Ctrler
>
void skeleton_motzkin_burger_min
(
    A& a, F& f, Q& q, E& e,
    Ctrler ctrler
);



template <typename A, typename F, typename Q, typename E>
inline void skeleton_motzkin_burger_min (A& a, F& f, Q& q, E& e)
{
    return skeleton_motzkin_burger_min
        (a, f, q, e, ctrl::skeleton_motzkin_burger_idler());
}



template
<
    typename A, typename F,
    typename Q, typename E,
    typename Ctrler
>
void skeleton_motzkin_burger_max
(
    A& a, F& f, Q& q, E& e,
    Ctrler ctrler
);



template <typename A, typename F, typename Q, typename E>
inline void skeleton_motzkin_burger_max (A& a, F& f, Q& q, E& e)
{
    return skeleton_motzkin_burger_max
        (a, f, q, e, ctrl::skeleton_motzkin_burger_idler());
}

// +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

template
<
    typename A, typename F,
    typename Q, typename E
>
void integer_conv_motzkin_burger (A& a, F& f, Q& q, E& e);  // Don't use!


} // namesapce Arageli


#ifdef ARAGELI_INCLUDE_CPP_WITH_EXPORT_TEMPLATE
    #define ARAGELI_INCLUDE_CPP_WITH_EXPORT_TEMPLATE_MOTZKIN_BURGER
    #include "motzkin_burger.cpp"
    #undef  ARAGELI_INCLUDE_CPP_WITH_EXPORT_TEMPLATE_MOTZKIN_BURGER
#endif

#endif

---