// Created on: 1991-03-28
// Created by: Jacques GOUSSARD
// Copyright (c) 1991-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.

#ifndef _IntWalk_TheInt2S_HeaderFile
#define _IntWalk_TheInt2S_HeaderFile

#include "Adaptor3d_Surface.hxx"
#include "IntSurf_PntOn2S.hxx"
#include "gp_Dir.hxx"
#include "gp_Dir2d.hxx"
#include "IntWalk_TheFunctionOfTheInt2S.hxx"
#include "Standard_Real.hxx"
#include "NCollection_Array1.hxx"
#include "IntImp_ConstIsoparametric.hxx"

class Adaptor3d_HSurfaceTool;
class math_FunctionSetRoot;

class IntWalk_TheInt2S
{
public:
  DEFINE_STANDARD_ALLOC

  //! compute the solution point with the close point
  Standard_EXPORT IntWalk_TheInt2S(const NCollection_Array1<double>&     Param,
                                   const occ::handle<Adaptor3d_Surface>& S1,
                                   const occ::handle<Adaptor3d_Surface>& S2,
                                   const double                          TolTangency);

  //! initialize the parameters to compute the solution point
  //! it 's possible to write to optimize:
  //! IntImp_Int2S inter(S1,S2,Func,TolTangency);
  //! math_FunctionSetRoot rsnld(inter.Function());
  //! while ...{
  //! Param(1)=...
  //! Param(2)=...
  //! param(3)=...
  //! inter.Perform(Param,rsnld);
  //! }
  Standard_EXPORT IntWalk_TheInt2S(const occ::handle<Adaptor3d_Surface>& S1,
                                   const occ::handle<Adaptor3d_Surface>& S2,
                                   const double                          TolTangency);

  //! returns the best constant isoparametric to find
  //! the next intersection's point +stores the solution
  //! point (the solution point is found with the close point
  //! to intersect the isoparametric with the other patch;
  //! the choice of the isoparametic is calculated)
  Standard_EXPORT IntImp_ConstIsoparametric Perform(const NCollection_Array1<double>& Param,
                                                    math_FunctionSetRoot&             Rsnld);

  //! returns the best constant isoparametric to find
  //! the next intersection's point +stores the solution
  //! point (the solution point is found with the close point
  //! to intersect the isoparametric with the other patch;
  //! the choice of the isoparametic is given by ChoixIso)
  Standard_EXPORT IntImp_ConstIsoparametric Perform(const NCollection_Array1<double>& Param,
                                                    math_FunctionSetRoot&             Rsnld,
                                                    const IntImp_ConstIsoparametric   ChoixIso);

  //! Returns TRUE if the creation completed without failure.
  bool IsDone() const;

  //! Returns TRUE when there is no solution to the problem.
  bool IsEmpty() const;

  //! Returns the intersection point.
  const IntSurf_PntOn2S& Point() const;

  //! Returns True if the surfaces are tangent at the
  //! intersection point.
  bool IsTangent() const;

  //! Returns the tangent at the intersection line.
  const gp_Dir& Direction() const;

  //! Returns the tangent at the intersection line in the
  //! parametric space of the first surface.
  const gp_Dir2d& DirectionOnS1() const;

  //! Returns the tangent at the intersection line in the
  //! parametric space of the second surface.
  const gp_Dir2d& DirectionOnS2() const;

  //! return the math function which
  //! is used to compute the intersection
  IntWalk_TheFunctionOfTheInt2S& Function();

  //! return the intersection point which is
  //! enable for changing.
  IntSurf_PntOn2S& ChangePoint();

private:
  bool                          done;
  bool                          empty;
  IntSurf_PntOn2S               pint;
  bool                          tangent;
  gp_Dir                        d3d;
  gp_Dir2d                      d2d1;
  gp_Dir2d                      d2d2;
  IntWalk_TheFunctionOfTheInt2S myZerParFunc;
  double                        tol;
  double                        ua0;
  double                        va0;
  double                        ua1;
  double                        va1;
  double                        ub0;
  double                        vb0;
  double                        ub1;
  double                        vb1;
  double                        ures1;
  double                        ures2;
  double                        vres1;
  double                        vres2;
};

#define ThePSurface occ::handle<Adaptor3d_Surface>
#define ThePSurface_hxx <Adaptor3d_Surface.hxx>
#define ThePSurfaceTool Adaptor3d_HSurfaceTool
#define ThePSurfaceTool_hxx <Adaptor3d_HSurfaceTool.hxx>
#define IntImp_TheFunction IntWalk_TheFunctionOfTheInt2S
#define IntImp_TheFunction_hxx <IntWalk_TheFunctionOfTheInt2S.hxx>
#define IntImp_Int2S IntWalk_TheInt2S
#define IntImp_Int2S_hxx <IntWalk_TheInt2S.hxx>

#include "IntImp_Int2S.lxx"

#undef ThePSurface
#undef ThePSurface_hxx
#undef ThePSurfaceTool
#undef ThePSurfaceTool_hxx
#undef IntImp_TheFunction
#undef IntImp_TheFunction_hxx
#undef IntImp_Int2S
#undef IntImp_Int2S_hxx

#endif // _IntWalk_TheInt2S_HeaderFile