/Users/mourrain/Devel/mmx/shape/include/shape/surface_rational.hpp

Go to the documentation of this file.

/*****************************************************************************
 * M a t h e m a g i x
 *****************************************************************************
 * RationalSurface
 * 2009-05-08
 * Bernard Mourrain & Julien Wintz
 *****************************************************************************
 *               Copyright (C) 2006 INRIA Sophia-Antipolis
 *****************************************************************************
 * Comments :
 ****************************************************************************/
# ifndef shape_surface_rational_hpp
# define shape_surface_rational_hppx

# include <realroot/ring_tens.hpp>
# include <shape/surface_parametric.hpp>

# define TMPL_DEF template<class K, int N=3>
# define TMPL  template<class K, int N>
# define ParametricSurface surface_parametric<K>
# define SELF   surface_rational<K,N>
//====================================================================
namespace mmx {
namespace shape {

TMPL struct surface_rational;
template<class K> struct surface_rational_def
  :public surface_def<K> {
  typedef surface_rational<K> RationalSurface;
  typedef polynom< ring<double, Bernstein> > Polynomial;
};

TMPL_DEF
class surface_rational : public ParametricSurface {
public:
  typedef typename surface_rational_def<K>::Scalar    Scalar;
  typedef typename surface_rational_def<K>::Point     Point;
  typedef bounding_box<K> BoundingBox;
  typedef typename surface_rational_def<K>::polynomial_type Polynomial;

  surface_rational(void) ;
  surface_rational(const BoundingBox &) ;
   
  surface_rational(const Polynomial& X, const Polynomial& Y, const Polynomial& Z, const Polynomial& W=1): 
    ParametricSurface(),m_umin(0), m_umax(1), m_vmin(0), m_vmax(1), m_w(W) {
    m_p[0]=X; m_p[1]=Y; m_p[2]=Z;
  }
  Point* eval(Scalar u, Scalar v) const ;
  void   eval(Point&, Scalar u, Scalar v) const ;
  Point* operator() (const Scalar& u, const Scalar& v) const ;


  Scalar umin(void) const { return m_umin ; }
  Scalar umax(void) const { return m_umax ; }
  Scalar vmin(void) const { return m_vmin ; }
  Scalar vmax(void) const { return m_vmax ; }

  void get_range(double & umin, double & umax, double & vmin, double & vmax ) const;
  void set_range(Scalar umin, Scalar umax, Scalar vmin, Scalar vmax) ;
  

  Polynomial denominator()    const { return m_w; }
  Polynomial numerator(int i) const { return m_p[i]; }


private:
  Scalar m_umin, m_umax, m_vmin, m_vmax;
    
  Polynomial m_w, m_p[N];

} ;


TMPL typename SELF::Point* 
SELF::eval(Scalar u, Scalar v) const {
  Scalar w = m_w(u,v);
  Point* p = new Point(m_p[0](u,v)/w, m_p[1](u,v)/w, m_p[2](u,v)/w);
  return p ;
}

TMPL void 
SELF::eval(Point& p, Scalar u, Scalar v) const {
  Scalar w=m_w(u,v);
  p.setx(m_p[0](u,v)/w);
  p.sety(m_p[1](u,v)/w);
  p.setz(m_p[2](u,v)/w);
}

TMPL typename SELF::Point* 
SELF::operator() (const Scalar& u, const Scalar& v) const {
  return this->eval(u,v); 
}

TMPL void 
SELF::set_range(Scalar umin, Scalar umax, Scalar vmin, Scalar vmax) {
  m_umin=umin;                                           
  m_umax=umax;                                           
  m_vmin=vmin;                                           
  m_vmax=vmax;                                           
}
TMPL void 
SELF::get_range(Scalar& umin, Scalar& umax, Scalar& vmin, Scalar& vmax) const {
  umin=m_umin;                                           
  umax=m_umax;                                           
  vmin=m_vmin;                                           
  vmax=m_vmax;                                           
}

//--------------------------------------------------------------------
TMPL1 struct viewer;   
TMPL AXEL&
operator<<(AXEL& out, const SELF& c) {
  out<<"\n <surface type=\"rational\">\n";
  out<<"   <domain>"<< c.umin()<<" "<<c.umax()<<" "<<c.vmin()<<" "<<c.vmax()<<"</domain>\n";
  for(int i=0; i<N;i++){
    out<<"   <polynomial>";
    print(out,c.numerator(i),variables("t"));
    out<<"</polynomial>\n";
  }
  out<<"   <polynomial>";
  print(out,c.denominator(),variables("t"));
  out<<"</polynomial>\n";
  out<<" </surface>\n";
  return out;
}
//====================================================================
} // namespace shape
//====================================================================
} // namespace mmx
//====================================================================
# undef TMPL
# undef TMPL1
# undef TMPL_DEF
# undef Scalar
# undef Point
# undef SELF
# undef ParametricSurface
# undef AXEL
# endif // RATIONALSURFACE_H