include/shape/algebraic_surface.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * AlgebraicSurface
 * 2008-03-28
 * Julien Wintz
 *****************************************************************************
 *               Copyright (C) 2006 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef shape_algebraic_surface_hpp
# define shape_algebraic_surface_hpp

# include <realroot/GMP.hpp>
# include <shape/algebraic_set.hpp>
# include <shape/surface.hpp>
# define TMPL   template<class C, class V>
# define TMPL1  template<class V>
# define SELF   algebraic_surface<C,V>
# define Viewer viewer<axel,V>
//====================================================================
namespace mmx {
namespace shape {

// TMPL struct algebraic_surface;

TMPL class algebraic_surface;

//--------------------------------------------------------------------
struct algebraic_surface_def {};

// template<> struct use<algebraic_surface_def>
//   : public use<algebraic_set_def> {
//   typedef algebraic_surface<double,double> AlgebraicSurface;
// };
//--------------------------------------------------------------------
template<class C, class V=default_env>
class algebraic_surface : public surface<REF_OF(V)>
{
public:
  typedef REF_OF(V) Ref;
  typedef typename algebraic_set<C,Ref>::Polynomial Polynomial;
  //typedef polynomial<C, with<Sparse,DegRevLex> >    Polynomial;
  typedef bounding_box<C,Ref>                       BoundingBox;
  typedef surface<Ref>                              Surface;

  algebraic_surface(void) {} ;

  algebraic_surface(char *) ;
  algebraic_surface(const Polynomial &) ;

  ~algebraic_surface(void) {};

  Polynomial set_equation(const Polynomial& p)  { return this->m_polynomial=p; }

  Polynomial equation() const  { return this->m_polynomial ; }



private:
    Polynomial m_polynomial ;
} ;

TMPL
SELF::algebraic_surface(char * s): Surface()
{
  this->set_equation(Polynomial(s));
}

TMPL
SELF::algebraic_surface(const Polynomial& p) : Surface()
{
  this->set_equation(p);
}


//--------------------------------------------------------------------
  TMPL struct viewer; struct axel;
TMPL Viewer& 
operator<<(Viewer& os, const SELF& s) {
  os<<"<surface type=\"implicit\" color=\""
    <<int(os.color.r*255)<<" "
    <<int(os.color.g*255)<<" "
    <<int(os.color.b*255)<<"\">\n";
  os<<"  <domain>"
    <<os(0,0)<<" "<<os(0,1)<<" "
    <<os(1,0)<<" "<<os(1,1)<<" "
    <<os(2,0)<<" "<<os(2,1)
    <<"</domain>\n";
  os<<"  <polynomial>";
  print_as_double(os,s.equation().rep(),variables("x y z"));
  os<<"</polynomial>\n";
  os<<"</surface>\n";
  return os;
}

//====================================================================
} ; // namespace shape
//====================================================================
} ; // namespace mmx
//====================================================================
# undef TMPL 
# undef TMPL1 
# undef Surface 
# undef SELF 
# undef Viewer
# endif // IMPLICITSURFACE_H

---