ParamCurveInt.h

//===========================================================================
// GoTools - SINTEF Geometry Tools
//
// GoTools module: Intersections, version 1.0
//
// Copyright (C) 2000-2007 SINTEF ICT, Applied Mathematics, Norway.
//
// This program is free software; you can redistribute it and/or          
// modify it under the terms of the GNU General Public License            
// as published by the Free Software Foundation version 2 of the License. 
//
// This program is distributed in the hope that it will be useful,        
// but WITHOUT ANY WARRANTY; without even the implied warranty of         
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the          
// GNU General Public License for more details.                           
//
// You should have received a copy of the GNU General Public License      
// along with this program; if not, write to the Free Software            
// Foundation, Inc.,                                                      
// 59 Temple Place - Suite 330,                                           
// Boston, MA  02111-1307, USA.                                           
//
// Contact information: E-mail: tor.dokken@sintef.no                      
// SINTEF ICT, Department of Applied Mathematics,                         
// P.O. Box 124 Blindern,                                                 
// 0314 Oslo, Norway.                                                     
//
// Other licenses are also available for this software, notably licenses
// for:
// - Building commercial software.                                        
// - Building software whose source code you wish to keep private.        
//===========================================================================
#ifndef _PARAMCURVEINT_H
#define _PARAMCURVEINT_H

#include "ParamGeomInt.h"
#include "ParamCurve.h"


namespace Go {



class RotatedBox;



class ParamCurveInt : public ParamGeomInt {
public:
    explicit ParamCurveInt(boost::shared_ptr<ParamCurve> curve,
                           ParamGeomInt* parent = 0); 

    virtual ~ParamCurveInt() {};

    virtual void point(Point& res, const double *par) const
    { curve_->point(res, par[0]); }

    virtual void point(std::vector<Point>& res, 
                       const double* par, 
                       int der,
                       const bool* from_right = 0,
                       double resolution = 1.0e-12) const 
    {
        (from_right == 0)
            ? curve_->point(res, *par, der, true)
            : curve_->point(res, *par, der, from_right[0]);
    }

    virtual ParamCurveInt* getParamCurveInt();

    virtual ParamSurfaceInt* getParamSurfaceInt()
    { return NULL; }

    boost::shared_ptr<ParamCurve> getParamCurve();

    boost::shared_ptr<const ParamCurve> getParamCurve() const;

    boost::shared_ptr<ParamCurve> 
    getParentParamCurve(double& start, double& end);

    boost::shared_ptr<ParamCurve> getParentParamCurve();

    virtual boost::shared_ptr<ParamCurveInt> 
    makeIntObject(boost::shared_ptr<ParamCurve> curve);
    
    virtual int numParams() const;

    virtual void getLengthAndWiggle(double *length, double *wiggle);

    virtual bool hasInnerKnots(int pardir) const;

    virtual bool hasCriticalVals(int pardir) const;

    virtual void setCriticalVal(int pardir, double par);

    virtual bool hasCriticalValsOrKnots(int pardir) const;

    virtual bool canDivide(int pardir);

    virtual std::vector<double> getCriticalVals(int pardir) const;

    virtual std::vector<double> getInnerKnotVals(int pardir,
                                                 bool sort = false) const;

    virtual std::vector<double> getCriticalValsAndKnots(int pardir) const;

    virtual int getMeshSize(int dir);

    virtual double paramFromMesh(int dir, int idx);

    virtual std::vector<double>::iterator getMesh();

    virtual double startParam(int pardir) const
    { return curve_->startparam(); }

    virtual double endParam(int pardir) const
    { return curve_->endparam(); }

    virtual bool boundaryPoint(const double* par, double eps) const; 

    virtual void
    subdivide(int pardir, double par, 
              std::vector<boost::shared_ptr<ParamGeomInt> >& subdiv_objs,
              std::vector<boost::shared_ptr<ParamGeomInt> >& bd_objs);

    virtual CompositeBox compositeBox() const;

    virtual DirectionCone directionCone() const;

    virtual void 
    getBoundaryObjects(std::vector<boost::
                       shared_ptr<BoundaryGeomInt> >& bd_objs);

    virtual int checkPeriodicity(int pardir = 0) const;

    virtual int dimension() const
    { return dim_; }

    double startparam() const 
    { return curve_->startparam(); }

    double endparam() const
    { return curve_->endparam(); }

    // Functions used from coincidence checking

    void assureInRange(double& t);
 
    virtual int knotIntervalFuzzy(double& t, double tol) const; 
    
    virtual double nextSegmentVal(double par,
                                  bool forward, double tol) const;

    virtual bool isDegenerate(double epsge, int dir, double *par);

    virtual bool isSpline();

    virtual const SplineCurve* getSpline()
        { return 0; }

    virtual double getOptimizedConeAngle(Point& axis1, Point& axis2);

    virtual RotatedBox getRotatedBox(std::vector<Point>& axis) const;

    void axisFromEndpts(Point& axis) const;

protected:
    // Data members
    boost::shared_ptr<ParamCurve> curve_;

    int dim_; // Space dimension.

    std::vector<std::pair<double, int> > segment_;

    // Approximating polygon
    mutable std::vector<double> mesh_;

    mutable bool lw_set_;
    mutable double length_, wiggle_;

private:
    // Make approximating polygon
    void makeMesh(int size) const;

};


} // namespace Go


#endif // _PARAMCURVEINT_H

---