OCCT

// Created on: 1993-08-24
// Created by: Bruno DUMORTIER
// 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 <gp_Circ.hxx>
#include <gp_Elips.hxx>
#include <gp_Hypr.hxx>
#include <gp_Parab.hxx>
#include <gp_Vec.hxx>
#include <ProjLib_Plane.hxx>

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

ProjLib_Plane::ProjLib_Plane() = default;

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl)
{
  Init(Pl);
}

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl, const gp_Lin& L)
{
  Init(Pl);
  Project(L);
}

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl, const gp_Circ& C)
{
  Init(Pl);
  Project(C);
}

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl, const gp_Elips& E)
{
  Init(Pl);
  Project(E);
}

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl, const gp_Parab& P)
{
  Init(Pl);
  Project(P);
}

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

ProjLib_Plane::ProjLib_Plane(const gp_Pln& Pl, const gp_Hypr& H)
{
  Init(Pl);
  Project(H);
}

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

void ProjLib_Plane::Init(const gp_Pln& Pl)
{
  myType       = GeomAbs_OtherCurve;
  isDone       = false;
  myIsPeriodic = false;
  myPlane      = Pl;
}

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

static gp_Pnt2d EvalPnt2d(const gp_Pnt& P, const gp_Pln& Pl)
{
  gp_Vec OP(Pl.Location(), P);
  return gp_Pnt2d(OP.Dot(gp_Vec(Pl.Position().XDirection())),
                  OP.Dot(gp_Vec(Pl.Position().YDirection())));
}

static gp_Dir2d EvalDir2d(const gp_Dir& D, const gp_Pln& Pl)
{
  return gp_Dir2d(D.Dot(Pl.Position().XDirection()), D.Dot(Pl.Position().YDirection()));
}

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

void ProjLib_Plane::Project(const gp_Lin& L)
{
  myType = GeomAbs_Line;
  myLin  = gp_Lin2d(EvalPnt2d(L.Location(), myPlane), EvalDir2d(L.Direction(), myPlane));
  isDone = true;
}

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

void ProjLib_Plane::Project(const gp_Circ& C)
{

  myType = GeomAbs_Circle;

  gp_Pnt2d P2d = EvalPnt2d(C.Location(), myPlane);
  gp_Dir2d X2d = EvalDir2d(C.Position().XDirection(), myPlane);
  gp_Dir2d Y2d = EvalDir2d(C.Position().YDirection(), myPlane);
  gp_Ax22d Ax(P2d, X2d, Y2d);

  myCirc       = gp_Circ2d(Ax, C.Radius());
  myIsPeriodic = true;
  isDone       = true;
}

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

void ProjLib_Plane::Project(const gp_Elips& E)
{
  myType = GeomAbs_Ellipse;

  gp_Pnt2d P2d = EvalPnt2d(E.Location(), myPlane);
  gp_Dir2d X2d = EvalDir2d(E.Position().XDirection(), myPlane);
  gp_Dir2d Y2d = EvalDir2d(E.Position().YDirection(), myPlane);
  gp_Ax22d Ax(P2d, X2d, Y2d);

  myElips      = gp_Elips2d(Ax, E.MajorRadius(), E.MinorRadius());
  myIsPeriodic = true;
  isDone       = true;
}

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

void ProjLib_Plane::Project(const gp_Parab& P)
{
  myType = GeomAbs_Parabola;

  gp_Pnt2d P2d = EvalPnt2d(P.Location(), myPlane);
  gp_Dir2d X2d = EvalDir2d(P.Position().XDirection(), myPlane);
  gp_Dir2d Y2d = EvalDir2d(P.Position().YDirection(), myPlane);
  gp_Ax22d Ax(P2d, X2d, Y2d);

  myParab = gp_Parab2d(Ax, P.Focal());
  isDone  = true;
}

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

void ProjLib_Plane::Project(const gp_Hypr& H)
{
  myType = GeomAbs_Hyperbola;

  gp_Pnt2d P2d = EvalPnt2d(H.Location(), myPlane);
  gp_Dir2d X2d = EvalDir2d(H.Position().XDirection(), myPlane);
  gp_Dir2d Y2d = EvalDir2d(H.Position().YDirection(), myPlane);
  gp_Ax22d Ax(P2d, X2d, Y2d);

  myHypr = gp_Hypr2d(Ax, H.MajorRadius(), H.MinorRadius());
  isDone = true;
}