home *** CD-ROM | disk | FTP | other *** search
/ ftp.hitl.washington.edu / ftp.hitl.washington.edu.tar / ftp.hitl.washington.edu / pub / people / habib / kodak / Draw_Vrml_new.cc < prev    next >
C/C++ Source or Header  |  2000-04-18  |  61KB  |  1,849 lines
  1.     /*===========================================================================*/
  2.     /* Includes.                                                                 */
  3.     /*===========================================================================*/
  4. //    #include <windows.h> //WinUnix
  5.     #include <GL/gl.h>
  6.     #include <stdio.h>
  7.     #include <iostream.h>
  8.     #include <stdlib.h>
  9.     #include <math.h>
  10.  
  11. #include "trackball.h"
  12. static unsigned int test[10];
  13.     struct point {float x,y,z;};
  14.     #include "f-gen.h"
  15. #define PI 3.14159265
  16.  
  17. int Counter_Id=0;
  18.  
  19. // Flags
  20. int ADJUST_TO_ZERO = 0; // Flattenes out the city to zero
  21. float TransX, TransY, TransZ; // Global Transformation.
  22. //    #include <malloc.h>
  23.  
  24. #define DimX 400 // Change them in interface.h.
  25. #define DimY 150
  26. #define DimZ 400
  27.  
  28. double curquat_vrml[4];
  29. double lastquat_vrml[4];
  30. double m_vrml[4][4];
  31.  
  32. struct Virtual_Espace{
  33.     char  FULL;
  34.     //struct point *N;
  35. } espace[DimX][DimY][DimZ];
  36. float VoxelX,VoxelY, VoxelZ; //It's the dimension of a voxel in real space.
  37. struct espace_indice{
  38.     float Xmin, Xmax, Ymin, Ymax, Zmin, Zmax;
  39.     int   Imin, Imax, Jmin, Jmax, Kmin, Kmax;
  40. } EI;
  41.     typedef enum { p_Transform, p_children, p_Material, p_ImageTexture, p_IndexedFaceSet, p_Coordinate, p_point,
  42.         p_TextureCoordinate, p_coordIndex, p_TimeSensor, p_texCoordIndex, p_Appearance,p_texCoord,p_Shape,
  43.     p_nothing} 
  44.     Type_AccTag;
  45.     Type_AccTag AccTag;
  46. struct Type_Pipe { int pipe[1000];
  47.     int indice; } Acc; // Keeps track of which accolade has just been closed when going into the vrml structure.
  48.  
  49. #define SFNode  struct
  50. #define SFBool char
  51.  
  52. extern "C" unsigned char * read_JPEG_file (char *, char *, unsigned char *, int *, int *);
  53.  
  54. typedef enum { EmptyDracula,machi, Top_Scene ,n_rotation, n_Transform, n_Shape, n_AccO, n_AccF, n_translation, n_scale,
  55.        n_children, n_Appearance, n_Material, n_shininess, n_transparency, n_diffuseColor, n_specularColor,
  56.        n_geometry, n_point, n_coordIndex, f_solid, f_colorIndex, f_color, f_texCoordIndex,
  57.        f_TextureCoordinate, f_texCoord, f_Coordinate, f_texture, nothing, n_material,nappearance,f_coord,
  58. n_USE, n_DEF, n_IndexedFaceSet, nccw, n_solid, n_TRUE, n_FALSE,n_CroO, n_CroF, n_appearance,n_ccw,
  59. n_TimeSensor, n_normal, n_nonexist} Node_type;
  60. Node_type node_type, previous_node_type;
  61.  
  62. FILE *file;
  63. int nbr_vertex = 0;
  64.  
  65. // Everything in those defenition is conform to the VRML2 notation aside :
  66. //    1 Hab_...
  67. //    2    _ins which stands for instance 
  68. //  3 SFNode have been abolished.
  69.  
  70. double XmaxVrmlCity, YmaxVrmlCity, ZmaxVrmlCity;
  71. double XminVrmlCity, YminVrmlCity, ZminVrmlCity;
  72.  
  73. typedef GLfloat SFFloat;
  74. typedef int SFInt32 ;
  75. typedef struct
  76. {
  77.     SFFloat    *pElements;
  78.     int        NumElements;
  79. } MFFloat;
  80.  
  81. typedef struct
  82. {
  83.     SFInt32    *pElements;
  84.     int        num_el;
  85. } MFInt32;
  86.  
  87. #define  SFVec3f struct point 
  88.  
  89. typedef struct Hab_MFVec3f { SFVec3f *pt;
  90.                              int     num_el;
  91. } MFVec3f;
  92.  
  93. typedef struct TypeRotation {    float x;
  94.                                 float y;
  95.                                 float z;
  96.                                 float teta; 
  97. } Type_Rotation;
  98.  
  99. typedef struct Hab_SFColor {    float r;
  100.                                 float g;
  101.                                 float b;
  102. } SFColor;
  103.  
  104. typedef struct Hab_MFColor {SFColor *pt;
  105.                             int     col_el;
  106. } MFColor;
  107.  
  108.  
  109. typedef struct TypeMaterial {
  110.   /*exposedField*/ SFFloat ambientIntensity;  //0.2         // #á[0,1]
  111.   /*exposedField*/ SFColor diffuseColor;      //0.8 0.8 0.8 #á[0,1]
  112.   /*exposedField*/ SFColor emissiveColor;     //0 0 0       #á[0,1]
  113.   /*exposedField*/ float     shininess;       //0.2         #á[0,1]
  114.   /*exposedField*/ SFColor specularColor;     //0 0 0       #á[0,1]
  115.   /*exposedField*/ SFFloat transparency;      //0           #á[0,1]
  116. } Material;
  117.  
  118. typedef struct TypeTexture {
  119.     unsigned char *image;
  120.     int im_width[1];
  121.     int im_height[1];
  122.     unsigned int tex_name[10];
  123.     char tex_blk[40];
  124. } Type_Texture;
  125.  
  126. typedef struct TypeAppearance { Material material;
  127.                             
  128.                                 //TextrueTransform textureTransform;
  129. } Type_Appearance;
  130.                             
  131. typedef struct Hab_Pixel {
  132.     float x,y;
  133. } Pixel;
  134.  
  135.  
  136. typedef struct TypeIndexedFaceSet { 
  137.   /*exposedField*/  MFColor     color;             
  138.   /*eventIn*/       MFInt32     colorIndex;
  139.  
  140.   /*exposedField*/  MFVec3f     coord;             //NULL
  141.   /*eventIn*/       MFInt32     coordIndex;
  142.  
  143.   /*exposedField*/  Pixel       *texCoord;          //NULL
  144.   /*eventIn*/       MFInt32     texCoordIndex;
  145.  
  146.                     MFVec3f     normal;
  147.                     MFInt32     normalIndex;
  148.  
  149.   /*field*/         SFBool      ccw;               //TRUE
  150.   /*field*/         SFBool      solid;             //TRUE
  151.  
  152.                     int         num_el;
  153. } Type_IndexedFaceSet;
  154.  
  155. typedef struct TypeGeometry {
  156.     Type_IndexedFaceSet IndexedFaceSet;
  157. } Type_Geometry;
  158.  
  159. typedef struct TypeShape {
  160.     //Type_Transform *FatherTransform;
  161.     Type_Appearance Appearance;
  162.     Type_Texture    Texture;
  163.     Type_Geometry   Geometry;
  164.     int Acc;
  165.     int Cro;
  166.     int Id;
  167. } Type_Shape;
  168. Type_Shape *CNS;
  169. typedef struct TypeTransform { 
  170.  
  171.     SFVec3f translation;
  172.     SFVec3f scale;
  173.     int FlagScale; // says if we have a scale or not.
  174.     Type_Rotation rotation;
  175.  
  176.     struct TypeTransform *FatherTransform;
  177.  
  178.     struct TypeTransform **ChildrenTransform;
  179.     int    numChildTrans; // this is the total number of sons.
  180.     int    indChildTrans; // this is the actual indice of sons.
  181.  
  182.     Type_Shape     **ChildrenShape;
  183.     int    numChildShape; // this is the total number of sons.
  184.     int    indChildShape; // this is the actual indice of sons.
  185.     int    FlagChildren;
  186.  
  187.     Type_Shape      *Shape;
  188.     int ShapeLeave;      // To say if This node has a Shape Leave or not.
  189.  
  190.     Node_type node_type;
  191.     int Acc;
  192.     int Cro;
  193.                 } Type_Transform;
  194.  
  195. Type_Transform *CN, *old_node, *CNV, *CNC, *CNSe, *CNVo[6][6],*CNSphere;
  196.                                            // CNC is for the city CNV is for the Vehicle. CN is Current Node.
  197.                                            // CNSe for the green Sensors.
  198. struct point ScaleNormals;
  199.  
  200. struct USE_DEF_HASH_TAB
  201. {
  202.     char blk_name[100];
  203.     char tex_blk_name[100];
  204.     MFVec3f coord;
  205.     Pixel  *texCoord;
  206.     int index;
  207.     int texIndex;
  208. } UseDef[5000],iz;
  209.  
  210. inline void wordd_to_node(char *wordd)
  211. {
  212.     node_type = nothing;
  213.     if (strcmp (wordd ,"Transform") ==0)            node_type = n_Transform    ; else
  214.     if (strcmp (wordd , "Shape") ==0)               node_type = n_Shape        ; else
  215.     if (strcmp (wordd , "{")==0)                    node_type = n_AccO   ; else
  216.     if (strcmp (wordd , "}")==0)                    node_type = n_AccF    ; else
  217.     if (strcmp (wordd , "[")==0)                    node_type = n_CroO   ; else
  218.     if (strcmp (wordd , "]")==0)                    node_type = n_CroF    ; else
  219.     if (strcmp (wordd , "translation")==0)          node_type = n_translation  ; else
  220.     if (strcmp (wordd , "scale")==0)                node_type = n_scale  ; else
  221.     if (strcmp (wordd , "children")==0)             node_type = n_children     ; else
  222.     if (strcmp (wordd , "Appearance")==0)           node_type = n_Appearance   ; else
  223.     if (strcmp (wordd , "appearance")==0)           node_type = n_appearance   ; else
  224.     if (strcmp (wordd , "Material"  )==0)           node_type = n_Material     ; else
  225.     if (strcmp (wordd , "material"  )==0)           node_type = n_material     ; else
  226.     if (strcmp (wordd , "shininess" )==0)           node_type = n_shininess    ; else
  227.     if (strcmp (wordd , "transparency" )==0)        node_type = n_transparency ; else
  228.     if (strcmp (wordd , "diffuseColor" )==0)        node_type = n_diffuseColor    ; else
  229.     if (strcmp (wordd , "specularColor" )==0)       node_type = n_specularColor    ; else
  230.     if (strcmp (wordd , "geometry" )==0)            node_type = n_geometry         ; else
  231.     if (strcmp (wordd , "point" ) ==0)              node_type = n_point            ; else
  232.     if (strcmp (wordd , "normal" ) ==0)             node_type = n_normal           ; else
  233.     if (strcmp (wordd , "coordIndex" )==0)          node_type = n_coordIndex       ; else
  234.     if (strcmp (wordd , "texCoordIndex" )==0)       node_type = f_texCoordIndex    ; else
  235.     if (strcmp (wordd , "solid" )==0)               node_type = f_solid            ; else
  236.     if (strcmp (wordd , "colorIndex" )==0)          node_type = f_colorIndex       ; else
  237.     if (strcmp (wordd , "color" )==0)               node_type = f_color            ; else
  238.     if (strcmp (wordd , "Coordinate" )==0)          node_type = f_Coordinate       ; else
  239.     if (strcmp (wordd , "coord" )==0)               node_type = f_coord            ; else
  240.     if (strcmp (wordd , "TimeSensor" )==0)          node_type = n_TimeSensor       ; else
  241.     if (strcmp (wordd , "texCoord"   )==0)          node_type = f_texCoord         ; else
  242.     if (strcmp (wordd , "texture"    )==0)          node_type = f_texture          ; else
  243.     if (strcmp (wordd , "rotation"    )==0)          node_type =n_rotation         ; else
  244.     if (strcmp (wordd , "TextureCoordinate" )==0)   node_type = f_TextureCoordinate; else
  245.     if (strcmp (wordd , "USE" )==0)   node_type = n_USE; else
  246.     if (strcmp (wordd , "DEF" )==0)   node_type = n_DEF; else
  247.     if (strcmp (wordd , "USE" )==0)   node_type = n_USE; else
  248.     if (strcmp (wordd , "IndexedFaceSet" )==0)   node_type = n_IndexedFaceSet; else
  249.     if (strcmp (wordd , "ccw" )==0)   node_type = n_ccw; else
  250.     if (strcmp (wordd , "solid" )==0)   node_type = n_solid; else
  251.     if (strcmp (wordd , "TRUE" )==0)   node_type = n_TRUE; else
  252.     if (strcmp (wordd , "FALSE" )==0)   node_type = n_FALSE; 
  253.  
  254.     
  255. }
  256.  
  257. // reads next float from wrl Ascii file
  258. float LisFlWrl()
  259. {
  260.     float f;
  261.     fscanf (file,"%f",&f);
  262.     return f;
  263. }
  264. // reads next int from wrl Ascii file
  265. int LisIntWrl()
  266. {
  267.     int f;
  268.     fscanf (file,"%d",&f);
  269.     return f;
  270. }
  271.  
  272. void initialize_Transform_node()
  273. {
  274.   CN->node_type = n_Transform;
  275.     //CN->node_type = n_Transform;
  276.     CN->translation.x = 0;
  277.     CN->translation.y = 0;
  278.     CN->translation.z = 0;
  279.  
  280.     CN->rotation.x = 0;
  281.     CN->rotation.y = 0;
  282.     CN->rotation.z = 1;
  283.     CN->rotation.teta = 0;
  284.  
  285.     CN->numChildTrans = 0;
  286.     CN->indChildTrans = 0;
  287.     CN->numChildShape = 0;
  288.     CN->indChildShape = 0;
  289.  
  290.     CN->ShapeLeave = 0;
  291.     CN->Acc = 0;
  292.     CN->Cro = 0;
  293.     CN->FlagChildren = 0;
  294.     CN->FlagScale = 0;
  295. }
  296.  
  297. static counter3 = 0;
  298. inline void Create_Child_Transform_Node()
  299. {
  300.   //cout << "CN->numChildTrans = " << CN->numChildTrans << endl;
  301.     // construct Son
  302.     CN->numChildTrans++;
  303.     int numChildTrans = CN->numChildTrans;
  304.   counter3++;
  305.     // If this is the first son do a malloc for the Node ** CN->Son else do a realloc.
  306.     if (CN->numChildTrans == 1){
  307.         if ( (CN->ChildrenTransform = (Type_Transform**) 
  308.           malloc (sizeof (Type_Transform*)) ) ==NULL) 
  309.             cout << "malloc error 1 " << endl;
  310.     } else
  311.         if ( (CN->ChildrenTransform = (Type_Transform**) realloc (CN->ChildrenTransform,numChildTrans*sizeof (Type_Transform*)) ) ==NULL) 
  312.             cout << "realloc error 11 " << endl;
  313.     // allocate the Son itself
  314.     if ( (CN->ChildrenTransform[CN->numChildTrans-1] = (Type_Transform*) malloc (sizeof (Type_Transform )) ) ==NULL)
  315.             cout << "malloc error 2 " << endl;
  316.     // Set current Type_Transform to son.
  317.     old_node = CN;
  318.     CN = CN->ChildrenTransform[CN->numChildTrans-1];
  319.     CN->FatherTransform = new Type_Transform;
  320.     CN->FatherTransform = old_node;
  321.     CN->numChildTrans = 0;
  322.     CN->indChildTrans = 0;
  323.     CN->numChildShape = 0;
  324.     CN->indChildShape = 0;
  325. }
  326.  
  327. inline 
  328. void Create_Child_Shape_Node()
  329. {
  330.     // construct Son
  331.     CN->numChildShape++;
  332.     int numChildShape = CN->numChildShape;
  333.     // If this is the first son do a malloc for the Node ** CN->Son else do a realloc.
  334.     if (CN->numChildShape == 1){
  335.         if ( (CN->ChildrenShape = (Type_Shape**) malloc (sizeof (Type_Shape*)) ) ==NULL) 
  336.             cout << "malloc error 1df " << endl;
  337.     } else
  338.         if ( (CN->ChildrenShape = (Type_Shape**) realloc (CN->ChildrenShape,numChildShape*sizeof (Type_Shape*)) ) ==NULL) 
  339.             cout << "realloc error 1341 " << endl;
  340.     // allocate the Son itself
  341.     if ( (CN->ChildrenShape[CN->numChildShape-1] = (Type_Shape*) malloc (sizeof (Type_Shape )) ) ==NULL)
  342.             cout << "malloc error 2443 " << endl;
  343.     // Put CNS
  344.     CNS = CN->ChildrenShape[CN->numChildShape-1];
  345. }
  346.  
  347. void
  348. Create_One_Shape_Node()
  349. {
  350.         // construct Son
  351.     CN->ShapeLeave =1;
  352.     if ( (CN->Shape = (Type_Shape*) malloc (sizeof (Type_Shape )) ) ==NULL)
  353.             cout << "malloc error 2443 " << endl;
  354.     CNS = CN->Shape;
  355. }
  356.  
  357. void
  358. InitializeChildrenShape()
  359. {
  360.     CNS->Acc = 0;
  361.     CNS->Cro = 0;
  362.     CNS->Texture.tex_name[0] = 99999;
  363.     CNS->Appearance.material.transparency = 0;
  364.  
  365.     //cout << "Initializing Counter_Id Number  " << Counter_Id <<endl;
  366.     CNS->Id = Counter_Id;
  367.     Counter_Id++;
  368.  
  369.     CNS->Geometry.IndexedFaceSet.color.col_el = 0;
  370.     CNS->Geometry.IndexedFaceSet.colorIndex.num_el = 0;
  371.  
  372.     CNS->Geometry.IndexedFaceSet.coord.num_el=0; // to verify that each Geom has a coord field.
  373.     CNS->Geometry.IndexedFaceSet.coordIndex.num_el=0; 
  374.  
  375.     CNS->Geometry.IndexedFaceSet.texCoordIndex.num_el=0; 
  376.     CNS->Geometry.IndexedFaceSet.num_el = 0;
  377. }
  378.  
  379. void Saffir_the_ys()
  380. {
  381.     MFVec3f &pt = CNS->Geometry.IndexedFaceSet.coord; //by reference
  382.     for (int i = 0 ; i<pt.num_el; i++) 
  383.         pt.pt[i].y = 0;
  384. }
  385. // Adjusts the buildings to the zero level.
  386. void adjust_to_zero()
  387. {
  388.     MFVec3f &pt = CNS->Geometry.IndexedFaceSet.coord; //by reference
  389.     float lowesty = 999999;
  390.     for (int i = 0 ; i<pt.num_el; i++) 
  391.         if (pt.pt[i].y < lowesty) lowesty = pt.pt[i].y;
  392.  
  393.     for (i = 0 ; i<pt.num_el; i++) 
  394.         pt.pt[i].y = pt.pt[i].y - lowesty;
  395. }
  396.  
  397. void read_all_points_into_this_node()
  398. {
  399.     char wordd[90];
  400.     fscanf (file,"%s", &wordd); // reads the [
  401.     if (strcmp (wordd,"[") !=0) {
  402.         cout << "mistakdke! a [ expected in here!!" << endl;
  403.         exit (0);
  404.     }
  405.  
  406.     MFVec3f &pt = CNS->Geometry.IndexedFaceSet.coord; // by reference
  407.     pt.num_el = 1;
  408.     // malloc for the first element read
  409.     if ( (pt.pt = (SFVec3f *) malloc (pt.num_el * sizeof ( SFVec3f)) ) ==NULL)
  410.             cout << "malloc error 1236187 " << endl;
  411.     fscanf (file,"%f %f %f",&pt.pt[pt.num_el-1].x,
  412.                             &pt.pt[pt.num_el-1].y,
  413.                             &pt.pt[pt.num_el-1].z);
  414.     fscanf (file,"%s", &wordd);
  415.     if (strcmp (wordd ,",") ==0) goto eid;
  416.     else goto end;
  417.     // realloc for the rest of the elements
  418. eid:
  419.     // look for , ]
  420.     fscanf (file,"%s", &wordd);
  421.     if (strcmp (wordd ,"]") ==0) goto end;
  422.     pt.num_el++;
  423.     if ( (pt.pt = 
  424.     (SFVec3f *) realloc (pt.pt, pt.num_el * sizeof ( SFVec3f)) ) ==NULL)
  425.             cout << "realloc error 1234187 " << endl;
  426.     pt.pt[pt.num_el-1].x = atof(wordd);
  427.     fscanf (file,"%f %f ",
  428.                             &pt.pt[pt.num_el-1].y,
  429.                             &pt.pt[pt.num_el-1].z);
  430.     fscanf (file,"%s", &wordd);
  431.     if (strcmp (wordd ,",") ==0) goto eid;
  432. end:
  433.     //cout << "I read" << pt.num_el << " points!!" << endl;
  434.     return;
  435. }
  436.    
  437.  
  438. void read_all_normal_vectors_into_this_node()
  439. {
  440.     char wordd[90];
  441.     fscanf (file,"%s", &wordd); // reads the [
  442.     if (strcmp (wordd,"[") !=0) {
  443.         cout << "mistake! a [ expected in here!!" << endl;
  444.         exit (0);
  445.     }
  446.     MFVec3f &pt = CNS->Geometry.IndexedFaceSet.normal; // by reference
  447.     pt.num_el = 1;
  448.     // malloc for the first element read
  449.     if ( (pt.pt = (SFVec3f *) malloc (pt.num_el * sizeof ( SFVec3f)) ) ==NULL)
  450.             cout << "malloc error 1236187 " << endl;
  451.     fscanf (file,"%f %f %f",&pt.pt[pt.num_el-1].x,
  452.                             &pt.pt[pt.num_el-1].y,
  453.                             &pt.pt[pt.num_el-1].z);
  454.     fscanf (file,"%s", &wordd);
  455.     if (strcmp (wordd ,",") ==0) goto eid12;
  456.     else goto end12;
  457.     // realloc for the rest of the elements
  458. eid12:
  459.     // look for , ]
  460.     fscanf (file,"%s", &wordd);
  461.     if (strcmp (wordd ,"]") ==0) goto end12;
  462.     pt.num_el++;
  463.     if ( (pt.pt = 
  464.     (SFVec3f *) realloc (pt.pt, pt.num_el * sizeof ( SFVec3f)) ) ==NULL)
  465.             cout << "realloc error 1234187 " << endl;
  466.     pt.pt[pt.num_el-1].x = atof(wordd);
  467.     fscanf (file,"%f %f ",
  468.                             &pt.pt[pt.num_el-1].y,
  469.                             &pt.pt[pt.num_el-1].z);
  470.     fscanf (file,"%s", &wordd);
  471.     if (strcmp (wordd ,",") ==0) goto eid12;
  472. end12:
  473.     //cout << "I read" << pt.num_el << " points!!" << endl;
  474.     return;
  475. }
  476.  
  477.  
  478. void read_all_texture_points_into_this_node()
  479. {
  480.     char wordd[90];
  481.     fscanf (file,"%s", &wordd); // reads the [
  482.     if (strcmp (wordd,"[") !=0) {
  483.         cout << "mistake!! a [ expected in here!!" << endl;
  484.         exit (0);
  485.     }
  486.     Type_IndexedFaceSet &pt = CNS->Geometry.IndexedFaceSet; // by reference
  487.     // malloc for the first element read
  488.     if ( (pt.texCoord = (Pixel *) malloc (sizeof ( Pixel)) ) ==NULL)
  489.             cout << "malloc error 1236187 " << endl;
  490.     pt.num_el = 1;
  491.     fscanf (file,"%f %f ",&pt.texCoord[pt.num_el-1].x,
  492.                           &pt.texCoord[pt.num_el-1].y );
  493.     fscanf (file,"%s", &wordd);
  494.     if (strcmp (wordd ,",") ==0) goto eid11;
  495.     else goto end11;
  496.     // realloc for the rest of the elements
  497. eid11:
  498.     // look for , ]
  499.     fscanf (file,"%s", &wordd);
  500.     if (strcmp (wordd ,"]") ==0) goto end11;
  501.     pt.num_el++;
  502.     if ( (pt.texCoord = 
  503.     (Pixel *) realloc (pt.texCoord, pt.num_el * sizeof ( Pixel)) ) ==NULL)
  504.             cout << "realloc error 1234187 " << endl;
  505.     pt.texCoord[pt.num_el-1].x = atof(wordd);
  506.     fscanf (file,"%f ",
  507.                             &pt.texCoord[pt.num_el-1].y );
  508.     fscanf (file,"%s", &wordd);
  509.     if (wordd[0] == '#') fscanf (file,"%s", &wordd);
  510.     if (strcmp (wordd ,",") ==0) goto eid11;
  511. end11:
  512.      //cout << "I read" << pt.num_el << " points!!" << endl;
  513.     return;
  514. }
  515. void read_all_texCoordIndex_into_CN()
  516. {
  517.     char wordd[90];
  518.     MFInt32 &pt = CNS->Geometry.IndexedFaceSet.texCoordIndex; // by reference
  519.  
  520.     fscanf (file,"%s", &wordd); // reads the [
  521.     if (strcmp (wordd,"[") !=0) {
  522.         cout << "mistake!!! a [ expected in here!!" << endl;
  523.         exit (0);
  524.     }
  525.     fscanf (file,"%s", &wordd);
  526.     if (strcmp (wordd ,"]") ==0) {
  527.        pt.num_el = 0;
  528.        goto end2tex;
  529.     }
  530.     pt.num_el = 1;
  531.     // malloc the first el.
  532.     if ( (pt.pElements = (SFInt32 *) malloc (pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  533.             cout << "malloc error 1964187tex " << endl;
  534.     pt.pElements[pt.num_el-1] = atoi(wordd);
  535.     fscanf (file,"%s", &wordd);
  536.     if (strcmp (wordd ,",") ==0) goto eid2tex;
  537.     else goto end2tex;
  538.     // realloc the rest el.
  539. eid2tex:
  540.     // look for , ]
  541.     fscanf (file,"%s", &wordd);
  542.     if (strcmp (wordd ,"]") ==0) goto end2tex;
  543.     pt.num_el++;
  544.     if ( (pt.pElements = 
  545.     (SFInt32 *) realloc (pt.pElements, pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  546.             cout << "realloc error 1444187 " << endl;
  547.  
  548.     pt.pElements[pt.num_el-1] = atoi(wordd);
  549.                         
  550.     fscanf (file,"%s", &wordd);
  551.     if (strcmp (wordd ,",") ==0) goto eid2tex;
  552. end2tex:
  553.     //cout << "texture pt.num_el=" << pt.num_el << endl;
  554.     return;
  555. }
  556. void read_all_coordIndex_into_CN()
  557. {
  558.     char wordd[90];
  559.     MFInt32 &pt = CNS->Geometry.IndexedFaceSet.coordIndex; // by reference
  560.  
  561.     fscanf (file,"%s", &wordd); // reads the [
  562.     if (strcmp (wordd,"[") !=0) {
  563.         cout << "mistake!!!! a [ expected in here!!" << endl;
  564.         exit (0);
  565.     }
  566.     fscanf (file,"%s", &wordd);
  567.     if (strcmp (wordd ,"]") ==0) {
  568.        pt.num_el = 0;
  569.        goto end2;
  570.     }
  571.     pt.num_el = 1;
  572.     // malloc the first el.
  573.     if ( (pt.pElements = (SFInt32 *) malloc (pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  574.             cout << "malloc error 1964187 " << endl;
  575.     pt.pElements[pt.num_el-1] = atoi(wordd);
  576.     fscanf (file,"%s", &wordd);
  577.     if (strcmp (wordd ,",") ==0) goto eid2;
  578.     else goto end2;
  579.     // realloc the rest el.
  580. eid2:
  581.     // look for , ]
  582.     fscanf (file,"%s", &wordd);
  583.     if (strcmp (wordd ,"]") ==0) goto end2;
  584.     pt.num_el++;
  585.     if ( (pt.pElements = 
  586.     (SFInt32 *) realloc (pt.pElements, pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  587.             cout << "realloc error 1444187 " << endl;
  588.  
  589.     pt.pElements[pt.num_el-1] = atoi(wordd);
  590.                         
  591.     fscanf (file,"%s", &wordd);
  592.     if (strcmp (wordd ,",") ==0) goto eid2;
  593. end2:
  594.     if (pt.num_el == 1) {
  595.         cout <<endl;
  596.     }
  597.  
  598.   //cout << "pt.num_el=" << pt.num_el << endl;
  599.     // cout << "I read" << pt.num_el << " Indexpoints!!" << endl;
  600.     return;
  601. }
  602. void read_all_colorIndex_into_CN()
  603. {
  604.     char wordd[90];
  605.     fscanf (file,"%s", &wordd); // reads the [
  606.     if (strcmp (wordd,"[") !=0) {
  607.         cout << "mistakeeee! a [ expected in here!!" << endl;
  608.         exit (0);
  609.     }
  610.  
  611.     MFInt32 &pt = CNS->Geometry.IndexedFaceSet.colorIndex; // by reference
  612.     pt.num_el = 1;
  613.     // malloc the first el.
  614.     if ( (pt.pElements = (SFInt32 *) malloc (pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  615.             cout << "malloc error 1764187 " << endl;
  616.     fscanf (file,"%d" , &pt.pElements[pt.num_el-1]);
  617.     fscanf (file,"%s", &wordd);
  618.     if (strcmp (wordd ,",") ==0) goto eid3;
  619.     else goto end3;
  620.     // realloc the rest el.
  621. eid3:
  622.     // look for , ]
  623.     fscanf (file,"%s", &wordd);
  624.     if (strcmp (wordd ,"]") ==0) goto end3;
  625.     pt.num_el++;
  626.     if ( (pt.pElements = 
  627.     (SFInt32 *) realloc (pt.pElements, pt.num_el * sizeof ( SFInt32)) ) ==NULL)
  628.             cout <<  "realloc error 14443287 " << endl;
  629.  
  630.     pt.pElements[pt.num_el-1] = atoi(wordd);
  631.                         
  632.     fscanf (file,"%s", &wordd);
  633.     if (strcmp (wordd ,",") ==0) goto eid3;
  634. end3:
  635.     //cout << "I read" << pt.num_el << " ColorIndexpoints!!" << endl;
  636.     return;
  637. }
  638. void read_color_values()
  639. {
  640.     char wordd[90];
  641.     MFColor &pt = CNS->Geometry.IndexedFaceSet.color; // by reference
  642.     pt.col_el = 1;
  643.     // malloc the first el.
  644.     if ( (pt.pt = (SFColor *) malloc( sizeof (SFColor)) ) == NULL)
  645.         cout << " malloc error 1221122" << endl;
  646.     fscanf (file, "%f %f %f", &pt.pt[0].r,
  647.                               &pt.pt[0].g,
  648.                               &pt.pt[0].b);
  649.     fscanf (file,"%s", &wordd);
  650.     if (strcmp(wordd, ",")==0) goto eid4;
  651.     else goto end4;
  652. eid4:
  653.     // look for , ]
  654.     fscanf (file,"%s", &wordd);
  655.     if (strcmp (wordd ,"]") ==0) goto end4;
  656.     pt.col_el++;
  657.     if ( (pt.pt = (SFColor *) realloc (pt.pt, pt.col_el * sizeof(SFColor)) ) == NULL)
  658.         cout << "realloc error 199834992" << endl;
  659.     pt.pt[pt.col_el-1].r = atof(wordd);
  660.     fscanf (file, "%f %f ", 
  661.                               &pt.pt[pt.col_el-1].g,
  662.                               &pt.pt[pt.col_el-1].b);
  663.     fscanf (file,"%s", &wordd);
  664.     if (strcmp (wordd ,",") ==0) goto eid4;
  665. end4:
  666.     //cout << "I read" << pt.col_el << " Color_ppoints!!" << endl;
  667.     return;
  668. }
  669.  
  670. int Pass_To_Next_Trans()
  671. {
  672.  
  673. recom:
  674.   //cout << "numChild " << CN->numChildTrans << "indChild "<<CN->indChildTrans << endl;
  675.       //cout << "Node_type in Transform" << CN->node_type << endl;
  676.     if (CN->numChildTrans > CN->indChildTrans){
  677.         CN->indChildTrans++;
  678.         CN = CN->ChildrenTransform[CN->indChildTrans-1];
  679.         if (ADJUST_TO_ZERO) 
  680.             glTranslatef (CN->translation.x, 0, CN->translation.z);
  681.         else {    
  682.             glTranslatef (CN->translation.x, CN->translation.y, CN->translation.z);
  683.             glRotatef (CN->rotation.teta*180.0/PI, CN->rotation.x, CN->rotation.y, CN->rotation.z);
  684.             if (CN->FlagScale){
  685.                 glScalef (CN->scale.x, CN->scale.y, CN->scale.z);
  686.                 ScaleNormals.x = CN->scale.x * ScaleNormals.x;
  687.                 ScaleNormals.y = CN->scale.y * ScaleNormals.y;
  688.                 ScaleNormals.z = CN->scale.z * ScaleNormals.z;
  689.             }
  690.  
  691.  
  692.               // I didn't do the above Translatef / Rotatef in order to calculate 
  693.               // correctly the ?maxVrmlCity and do the espace calculation 
  694.               // correctly as well.
  695.             /* TransX = TransX + CN->translation.x;
  696.             TransY = TransY + CN->translation.y;
  697.             TransZ = TransZ + CN->translation.z;
  698.  
  699.             lastquat_vrml[3] = cos(CN->rotation.teta/2.0);
  700.             lastquat_vrml[0] = CN->rotation.x*sin(CN->rotation.teta/2.0);
  701.             lastquat_vrml[1] = CN->rotation.y*sin(CN->rotation.teta/2.0);
  702.             lastquat_vrml[2] = CN->rotation.z*sin(CN->rotation.teta/2.0);
  703.             add_quats(lastquat_vrml, curquat_vrml, curquat_vrml);
  704.             build_rotmatrix(m_vrml, curquat_vrml); */
  705.         }
  706.         return 1;
  707.     }else { // if going up the tree.
  708.         CN->indChildTrans = 0; // initializing for another evantual pass
  709.         if (ADJUST_TO_ZERO)
  710.             glTranslatef (-CN->translation.x,0 , -CN->translation.z);
  711.         else{
  712.             if (CN->FlagScale){
  713.                 glScalef (1/CN->scale.x, 1/CN->scale.y, 1/CN->scale.z);
  714.                 ScaleNormals.x = ScaleNormals.x / CN->scale.x ;
  715.                 ScaleNormals.y = ScaleNormals.y / CN->scale.y ;
  716.                 ScaleNormals.z = ScaleNormals.z / CN->scale.z ;
  717.             }
  718.             glRotatef (-CN->rotation.teta*180.0/PI, CN->rotation.x, CN->rotation.y, CN->rotation.z);
  719.             glTranslatef (-CN->translation.x, -CN->translation.y, -CN->translation.z);
  720.  
  721.             /*
  722.             lastquat_vrml[3] = cos(CN->rotation.teta/2.0);
  723.             lastquat_vrml[0] = -CN->rotation.x*sin(CN->rotation.teta/2.0);
  724.             lastquat_vrml[1] = -CN->rotation.y*sin(CN->rotation.teta/2.0);
  725.             lastquat_vrml[2] = -CN->rotation.z*sin(CN->rotation.teta/2.0);
  726.             add_quats(lastquat_vrml, curquat_vrml, curquat_vrml);
  727.             build_rotmatrix(m_vrml, curquat_vrml);
  728.  
  729.             TransX = TransX - CN->translation.x;
  730.             TransY = TransY - CN->translation.y;
  731.             TransZ = TransZ - CN->translation.z; */
  732.  
  733.         }
  734.  
  735.         CN = CN->FatherTransform;
  736.         if ((CN->node_type == Top_Scene)&&(CN->indChildTrans==CN->numChildTrans)){
  737.             CN->indChildTrans = 0;
  738.       //cout << "Node_type in returning 0" << CN->node_type << endl;
  739.       //cout << " I am returning 0 to the calling function " << endl;
  740.             return 0;
  741.         }
  742.     }
  743.     goto recom;
  744. }
  745.  
  746. void 
  747. Fill_Virtual_Espace_With_Voxels(){
  748.     // The idea is to calculate the intersection of all the triangles 
  749.     // with the voxels in the 3D space, and put the normal value in the
  750.     // corresponding voxel.
  751.     //Geometry
  752.     MFInt32 *ind       = &CNS->Geometry.IndexedFaceSet.coordIndex;
  753.     MFInt32 *indtex    = &CNS->Geometry.IndexedFaceSet.texCoordIndex;
  754.     MFVec3f *pt     = &CNS->Geometry.IndexedFaceSet.coord;
  755.     MFColor &col    =  CNS->Geometry.IndexedFaceSet.color;
  756.     MFInt32    &col_ind=  CNS->Geometry.IndexedFaceSet.colorIndex;
  757.     SFInt32 j,j0,j1,j2, t, t0, t1, t2,k,ii;
  758.     int c=0;
  759.     for (int i = 0;i<ind->num_el; i++) {
  760.         // count the number of vertex.
  761.         int num_vertx = 0;
  762.         while (ind->pElements[i+num_vertx]!= -1) {
  763.             num_vertx++; // q now is the number of vertices.
  764.         }
  765.         // Calculate the normal to the face
  766.         j0=ind->pElements[i  ];
  767.         j1=ind->pElements[i+1];
  768.         j2=ind->pElements[i+2];
  769.         SFVec3f norm, vect0, vect1, vect2, vect3, ScanPt1, Pas, FPt;
  770.         vect1 = sous_vect_from_vect (pt->pt[j1], pt->pt[j0]);
  771.         vect2 = sous_vect_from_vect (pt->pt[j1], pt->pt[j2]);
  772.         norm = vectoriel (vect1 , vect2);
  773.         unitaire (norm); // makes it 1 => norm.x norm.y norm.z 
  774.         // glNormal3f(norm.x, norm.y, norm.z);
  775.             // if 3 vertex
  776.         if (num_vertx==3) {
  777.             float pas , pas1, pas2; // vect subdivision
  778.             int   pasi1, pasi2;
  779.             pas = VoxelX; // dimension of one voxel.
  780.             if (pas > VoxelY) pas = VoxelY;
  781.             if (pas > VoxelZ) pas = VoxelZ;
  782.             pas = pas /3.0;
  783.             pasi2 = longueurv(vect2) / pas;
  784.             vect2 = div_vect_par_nbr (vect2, pasi2);
  785.             // remplit espace (voxel) en scannant chaque triangle.
  786.             for (int zik=0;zik<pasi2;zik++){
  787.                 struct point pt0, pt1, pt2;
  788.                 pt0.x = pt->pt[j0].x + TransX;
  789.                 pt0.y = pt->pt[j0].y + TransY;
  790.                 pt0.z = pt->pt[j0].z + TransZ;
  791.                 pt1.x = pt->pt[j1].x + TransX;
  792.                 pt1.y = pt->pt[j1].y + TransY;
  793.                 pt1.z = pt->pt[j1].z + TransZ;
  794.                 pt2.x = pt->pt[j2].x + TransX;
  795.                 pt2.y = pt->pt[j2].y + TransY;
  796.                 pt2.z = pt->pt[j2].z + TransZ;
  797.                 ScanPt1 = add_vect_to_pt (pt2, mul_vect_par_nbr(vect2,zik));
  798.                 vect0   = sous_vect_from_vect(ScanPt1,pt0);
  799.                 pasi1   = longueurv (vect0) / pas;
  800.                 Pas = div_vect_par_nbr (vect0 ,pasi1);
  801.                 for (int ayr=0;ayr<pasi1;ayr++){
  802.                     FPt = add_vect_to_pt(mul_vect_par_nbr(Pas,ayr), pt0);
  803.                    ii = (FPt.x - XminVrmlCity) * DimX / (XmaxVrmlCity - XminVrmlCity);
  804.                     j = (FPt.y - YminVrmlCity) * DimY / (YmaxVrmlCity - YminVrmlCity);
  805.                     k = (FPt.z - ZminVrmlCity) * DimZ / (ZmaxVrmlCity - ZminVrmlCity);
  806.                     if (espace[ii][j][k].FULL !=1){
  807.                         espace[ii][j][k].FULL = 1;
  808.                         /*if ( (espace[ii][j][k].N = (struct point*) malloc (sizeof (struct point)) ) ==NULL) 
  809.                                 cout << "malloc error espace " << endl;
  810.                         espace[ii][j][k].N[0].x = norm.x;
  811.                         espace[ii][j][k].N[0].y = norm.y;
  812.                         espace[ii][j][k].N[0].z = norm.z; */
  813.                     }
  814.                 }
  815.             }
  816.             i=i+3;
  817.         } // if (num_vertx==3
  818.         if (num_vertx==4) {
  819.             for (int marra = 0; marra<num_vertx; marra++) {
  820.                 j=ind->pElements[i];                
  821.                 glVertex3f (pt->pt[j].x,
  822.                             pt->pt[j].y,
  823.                             pt->pt[j].z);
  824.                     i++;
  825.             }
  826.             glEnd();
  827.         } // if (num_vertx>4
  828.         if (num_vertx>4) {
  829.             glBegin (GL_POLYGON);
  830.             for (int marra = 0; marra<num_vertx; marra++) {
  831.                 j=ind->pElements[i];
  832.                 glVertex3f (pt->pt[j].x,
  833.                             pt->pt[j].y,
  834.                             pt->pt[j].z);
  835.                 i++;
  836.             }
  837.         } // if (num_vertx>4
  838.     // i++ to skip the -1 and position on the next valid input. ( no not necessary since i++ in for...) 
  839.         c++;
  840.     } // for (int i = 0; i<ind->num_el; i++)
  841.  
  842. }    
  843.  
  844.  
  845. //Places point in the espace, for collision detection, This function works even if the 
  846. // espace is empty (no pts yet).
  847. inline 
  848. void 
  849. Virtual_Espace() {
  850.     TransX= TransY= TransZ=0; // Global Transformation.
  851.     VoxelX = (XmaxVrmlCity - XminVrmlCity)/DimX;
  852.     VoxelY = (YmaxVrmlCity - YminVrmlCity)/DimY;
  853.     VoxelZ = (ZmaxVrmlCity - ZminVrmlCity)/DimZ;
  854.     cout << "VoxelX = " << VoxelX << "  -  VoxelY = " << VoxelY << " - VoxelZ = " << VoxelZ << endl;
  855.     int i,j,k,toto1,toto2;
  856.     toto1 = toto2 = 0;
  857.     //Initialize the freeken space.
  858.     for ( i= 0; i<DimX; i++)
  859.         for ( j=0; j<DimY; j++)
  860.             for ( k=0; k<DimZ; k++)
  861.                 espace[i][j][k].FULL = 0;
  862.     /* Calculate the intersection with the surfaces with the city.
  863.     while (CN->node_type != Top_Scene)
  864.         CN = CN->FatherTransform;
  865.     while (Pass_To_Next_Trans()){
  866.      if (CN->ShapeLeave) {
  867.       CNS = CN->Shape;
  868.       Fill_Virtual_Espace_With_Voxels();
  869.      }
  870.         else for (CN->indChildShape=0; CN->indChildShape<CN->numChildShape; CN->indChildShape++) {
  871.             CNS = CN->ChildrenShape[CN->indChildShape];
  872.             Fill_Virtual_Espace_With_Voxels();
  873.         }
  874.     }// while
  875.  
  876.     // writing espace to file.
  877.     FILE *file_espace;
  878.     if ((file_espace = fopen ("espace","wb")) == NULL)
  879.         printf ("\n pas de chance l'ami ton fichier est corrompu \n");
  880.     for (i=0;i<DimX;i++)
  881.         for (j=0;j<DimY;j++)
  882.             for (k=0;k<DimZ;k++){
  883.                 fwrite (&espace[i][j][k].FULL, sizeof(char),1,file_espace);
  884.                 //if (espace[i][j][k].FULL==1){
  885.                 //    toto1++;
  886.                 //    fwrite (&espace[i][j][k].N[0].x, sizeof(float), 1, file_espace);
  887.                 //    fwrite (&espace[i][j][k].N[0].y, sizeof(float), 1, file_espace);
  888.                 //    fwrite (&espace[i][j][k].N[0].z, sizeof(float), 1, file_espace);
  889.                 //}
  890.             } */
  891.  
  892.     // reading espace to file.
  893.     FILE *file_espace;
  894.     if ((file_espace = fopen ("espace","rb")) == NULL)
  895.         printf ("\n pas de chance l'ami ton fichier est corrompu 185739945\n");
  896.     for (i=0;i<DimX;i++)
  897.         for (j=0;j<DimY;j++)
  898.             for (k=0;k<DimZ;k++){
  899.                 fread (&espace[i][j][k].FULL, sizeof(char),1,file_espace);
  900.                 //if (espace[i][j][k].FULL==1){
  901.                 //    toto2++;
  902.                 //    if ( (espace[i][j][k].N = (struct point*) malloc (sizeof (struct point)) ) ==NULL) 
  903.                 //        cout << "malloc error espaceeee " << endl;
  904.                 //    fread (&espace[i][j][k].N[0].x, sizeof(float), 1, file_espace);
  905.                 //    fread (&espace[i][j][k].N[0].y, sizeof(float), 1, file_espace);
  906.                 //    fread (&espace[i][j][k].N[0].z, sizeof(float), 1, file_espace);
  907.                 //}
  908.             } 
  909.             
  910.     // close the file, flush the buffer.
  911.     fclose (file_espace);
  912.  
  913.     // Fill the empty bottom...
  914.     for ( i= 0; i<DimX; i++)
  915.         for ( k=0; k<DimZ; k++)
  916.             for ( j=0; j<DimY; j++)
  917.                 if (espace[i][j][k].FULL == 1){
  918.                     for (int jj = 0; jj<=j+3; jj++){//The j+3 is due to the fact that I want the earth to be located 3 voxels above the real earth.
  919.                         espace[i][jj][k].FULL = j-jj+1+3; //The +3 is due to the fact that I want the earth to be located 3 voxels above the real earth.
  920.                         //espace[i][jj][k].N = espace[i][j][k].N;
  921.                     }
  922.                     j=DimY;
  923.                 } 
  924.     extern float HauteurEye;
  925.     HauteurEye = 2*VoxelY;
  926.     cout << "HauteurEye = " << HauteurEye << endl;
  927. }
  928.  
  929. inline 
  930. void 
  931. Draw_Shape_Node(int listnum){
  932.         extern double TOP, BOTTOM, LEFT, RIGHT;
  933.         GLfloat v[4],v1[4];
  934.         v[0] = CNS->Appearance.material.diffuseColor.r;
  935.         v[1] = CNS->Appearance.material.diffuseColor.g; 
  936.         v[2] = CNS->Appearance.material.diffuseColor.b;
  937.         v[3] = CNS->Appearance.material.transparency;
  938.     //cout << v[0] << "  vs  " << v[1] << endl;
  939.     //cout << v[2] << "  vs  " << v[3] << endl;
  940.         v1[0] = CNS->Appearance.material.specularColor.r;
  941.         v1[1] = CNS->Appearance.material.specularColor.g; 
  942.         v1[2] = CNS->Appearance.material.specularColor.b;
  943.         v1[3] = CNS->Appearance.material.transparency;
  944.         //cout << v[0] << " " << v[1] << " " << v[2] << " " << v[3] << endl;
  945.         
  946.         //glColorMaterial (GL_FRONT_AND_BACK, GL_DIFFUSE);
  947.         //glMateriali (GL_FRONT, GL_SHININESS,CNS->Appearance.material.shininess);
  948.         //glColorMaterial (GL_FRONT, GL_SPECULAR);
  949.         if (CNS->Geometry.IndexedFaceSet.color.col_el == 0)
  950.          glDisable (GL_COLOR_MATERIAL);
  951.         else 
  952.             glEnable (GL_COLOR_MATERIAL);
  953.         glMaterialfv (GL_FRONT_AND_BACK, GL_DIFFUSE,v); 
  954.         glMaterialfv (GL_FRONT_AND_BACK, GL_SPECULAR,v1); 
  955.         //glMaterialfv (GL_BACK , GL_DIFFUSE,v1); 
  956.         
  957.         //Geometry
  958.         MFInt32 *ind       = &CNS->Geometry.IndexedFaceSet.coordIndex;
  959.         MFInt32 *indtex    = &CNS->Geometry.IndexedFaceSet.texCoordIndex;
  960.         MFVec3f *pt     = &CNS->Geometry.IndexedFaceSet.coord;
  961.         MFColor &col    =  CNS->Geometry.IndexedFaceSet.color;
  962.         MFInt32    &col_ind=  CNS->Geometry.IndexedFaceSet.colorIndex;
  963.         SFInt32 j,j0,j1,j2, t, t0, t1, t2;
  964.         
  965.         int c=0;
  966.     //cout << "ind->num_el " << ind->num_el << " Of Counter_Id " << CNS->Id << endl;
  967.  
  968.         for (int i = 0;i<ind->num_el; i++) {
  969.             // count the number of vertex.
  970.             int num_vertx = 0;
  971.             while (ind->pElements[i+num_vertx]!= -1) {
  972.         //cout << "ind->pElements[" << i+num_vertx << "]="<< ind->pElements[i+num_vertx] << endl;
  973.                 num_vertx++; // q now is the number of vertices.
  974.             }
  975.             // Sets the normal for light calculation (see journals Dec 9th 1998).
  976.             j0=ind->pElements[i  ];
  977.       //cout << "j0=" << j0;
  978.             j1=ind->pElements[i+1];
  979.       //cout << "j0=" << j1;
  980.             j2=ind->pElements[i+2];
  981.       //cout << "j0=" << j2;
  982.             SFVec3f norm, vect1, vect2;
  983.             vect1 = sous_vect_from_vect (pt->pt[j1], pt->pt[j0]);
  984.             vect2 = sous_vect_from_vect (pt->pt[j1], pt->pt[j2]);
  985.             norm = vectoriel (vect1 , vect2);
  986.             //if (norm.y == 0) norm.y = 0.01; // For collision detection.
  987.             unitaire (norm); // makes it normal.
  988.             // I mult with ScaleNormals to correct for the light calculations.
  989.             glNormal3f( ScaleNormals.x * norm.x, ScaleNormals.y * norm.y,  ScaleNormals.z * norm.z);
  990.  
  991.             // Binds the texture if one exists
  992.             if (CNS->Texture.tex_name[0] !=99999) {
  993.                 glEnable(GL_TEXTURE_2D);
  994.                 glBindTexture(GL_TEXTURE_2D, CNS->Texture.tex_name[0]);
  995.             }else
  996.             glDisable(GL_TEXTURE_2D);
  997.             //glTexEnvf(GL_TEXTURE_ENV, GL_TEXTURE_ENV_MODE, GL_DECAL);
  998.             // if 3 vertex
  999.             if (num_vertx==3) {
  1000.                 glPointSize (3);
  1001.                 glBegin (GL_TRIANGLES);
  1002.                 if (CNS->Geometry.IndexedFaceSet.solid == GL_TRUE) {
  1003.                     //glPolygonMode(GL_FRONT,GL_LINE);
  1004.                     glFrontFace(GL_CCW);
  1005.           cout << "culling" << endl;
  1006.                     glCullFace (GL_FRONT_AND_BACK);// The culling must be within the glBegin
  1007.                     glEnable(GL_CULL_FACE);// if we want a triangle rendered when outside
  1008.                 }                      // of window boundary.
  1009.                 for (int marra = 0; marra<num_vertx; marra++) {
  1010.                     struct point mempt;
  1011.                     j=ind->pElements[i];
  1012.                     mempt.x = pt->pt[j].x;
  1013.                     mempt.y = pt->pt[j].y;
  1014.                     mempt.z = pt->pt[j].z;
  1015.           //cout << "mempt.xyz = " << mempt.x << " " << mempt.y << " " <<mempt.z << endl;
  1016.                     // pt->pt[j] = mult_matrix_vector (m_vrml, pt->pt[j]); // Rotate.
  1017.                     //pt->pt[j].x = pt->pt[j].x + TransX; // I did this instead of Translatef
  1018.                     //pt->pt[j].y = pt->pt[j].y + TransY; // In order to calculate ?maxVrmlCity
  1019.                     //pt->pt[j].z = pt->pt[j].z + TransZ; // correctly and do the espace thing.
  1020.                     if (col_ind.num_el > c)
  1021.                     glColor3f (col.pt[col_ind.pElements[c]].r,
  1022.                                col.pt[col_ind.pElements[c]].g,
  1023.                                col.pt[col_ind.pElements[c]].b); 
  1024.           //if there is a texture, and if the indexed 3Dpoint has a correspondant 2D texture pt.
  1025.                     if ((CNS->Texture.tex_name[0] !=99999) && (i<indtex->num_el) ) {
  1026.                         //glBindTexture(GL_TEXTURE_2D, CNB->Geometry_ins.texture.tex_name[0]);
  1027.                         t=indtex->pElements[i];
  1028.                         glTexCoord2f(CNS->Geometry.IndexedFaceSet.texCoord[t].x , CNS->Geometry.IndexedFaceSet.texCoord[t].y);
  1029.                     }
  1030.                     glVertex3f (pt->pt[j].x,
  1031.                                 pt->pt[j].y,
  1032.                                 pt->pt[j].z);
  1033.                     nbr_vertex++;
  1034.                     if (listnum == 9) if (TOP    < pt->pt[j].y) TOP    = pt->pt[j].y;
  1035.                     if (listnum == 9) if (BOTTOM > pt->pt[j].y) BOTTOM = pt->pt[j].y;
  1036.                     if (listnum == 9) if (RIGHT  < pt->pt[j].x) RIGHT  = pt->pt[j].x;
  1037.                     if (listnum == 9) if (LEFT   > pt->pt[j].x) LEFT   = pt->pt[j].x;
  1038.                     if (listnum == 9) if (ZmaxVrmlCity   < pt->pt[j].z) ZmaxVrmlCity = pt->pt[j].z;
  1039.                     if (listnum == 9) if (ZminVrmlCity   > pt->pt[j].z) ZminVrmlCity = pt->pt[j].z;
  1040.                     pt->pt[j].x = mempt.x;  // For an eventual second pass.
  1041.                     pt->pt[j].y = mempt.y;
  1042.                     pt->pt[j].z = mempt.z;
  1043.                     i++;
  1044.                 }
  1045.                 glEnd();
  1046.             } // if (num_vertx==4
  1047.             if (num_vertx==4) {
  1048.                 glBegin (GL_QUADS);
  1049.                 if (CNS->Geometry.IndexedFaceSet.solid == GL_FALSE) {
  1050.                     glCullFace (GL_FRONT_AND_BACK);// The culling must be within the glBegin
  1051.                     glEnable(GL_CULL_FACE);// if we want a triangle rendered when outside
  1052.                 }                          // of window boundary.
  1053.                 for (int marra = 0; marra<num_vertx; marra++) {
  1054.                     j=ind->pElements[i];
  1055.                     pt->pt[j].x = pt->pt[j].x + TransX;
  1056.                     pt->pt[j].y = pt->pt[j].y + TransY;
  1057.                     pt->pt[j].z = pt->pt[j].z + TransZ;
  1058.                     if (col_ind.num_el > c)
  1059.                         glColor3f (col.pt[col_ind.pElements[c]].r,
  1060.                                    col.pt[col_ind.pElements[c]].g,
  1061.                                    col.pt[col_ind.pElements[c]].b);
  1062.                                                    
  1063.                     if (CNS->Texture.tex_name[0] !=99999) {
  1064.                         t=indtex->pElements[i];
  1065.                         //glBindTexture(GL_TEXTURE_2D, CNB->Geometry_ins.texture.tex_name[0]);
  1066.                         glTexCoord2f(CNS->Geometry.IndexedFaceSet.texCoord[t].x , CNS->Geometry.IndexedFaceSet.texCoord[t].y);
  1067.                     }
  1068.                     glVertex3f (pt->pt[j].x,
  1069.                                 pt->pt[j].y,
  1070.                                 pt->pt[j].z);
  1071.                     if (listnum == 9) if (TOP    < pt->pt[j].y) TOP    = pt->pt[j].y;
  1072.                     if (listnum == 9) if (BOTTOM > pt->pt[j].y) BOTTOM = pt->pt[j].y;
  1073.                     if (listnum == 9) if (RIGHT  < pt->pt[j].x) RIGHT  = pt->pt[j].x;
  1074.                     if (listnum == 9) if (LEFT   > pt->pt[j].x) LEFT   = pt->pt[j].x;
  1075.                     if (listnum == 9) if (ZmaxVrmlCity   < pt->pt[j].z) ZmaxVrmlCity = pt->pt[j].z;
  1076.                     if (listnum == 9) if (ZminVrmlCity   > pt->pt[j].z) ZminVrmlCity = pt->pt[j].z;
  1077.                     i++;
  1078.                 }
  1079.                 glEnd();
  1080.             } // if (num_vertx>4
  1081.             if (num_vertx>4) {
  1082.                 glBegin (GL_POLYGON);
  1083.                 if (CNS->Geometry.IndexedFaceSet.solid == GL_FALSE) {
  1084.                     glCullFace (GL_FRONT_AND_BACK);// The culling must be within the glBegin
  1085.                     glEnable(GL_CULL_FACE);// if we want a triangle rendered when outside
  1086.                 }                          // of window boundary.
  1087.                 for (int marra = 0; marra<num_vertx; marra++) {
  1088.                     j=ind->pElements[i];
  1089.                     pt->pt[j].x = pt->pt[j].x + TransX;
  1090.                     pt->pt[j].y = pt->pt[j].y + TransY;
  1091.                     pt->pt[j].z = pt->pt[j].z + TransZ;
  1092.                     
  1093.                     if (col_ind.num_el > c)
  1094.                     glColor3f (col.pt[col_ind.pElements[c]].r,
  1095.                                col.pt[col_ind.pElements[c]].g,
  1096.                                col.pt[col_ind.pElements[c]].b); 
  1097.  
  1098.                     if (CNS->Texture.tex_name[0] !=99999) {
  1099.                         t=indtex->pElements[i];
  1100.                         //glBindTexture(GL_TEXTURE_2D, CNB->Geometry_ins.texture.tex_name[0]);
  1101.                     glTexCoord2f(CNS->Geometry.IndexedFaceSet.texCoord[t].x , CNS->Geometry.IndexedFaceSet.texCoord[t].y);
  1102.                     }
  1103.                     glVertex3f (pt->pt[j].x,
  1104.                                 pt->pt[j].y,
  1105.                                 pt->pt[j].z);
  1106.                     if (listnum == 9) if (TOP    < pt->pt[j].y) TOP    = pt->pt[j].y;
  1107.                     if (listnum == 9) if (BOTTOM > pt->pt[j].y) BOTTOM = pt->pt[j].y;
  1108.                     if (listnum == 9) if (RIGHT  < pt->pt[j].x) RIGHT  = pt->pt[j].x;
  1109.                     if (listnum == 9) if (LEFT   > pt->pt[j].x) LEFT   = pt->pt[j].x;
  1110.                     if (listnum == 9) if (ZmaxVrmlCity   < pt->pt[j].z) ZmaxVrmlCity = pt->pt[j].z;
  1111.                     if (listnum == 9) if (ZminVrmlCity   > pt->pt[j].z) ZminVrmlCity = pt->pt[j].z;
  1112.                     i++;
  1113.                 }
  1114.             glEnd();
  1115.             } // if (num_vertx==3
  1116.         // i++ to skip the -1 and position on the next valid input. ( no not necessary since i++ in for...) 
  1117.             c++;
  1118.         } // for (int i = 0; i<ind->num_el; i++)
  1119.  
  1120. }    
  1121. void Draw_Vrml(Type_Transform *CNlocal, int listnum) // CNlocal is either CNC or CNV or...
  1122. {
  1123.     extern double TOP, BOTTOM, LEFT, RIGHT;
  1124.     CN = CNlocal;
  1125.     // calculate the correct Bounding box
  1126.     if (listnum == 9) TOP = -9999999;
  1127.     if (listnum == 9) BOTTOM = 9999999; 
  1128.  
  1129.     if (listnum == 9) LEFT = 9999999;
  1130.     if (listnum == 9) RIGHT =-9999999; 
  1131.  
  1132.     if (listnum == 9) ZmaxVrmlCity = -9999999;
  1133.     if (listnum == 9) ZminVrmlCity =  9999999;
  1134.     TransX= TransY= TransZ=0; // Global Transformation.
  1135.     trackball(curquat_vrml, 0.0, 0.0, 0.0,0.0);
  1136.     nbr_vertex = 0;
  1137.     glEnable (GL_TEXTURE_2D);
  1138.     // find beginning of nodes
  1139.   int counter = 0;
  1140.     while (CN->node_type != Top_Scene)
  1141.         CN = CN->FatherTransform;
  1142.     // compile the list.
  1143.     glNewList(listnum, GL_COMPILE_AND_EXECUTE);
  1144.     while (Pass_To_Next_Trans()){
  1145.     counter ++;
  1146.         if (CN->ShapeLeave) {
  1147.             CNS = CN->Shape;
  1148.             Draw_Shape_Node(listnum);
  1149.         }
  1150.         else for (CN->indChildShape=0; CN->indChildShape<CN->numChildShape; CN->indChildShape++) {
  1151.             CNS = CN->ChildrenShape[CN->indChildShape];
  1152.             Draw_Shape_Node(listnum);
  1153.         }
  1154.     }// while
  1155.   //cout << "number of Transform nodes" << counter << endl;
  1156.     glDisable(GL_TEXTURE_2D);
  1157.     glEndList();
  1158.     cout << "nbr_vertex = " << nbr_vertex << endl;
  1159.     cout << "top=" << TOP << "  bottom=" << BOTTOM << "RIGHT=" << RIGHT <<" LEFT=" << LEFT << endl;
  1160.     if (listnum == 9) XmaxVrmlCity = RIGHT;
  1161.     if (listnum == 9) YmaxVrmlCity = TOP;
  1162.     if (listnum == 9) XminVrmlCity = LEFT;
  1163.     if (listnum == 9) YminVrmlCity = BOTTOM;
  1164.     //if (listnum == 9) Virtual_Espace();
  1165. }
  1166.  
  1167. /* void Draw_Vrml_Vehicle()
  1168. {
  1169.     CN = CNV;
  1170.     //Draw_Vrml(19);
  1171.     glNewList(19, GL_COMPILE_AND_EXECUTE);
  1172.     glBegin(GL_TRIANGLES);
  1173.     glVertex3f (0,0,0);
  1174.     glVertex3f (0.1,0.1,0.1);
  1175.     glVertex3f (0.1, 0, 0.2);
  1176.     glEnd();
  1177.     glEndList(); 
  1178.  
  1179. }
  1180. void Draw_Vrml_City()
  1181. {
  1182.     CN = CNC;
  1183.     Draw_Vrml(9);
  1184.     CN = CNC;
  1185.     Virtual_Espace();
  1186. } */
  1187.  
  1188. Test_For_Coord_In_Geometry_node()
  1189. {    /*if ((CN->node_type == n_geometry) &&
  1190.          (CN->Geometry_ins.IndexedFaceSet_Ins.coord.point.num_el == 0)) 
  1191.     CN->Geometry.IndexedFaceSet.coord = CN->Father->Geometry_ins.IndexedFaceSet_Ins.coord;*/
  1192.     //cout << "forget me" << endl;
  1193. }
  1194. void
  1195. inline
  1196. OuvreAccolade()
  1197. {
  1198.     CN->Acc++; // Not necessary yet.
  1199.     AccTag = p_nothing;          // Keeps track of the Accolades
  1200.     Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1201.     Acc.indice = Acc.indice +1;
  1202. }
  1203. void 
  1204. inline
  1205. OuvreCrochet()
  1206. {
  1207.     CN->Cro++; // not necessary yet
  1208.     AccTag = p_nothing;          // Keeps track of the Accolades
  1209.     Acc.indice = Acc.indice +1;
  1210.     Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1211. }
  1212.  
  1213. void
  1214. inline
  1215. FermeCrochet()
  1216. {
  1217.     CN->Cro--; // not necessary yet
  1218.     Acc.indice = Acc.indice -1;
  1219.     if (Acc.indice <0) {
  1220.         cout << "Acc.indice < 0 search for me..." << endl;
  1221.         getchar();
  1222.         exit(1);
  1223.     }
  1224. }
  1225.  
  1226.  
  1227. void
  1228. inline 
  1229. read_RAW_file (char *directory, unsigned char *image, int *im_width, int *im_height)
  1230. {
  1231.   FILE * raw;
  1232.   im_width [0] = 8;
  1233.   im_height[0] = 8;
  1234.   if ((raw = fopen (directory, "rb"))==NULL){
  1235.     cout << "Error reading the raw file" << endl;
  1236.     cout << "width = " << im_width[0] << "   height=" << im_height[0]<< endl;
  1237.     // But still allocates a 64 byte image, just so that the program doesn't crash...
  1238.     if ((CNS->Texture.image = (unsigned char *)malloc (im_width[0]*im_height[0]*3) ) == NULL){
  1239.       cout << "Malloc error here 938488uf74" << endl;
  1240.       exit (1);
  1241.     }
  1242.     // fclose (raw);
  1243.     return ;
  1244.   }
  1245.   fread (im_width , sizeof(int), 1,raw);
  1246.   fread (im_height, sizeof(int), 1,raw);
  1247.   //cout << "width = " << im_width[0] << "   height=" << im_height[0]<< endl;
  1248.   if ((CNS->Texture.image = (unsigned char *)malloc (im_width[0]*im_height[0]*3) ) == NULL){
  1249.         cout << "Malloc error here 938488uf74" << endl;
  1250.         exit (1);
  1251.   }
  1252.   fread (CNS->Texture.image, sizeof(unsigned char), im_width[0]*im_height[0]*3,raw);
  1253.   //cout << "image read " << endl;
  1254.   fclose (raw); /*
  1255.   if ((CNS->Texture.image = (unsigned char *)malloc (128*256*3) ) == NULL){
  1256.         cout << "Malloc error here 9384793874" << endl;
  1257.         exit (1);
  1258.         }  
  1259.       for (int muk = 0 ; muk <64*64*3; muk++)
  1260.         CNS->Texture.image[muk] = 230;
  1261.         CNS->Texture.im_width[0]= 128; CNS->Texture.im_height[0]=256; */
  1262. }
  1263.  
  1264. Type_Transform *
  1265. Read_Vrml_File(char *vrml_filename)
  1266. {
  1267.     extern char *texture_directory;
  1268.     // Initialize ScaleNormals to 1.
  1269.     ScaleNormals.x = 1;
  1270.     ScaleNormals.y = 1;
  1271.     ScaleNormals.z = 1;
  1272.  
  1273.     int size,i;
  1274.     char wordd[90];
  1275.     Acc.indice = 0;
  1276.     iz.index = -1;
  1277.     iz.texIndex = -1;
  1278.  
  1279.     //glGenTextures(10,test);
  1280.  
  1281.     // Open vrml file for reading
  1282.     if ((file = fopen (vrml_filename,"r")) == NULL)
  1283.         printf ("\n pas de chance l'ami ton fichier est corrompu 89472229099\n");
  1284.   //cout << vrml_filename << endl;
  1285.  
  1286.     // initialising the node top_scene.
  1287.     CN = new Type_Transform;
  1288.     initialize_Transform_node();
  1289.     CN->node_type = Top_Scene;
  1290.     previous_node_type = Top_Scene;
  1291.   int counter2 = 0;
  1292.     while (fscanf (file,"%s",&wordd) != EOF) {
  1293.         // cout << wordd << endl;
  1294.         wordd_to_node(wordd); // result in node
  1295.         switch (node_type) {
  1296.         case n_CroO:
  1297.             OuvreCrochet();
  1298.             break;
  1299.         case n_CroF:
  1300.             FermeCrochet();
  1301.             break;
  1302.         case n_AccO:
  1303.             OuvreAccolade();
  1304.             break;
  1305.         case n_AccF:
  1306.             CN->Acc--; // Not necessary yet.
  1307.             Acc.indice = Acc.indice -1;
  1308.             if (Acc.indice >=1000){
  1309.                 cout << "Change your indice, exit 38959988477683" << endl;
  1310.                 exit(1);
  1311.             }
  1312.             if (Acc.pipe[Acc.indice] == p_Transform)
  1313.                 CN = CN->FatherTransform;
  1314.             if (Acc.indice <0) {
  1315.                 cout << "Acc.indice < 0 search for meeeee 135335..." << endl;
  1316.                 getchar();
  1317.                 exit(1);
  1318.             }
  1319.             break;
  1320.         case n_IndexedFaceSet:
  1321.             AccTag = p_IndexedFaceSet;          // Keeps track of the Accolades
  1322.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1323.       if (Acc.indice >=1000){
  1324.         cout << "Change your indice, exit 38959988477683" << endl;
  1325.         exit(1);
  1326.       }
  1327.             Acc.indice++;
  1328.             fscanf (file,"%s",&wordd);
  1329.             if (wordd[0] != '{') {
  1330.                 cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1331.                 getchar();
  1332.             }
  1333.             break;
  1334.         case n_Transform : 
  1335.             //Test_For_Coord_In_Geometry_node();
  1336.       //cout << "Create_Child_Transform_Node" << endl;
  1337.             Create_Child_Transform_Node();
  1338.             initialize_Transform_node();
  1339.             //cout << "CN-> " << CN->node_type << endl;
  1340.             AccTag = p_Transform;          // Keeps track of the Accolades
  1341.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1342.       if (Acc.indice >=1000){
  1343.         cout << "Change your indice, exit 38959988477683" << endl;
  1344.         exit(1);
  1345.       }
  1346.             Acc.indice++;
  1347.       counter2 ++;
  1348.             fscanf (file,"%s",&wordd);
  1349.             if (wordd[0] != '{') {
  1350.                 cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1351.                 getchar();
  1352.             }
  1353.             node_type = nothing;
  1354.             break;
  1355.         case n_Appearance:
  1356.             //Test_For_Coord_In_Geometry_node();
  1357.             //Create_Son_Node();
  1358.             //CN->node_type = n_Appearance;
  1359.             //previous_node_type = n_Appearance;
  1360.             //CN->Geometry_ins.texture.tex_name[0] = 99999;
  1361.             AccTag = p_Appearance;          // Keeps track of the Accolades
  1362.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1363.       if (Acc.indice >=1000){
  1364.         cout << "Change your indice, exit 38959988477683" << endl;
  1365.         exit(1);
  1366.       }
  1367.             Acc.indice++;
  1368.             fscanf (file,"%s",&wordd);
  1369.             if (wordd[0] != '{') {
  1370.                 cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1371.                 getchar();
  1372.             }
  1373.             node_type = nothing;
  1374.             break;
  1375.         case n_geometry:
  1376.             //if (previous_node_type == n_Appearance) 
  1377.             //    Create_Brother_Node();
  1378.             //else Create_Son_Node();
  1379.             //CN->node_type = n_geometry;
  1380.             // strcpy (CN->Geometry_ins.texture.tex_blk , "nothing");
  1381.             previous_node_type = n_geometry; 
  1382.             node_type = nothing;
  1383.             break;
  1384.         case n_children:
  1385.             CN->FlagChildren = 1;
  1386.             AccTag = p_children;          // Keeps track of the Accolades
  1387.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1388.       if (Acc.indice >=1000){
  1389.         cout << "Change your indice, exit 378838789583" << endl;
  1390.         exit(1);
  1391.       }
  1392.             Acc.indice++;
  1393.             fscanf (file,"%s",&wordd);
  1394.             if (wordd[0] != '[') {
  1395.                 cout << "GetChar: An '[' is expected here!!! look for me 193843" << endl;
  1396.                 getchar();
  1397.             }
  1398.             break; 
  1399.         case n_Shape:
  1400.             if (CN->FlagChildren){
  1401.                 Create_Child_Shape_Node();
  1402.                 InitializeChildrenShape();
  1403.             } else {
  1404.                 Create_One_Shape_Node();
  1405.                 InitializeChildrenShape();
  1406.             }
  1407.             AccTag = p_Shape;              // Keeps track of the Accolades
  1408.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1409.       if (Acc.indice >=1000){
  1410.         cout << "Change your indice, exit 879384733" << endl;
  1411.         exit(1);
  1412.       }
  1413.             Acc.indice++;
  1414.             fscanf (file,"%s",&wordd);
  1415.             if (wordd[0] != '{') {
  1416.                 cout << "GetChar: An '{' is expected here!!! look for me 9593962" << endl;
  1417.                 getchar();
  1418.             }
  1419.             node_type = nothing;
  1420.             break; 
  1421.         case n_Material:
  1422.             //Create_Son_node();
  1423.             //CN->node_type = n_Material;
  1424.             AccTag = p_Material;          // Keeps track of the Accolades
  1425.             Acc.pipe[Acc.indice] = AccTag; // Keeps track of the Accolades
  1426.       if (Acc.indice >=1000){
  1427.         cout << "Change your indice, exit 98374983" << endl;
  1428.         exit(1);
  1429.       }
  1430.             Acc.indice++;
  1431.             fscanf (file,"%s",&wordd);
  1432.             if (wordd[0] != '{') {
  1433.                 cout << "GetChar: An '{' is expected here!!! look for me 991995922" << endl;
  1434.                 getchar();
  1435.             }
  1436.             break; 
  1437.         case n_translation:
  1438.             // the CN is a Transform.
  1439.             CN->translation.x = LisFlWrl();
  1440.             CN->translation.y = LisFlWrl();
  1441.             CN->translation.z = LisFlWrl();
  1442.             node_type = nothing;
  1443.             break; 
  1444.         case n_scale:
  1445.             // the CN is a Transform.
  1446.             CN->FlagScale = 1;
  1447.             CN->scale.x = LisFlWrl();
  1448.             CN->scale.y = LisFlWrl();
  1449.             CN->scale.z = LisFlWrl();
  1450.             node_type = nothing;
  1451.             break; 
  1452.         case n_rotation:
  1453.             // the CN is a Transform.
  1454.             CN->rotation.x = LisFlWrl();
  1455.             CN->rotation.y = LisFlWrl();
  1456.             CN->rotation.z = LisFlWrl();
  1457.             CN->rotation.teta = LisFlWrl();
  1458.             node_type = nothing;
  1459.             break; 
  1460.         case n_diffuseColor:
  1461.             CNS->Appearance.material.diffuseColor.r = LisFlWrl(); 
  1462.             CNS->Appearance.material.diffuseColor.g = LisFlWrl(); 
  1463.             CNS->Appearance.material.diffuseColor.b = LisFlWrl(); 
  1464.             node_type = nothing;
  1465.             break;
  1466.         case n_specularColor:
  1467.             CNS->Appearance.material.specularColor.r = LisFlWrl(); 
  1468.             CNS->Appearance.material.specularColor.g = LisFlWrl(); 
  1469.             CNS->Appearance.material.specularColor.b = LisFlWrl(); 
  1470.             node_type = nothing;
  1471.             break;
  1472.  
  1473.         case n_shininess:
  1474.             CNS->Appearance.material.shininess = LisFlWrl();
  1475.             node_type = nothing;
  1476.             break;
  1477.         case n_transparency:
  1478.             CNS->Appearance.material.transparency = LisFlWrl();
  1479.             node_type = nothing;
  1480.             break;
  1481.         case f_solid:
  1482.             fscanf (file,"%s",&wordd);
  1483.             if ((strcmp (wordd,"FALSE"))==0)
  1484.                 CNS->Geometry.IndexedFaceSet.solid = GL_FALSE;
  1485.             else
  1486.                 CNS->Geometry.IndexedFaceSet.solid = GL_TRUE;
  1487.             break;
  1488.         case n_TimeSensor:
  1489.                 AccTag = p_TimeSensor;          // Keeps track of the Accolades
  1490.                 Acc.pipe[Acc.indice] = AccTag;  // Keeps track of the Accolades
  1491.       if (Acc.indice >=1000){
  1492.         cout << "Change your indice, exit 3398439899903" << endl;
  1493.         exit(1);
  1494.       }
  1495.                 Acc.indice++;
  1496.             fscanf (file,"%s",&wordd);
  1497.             if (wordd[0] != '{') {
  1498.                 cout << "GetChar: An '{' is expected here!!! look for me 159383395" << endl;
  1499.                 getchar();
  1500.             }
  1501.                 break;
  1502.         case f_coord:
  1503.             fscanf (file,"%s",&wordd);
  1504.             if (strcmp(wordd,"DEF") == 0){
  1505.                 // Put Blk_name in hash table
  1506.                 iz.index++;
  1507.                 fscanf (file, "%s", &UseDef[iz.index].blk_name);
  1508.                 //UseDef[iz.index].coord = CNS->Geometry.IndexedFaceSet.coord; after reading the pts bozzo!!
  1509.                 // Reads p_Coordinate
  1510.                 fscanf (file,"%s",&wordd);
  1511.                 fscanf (file,"%s",&wordd);
  1512.                 if (wordd[0] != '{') {
  1513.                     cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1514.                     getchar();
  1515.                 }
  1516.                 AccTag = p_Coordinate;          // Keeps track of the Accolades
  1517.                 Acc.pipe[Acc.indice] = AccTag;  // Keeps track of the Accolades
  1518.       if (Acc.indice >=1000){
  1519.         cout << "Change your indice, exit 3899932878897213" << endl;
  1520.         exit(1);
  1521.       }
  1522.                 Acc.indice++;
  1523.                 fscanf (file,"%s",&wordd);
  1524.                 wordd_to_node(wordd); 
  1525.                 if (node_type == n_point)
  1526.                     read_all_points_into_this_node();
  1527.                 if (ADJUST_TO_ZERO) {
  1528.                     adjust_to_zero();
  1529.                     if ((UseDef[iz.index].blk_name[2] == 'R') && (UseDef[iz.index].blk_name[3] == 'O') &&
  1530.                         (UseDef[iz.index].blk_name[4] == 'A') && (UseDef[iz.index].blk_name[5] == 'D') )
  1531.                         Saffir_the_ys();
  1532.                     if ((UseDef[iz.index].blk_name[2] == 'C') && (UseDef[iz.index].blk_name[3] == 'U') &&
  1533.                         (UseDef[iz.index].blk_name[4] == 'R') && (UseDef[iz.index].blk_name[5] == 'B') )
  1534.                         Saffir_the_ys();
  1535.                     if ((UseDef[iz.index].blk_name[2] == 'T') && (UseDef[iz.index].blk_name[3] == 'O') &&
  1536.                         (UseDef[iz.index].blk_name[4] == 'P') && (UseDef[iz.index].blk_name[5] == 'O') )
  1537.                         Saffir_the_ys();
  1538.                 }
  1539.                 UseDef[iz.index].coord = CNS->Geometry.IndexedFaceSet.coord;
  1540.                 node_type = nothing;
  1541.             }
  1542.             // If Use coord of another blk, just put it in!!
  1543.             if (strcmp(wordd,"USE") == 0) {
  1544.                 fscanf (file, "%s", &wordd);
  1545.                 for (int ozi=0;ozi<=iz.index; ozi++)
  1546.                     if (strcmp(wordd,UseDef[ozi].blk_name) == 0)
  1547.                         CNS->Geometry.IndexedFaceSet.coord = UseDef[ozi].coord;
  1548.             }
  1549.             break;
  1550.         case n_normal:
  1551.             fscanf (file, "%s", &wordd); //reads Normal
  1552.             if (strcmp (wordd, "Normal") != 0) {
  1553.                 cout << "What could be in front of normal else then Normal?? " << endl;
  1554.                 exit (0);
  1555.                 getchar();
  1556.             }
  1557.  
  1558.             fscanf (file, "%s", &wordd); //reads {
  1559.             if (strcmp (wordd, "{") != 0) {
  1560.                 cout << "What could be in front of normal Normal else then {?? " << endl;
  1561.                 exit (0);
  1562.                 getchar();
  1563.             }
  1564.             else 
  1565.                 OuvreAccolade();
  1566.                 
  1567.             fscanf (file, "%s", &wordd); //reads vector
  1568.             if (strcmp (wordd, "vector") != 0) {
  1569.                 cout << "What could be in front of normal Normal { else then vector?? " << endl;
  1570.                 exit (0);
  1571.                 getchar();
  1572.             }
  1573.  
  1574.             read_all_normal_vectors_into_this_node();
  1575.  
  1576.  
  1577.             break;
  1578. #if 1
  1579.         case f_texture:
  1580.             fscanf (file,"%s",&wordd); // Reads ImageTexture
  1581.             fscanf (file,"%s",&wordd); // Reads {
  1582.             fscanf (file,"%s",&wordd); // Reads url
  1583.             fscanf (file,"%s",&wordd); // Reads "image.jpg"
  1584.             char directory [200];
  1585.             strcpy (directory, texture_directory);
  1586.             //strcpy (directory, "/home/grof/MRcode/ver.5.1/"); WinUnix
  1587.             //strcpy (directory, "D:\\Arch98_Class\\VirtCity\\"); // WinUnix
  1588.             //strcpy (directory, "D:\\Arch98_Class\\VirtCity\\"); // WinUnix
  1589.             for (i = 1; i<89; i++){
  1590.                 if (wordd[i] == '"') wordd[i] = wordd[i+1];
  1591.                 wordd[i-1] = wordd[i];
  1592.             }
  1593.             // cout << wordd  ;
  1594.             strcat (directory, wordd);
  1595.         /* read_RAW_file (directory, 
  1596.                             CNS->Texture.image, 
  1597.                             CNS->Texture.im_width, 
  1598.                             CNS->Texture.im_height); */
  1599.  
  1600.       CNS->Texture.image=
  1601.                 read_JPEG_file (wordd, directory, 
  1602.                             CNS->Texture.image, 
  1603.                             CNS->Texture.im_width, 
  1604.                             CNS->Texture.im_height); 
  1605.  
  1606.       // cout << CNS->Texture.im_width[0] << "  " << CNS->Texture.im_height[0] << endl;
  1607.  
  1608.       /* if ((CNS->Texture.image = (unsigned char *)malloc (128*256*3) ) == NULL){
  1609.         cout << "Malloc error here 9384793874" << endl;
  1610.         exit (1);
  1611.         }  
  1612.       for (int muk = 0 ; muk <64*64*3; muk++)
  1613.         CNS->Texture.image[muk] = 230;
  1614.         CNS->Texture.im_width[0]= 128; CNS->Texture.im_height[0]=256; */
  1615.       //cout << wordd << endl;
  1616.       //cout << CNS->Texture.im_width[0] << "  " << CNS->Texture.im_height[0] << endl; 
  1617.             glPixelStorei(GL_UNPACK_ALIGNMENT, 1);
  1618.             glClearColor (0.0, 0.0, 0.0, 0.0);
  1619.             glShadeModel(GL_FLAT);
  1620.             glEnable(GL_DEPTH_TEST);
  1621.             glGenTextures(1,CNS->Texture.tex_name);
  1622.             glBindTexture(GL_TEXTURE_2D, CNS->Texture.tex_name[0]);
  1623.             glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_S, GL_REPEAT); //GL_CLAMP); //GL_REPEAT);
  1624.             glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_WRAP_T, GL_REPEAT); //GL_CLAMP); //GL_REPEAT);
  1625.             glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MAG_FILTER, GL_NEAREST);
  1626.             glTexParameteri(GL_TEXTURE_2D, GL_TEXTURE_MIN_FILTER, GL_NEAREST);
  1627.             glTexImage2D(GL_TEXTURE_2D,           // target
  1628.                     0,                       // level
  1629.                     GL_RGB,                  // internal format
  1630.                     CNS->Texture.im_width[0], CNS->Texture.im_height[0],   
  1631.                     0,                       // border
  1632.                     GL_RGB,                  // format
  1633.                     GL_UNSIGNED_BYTE,        // type
  1634.                     CNS->Texture.image);
  1635.             fscanf (file,"%s",&wordd); // Reads }
  1636.             if (wordd[0] != '}') { cout << "error reading }..!! getchar()" << endl;
  1637.                             getchar();
  1638.             }
  1639.         break;
  1640.         case f_texCoord:
  1641.             fscanf (file, "%s", &wordd);
  1642.             if (strcmp (wordd, "DEF") == 0) {
  1643.                 iz.texIndex++;
  1644.                 fscanf (file, "%s", &UseDef[iz.texIndex].tex_blk_name);
  1645.                 strcpy(CNS->Texture.tex_blk,UseDef[iz.texIndex].tex_blk_name);
  1646.                 fscanf (file, "%s", &wordd);
  1647.                 if (strcmp (wordd, "TextureCoordinate") !=0){
  1648.                     cout << "Getchar! Que peut il y avoir apres texCoord qui ne soit pas DEF, USE ou TextureCoordinate?" << endl;
  1649.                     getchar();
  1650.                 }
  1651.                 AccTag = p_TextureCoordinate ;  // Keeps track of the Accolades
  1652.                 Acc.pipe[Acc.indice] = AccTag;  // Keeps track of the Accolades        
  1653.                 if (Acc.indice >=1000){
  1654.         cout << "Change your indice, exit 9899399399587778" << endl;
  1655.         exit(1);
  1656.       }
  1657.                 Acc.indice++;
  1658.                 fscanf (file,"%s",&wordd);
  1659.                 if (wordd[0] != '{') {
  1660.                     cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1661.                     getchar();
  1662.                 }
  1663.                 fscanf (file,"%s",&wordd);
  1664.                 wordd_to_node(wordd); 
  1665.                 if (node_type == n_point)
  1666.                     read_all_texture_points_into_this_node();
  1667.                 else {
  1668.                     cout << "Getchar! Que peut il y avoir ici qui ne soit pas point???" << endl;
  1669.                     getchar();
  1670.                 }
  1671.                 UseDef[iz.texIndex].texCoord = CNS->Geometry.IndexedFaceSet.texCoord;
  1672.             }else
  1673.             if (strcmp (wordd, "USE") == 0){
  1674.                 //Type_Transform *CN_UP = CN;
  1675.                 fscanf (file,"%s", &wordd);
  1676.                 for (int ozi=0;ozi<=iz.texIndex; ozi++)
  1677.                     if (strcmp(wordd,UseDef[ozi].tex_blk_name) == 0)
  1678.                         CNS->Geometry.IndexedFaceSet.texCoord = UseDef[ozi].texCoord;
  1679.                 //while (strcmp (CN_UP->Geometry_ins.texture.tex_blk , wordd) != 0)
  1680.                 //    CN_UP = CN_UP->Father;
  1681.                 //CN->Geometry_ins.IndexedFaceSet_Ins.texCoord = CN_UP->Geometry_ins.IndexedFaceSet_Ins.texCoord;
  1682.             }else
  1683.             if (strcmp (wordd, "TextureCoordinate") == 0) {
  1684.                 fscanf (file,"%s",&wordd);
  1685.                 if (wordd[0] != '{') {
  1686.                     cout << "GetChar: An '{' is expected here!!! look for me 159395" << endl;
  1687.                     getchar();
  1688.                 }
  1689.                 AccTag = p_TextureCoordinate ;  // Keeps track of the Accolades
  1690.                 Acc.pipe[Acc.indice] = AccTag;  // Keeps track of the Accolades
  1691.       if (Acc.indice >=1000){
  1692.         cout << "Change your indice, exit 99887763325621" << endl;
  1693.         exit(1);
  1694.       }
  1695.                 Acc.indice++;
  1696.                 fscanf (file,"%s",&wordd);
  1697.                 wordd_to_node(wordd); 
  1698.                 if (node_type == n_point)
  1699.                     read_all_texture_points_into_this_node();
  1700.                 else {
  1701.                     cout << "Getchar! Que peut il y avoir ici qui ne soit pas point???" << endl;
  1702.                     getchar();
  1703.                 }
  1704.             }
  1705.                 node_type = nothing;
  1706.             break;
  1707.         case f_texCoordIndex:
  1708.             read_all_texCoordIndex_into_CN();
  1709.             break;
  1710. #endif
  1711.  
  1712.         case n_coordIndex:
  1713.       //cout << "reading coordIndex " << endl;
  1714.             read_all_coordIndex_into_CN();
  1715.             break;
  1716.         case f_colorIndex:
  1717.              read_all_colorIndex_into_CN();
  1718.              break;
  1719.         case f_color:
  1720.             fscanf (file,"%s",&wordd);
  1721.             if ((strcmp (wordd,"["))==0)
  1722.                 read_color_values();
  1723.             break;
  1724.         }
  1725.     }
  1726.     // Maybe the last node is a Geometry node and has the coord field not set yet.
  1727.     //if (CN->node_type == n_geometry) Test_For_Coord_In_Geometry_node();
  1728.  
  1729.   //cout << "After reading the vrml find there was " << counter2 << "transform nodes" << endl;
  1730.   //cout << "Number of nodes transform Created" << counter3 << endl;
  1731.     return (CN);
  1732.  
  1733. }
  1734.  
  1735. void test_bassit()
  1736. {
  1737.     float xmin = -500;
  1738.     float xmax =  500;
  1739.     float ymin = -500;
  1740.     float ymax =  500;
  1741.     float zmin = -500;
  1742.     float zmax =  500;
  1743.  
  1744. /*    glPopMatrix();
  1745.    glPushMatrix(); */
  1746.  
  1747.     glBegin(GL_LINES);
  1748.     glColor3f (0.45,0.45,0.45);
  1749.     glVertex3f (xmin,ymax,zmax);
  1750.     glVertex3f (xmax,ymax,zmax);
  1751.     glEnd();
  1752.  
  1753.     glBegin(GL_LINES);
  1754.     glColor3f (0.45,0.45,0.45);
  1755.     glVertex3f (xmax,ymax,zmax);
  1756.     glVertex3f (xmax,ymin,zmax);
  1757.     glEnd();
  1758.  
  1759.     glBegin(GL_LINES);
  1760.     glColor3f (0.45,0.45,0.45);
  1761.     // glColor3f (0.5,0,0);
  1762.     glVertex3f (xmax,ymin,zmax);
  1763.     glVertex3f (xmin,ymin,zmax);
  1764.     glEnd();
  1765.  
  1766.     glBegin(GL_LINES);
  1767.     glColor3f (0.45,0.45,0.45);
  1768.     // glColor3f (0.0,0.5,0.0);
  1769.     glVertex3f (xmin,ymin,zmax);
  1770.     glVertex3f (xmin,ymax,zmax);
  1771.     glEnd();
  1772.  
  1773. /*    glBegin(GL_LINES);
  1774.     glColor3f (0.45,0.45,0.45);
  1775.     // glColor3f (0,0,0.8);
  1776.     glVertex3f (xmin,ymax,zmax);
  1777.     glVertex3f (xmin,ymax,zmin);
  1778.     glEnd();
  1779.  
  1780.     glBegin(GL_LINES);
  1781.     glColor3f (0.45,0.45,0.45);
  1782.     //glColor3f (0,0,0.8);
  1783.     glVertex3f (xmax,ymax,zmax);
  1784.     glVertex3f (xmax,ymax,zmin);
  1785.     glEnd();
  1786.  
  1787.     glBegin(GL_LINES);
  1788.     glColor3f (0,0,0.8);
  1789.     glVertex3f (xmax,ymin,zmax);
  1790.     glVertex3f (xmax,ymin,zmin);
  1791.     glEnd();
  1792.  
  1793.     glBegin(GL_LINES);
  1794.     glColor3f (0,0,0.8);
  1795.     glVertex3f (xmin,ymin,zmax);
  1796.     glVertex3f (xmin,ymin,zmin);
  1797.     glEnd();*/
  1798.  
  1799.     glTranslatef(0.0,0.0,-1000.0);
  1800.  
  1801.     glBegin(GL_LINES);
  1802.     glColor3f (0.0,0.5,0.0);
  1803.     glVertex3f (xmin,ymax,zmax);
  1804.     glVertex3f (xmin,ymin,zmax);
  1805.     glEnd();
  1806.  
  1807.     glBegin(GL_LINES);
  1808.     glColor3f (0.5,0,0);
  1809.     glVertex3f (xmin,ymin,zmax);
  1810.     glVertex3f (xmax,ymin,zmax);
  1811.     glEnd();
  1812.  
  1813.     glBegin(GL_LINES);
  1814.     glColor3f (0.0,0.5,0.0);
  1815.     glVertex3f (xmax,ymin,zmax);
  1816.     glVertex3f (xmax,ymax,zmax);
  1817.     glEnd();
  1818.  
  1819.     glBegin(GL_LINES);
  1820.     glColor3f (0.5,0,0);
  1821.     glVertex3f (xmax,ymax,zmax);
  1822.     glVertex3f (xmin,ymax,zmax);
  1823.     glEnd();
  1824. }
  1825.  
  1826. void ReadBufferSphere(Type_Transform *CNSphere)
  1827. {
  1828.     // must be in this file for some puttain unknown reason!!!
  1829.    glReadBuffer(GL_BACK);    
  1830.    glReadPixels (0,0, CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.im_width[0],
  1831.                       CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.im_height[0], 
  1832.                       GL_RGB, GL_UNSIGNED_BYTE,
  1833.                       CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.image); // image[0][0]); 
  1834.    glBindTexture (GL_TEXTURE_2D, CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.tex_name[0]);
  1835.     glTexImage2D(  GL_TEXTURE_2D,    // target
  1836.                      0,                // level
  1837.                      //0, 0,                              // xoffset, yoffset
  1838.                      GL_RGB, // format
  1839.                      CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.im_width[0], 
  1840.                      CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.im_height[0],             // Width, Height
  1841.                      //im_width[0], im_height[0],
  1842.                      0, // border
  1843.                      GL_RGB,           // format
  1844.                      GL_UNSIGNED_BYTE, // type
  1845.                      //image);
  1846.                      CNSphere->ChildrenTransform[0]->ChildrenShape[0]->Texture.image);        // GLvoid *pixels */ 
  1847.  
  1848. }
  1849.