src/gp/gp_Ax3.cxx

   1 // Created on: 1993-08-03
   2 // Created by: Laurent BOURESCHE
   3 // Copyright (c) 1993-1999 Matra Datavision
   4 // Copyright (c) 1999-2014 OPEN CASCADE SAS
   5 //
   6 // This file is part of Open CASCADE Technology software library.
   7 //
   8 // This library is free software; you can redistribute it and/or modify it under
   9 // the terms of the GNU Lesser General Public License version 2.1 as published
  10 // by the Free Software Foundation, with special exception defined in the file
  11 // OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
  12 // distribution for complete text of the license and disclaimer of any warranty.
  13 //
  14 // Alternatively, this file may be used under the terms of Open CASCADE
  15 // commercial license or contractual agreement.
  16
  17
  18 #include <gp_Ax1.hxx>
  19 #include <gp_Ax2.hxx>
  20 #include <gp_Ax3.hxx>
  21 #include <gp_Dir.hxx>
  22 #include <gp_Pnt.hxx>
  23 #include <gp_Trsf.hxx>
  24 #include <gp_Vec.hxx>
  25 #include <Standard_ConstructionError.hxx>
  26
  27 //=======================================================================
  28 //function : gp_Ax3
  29 //purpose  :
  30 //=======================================================================
  31 gp_Ax3::gp_Ax3 (const gp_Pnt& P,
  32                 const gp_Dir& V) : axis(P,V)
  33 {
  34   Standard_Real A = V.X();
  35   Standard_Real B = V.Y();
  36   Standard_Real C = V.Z();
  37   Standard_Real Aabs = A;
  38   if (Aabs < 0) Aabs = - Aabs;
  39   Standard_Real Babs = B;
  40   if (Babs < 0) Babs = - Babs;
  41   Standard_Real Cabs = C;
  42   if (Cabs < 0) Cabs = - Cabs;
  43   gp_Dir D;
  44
  45   //  pour determiner l axe X :
  46   //  on dit que le produit scalaire Vx.V = 0.
  47   //  et on recherche le max(A,B,C) pour faire la division.
  48   //  l une des coordonnees du vecteur est nulle.
  49
  50   if     ( Babs <= Aabs && Babs <= Cabs) {
  51     if (Aabs > Cabs) D.SetCoord(-C,0., A);
  52     else             D.SetCoord( C,0.,-A);
  53   }
  54   else if( Aabs <= Babs && Aabs <= Cabs) {
  55     if (Babs > Cabs) D.SetCoord(0.,-C, B);
  56     else             D.SetCoord(0., C,-B);
  57   }
  58   else {
  59     if (Aabs > Babs) D.SetCoord(-B, A,0.);
  60     else             D.SetCoord( B,-A,0.);
  61   }
  62   vxdir = D;
  63   vydir = V.Crossed(vxdir);
  64 }
  65
  66 void  gp_Ax3::Mirror(const gp_Pnt& P)
  67 {
  68   axis.Mirror (P);
  69   vxdir.Reverse ();
  70   vydir.Reverse ();
  71 }
  72
  73 gp_Ax3  gp_Ax3::Mirrored(const gp_Pnt& P)const
  74 {
  75   gp_Ax3 Temp = *this;
  76   Temp.Mirror (P);
  77   return Temp;
  78 }
  79
  80 void  gp_Ax3::Mirror(const gp_Ax1& A1)
  81 {
  82   vydir.Mirror (A1);
  83   vxdir.Mirror (A1);
  84   axis.Mirror (A1);
  85 }
  86
  87 gp_Ax3  gp_Ax3::Mirrored(const gp_Ax1& A1)const
  88 {
  89   gp_Ax3 Temp = *this;
  90   Temp.Mirror (A1);
  91   return Temp;
  92 }
  93
  94 void  gp_Ax3::Mirror(const gp_Ax2& A2)
  95 {
  96   vydir.Mirror (A2);
  97   vxdir.Mirror (A2);
  98   axis.Mirror (A2);
  99 }
 100
 101 gp_Ax3  gp_Ax3::Mirrored(const gp_Ax2& A2)const
 102 {
 103   gp_Ax3 Temp = *this;
 104   Temp.Mirror (A2);
 105   return Temp;
 106 }
 107