src/IGESGeom/IGESGeom_Plane.cxx

   1 // Created by: CKY / Contract Toubro-Larsen
   2 // Copyright (c) 1993-1999 Matra Datavision
   3 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   4 //
   5 // This file is part of Open CASCADE Technology software library.
   6 //
   7 // This library is free software; you can redistribute it and/or modify it under
   8 // the terms of the GNU Lesser General Public License version 2.1 as published
   9 // by the Free Software Foundation, with special exception defined in the file
  10 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  11 // distribution for complete text of the license and disclaimer of any warranty.
  12 //
  13 // Alternatively, this file may be used under the terms of Open CASCADE
  14 // commercial license or contractual agreement.
  15
  16 //--------------------------------------------------------------------
  17 //--------------------------------------------------------------------
  18
  19 #include <gp_GTrsf.hxx>
  20 #include <gp_Pnt.hxx>
  21 #include <gp_XYZ.hxx>
  22 #include <IGESData_IGESEntity.hxx>
  23 #include <IGESGeom_Plane.hxx>
  24 #include <Standard_Type.hxx>
  25
  26 IMPLEMENT_STANDARD_RTTIEXT(IGESGeom_Plane,IGESData_IGESEntity)
  27
  28 IGESGeom_Plane::IGESGeom_Plane ()     {  }
  29
  30
  31     void IGESGeom_Plane::Init
  32   (const Standard_Real A, const Standard_Real B,
  33    const Standard_Real C, const Standard_Real D,
  34    const Handle(IGESData_IGESEntity)& aCurve,
  35    const gp_XYZ& attach, const Standard_Real aSize)
  36 {
  37   theA = A;
  38   theB = B;
  39   theC = C;
  40   theD = D;
  41   theCurve  = aCurve;
  42   theAttach = attach;
  43   theSize   = aSize;
  44   InitTypeAndForm(108,FormNumber());
  45 // FormNumber : 0 No Curve. +1 Bound. -1 Hole
  46 }
  47
  48     void  IGESGeom_Plane::SetFormNumber (const Standard_Integer form)
  49 {
  50   Standard_Integer fn = 0;
  51   if (form < 0) fn = -1;
  52   if (form > 0) fn = 1;
  53   InitTypeAndForm(108,fn);
  54 }
  55
  56     void IGESGeom_Plane::Equation
  57   (Standard_Real& A, Standard_Real& B, Standard_Real& C, Standard_Real& D) const
  58 {
  59   A = theA;
  60   B = theB;
  61   C = theC;
  62   D = theD;
  63 }
  64
  65     Standard_Boolean IGESGeom_Plane::HasBoundingCurve () const
  66 {
  67   return (!theCurve.IsNull());
  68 }
  69
  70     Standard_Boolean IGESGeom_Plane::HasBoundingCurveHole () const
  71 {
  72   return ((FormNumber() == -1) && (!theCurve.IsNull()));
  73 }
  74
  75     Handle(IGESData_IGESEntity) IGESGeom_Plane::BoundingCurve () const
  76 {
  77   return theCurve;
  78 }
  79
  80     Standard_Boolean IGESGeom_Plane::HasSymbolAttach () const
  81 {
  82   return (theSize > 0);
  83 }
  84
  85     gp_Pnt IGESGeom_Plane::SymbolAttach () const
  86 {
  87   gp_Pnt attach(theAttach);
  88   return attach;
  89 }
  90
  91     gp_Pnt IGESGeom_Plane::TransformedSymbolAttach () const
  92 {
  93   if (theSize > 0 && HasTransf())
  94     {
  95       gp_XYZ Symbol = theAttach;
  96       Location().Transforms(Symbol);
  97       return gp_Pnt(Symbol);
  98     }
  99   else return gp_Pnt(0, 0, 0);
 100 }
 101
 102     Standard_Real IGESGeom_Plane::SymbolSize () const
 103 {
 104   return theSize;
 105 }
 106
 107     void  IGESGeom_Plane::TransformedEquation
 108   (Standard_Real& A, Standard_Real& B, Standard_Real& C, Standard_Real& D) const
 109 {
 110   //eqn of plane AX + BY + CZ = D
 111
 112   Standard_Real x1,y1,z1,x2,y2,z2,x3,y3,z3;
 113
 114   //case 1 intersection of the plane with the XY plane.
 115   x1 = 0.0;
 116   y1 = 0.0;
 117   z1 = theD / theC;
 118
 119   //case 2 intersection of the plane with the XZ plane.
 120   x2 = 0.0;
 121   y2 = theD / theB;
 122   z2 = 0.0;
 123
 124   //case 3 intersection of the plane with the YZ plane.
 125   x3 = theD / theA;
 126   y3 = 0.0;
 127   z3 = 0.0;
 128   gp_XYZ P1(x1,y1,z1);
 129   gp_XYZ P2(x2,y2,z2);
 130   gp_XYZ P3(x3,y3,z3);
 131   Location().Transforms(P1);
 132   Location().Transforms(P2);
 133   Location().Transforms(P3);
 134   x1 = P1.X();
 135   y1 = P1.Y();
 136   z1 = P1.Z();
 137   x2 = P2.X();
 138   y2 = P2.Y();
 139   z2 = P2.Z();
 140   x3 = P3.X();
 141   y3 = P3.Y();
 142   z3 = P3.Z();
 143
 144 /*
 145   General eqn of plane can also be written as
 146   a(x1 -x2) + b(y1 - y2) + c(z1 - z2) = 0
 147   a(x3 - x2) + b(y3 - y2) + c(z3 -z2) = 0
 148   Applying Cramer's Rule :
 149   a                       b                       c
 150   -------------   =   ---------------     =    ---------------   =   k
 151   |y3-y2  z3-z2|      |z3-z2   x3-x2|          |x3-x2   y3-y2|
 152   |y1-y2  z1-z2|      |z1-z2   x1-x2|          |x1-x2   y1-y2|
 153
 154   .
 155   . .  a = c1*k , b = c2*k , c = c3*k
 156   hence c1(x - x2) + c2(y - y2) + c3(z - z2) = 0
 157
 158 */
 159   Standard_Real c1,c2,c3;
 160
 161   c1 = (y1*(-z3 + z2) + y2*(-z1 + z3) + y3*(z1 - z2));
 162   c2 = (x1*(z3 - z2) + x2*(-z3 + z1) + x3*(-z1 + z2));
 163   c3 = (x1*(-y3 + y2) + x2*(-y1 + y3) + x3*(y1 - y2));
 164
 165   A = c1;
 166   B = c2;
 167   C = c3;
 168   D = c1*x2 + c2*y2 + c3*z3;
 169 }