cell_uv_bernstein.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * cell for univariate bernstein equation
 * 2008-05-16
 *****************************************************************************
 *               Copyright (C) 2010 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef realroot_cell_bernstein_hpp
# define realroot_cell_bernstein_hpp

# include <realroot/Seq.hpp>
# include <realroot/sign_variation.hpp>
# include <realroot/Interval.hpp>
# include <realroot/polynomial_bernstein.hpp>
# include <realroot/polynomial_sparse.hpp>

# undef C
# define TMPL template<typename C>
# define SELF cell_uv_bernstein<C>

//====================================================================
namespace mmx {

 
TMPL 
class cell_uv_bernstein {
  
public:

  typedef C Scalar;
  typedef polynomial<C, with<Bernstein> > Polynomial;

  cell_uv_bernstein(void) ;
  cell_uv_bernstein(const Polynomial& E): m_equs(E), m_doms(Interval<C>(0,1)) {}
  cell_uv_bernstein(const Polynomial& E, const Interval<double>& D): m_equs(E), m_doms(D) {}

  template<class POL>
  cell_uv_bernstein(const POL& E, const Interval<double>& D) : m_doms(D) {
    let::assign(m_equs,E,D);
  }

  ~cell_uv_bernstein(void) ;

  Polynomial  equation() const {return m_equs;}
  Polynomial& equation()       {return m_equs;}

  Interval<double>   domain() const {return m_doms;}
  Interval<double> & domain()       {return m_doms;}

  C size(void) const { return m_doms.size(); }

private:
  Polynomial         m_equs;
  Interval<double>   m_doms;
};

//--------------------------------------------------------------------
TMPL
SELF::cell_uv_bernstein(void)  {}

TMPL
SELF::~cell_uv_bernstein(void) {}

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