src/ChFi3d/ChFi3d_SearchSing.cxx

   1 // Created on: 1997-03-28
   2 // Created by: Philippe MANGIN
   3 // Copyright (c) 1997-1999 Matra Datavision
   4 // Copyright (c) 1999-2012 OPEN CASCADE SAS
   5 //
   6 // The content of this file is subject to the Open CASCADE Technology Public
   7 // License Version 6.5 (the "License"). You may not use the content of this file
   8 // except in compliance with the License. Please obtain a copy of the License
   9 // at http://www.opencascade.org and read it completely before using this file.
  10 //
  11 // The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
  12 // main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
  13 //
  14 // The Original Code and all software distributed under the License is
  15 // distributed on an "AS IS" basis, without warranty of any kind, and the
  16 // Initial Developer hereby disclaims all such warranties, including without
  17 // limitation, any warranties of merchantability, fitness for a particular
  18 // purpose or non-infringement. Please see the License for the specific terms
  19 // and conditions governing the rights and limitations under the License.
  20
  21
  22
  23 #include <ChFi3d_SearchSing.ixx>
  24
  25 #include <gp_Pnt.hxx>
  26 #include <gp_Vec.hxx>
  27
  28 ChFi3d_SearchSing::ChFi3d_SearchSing(const Handle(Geom_Curve)& C1,
  29                                      const Handle(Geom_Curve)& C2)
  30 {
  31   myC1 = C1;
  32   myC2 = C2;
  33 }
  34
  35
  36
  37 Standard_Boolean ChFi3d_SearchSing::Value(const Standard_Real X,
  38                                           Standard_Real& F)
  39 {
  40   gp_Pnt P1, P2;
  41   gp_Vec V1, V2;
  42   myC1->D1(X, P1, V1);
  43   myC2->D1(X, P2, V2);
  44   gp_Vec V(P1,P2);
  45   F = V * (V2-V1);
  46   return Standard_True;
  47 }
  48
  49 Standard_Boolean ChFi3d_SearchSing::Derivative(const Standard_Real X,
  50                                                Standard_Real& D )
  51 {
  52   gp_Pnt P1, P2;
  53   gp_Vec V1, V2, W1, W2;
  54   myC1->D2(X, P1, V1, W1);
  55   myC2->D2(X, P2, V2, W2);
  56   gp_Vec V(P1,P2), VPrim;
  57   VPrim = V2 -V1;
  58   D = VPrim.SquareMagnitude() + (V * (W2-W1));
  59   return Standard_True;
  60 }
  61
  62 Standard_Boolean ChFi3d_SearchSing::Values(const Standard_Real X,
  63                                            Standard_Real& F,
  64                                            Standard_Real& D )
  65 {
  66   gp_Pnt P1, P2;
  67   gp_Vec V1, V2, W1, W2;
  68   myC1->D2(X, P1, V1, W1);
  69   myC2->D2(X, P2, V2, W2);
  70   gp_Vec V(P1,P2), VPrim;
  71   VPrim = V2 -V1;
  72   F = V * VPrim;
  73   D = VPrim.SquareMagnitude() + (V * (W2-W1));
  74   return Standard_True;
  75 }