subdivisor.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * Subdivisor
 * 2010-05-16
 *****************************************************************************
 *               Copyright (C) 2010 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef realroot_solver_subdivisor_hpp
# define realroot_solver_subdivisor_hpp

# include <stack>
# include <realroot/Seq.hpp>

# define TMPL template<class CELL, class V>
# define SELF subdivisor<CELL, V>
# define Stack std::stack<Cell*>
//====================================================================
namespace mmx {

TMPL 
class subdivisor {
  
public:
  typedef CELL Cell;
  subdivisor(Cell* c) ;
  ~subdivisor(void) ;
  
  void run(void) ;
  
  Seq<Cell*> solutions() {return m_sols;}
private:
  Stack*  m_stack ;
  Seq<Cell*> m_sols;
};

//--------------------------------------------------------------------
TMPL
SELF::subdivisor(CELL * object) {
  m_stack = new Stack ;
  m_stack -> push(object);
}

TMPL
SELF::~subdivisor(void) {
  delete m_stack ;
}

TMPL void 
SELF::run() {

  Cell* cl=NULL;  
  while( !m_stack->empty() ) {
    
    cl = m_stack->top();    m_stack->pop();

    if(V::reduce(cl)) {
      V::subdivide(cl,m_stack) ;
      delete cl;
    } else if(V::regular(cl)) {
      m_sols<<cl;
    } else
      delete cl;

    
  }
}
//====================================================================
} // namespace mmx
//====================================================================
# undef TMPL
# undef SELF
# undef Stack
# endif