0024166: Unable to create file with "Save" menu of voxeldemo Qt sample

[occt.git] / src / gp / gp_Sphere.cxx

// Copyright (c) 1995-1999 Matra Datavision
// Copyright (c) 1999-2012 OPEN CASCADE SAS
//
// The content of this file is subject to the Open CASCADE Technology Public
// License Version 6.5 (the "License"). You may not use the content of this file
// except in compliance with the License. Please obtain a copy of the License
// at http://www.opencascade.org and read it completely before using this file.
//
// The Initial Developer of the Original Code is Open CASCADE S.A.S., having its
// main offices at: 1, place des Freres Montgolfier, 78280 Guyancourt, France.
//
// The Original Code and all software distributed under the License is
// distributed on an "AS IS" basis, without warranty of any kind, and the
// Initial Developer hereby disclaims all such warranties, including without
// limitation, any warranties of merchantability, fitness for a particular
// purpose or non-infringement. Please see the License for the specific terms
// and conditions governing the rights and limitations under the License.



#include <gp_Sphere.ixx>
#include <Standard_ConstructionError.hxx>

void gp_Sphere::Coefficients
(Standard_Real& A1, Standard_Real& A2, Standard_Real& A3, 
 Standard_Real& B1, Standard_Real& B2, Standard_Real& B3, 
 Standard_Real& C1, Standard_Real& C2, Standard_Real& C3,
 Standard_Real& D) const
{
  // Dans le repere local de la sphere :
  // X*X + Y*Y + Z*Z - radius * radius = 0
  gp_Trsf T;
  T.SetTransformation (pos);
  Standard_Real T11 = T.Value (1, 1);
  Standard_Real T12 = T.Value (1, 2);
  Standard_Real T13 = T.Value (1, 3);
  Standard_Real T14 = T.Value (1, 4);
  Standard_Real T21 = T.Value (2, 1);
  Standard_Real T22 = T.Value (2, 2);
  Standard_Real T23 = T.Value (2, 3);
  Standard_Real T24 = T.Value (2, 4);
  Standard_Real T31 = T.Value (3, 1);
  Standard_Real T32 = T.Value (3, 2);
  Standard_Real T33 = T.Value (3, 3);
  Standard_Real T34 = T.Value (3, 4);
  A1 = T11 * T11 + T21 * T21 + T31 * T31;
  A2 = T12 * T12 + T22 * T22 + T32 * T32;
  A3 = T13 * T13 + T23 * T23 + T33 * T33;
  B1 = T11 * T12 + T21 * T22 + T31 * T32;
  B2 = T11 * T13 + T21 * T23 + T31 * T33;
  B3 = T12 * T13 + T22 * T23 + T32 * T33;
  C1 = T11 * T14 + T21 * T24 + T31 * T34;
  C2 = T12 * T14 + T22 * T24 + T32 * T34;
  C3 = T13 * T14 + T23 * T24 + T33 * T34;
  D  = T14 * T14 + T24 * T24 + T34 * T34 - radius * radius;
}

void gp_Sphere::Mirror (const gp_Pnt& P)
{ pos.Mirror (P); }

gp_Sphere gp_Sphere::Mirrored (const gp_Pnt& P) const
{
  gp_Sphere C = *this;
  C.pos.Mirror (P);
  return C;
}

void gp_Sphere::Mirror (const gp_Ax1& A1)
{ pos.Mirror (A1); }

gp_Sphere gp_Sphere::Mirrored (const gp_Ax1& A1) const
{
  gp_Sphere C = *this;
  C.pos.Mirror (A1);
  return C;
}

void gp_Sphere::Mirror (const gp_Ax2& A2)
{ pos.Mirror (A2); }

gp_Sphere gp_Sphere::Mirrored (const gp_Ax2& A2) const
{
  gp_Sphere C = *this;
  C.pos.Mirror (A2);
  return C;
}