closestPtCurves2D.h

Go to the documentation of this file.

//===========================================================================
// 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 _CLOSESTPTCURVES2D_H
#define _CLOSESTPTCURVES2D_H

#include "ParamCurve.h"

namespace Go {


void closestPtCurves2D(ParamCurve* cv1, ParamCurve* cv2, double aepsge,
                       double tmin1, double tmax1, double tmin2, double tmax2,
                       double seed1, double seed2, int method, bool quick, 
                       double& par1, double& par2, double& dist, 
                       Point& ptc1, Point& ptc2, int& istat);


// Anonymous namespace
namespace {

void insideParamDomain2D(Point& gd, const Point& acoef, double astart1, 
                         double aend1, double astart2,double aend2, int& corr);


void nextStep2D(double& dist, Point& diff, Point& delta, double& det,
                int& kstat, std::vector<Point>& eval1, 
                std::vector<Point>& eval2, int method);


void makeEqSys1(double A[], double b[], const std::vector<Point>& eval1,
                const std::vector<Point>&eval2, const Point& diff);


void makeEqSys2(double A[], double b[], const std::vector<Point>& eval1,
                const std::vector<Point>&eval2, const Point& diff);


void makeEqSys3(double A[], double b[], const std::vector<Point>& eval1,
                const std::vector<Point>&eval2, const Point& diff);


int localPretop(double dist,const Point& diff,
                const std::vector<Point>& eval1,
                const std::vector<Point>& eval2);


void singular(ParamCurve* cv1, ParamCurve* cv2,Point& par_val,double& dist,
              int quick,double aepsge,double delta,const Point& diff,
              const std::vector<Point>&eval1,const std::vector<Point>&eval2,
              double astart1,double astart2,double aend1,double aend2);


void secant2D(ParamCurve *pcurve1,ParamCurve *pcurve2,Point& par_val,
              double& dist,int& jstat, double delta,double aepsge,
              double astart1,double astart2,double aend1,double aend2);


void setReturnValues(const Point& par_val,ParamCurve* cv1,ParamCurve* cv2,
                     double aepsge,  int& jstat,double& par1,double& par2,
                     double& dist,Point& ptc1,Point& ptc2);


void newPointEval(ParamCurve *pcurve1,double par1,ParamCurve *pcurve2,
                  double astart2,double aend2,double aepsge,
                  double& par2,double& y,double& dist,int& jstat);

} // Anonymous namespace

 
} // namespace Go

#endif // _CLOSESTPTCURVES2D_H

---