Changeset 834 for GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/leaves

Timestamp:

04/26/06 19:02:47 (19 years ago)

Author:

gumbau

Message:

GeoLodTreeLibrary? leaves load from older leaf file format

Location:

GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/leaves

Files:

• Arbol.cpp (modified) (1 diff)
• Hoja.cpp (modified) (6 diffs)

GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/leaves/Arbol.cpp

@@ -522 +522 @@
         {
                 j=i;
-                while (Hojas[j].padre != -1) j= Hojas[j].padre;
+                while (Hojas[j].padre > -1) j= Hojas[j].padre;
                 Hojas[i].raiz = j;
 

GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/leaves/Hoja.cpp

@@ -1 @@
-#include <GL\glut.h>
+#include <math.h>
 #include "Hoja.h"
 
@@ -7 +7 @@
 //--------------------------------------------------------------------------------------------------------------------------------
 Hoja::Hoja(void) 
-{
-        
-        Vert_Hoja[0] = Vert_Hoja[1] = Vert_Hoja[2] = Vert_Hoja[3] = 0;
+{       
+        Vert_Hoja[0] = Vert_Hoja[1] = Vert_Hoja[2] = Vert_Hoja[3] =0; 
+        Centro[0] = Centro[1] = Centro[2] = 0;
         Normal[0] = Normal[1] = Normal[2] = 0;
+        hoja_cerca=-1;
+        Cuantas_hojas = 1;
+        dist = coplanar = criterio = hoja_cop = hoja_crit =-1;
+        existe = true;
         padre=-1;
-        Cuantas_hojas = 1;
         hijoi = hijod = -1;
         raiz= visible= -1;
-//      prev = next = -1;
-        
-                
 }
 
@@ -27 +27 @@
 Hoja::Hoja (const Hoja& aHoja)
 {
-
-
         Cuantas_hojas = aHoja.Cuantas_hojas;;
+        hoja_cerca = aHoja.hoja_cerca;
+        hoja_crit = aHoja.hoja_crit;
+        dist = aHoja.dist;
+        existe = aHoja.existe;
+        coplanar = aHoja.coplanar;
+        hoja_cop = aHoja.hoja_cop;
+        criterio = aHoja.criterio;
+        for ( int i=0;i<3;i++){
+                Centro[i] = aHoja.Centro[i];
+                Normal[i] = aHoja.Normal[i];
+        }
+        for (i = 0L; i < 4; i++)
+                Vert_Hoja[i] = aHoja.Vert_Hoja[i]; 
+        
         padre = aHoja.padre;
         hijoi = aHoja.hijoi;
@@ -35 +47 @@
         visible = aHoja.visible;
         raiz = aHoja.raiz;
-        
-        for (int i = 0L; i < 4; i++)
-                Vert_Hoja[i] = aHoja.Vert_Hoja[i]; 
-        for ( i = 0L; i < 3; i++)
-                Normal[i] = aHoja.Normal[i]; 
-        
-
 }
 
@@ -47 +52 @@
 
 //--------------------------------------------------------------------------------------------------------------------------------
-// Destructor. We must deallocate the memory allocated for pointers to Vert_Hoja and edges
+// Destructor. We must deallocate the memory allocated for pointers to vertices and edges
 //--------------------------------------------------------------------------------------------------------------------------------
 Hoja::~Hoja (void)
@@ -54 +59 @@
 
 
+//--------------------------------------------------------------------------------------------------------------------------------
+//  CALCULA LA DISTANCIA ENTRE HOJAS
+//--------------------------------------------------------------------------------------------------------------------------------
+
+float Hoja::Distancia (Hoja& Hoja1)
+{
+        float dist =0;
+        float x1,y1,z1;
+
+        x1 = Hoja1.Centro[0]; y1 = Hoja1.Centro[1] ; z1 = Hoja1.Centro[2];
+
+
+        //DISTANCIA ENTRE CENTROS
+
+        dist = ((Centro[0]-x1)*(Centro[0]-x1)) + ((Centro[1]-y1)*(Centro[1]-y1)) + ((Centro[2]-z1)*(Centro[2]-z1));
+
+        return (dist);
+}
+
+
+//--------------------------------------------------------------------------------------------------------------------------------
+//  CALCULA LA COPLANARIDAD ENTRE HOJAS
+//--------------------------------------------------------------------------------------------------------------------------------
+
+float Hoja::Coplanaridad (Hoja& Hoja1)
+{
+        float cop =0;
+        float modulo1, modulo2;
+        float x1,y1,z1, nx1, ny1, nz1;
+        float nx, ny, nz;
+
+        //hoja pasada como parametro, normalizo las componentes
+        x1 = Hoja1.Normal[0]; y1 = Hoja1.Normal[1] ; z1 = Hoja1.Normal[2];
+        modulo1 = sqrt ( (x1*x1) + (y1*y1) + (z1*z1));
+        nx1 = x1 / modulo1;
+        ny1 = y1 / modulo1;
+        nz1 = z1 / modulo1;
+
+        // hoja desde la que llamo
+
+        modulo2 = sqrt ( (Normal[0]*Normal[0]) + (Normal[1]*Normal[1]) + (Normal[2]*Normal[2]));
+        nx = Normal[0] / modulo2;
+        ny = Normal[1] / modulo2;
+        nz = Normal[2] / modulo2;
+
+        // producto escalar : si es proximo a 0, perpendiculares, a 1, coplanares
+
+        cop = (nx1*nx) + (ny1*ny) + (nz1*nz);
+
+
+        return (fabs(cop));
+
+}
+