OCCT

// Created on: 1993-06-07
// Created by: Jacques GOUSSARD
// Copyright (c) 1993-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 <IntPatch_ArcFunction.hxx>
#include <IntPatch_HInterTool.hxx>
#include <IntSurf_Quadric.hxx>

IntPatch_ArcFunction::IntPatch_ArcFunction() = default;

bool IntPatch_ArcFunction::Value(const double X, double& F)
{
  gp_Pnt2d p2d(myArc->Value(X));
  mySurf->D0(p2d.X(), p2d.Y(), ptsol);
  F = myQuad.Distance(ptsol);
  return true;
}

bool IntPatch_ArcFunction::Derivative(const double X, double& D)
{
  gp_Pnt2d p2d;
  gp_Vec2d d2d;
  gp_Vec   v, d1u, d1v;
  myArc->D1(X, p2d, d2d);
  mySurf->D1(p2d.X(), p2d.Y(), ptsol, d1u, d1v);
  v.SetLinearForm(d2d.X(), d1u, d2d.Y(), d1v);
  D = v.Dot(myQuad.Gradient(ptsol));
  return true;
}

bool IntPatch_ArcFunction::Values(const double X, double& F, double& D)
{
  gp_Pnt2d p2d;
  gp_Vec2d d2d;
  gp_Vec   d1u, d1v;

  gp_Vec v1, v2;
  myArc->D1(X, p2d, d2d);
  mySurf->D1(p2d.X(), p2d.Y(), ptsol, d1u, d1v);
  v1.SetLinearForm(d2d.X(), d1u, d2d.Y(), d1v);

  myQuad.ValAndGrad(ptsol, F, v2);
  D = v1.Dot(v2);
  return true;
}

int IntPatch_ArcFunction::GetStateNumber()
{
  seqpt.Append(ptsol);
  return seqpt.Length();
}

int IntPatch_ArcFunction::NbSamples() const
{
  return std::max(std::max(IntPatch_HInterTool::NbSamplesU(mySurf, 0., 0.),
                           IntPatch_HInterTool::NbSamplesV(mySurf, 0., 0.)),
                  IntPatch_HInterTool::NbSamplesOnArc(myArc));
}