OCCT

// Created on: 1996-07-31
// Created by: Herve LOUESSARD
// Copyright (c) 1996-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 <LocalAnalysis.hxx>
#include <LocalAnalysis_SurfaceContinuity.hxx>
#include <LocalAnalysis_CurveContinuity.hxx>
#include <Geom2d_TrimmedCurve.hxx>
#include <DrawTrSurf.hxx>
#include <GeometryTest.hxx>
#include <GeomAbs_Shape.hxx>
#include <Draw_Appli.hxx>
#include <Draw.hxx>
#include <Draw_Interpretor.hxx>
#include <Message.hxx>

/*********************************************************************************/
void InitEpsSurf(double& epsnl,
                 double& epsdis,
                 double& epsangk1,
                 double& epsangk2,
                 double& epsangn1,
                 double& perce,
                 double& maxlen)
{
  epsnl    = 0.001;
  epsdis   = 0.001;
  epsangk1 = 0.001;
  epsangk2 = 0.001;
  epsangn1 = 0.001;
  perce    = 0.01;
  maxlen   = 10000;
}

/*********************************************************************************/
void InitEpsCurv(double& epsnl,
                 double& epsdis,
                 double& epsangk1,
                 double& epsangk2,
                 double& epsangg1,
                 double& epsangg2,
                 double& percent,
                 double& maxlen)
{
  epsnl    = 0.001;
  epsdis   = 0.001;
  epsangk1 = 0.001;
  epsangk2 = 0.001;
  epsangg1 = 0.001;
  epsangg2 = 0.001;
  percent  = 0.01;
  maxlen   = 10000;
}

/*********************************************************************************/

static int surfaceCcontinuity(Draw_Interpretor& di, int n, const char** a)
{
  double        U1, U2, V1, V2, u1, u2, v1, v2;
  GeomAbs_Shape TypeCont;
  int           ord;
  bool          b1, b2, b3, b4;
  double        epsnl, epsC0, epsC1, epsC2, epsG1, maxlen, perce;
  if (n < 8)
  {
    return 1;
  }

  InitEpsSurf(epsnl, epsC0, epsC1, epsC2, epsG1, perce, maxlen);

  occ::handle<Geom_Surface> surf1 = DrawTrSurf::GetSurface(a[2]);
  if (surf1.IsNull())
  {
    return 1;
  }

  occ::handle<Geom_Surface> surf2 = DrawTrSurf::GetSurface(a[5]);
  if (surf2.IsNull())
  {
    return 1;
  }

  U1 = Draw::Atof(a[3]);
  U2 = Draw::Atof(a[6]);
  V1 = Draw::Atof(a[4]);
  V2 = Draw::Atof(a[7]);

  ord = Draw::Atoi(a[1]);

  surf1->Bounds(u1, u2, v1, v2);
  b1 = ((((U1 >= u1) && (U1 <= u2)) || ((U1 <= u1) && (U1 >= u2))));
  b3 = ((((V1 >= v1) && (V1 <= v2)) || ((V1 <= v1) && (V1 >= v2))));

  surf2->Bounds(u1, u2, v1, v2);
  b2 = ((((U2 >= u1) && (U2 <= u2)) || ((U2 <= u1) && (U2 >= u2))));
  b4 = ((((V2 >= v1) && (V2 <= v2)) || ((V2 <= v1) && (V2 >= v2))));

  if (!((b1 && b2) && (b3 && b4)))
  {
    return 1;
  }

  switch (ord)
  {
    case 0: {
      TypeCont = GeomAbs_C0;
      switch (n)
      {
        case 10:
          epsC0 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsnl = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 8 - 10 parameters";
          return 1;
      }
    }
    break;

    case 1: {
      TypeCont = GeomAbs_C1;
      switch (n)
      {
        case 11:
          epsC1 = Draw::Atof(a[10]);
          [[fallthrough]];
        case 10:
          epsC0 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsnl = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 8 - 11 parameters";
          return 1;
      }
    }
    break;

    case 2: {
      TypeCont = GeomAbs_C2;
      switch (n)
      {
        case 12:
          epsC2 = Draw::Atof(a[11]);
          [[fallthrough]];
        case 11:
          epsC1 = Draw::Atof(a[10]);
          [[fallthrough]];
        case 10:
          epsC0 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsnl = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 8 - 12 parameters";
          return 1;
      }
    }
    break;
    default:
      Message::SendFail() << "Error: invalid value of parameter 1 (" << a[1]
                          << "): should be 0, 1, or 2";
      return 1;
  }

  LocalAnalysis_SurfaceContinuity
    RES(surf1, U1, V1, surf2, U2, V2, TypeCont, epsnl, epsC0, epsC1, epsC2, epsG1, perce, maxlen);
  // LocalAnalysis::Dump(RES,std::cout);
  Standard_SStream aSStream;
  LocalAnalysis::Dump(RES, aSStream);
  di << aSStream;
  return 0;
}

/*********************************************************************************/

static int surfaceGcontinuity(Draw_Interpretor& di, int n, const char** a)
{
  double        U1, U2, V1, V2, u1, u2, v1, v2;
  GeomAbs_Shape TypeCont;
  int           ord;
  bool          b1, b2, b3, b4;
  double        epsnl, epsC0, epsC1, epsC2, epsG1, maxlen, perce;

  if (n < 8)
  {
    return 1;
  }
  InitEpsSurf(epsnl, epsC0, epsC1, epsC2, epsG1, perce, maxlen);

  occ::handle<Geom_Surface> surf1 = DrawTrSurf::GetSurface(a[2]);
  if (surf1.IsNull())
  {
    return 1;
  }

  occ::handle<Geom_Surface> surf2 = DrawTrSurf::GetSurface(a[5]);
  if (surf2.IsNull())
  {
    return 1;
  }

  U1 = Draw::Atof(a[3]);
  U2 = Draw::Atof(a[6]);
  V1 = Draw::Atof(a[4]);
  V2 = Draw::Atof(a[7]);

  ord = Draw::Atoi(a[1]);

  surf1->Bounds(u1, u2, v1, v2);
  b1 = ((((U1 >= u1) && (U1 <= u2)) || ((U1 <= u1) && (U1 >= u2))));
  b3 = ((((V1 >= v1) && (V1 <= v2)) || ((V1 <= v1) && (V1 >= v2))));

  surf2->Bounds(u1, u2, v1, v2);
  b2 = ((((U2 >= u1) && (U2 <= u2)) || ((U2 <= u1) && (U2 >= u2))));
  b4 = ((((V2 >= v1) && (V2 <= v2)) || ((V2 <= v1) && (V2 >= v2))));

  if (!((b1 && b2) && (b3 && b4)))
  {
    return 1;
  }

  switch (ord)
  {
    case 1: {
      TypeCont = GeomAbs_G1;
      switch (n)
      {
        case 11:
          epsG1 = Draw::Atof(a[10]);
          [[fallthrough]];
        case 10:
          epsC0 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsnl = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 8 - 11 parameters";
          return 1;
      }
    }
    break;

    case 2: {
      TypeCont = GeomAbs_G2;
      switch (n)
      {
        case 13:
          maxlen = Draw::Atof(a[12]);
          [[fallthrough]];
        case 12:
          perce = Draw::Atof(a[11]);
          [[fallthrough]];
        case 11:
          epsG1 = Draw::Atof(a[10]);
          [[fallthrough]];
        case 10:
          epsC0 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsnl = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 8 - 13 parameters";
          return 1;
      }
    }
    break;
    default:
      Message::SendFail() << "Error: invalid value of parameter 1 (" << a[1]
                          << "): should be 1 or 2";
      return 1;
  }

  LocalAnalysis_SurfaceContinuity
    RES(surf1, U1, V1, surf2, U2, V2, TypeCont, epsnl, epsC0, epsC1, epsC2, epsG1, perce, maxlen);
  // LocalAnalysis::Dump(RES,std::cout);
  Standard_SStream aSStream;
  LocalAnalysis::Dump(RES, aSStream);
  di << aSStream;
  return 0;
}

/*********************************************************************************/

static int curveGcontinuity(Draw_Interpretor& di, int n, const char** a)
{
  double        U1, U2, u1, u2;
  GeomAbs_Shape TypeCont;
  int           ord;
  bool          b1, b2;
  double        epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, maxlen, percent;
  if (n < 6)
  {
    return 1;
  }

  occ::handle<Geom_Curve> curv1 = DrawTrSurf::GetCurve(a[2]);
  if (curv1.IsNull())
  {
    return 1;
  }
  occ::handle<Geom_Curve> curv2 = DrawTrSurf::GetCurve(a[4]);
  if (curv2.IsNull())
  {
    return 1;
  }

  U1 = Draw::Atof(a[3]);
  U2 = Draw::Atof(a[5]);

  ord = Draw::Atoi(a[1]);

  u1 = curv1->FirstParameter();
  u2 = curv1->LastParameter();
  b1 = ((((U1 >= u1) && (U1 <= u2)) || ((U1 <= u1) && (U1 >= u2))));

  u1 = curv2->FirstParameter();
  u2 = curv2->LastParameter();
  b2 = ((((U2 >= u1) && (U2 <= u2)) || ((U2 <= u1) && (U2 >= u2))));

  if (!(b1 && b2))
  {
    return 1;
  }

  InitEpsCurv(epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, percent, maxlen);
  switch (ord)
  {
    case 1: {
      TypeCont = GeomAbs_G1;
      switch (n)
      {
        case 9:
          epsG1 = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          epsC0 = Draw::Atof(a[7]);
          [[fallthrough]];
        case 7:
          epsnl = Draw::Atof(a[6]);
          [[fallthrough]];
        case 6:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 6 - 8 parameters";
          return 1;
      }
    }
    break;
    case 2: {
      TypeCont = GeomAbs_G2;
      switch (n)
      {
        case 12:
          maxlen = Draw::Atof(a[11]);
          [[fallthrough]];
        case 11:
          percent = Draw::Atof(a[10]);
          [[fallthrough]];
        case 10:
          epsG2 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsG1 = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          epsC0 = Draw::Atof(a[7]);
          [[fallthrough]];
        case 7:
          epsnl = Draw::Atof(a[6]);
          [[fallthrough]];
        case 6:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 6 - 12 parameters";
          return 1;
      }
    }
    break;

    default:
      Message::SendFail() << "Error: invalid value of parameter 1 (" << a[1]
                          << "): should be 1 or 2";
      return 1;
  }

  LocalAnalysis_CurveContinuity
    RES(curv1, U1, curv2, U2, TypeCont, epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, percent, maxlen);
  // LocalAnalysis::Dump(RES,std::cout);
  Standard_SStream aSStream;
  LocalAnalysis::Dump(RES, aSStream);
  di << aSStream;
  return 0;
}

/*********************************************************************************/

static int curveCcontinuity(Draw_Interpretor& di, int n, const char** a)
{
  double        U1, U2, u1, u2;
  GeomAbs_Shape TypeCont;
  int           ord;
  bool          b1, b2;
  double        epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, maxlen, percent;
  if (n < 6)
  {
    return 1;
  }

  occ::handle<Geom_Curve> curv1 = DrawTrSurf::GetCurve(a[2]);
  if (curv1.IsNull())
  {
    return 1;
  }
  occ::handle<Geom_Curve> curv2 = DrawTrSurf::GetCurve(a[4]);
  if (curv2.IsNull())
  {
    return 1;
  }

  U1 = Draw::Atof(a[3]);
  U2 = Draw::Atof(a[5]);

  ord = Draw::Atoi(a[1]);

  u1 = curv1->FirstParameter();
  u2 = curv1->LastParameter();
  b1 = ((((U1 >= u1) && (U1 <= u2)) || ((U1 <= u1) && (U1 >= u2))));

  u1 = curv2->FirstParameter();
  u2 = curv2->LastParameter();
  b2 = ((((U2 >= u1) && (U2 <= u2)) || ((U2 <= u1) && (U2 >= u2))));

  if (!(b1 && b2))
  {
    return 1;
  }

  InitEpsCurv(epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, percent, maxlen);
  switch (ord)
  {
    case 0: {
      TypeCont = GeomAbs_C0;
      switch (n)
      {
        case 8:
          epsC0 = Draw::Atof(a[7]);
          [[fallthrough]];
        case 7:
          epsnl = Draw::Atof(a[6]);
          [[fallthrough]];
        case 6:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 6 - 8 parameters";
          return 1;
      }
    }
    break;
    case 1: {
      TypeCont = GeomAbs_C1;
      switch (n)
      {
        case 9:
          epsC1 = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          epsC0 = Draw::Atof(a[7]);
          [[fallthrough]];
        case 7:
          epsnl = Draw::Atof(a[6]);
          [[fallthrough]];
        case 6:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 6 - 9 parameters";
          return 1;
      }
    }
    break;

    case 2: {
      TypeCont = GeomAbs_C2;
      switch (n)
      {
        case 10:
          epsC2 = Draw::Atof(a[9]);
          [[fallthrough]];
        case 9:
          epsC1 = Draw::Atof(a[8]);
          [[fallthrough]];
        case 8:
          epsC0 = Draw::Atof(a[7]);
          [[fallthrough]];
        case 7:
          epsnl = Draw::Atof(a[6]);
          [[fallthrough]];
        case 6:
          break;
        default:
          Message::SendFail()
            << "Error: invalid number of arguments: expected to get 6 - 10 parameters";
          return 1;
      }
    }
    break;

    default:
      Message::SendFail() << "Error: invalid value of parameter 1 (" << a[1]
                          << "): should be 0, 1, or 2";
      return 1;
  }

  LocalAnalysis_CurveContinuity
    RES(curv1, U1, curv2, U2, TypeCont, epsnl, epsC0, epsC1, epsC2, epsG1, epsG2, percent, maxlen);
  // LocalAnalysis::Dump(RES,std::cout);
  Standard_SStream aSStream;
  LocalAnalysis::Dump(RES, aSStream);
  di << aSStream;
  return 0;
}

/***************************************************************************/
void GeometryTest::ContinuityCommands(Draw_Interpretor& theCommands)
{
  static bool loaded = false;
  if (loaded)
  {
    return;
  }
  loaded = true;

  DrawTrSurf::BasicCommands(theCommands);
  const char* g;

  g = "GEOMETRY curves and surfaces continuity analysis ";

  theCommands.Add("surfaceCcontinuity",
                  " surfaceCcontinuity   order surf1 parU1 parV1 surf2 parU2 parV2  [eps_nul[ "
                  "epsC0 [epsC1 [epsC2]]]]",
                  __FILE__,
                  surfaceCcontinuity,
                  g);

  theCommands.Add("surfaceGcontinuity",
                  " surfaceGcontinuity   order surf1 parU1 parV1 surf2 parU2 parV2  [eps_nul[ "
                  "epsG0 [epsG1[percent [maxlen]]]]",
                  __FILE__,
                  surfaceGcontinuity,
                  g);

  theCommands.Add(
    "curveCcontinuity",
    "curveCcontinuity   order curv1 u1  curv2  u2   [epsnul  [epsC0 [epsC1  [epsC2 ]]]]  ",
    __FILE__,
    curveCcontinuity,
    g);

  theCommands.Add("curveGcontinuity",
                  "curveGcontinuity   order  curv1 u1  curv2  u2   [epsnul  [epsG0  [epsG1  [epsG2 "
                  "[percent  [maxlen ]]]]]] ",
                  __FILE__,
                  curveGcontinuity,
                  g);
}