summaryrefslogtreecommitdiff
path: root/src/GccAna/GccAna_LinPnt2dBisec.cxx
blob: 90ffb7fdbca93f04964e97191e3d1f0e6cd3e6d5 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
// File:	GccAna_LinPnt2dBisec.cxx
// Created:	Fri Oct  4 09:53:31 1991
// Author:	Remi GILET
//		<reg@phobox>


//=========================================================================
//   CREATION DE LA BISSECTRICE ENTRE UNE DROITE ET UN POINTS.             +
//=========================================================================

#include <GccAna_LinPnt2dBisec.ixx>

#include <gp_XY.hxx>
#include <gp_Dir2d.hxx>
#include <gp_Ax2d.hxx>
#include <GccInt_BParab.hxx>
#include <GccInt_BLine.hxx>
#include <Standard_ConstructionError.hxx>
#include <StdFail_NotDone.hxx>
#include <gp.hxx>

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

GccAna_LinPnt2dBisec::
   GccAna_LinPnt2dBisec (const gp_Lin2d& Line1 ,
		         const gp_Pnt2d& Point2) {

   WellDone = Standard_False;

   Standard_Real xdir = Line1.Direction().X();
   Standard_Real ydir = Line1.Direction().Y();
   Standard_Real xloc = Line1.Location().X();
   Standard_Real yloc = Line1.Location().Y();
   Standard_Real dist = Line1.Distance(Point2);
//   if ( dist > gp::Resolution()) {
   if ( dist > 1.e-10)
     {
       Standard_Real xpoint2 = Point2.X();
       Standard_Real ypoint2 = Point2.Y();
       if ((-ydir*(xpoint2-xloc)+xdir*(ypoint2-yloc)) > 0.0)
	 {
	   gp_Ax2d axeparab(gp_Pnt2d(Point2.XY()-dist/2.*gp_XY(-ydir,xdir)),
			    gp_Dir2d(-ydir,xdir));
	   gp_Parab2d bislinpnt(axeparab,dist/2.0);
	   bissol = new GccInt_BParab(bislinpnt);
//         =====================================
	 }
       else
	 {
	   gp_Ax2d axeparab(gp_Pnt2d(Point2.XY()+dist/2.*gp_XY(-ydir,xdir)),
			    gp_Dir2d(ydir,-xdir));
	   gp_Parab2d bislinpnt(axeparab,dist/2.0);
	   bissol = new GccInt_BParab(bislinpnt);
//         =====================================
	 }
       WellDone = Standard_True;
     }
   else
     {
       gp_Lin2d bislinpnt(Point2,gp_Dir2d(-ydir,xdir));
       bissol = new GccInt_BLine(bislinpnt);
//     ====================================
       WellDone = Standard_True;
     }
 }

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

Standard_Boolean GccAna_LinPnt2dBisec::
   IsDone () const { return WellDone; }

Handle(GccInt_Bisec) GccAna_LinPnt2dBisec::
   ThisSolution () const 
{
  if (!WellDone)
    StdFail_NotDone::Raise();

  return bissol;
}