Changeset 544 for trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp

Timestamp:

01/16/06 19:33:48 (18 years ago)

Author:

mattausch

Message:

 

File:

• trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp (modified) (2 diffs)

trunk/VUT/GtpVisibilityPreprocessor/src/ViewCellBsp.cpp

@@ -9 +9 @@
 #include "AxisAlignedBox3.h"
 #include "Triangle3.h"
+#include "Tetrahedron3.h"
 
 #include <stack>
@@ -2549 +2550 @@
 
 
+float BspNodeGeometry::GetVolume() const
+{
+        //-- compute volume using tetrahedralization of the geometry
+        //   and adding the volume of the single tetrahedrons
+        float volume = 0;
+        const float f = 1.0f / 6.0f;
+
+        PolygonContainer::const_iterator pit, pit_end = mPolys.end();
+
+        for (pit = mPolys.begin(); pit != pit_end; ++ pit)
+        {
+                Polygon3 *poly = *pit;
+                const Vector3 v = poly->mVertices[0];
+
+                for (int i = 1; i < (int)poly->mVertices.size() - 1; ++ i)
+                {
+                        volume += f * (DotProd(v, 
+                                CrossProd(poly->mVertices[i], poly->mVertices[i + 1])));
+                }
+        }
+
+        Debug << "volume: " << volume << endl;
+        return volume;
+}
+
+
+Vector3 BspNodeGeometry::ComputeMassCenter() const
+{
+        int n = 0;
+
+        Vector3 center(0,0,0);
+
+        PolygonContainer::const_iterator pit, pit_end = mPolys.end();
+
+        for (pit = mPolys.begin(); pit != pit_end; ++ pit)
+        {
+                Polygon3 *poly = *pit;
+                
+                VertexContainer::const_iterator vit, vit_end = poly->mVertices.end();
+
+                for(vit = poly->mVertices.begin(); vit != vit_end; ++ vit)
+                {
+                        center += *vit;
+                        ++ n;
+                }
+        }
+
+        return center / (float)n;
+}
+
+
 void BspNodeGeometry::AddToMesh(Mesh &mesh)
 {