src/CPnts/CPnts_MyRootFunction.cxx

   1 // Copyright (c) 1995-1999 Matra Datavision
   2 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   3 //
   4 // This file is part of Open CASCADE Technology software library.
   5 //
   6 // This library is free software; you can redistribute it and/or modify it under
   7 // the terms of the GNU Lesser General Public License version 2.1 as published
   8 // by the Free Software Foundation, with special exception defined in the file
   9 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  10 // distribution for complete text of the license and disclaimer of any warranty.
  11 //
  12 // Alternatively, this file may be used under the terms of Open CASCADE
  13 // commercial license or contractual agreement.
  14
  15
  16 #include <CPnts_MyRootFunction.hxx>
  17 #include <math_GaussSingleIntegration.hxx>
  18 #include <Standard_DomainError.hxx>
  19
  20 void CPnts_MyRootFunction::Init(const CPnts_RealFunction& F,
  21                                 const Standard_Address D,
  22                                 const Standard_Integer Order)
  23 {
  24   myFunction.Init(F,D);
  25   myOrder = Order;
  26 }
  27
  28 void CPnts_MyRootFunction::Init(const Standard_Real X0,
  29                                 const Standard_Real L)
  30 {
  31   myX0 = X0;
  32   myL = L;
  33   myTol = -1; //to supress the tolerance
  34 }
  35
  36 void CPnts_MyRootFunction::Init(const Standard_Real X0,
  37                                 const Standard_Real L,
  38                                 const Standard_Real Tol)
  39 {
  40   myX0 = X0;
  41   myL = L;
  42   myTol = Tol;
  43 }
  44
  45 Standard_Boolean CPnts_MyRootFunction::Value(const Standard_Real X,
  46                                              Standard_Real& F)
  47 {
  48   math_GaussSingleIntegration Length;
  49
  50   if (myTol <= 0) Length = math_GaussSingleIntegration(myFunction, myX0, X, myOrder);
  51   else Length = math_GaussSingleIntegration(myFunction, myX0, X, myOrder, myTol);
  52
  53   if (Length.IsDone()){
  54     F= Length.Value() - myL;
  55     return Standard_True;
  56   }
  57   else {
  58     return Standard_False;
  59   }
  60 }
  61
  62 Standard_Boolean CPnts_MyRootFunction::Derivative(const Standard_Real X,
  63                                                   Standard_Real& Df)
  64 {
  65   return myFunction.Value(X,Df);
  66 }
  67
  68 Standard_Boolean CPnts_MyRootFunction::Values(const Standard_Real X,
  69                                               Standard_Real& F,
  70                                               Standard_Real& Df)
  71 {
  72   math_GaussSingleIntegration Length;
  73
  74   if (myTol <= 0) Length = math_GaussSingleIntegration(myFunction, myX0, X, myOrder);
  75   else Length = math_GaussSingleIntegration(myFunction, myX0, X, myOrder, myTol);
  76
  77   if (Length.IsDone()){
  78     F= Length.Value() - myL;
  79     return myFunction.Value(X,Df);
  80   }
  81   else {
  82     return Standard_False;
  83   }
  84 }