OCCT

// 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 <BRep_Tool.hxx>
#include <Geom2d_Curve.hxx>
#include <Geom_Curve.hxx>
#include <gp.hxx>
#include <ShapeAnalysis_Edge.hxx>
#include <ShapeAnalysis_TransferParameters.hxx>
#include <ShapeBuild_Edge.hxx>
#include <Standard_Type.hxx>
#include <TopoDS_Edge.hxx>
#include <TopoDS_Face.hxx>

IMPLEMENT_STANDARD_RTTIEXT(ShapeAnalysis_TransferParameters, Standard_Transient)

//=================================================================================================

ShapeAnalysis_TransferParameters::ShapeAnalysis_TransferParameters()
{
  myScale = 1.;
  myShift = 0.;
}

//=================================================================================================

ShapeAnalysis_TransferParameters::ShapeAnalysis_TransferParameters(const TopoDS_Edge& E,
                                                                   const TopoDS_Face& F)
{
  Init(E, F);
}

//=================================================================================================

void ShapeAnalysis_TransferParameters::Init(const TopoDS_Edge& E, const TopoDS_Face& F)
{
  myScale = 1.;
  myShift = 0.;
  double          l, f, l2d = 0.0, f2d = 0.0;
  TopLoc_Location L;
  myEdge = E;
  ShapeAnalysis_Edge      sae;
  occ::handle<Geom_Curve> curve3d; // = BRep_Tool::Curve (E,f,l);
  sae.Curve3d(E, curve3d, f, l, false);
  myFirst = f;
  myLast  = l;
  occ::handle<Geom2d_Curve> curve2d; // = BRep_Tool::CurveOnSurface (E, F, f2d,l2d);
  // ShapeAnalysis_Edge sae;
  if (!F.IsNull())
  { // process free edges
    sae.PCurve(E, F, curve2d, f2d, l2d, false);
  }
  myFirst2d = f2d;
  myLast2d  = l2d;
  myFace    = F;
  if (curve3d.IsNull() || curve2d.IsNull())
  {
    return;
  }

  double ln2d = l2d - f2d;
  double ln3d = l - f;
  myScale     = (ln3d <= gp::Resolution() ? 1. : ln2d / ln3d);
  myShift     = f2d - f * myScale;
}

//=================================================================================================

void ShapeAnalysis_TransferParameters::SetMaxTolerance(const double maxtol)
{
  myMaxTolerance = maxtol;
}

//=================================================================================================

occ::handle<NCollection_HSequence<double>> ShapeAnalysis_TransferParameters::Perform(
  const occ::handle<NCollection_HSequence<double>>& Params,
  const bool                                        To2d)
{
  occ::handle<NCollection_HSequence<double>> res = new NCollection_HSequence<double>;
  for (int i = 1; i <= Params->Length(); i++)
  {
    res->Append(Perform(Params->Value(i), To2d));
  }
  return res;
}

//=================================================================================================

double ShapeAnalysis_TransferParameters::Perform(const double Param, const bool To2d)
{
  double NewParam;
  if (To2d)
  {
    NewParam = myShift + Param * myScale;
  }
  else
  {
    NewParam = -myShift / myScale + Param * 1. / myScale;
  }
  return NewParam;
}

//=================================================================================================

void ShapeAnalysis_TransferParameters::TransferRange(TopoDS_Edge& newEdge,
                                                     const double prevPar,
                                                     const double currPar,
                                                     const bool   Is2d)
{
  ShapeBuild_Edge sbe;
  if (Is2d)
  {
    double span2d = myLast2d - myFirst2d;
    double tmp1, tmp2;
    if (prevPar > currPar)
    {
      tmp1 = currPar;
      tmp2 = prevPar;
    }
    else
    {
      tmp1 = prevPar;
      tmp2 = currPar;
    }
    double alpha = (tmp1 - myFirst2d) / span2d;
    double beta  = (tmp2 - myFirst2d) / span2d;
    sbe.CopyRanges(newEdge, myEdge, alpha, beta);
  }
  else
  {
    double alpha = (prevPar - myFirst) / (myLast - myFirst);
    double beta  = (currPar - myFirst) / (myLast - myFirst);
    sbe.CopyRanges(newEdge, myEdge, alpha, beta);
  }
}

//=================================================================================================

bool ShapeAnalysis_TransferParameters::IsSameRange() const
{
  return myShift == 0. && myScale == 1.;
}