ApproxCurve.h

//===========================================================================
// GoTools - SINTEF Geometry Tools version 1.1
//
// GoTools module: CORE
//
// 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 _GOAPPROXCURVE_H_
#define _GOAPPROXCURVE_H_

//   -----------------------------------------------------------------------
//      Interface file for class ApproxCurve
//   -----------------------------------------------------------------------
//
//       Approximate a set of points by a B-spline curve to
//       satisfy a given accuracy
//
//       Implementation of the member functions are given in the
//       following files:
//
//          1. ApproxCurve.C
//
//   -----------------------------------------------------------------------
//    Written by: Vibeke Skytt                           04-02
//   -----------------------------------------------------------------------

#include "SplineCurve.h"
#include <vector>
#include "SmoothCurve.h"

namespace Go
{

class ApproxCurve
{
public:
    ApproxCurve(const std::vector<double>& points, 
                  const std::vector<double>& parvals,
                  int dim, double aepsge);

    ApproxCurve(const std::vector<double>& points, 
                  const std::vector<double>& parvals,
                  int dim, double aepsge, int in, int ik);

    ApproxCurve(const std::vector<double>& points, 
                  const std::vector<double>& parvals, int dim,
                  double aepsge, int in, int ik,
                  std::vector<double>& knots);

    ~ApproxCurve();

    void setSmooth(double w);


    void unsetSmooth();


    void setEndPoints(const std::vector<Point>& start_point, 
                      const std::vector<Point>& end_point);

    boost::shared_ptr<SplineCurve> getApproxCurve(double& maxdist, 
                                                  double& avdist,
                                                  int max_iter = 5);
protected:
    ApproxCurve();

private:
  boost::shared_ptr<SplineCurve> curr_crv_;
  double maxdist_;
  double avdist_;
  double aepsge_;
  double smoothweight_;
  double smoothfac_;

  int dim_;
  std::vector<double> points_;
  std::vector<double> parvals_;
  std::vector<Point> start_pt_; // Pt, der. May be empty.
  std::vector<Point> end_pt_; // Pt, der. May be empty.

  void makeInitCurve();

  void makeInitCurve(int in, int ik, std::vector<double> knots);

  void makeInitCurve(int in, int ik);

  void adjustSmoothWeight();

  void makeSmoothCurve();

  void checkAccuracy(std::vector<double>& newknots, int uniform=1);

  int doApprox(int max_iter);

  void getConstraints(std::vector<sideConstraint>& pt_constraints,
                      std::vector<sideConstraint>& tangent_constraints);

};
}

#endif

---