#include <DrawTrSurf_Curve2d.ixx>
#include <Geom2d_Curve.hxx>
#include <Geom2d_OffsetCurve.hxx>
#include <Geom2d_Line.hxx>
#include <Geom2d_Parabola.hxx>
#include <Geom2d_Hyperbola.hxx>
#include <Geom2dAdaptor_Curve.hxx>
#include <Geom_Plane.hxx>
#include <GeomTools_Curve2dSet.hxx>
#include <Geom2dLProp_CLProps2d.hxx>
#include <gp.hxx>
#include <TColStd_Array1OfReal.hxx>
#include <Precision.hxx>

static Standard_Real DrawTrSurf_CurveLimit = 400;
//extern Standard_Boolean Draw_Bounds;
Standard_IMPORT Standard_Boolean Draw_Bounds;

//=======================================================================
//function : DrawTrSurf_Curve2d
//purpose  : 
//=======================================================================

DrawTrSurf_Curve2d::DrawTrSurf_Curve2d (const Handle(Geom2d_Curve)& C,
					const Standard_Boolean DispOrigin) :
       DrawTrSurf_Drawable (50)
{
  curv = C;
  look = Draw_vert;
  disporigin = DispOrigin ;
  dispcurvradius = Standard_False;
  radiusmax =  1.0e3;
  radiusratio = 0.1;
}


//=======================================================================
//function : DrawTrSurf_Curve2d
//purpose  : 
//=======================================================================

DrawTrSurf_Curve2d::DrawTrSurf_Curve2d (const Handle(Geom2d_Curve)& C,
					const Draw_Color& aColor,
					const Standard_Integer Discret,
					const Standard_Boolean DispOrigin,
					const Standard_Boolean DispCurvRadius,
					const Standard_Real  RadiusMax,
					const Standard_Real  RadiusRatio) :
       DrawTrSurf_Drawable (Discret)
{
  curv = C;
  look = aColor;
  disporigin = DispOrigin ;
  dispcurvradius = DispCurvRadius;
  radiusmax =  RadiusMax;
  radiusratio = RadiusRatio;
}


//=======================================================================
//function : DrawOn
//purpose  : 
//=======================================================================

void DrawTrSurf_Curve2d::DrawOn (Draw_Display& dis) const 
{
  
  Standard_Real First = curv->FirstParameter();
  Standard_Real Last  = curv->LastParameter();
  Standard_Boolean firstInf = Precision::IsNegativeInfinite(First);
  Standard_Boolean lastInf  = Precision::IsPositiveInfinite(Last);

  if (firstInf || lastInf) {
    gp_Pnt2d P1,P2;
    Standard_Real delta = 1;
    if (firstInf && lastInf) {
      do {
	delta *= 2;
	First = - delta;
	Last  =   delta;
	curv->D0(First,P1);
	curv->D0(Last,P2);
      } while (P1.Distance(P2) < DrawTrSurf_CurveLimit);
    }
    else if (firstInf) {
      curv->D0(Last,P2);
      do {
	delta *= 2;
	First = Last - delta;
	curv->D0(First,P1);
      } while (P1.Distance(P2) < DrawTrSurf_CurveLimit);
    }
    else if (lastInf) {
      curv->D0(First,P1);
      do {
	delta *= 2;
	Last = First + delta;
	curv->D0(Last,P2);
      } while (P1.Distance(P2) < DrawTrSurf_CurveLimit);
    }
  }    
  
  dis.SetColor (look);
  
  Geom2dAdaptor_Curve C2d(curv,First,Last);
  DrawCurve2dOn (C2d, dis);

// mark the origin
  if (disporigin) {
    Draw_Bounds = Standard_False;
    gp_Pnt2d p1,p2;
    gp_Vec2d v;
    C2d.D1(Last,p1,v);
    if (v.Magnitude() > gp::Resolution()) {
      Standard_Real L = 20 / dis.Zoom();
      Standard_Real H = 10 / dis.Zoom();
      gp_Dir2d d(v);
      p2.SetCoord(p1.X() - L*d.X() - H*d.Y(), p1.Y() - L*d.Y() + H*d.X());
      dis.MoveTo(p2);
      p2.SetCoord(p1.X() - L*d.X() + H*d.Y(), p1.Y() - L*d.Y() - H*d.X());
      dis.DrawTo(p1);
      dis.DrawTo(p2);
    }
    Draw_Bounds = Standard_True;
  }

// Draw the curvature Radius      
  if (dispcurvradius && (C2d.GetType() != GeomAbs_Line)) {
    Standard_Integer ii;
    Standard_Integer intrv, nbintv = C2d.NbIntervals(GeomAbs_CN);
    TColStd_Array1OfReal TI(1,nbintv+1);
    C2d.Intervals(TI,GeomAbs_CN);
    Standard_Real Resolution = 1.0e-9, Curvature;
    Geom2dLProp_CLProps2d LProp(curv, 2, Resolution);
    gp_Pnt2d P1, P2;    

    for (intrv = 1; intrv <= nbintv; intrv++) {
	Standard_Real t = TI(intrv);
	Standard_Real step = (TI(intrv+1) - t) / GetDiscretisation();
	Standard_Real LRad, ratio;
	for (ii = 1; ii <= GetDiscretisation(); ii++) {	 
	  LProp.SetParameter(t);
          if (LProp.IsTangentDefined()) {
	     Curvature = Abs(LProp.Curvature());
	     if ( Curvature >  Resolution) {
	       curv->D0(t, P1);
	       dis.MoveTo(P1);
	       LRad = 1./Curvature;
	       ratio = ( (  LRad > radiusmax) ? radiusmax/LRad : 1 );
	       ratio *= radiusratio;
	       LProp.CentreOfCurvature(P2);
	       gp_Vec2d V(P1, P2);
	       dis.DrawTo(P1.Translated(ratio*V));
	    }
	   }
	   t += step;
	}
      }
  }
}



//=======================================================================
//function : Copy
//purpose  : 
//=======================================================================
Handle(Draw_Drawable3D)  DrawTrSurf_Curve2d::Copy()const 
{
  Handle(DrawTrSurf_Curve2d) DC = new DrawTrSurf_Curve2d
    (Handle(Geom2d_Curve)::DownCast(curv->Copy()),
     look,
     GetDiscretisation());
     
  return DC;
}


//=======================================================================
//function : Dump
//purpose  : 
//=======================================================================

void  DrawTrSurf_Curve2d::Dump(Standard_OStream& S) const 
{
  GeomTools_Curve2dSet::PrintCurve2d(curv,S);
}


//=======================================================================
//function : Whatis
//purpose  : 
//=======================================================================

void  DrawTrSurf_Curve2d::Whatis(Draw_Interpretor& S)const 
{
  S << "2d curve";
}


//=======================================================================
//function : Is3D
//purpose  : 
//=======================================================================

Standard_Boolean DrawTrSurf_Curve2d::Is3D() const
{
  return Standard_False;
}
