Param2FunctionInt.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 _PARAM2FUNCTIONINT_H
#define _PARAM2FUNCTIONINT_H


#include "ParamFunctionInt.h"
#include "ParamSurface.h"


namespace Go {



class SplineSurface;



class Param2FunctionInt : public ParamFunctionInt {
public:
    explicit Param2FunctionInt(boost::shared_ptr<ParamSurface> surf);

    explicit Param2FunctionInt(boost::shared_ptr<ParamSurface> surf,
                               Param2FunctionInt *parent);

    virtual ~Param2FunctionInt() {};

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

    virtual void point(std::vector<Point>& pt, 
                       const double* tpar, 
                       int der,
                       const bool* from_right = 0,
                       double resolution = 1.0e-12) const
    {
        if (from_right == 0) { // User did not specify direction of
                               // differentiation
            surf_->point(pt, tpar[0], tpar[1], der, true, true);
        } else {
            surf_->point(pt, tpar[0], tpar[1], der,
                         from_right[0], from_right[1], resolution);
        } 
    }

    virtual Param2FunctionInt* getParam2FunctionInt();

    boost::shared_ptr<ParamSurface> getParamSurface();
    boost::shared_ptr<const ParamSurface> getParamSurface() const;

    boost::shared_ptr<ParamSurface>
    getParentParamSurface(RectDomain& domain);
    boost::shared_ptr<ParamSurface> getParentParamSurface();

    virtual boost::shared_ptr<Param2FunctionInt> 
    makeIntFunction(boost::shared_ptr<ParamSurface> surf);
    
    virtual int numParams() const;

    boost::shared_ptr<ParamCurve> 
    getIsoCurve(double param_start, double param_end,
                double isoval, bool pardir_is_u) const;

    boost::shared_ptr<ParamCurve>
        getConstantParameterCurve(int dir, double par);

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

    virtual bool hasInnerKnots(int pardir) const;

    virtual bool hasCriticalVals(int pardir) const;

    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();

    // @@sbr A trimmed sf need not have a startparam/endparam ...
    virtual double startParam(int pardir) const;

    virtual double endParam(int pardir) const;

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

    std::vector<double> getMima() const;

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

    // We want to know if there exists par direction in which the
    // function is monotone (guarantees at most one solution for a
    // continuous function).  We also would like to know in what
    // direction the function is monotone (will then move in
    // orthogonal direction).
    virtual bool monotone(Point& dir, double tol=1.0e-15) const;

    virtual CompositeBox compositeBox() const;

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

    int dimension()
    { return dim_; }
    
    double isolateDegPar(int dir, int deg_edge, double threshold);

    // Functions used from coincidence checking

    void assureInRange(int pardir, double& t);
 
    // Find the interval in which t lies, if t is within tol of an
    // existing knot, t will be changed to that knot.

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

    bool isDegenerate(double epsge, int pardir);

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

    boost::shared_ptr<ParamSurface> getSurface()
    { return surf_; }

    boost::shared_ptr<const ParamSurface> getSurface() const
    { return surf_;}

    void derivs(double u, double v,
                Point& deriv_u, Point& deriv_v) const
    {
        surf_->point(temp_point_array_, u, v, 1);
        deriv_u = temp_point_array_[1];
        deriv_v = temp_point_array_[2];
    }

protected:
    // Data members

    boost::shared_ptr<ParamSurface> surf_; // Our 2D-parametric 1D
                                           // function (R^2->R).

    int dim_; // Space dimension. 1.

    Param2FunctionInt *parentfunction_; // Probably undivided sf, to
                                        // avoid num error.

    std::vector<std::pair<double, int> > segment_[2];

    // pardir (0 || 1) is that opposite bdidx (i.e. dir of constant
    // parameter tpar along edge given by bdidx).
    Go::SplineCurve* extractBdCurve(const Go::SplineSurface&, int bdidx,
                                    int& pardir, double& tpar) const;

    // Generated data
    // @@sbr Possibly include more of the computed datas. For instance
    // structure conc monotonicity.
    bool bd_deg_[4];  // Degeneracy, sequence left, right, bottom, top
    double deg_tol_;
    RectDomain domain_;
    // Approximating mesh
    mutable int nmesh_[2];
    mutable std::vector<double> mesh_;

    mutable std::vector<Point> temp_point_array_;

};


} // namespace Go


#endif // _PARAM2FUNCTIONINT_H

---