0024484: sprops gives incorrect matrix of inertia and moments
[occt.git] / src / GProp / GProp_SGProps.gxx
old mode 100755 (executable)
new mode 100644 (file)
index 099a6c0..f0d812e
@@ -1,20 +1,16 @@
 // Copyright (c) 1995-1999 Matra Datavision
-// Copyright (c) 1999-2012 OPEN CASCADE SAS
+// Copyright (c) 1999-2014 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.
+// This file is part of Open CASCADE Technology software library.
 //
-// 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.
+// This library is free software; you can redistribute it and / or modify it
+// under the terms of the GNU Lesser General Public 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.
 //
-// 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.
+// Alternatively, this file may be used under the terms of Open CASCADE
+// commercial license or contractual agreement.
 
 #include <Standard_NotImplemented.hxx>
 #include <math_Vector.hxx>
@@ -522,12 +518,14 @@ static Standard_Real CCompute(Face&                  S,
                           LocIxy = FuncAdd(LocIxy, FuncMul(x, YdS));
                           LocIyz = FuncAdd(LocIyz, FuncMul(y, ZdS));
                           LocIxz = FuncAdd(LocIxz, FuncMul(x, ZdS));
-                          x = Precision::IsInfinite(x) ? Precision::Infinite() : x*x;
-                          y = Precision::IsInfinite(y) ? Precision::Infinite() : y*y;
-                          z = Precision::IsInfinite(z) ? Precision::Infinite() : z*z;
-                          LocIxx = FuncAdd(LocIxx, FuncAdd(YdS, ZdS));
-                          LocIyy = FuncAdd(LocIyy, FuncAdd(XdS, ZdS));
-                          LocIzz = FuncAdd(LocIzz, FuncAdd(XdS, YdS));
+
+                          const Standard_Real XXdS = FuncMul(x, XdS);
+                          const Standard_Real YYdS = FuncMul(y, YdS);
+                          const Standard_Real ZZdS = FuncMul(z, ZdS);
+
+                          LocIxx = FuncAdd(LocIxx, FuncAdd(YYdS, ZZdS));
+                          LocIyy = FuncAdd(LocIyy, FuncAdd(XXdS, ZZdS));
+                          LocIzz = FuncAdd(LocIzz, FuncAdd(XXdS, YYdS));
                         }//for: iU
                       }//for: iGU
 
@@ -792,12 +790,14 @@ static void Compute(Face& S, Domain& D, const gp_Pnt& loc, Standard_Real& dim, g
         LocIxy = FuncAdd(LocIxy, FuncMul(x, YdS));
         LocIyz = FuncAdd(LocIyz, FuncMul(y, ZdS));
         LocIxz = FuncAdd(LocIxz, FuncMul(x, ZdS));
-        x = Precision::IsInfinite(x) ? Precision::Infinite() : x*x;
-        y = Precision::IsInfinite(y) ? Precision::Infinite() : y*y;
-        z = Precision::IsInfinite(z) ? Precision::Infinite() : z*z;
-        LocIxx = FuncAdd(LocIxx, FuncAdd(YdS, ZdS));
-        LocIyy = FuncAdd(LocIyy, FuncAdd(XdS, ZdS));
-        LocIzz = FuncAdd(LocIzz, FuncAdd(XdS, YdS));
+
+        const Standard_Real XXdS = FuncMul(x, XdS);
+        const Standard_Real YYdS = FuncMul(y, YdS);
+        const Standard_Real ZZdS = FuncMul(z, ZdS);
+        
+        LocIxx = FuncAdd(LocIxx, FuncAdd(YYdS, ZZdS));
+        LocIyy = FuncAdd(LocIyy, FuncAdd(XXdS, ZZdS));
+        LocIzz = FuncAdd(LocIzz, FuncAdd(XXdS, YYdS));
       }
 
       CArea = FuncAdd(CArea, FuncMul(LocArea, ur));
@@ -923,12 +923,14 @@ static void Compute(const Face& S,
       Ixyi = FuncAdd(Ixyi, FuncMul(x, YdS));
       Iyzi = FuncAdd(Iyzi, FuncMul(y, ZdS));
       Ixzi = FuncAdd(Ixzi, FuncMul(x, ZdS));
-      x = Precision::IsInfinite(x) ? Precision::Infinite() : x*x;
-      y = Precision::IsInfinite(y) ? Precision::Infinite() : y*y;
-      z = Precision::IsInfinite(z) ? Precision::Infinite() : z*z;
-      Ixxi = FuncAdd(Ixxi, FuncAdd(YdS, ZdS));
-      Iyyi = FuncAdd(Iyyi, FuncAdd(XdS, ZdS));
-      Izzi = FuncAdd(Izzi, FuncAdd(XdS, YdS));
+
+      const Standard_Real XXdS = FuncMul(x, XdS);
+      const Standard_Real YYdS = FuncMul(y, YdS);
+      const Standard_Real ZZdS = FuncMul(z, ZdS);
+
+      Ixxi = FuncAdd(Ixxi, FuncAdd(YYdS, ZZdS));
+      Iyyi = FuncAdd(Iyyi, FuncAdd(XXdS, ZZdS));
+      Izzi = FuncAdd(Izzi, FuncAdd(XXdS, YYdS));
     }
 
     dim   = FuncAdd(dim, FuncMul(dsi, GaussWV (j)));