Changeset 540 for trunk/VUT/GtpVisibilityPreprocessor/src/Beam.cpp

Timestamp:

01/15/06 04:23:51 (19 years ago)

Author:

mattausch

Message:

 

File:

• trunk/VUT/GtpVisibilityPreprocessor/src/Beam.cpp (modified) (4 diffs)

trunk/VUT/GtpVisibilityPreprocessor/src/Beam.cpp

@@ -1 +1 @@
 #include "VssRay.h"
 #include "Beam.h"
-
+#include "Mesh.h"
+#include "Polygon3.h"
 
 
 void
-Beam::Construct(const  AxisAlignedBox3 &box, const  AxisAlignedBox3 &dBox)
+Beam::Construct(const AxisAlignedBox3 &box, const AxisAlignedBox3 &dBox)
 {
   // the frustum is defined by set of negative halfspace described by mPlanes
@@ -73 +74 @@
 
 
-void Beam::ComputeFrustum(float &left, float &right,
-                                                  float &bottom, float &top,
-                                                  float &near, float &far,
-                                                  const AxisAlignedBox3 &sceneBBox) const
+void Beam::ComputePerspectiveFrustum(float &left, float &right,
+                                                                         float &bottom, float &top,
+                                                                         float &near, float &far,
+                                                                         const AxisAlignedBox3 &sceneBBox) const
 {
         const float xDirRange = mDirBox.Max().x - mDirBox.Min().x;
         const float yDirRange = mDirBox.Max().y - mDirBox.Min().y;
+        //Debug << "xdir range: " << xDirRange << endl;
+        //Debug << "ydir range: " << yDirRange << endl;
 
         near = 0.1f; // NOTE: what is the best value for near and far plane?
         far = 2.0f * Magnitude(sceneBBox.Diagonal());
-
-        left = near / tan(xDirRange * 0.5f);
-        right = near / tan(xDirRange * 0.5f);
-
-        bottom = near / tan(yDirRange * 0.5f);
-        top = near / tan(yDirRange * 0.5f);
+        
+        right = fabs(near * tan(xDirRange * 0.5f));
+        left = -right;
+
+        top = fabs(near * tan(yDirRange * 0.5f));
+        bottom = -top;
+}
+
+
+void Beam::ComputeOrthoFrustum(float &left, float &right,
+                                                           float &bottom, float &top,
+                                                           float &near, float &far,
+                                                           const AxisAlignedBox3 &sceneBBox) const
+{
+        const int vAxis = GetMainDirection().DrivingAxis();
+        const int axis2 = (vAxis + 1) % 3;
+        const int axis3 = (vAxis + 2) % 3;
+
+        const float xDirRange = mBox.Max()[axis2] - mDirBox.Min()[axis2];
+        const float yDirRange = mDirBox.Max()[axis3] - mDirBox.Min()[axis3];
+
+        near = 0.1f; // NOTE: what is the best value for near and far plane?
+        far = 2.0f * Magnitude(sceneBBox.Diagonal());
+
+        right = xDirRange * 0.5f;
+        left = -right;
+
+        top = yDirRange * 0.5;
+        bottom = -top;
 }
 
@@ -94 +120 @@
 Vector3 Beam::GetMainDirection() const
 {
-        const Vector3 dCenter = mDirBox.Center();
-        return VssRay::GetDirection(dCenter.x, dCenter.y);
+        return - mPlanes[0].mNormal;
 }
 
@@ -123 +148 @@
 }
 
+
+void Beam::CreateMesh(const float zfar)
+{
+        if (mMesh)
+                return;
+
+        mMesh = new Mesh();
+        
+        // -- compute far plane
+
+        // box should not never remove part of beam polygons
+        Vector3 bmin = mBox.Min() - Vector3(zfar * 2.0);
+        Vector3 bmax = mBox.Max() + Vector3(zfar * 2.0);
+
+        AxisAlignedBox3 bbox(bmin, bmax);
+        Plane3 fplane;
+        fplane.mNormal = -mPlanes[0].mNormal;
+
+        // NOTE: beam far plane must not not be culled by gl far plane
+        fplane.mD = mPlanes[0].mD - zfar - 1.0f; 
+        mPlanes.push_back(fplane);
+
+        for (int i = 0; i < mPlanes.size(); ++ i)
+        {
+                Polygon3 *poly = bbox.CrossSection(mPlanes[i]);
+                
+                if (!poly->Valid(Limits::Small))
+                        DEL_PTR(poly);
+                
+                for (int j = 0; (j < mPlanes.size()) && poly; ++ j)
+                {
+                        if (j != i)
+                        {
+                                Polygon3 *front = new Polygon3();
+                                Polygon3 *back = new Polygon3();
+                                
+                                poly->Split(mPlanes[j], *front, *back, Limits::Small);
+                                DEL_PTR(poly);
+                                DEL_PTR(front);
+
+                                if (!back->Valid(Limits::Small))
+                                        DEL_PTR(back);
+                                poly = back;
+                        }
+                }
+        
+                if (poly)
+                {
+                        poly->AddToMesh(*mMesh);
+                }
+        }
+
+        // remove far plane
+        mPlanes.pop_back();
+}