OCCT

// Created on: 1997-03-28
// Created by: Philippe MANGIN
// Copyright (c) 1997-1999 Matra Datavision
// Copyright (c) 1999-2014 OPEN CASCADE SAS
//
// This file is part of Open CASCADE Technology software library.
//
// This library is free software; you can redistribute it and/or modify it under
// the terms of the GNU Lesser General Public License version 2.1 as published
// by the Free Software Foundation, with special exception defined in the file
// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
// distribution for complete text of the license and disclaimer of any warranty.
//
// Alternatively, this file may be used under the terms of Open CASCADE
// commercial license or contractual agreement.

#include <ChFi3d_SearchSing.hxx>
#include <Geom_Curve.hxx>
#include <gp_Pnt.hxx>
#include <gp_Vec.hxx>

ChFi3d_SearchSing::ChFi3d_SearchSing(const occ::handle<Geom_Curve>& C1,
                                     const occ::handle<Geom_Curve>& C2)
{
  myC1 = C1;
  myC2 = C2;
}

bool ChFi3d_SearchSing::Value(const double X, double& F)
{
  gp_Pnt P1, P2;
  gp_Vec V1, V2;
  myC1->D1(X, P1, V1);
  myC2->D1(X, P2, V2);
  gp_Vec V(P1, P2);
  F = V * (V2 - V1);
  return true;
}

bool ChFi3d_SearchSing::Derivative(const double X, double& D)
{
  gp_Pnt P1, P2;
  gp_Vec V1, V2, W1, W2;
  myC1->D2(X, P1, V1, W1);
  myC2->D2(X, P2, V2, W2);
  gp_Vec V(P1, P2), VPrim;
  VPrim = V2 - V1;
  D     = VPrim.SquareMagnitude() + (V * (W2 - W1));
  return true;
}

bool ChFi3d_SearchSing::Values(const double X, double& F, double& D)
{
  gp_Pnt P1, P2;
  gp_Vec V1, V2, W1, W2;
  myC1->D2(X, P1, V1, W1);
  myC2->D2(X, P2, V2, W2);
  gp_Vec V(P1, P2), VPrim;
  VPrim = V2 - V1;
  F     = V * VPrim;
  D     = VPrim.SquareMagnitude() + (V * (W2 - W1));
  return true;
}