Changeset 1007 for GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/vmi/src/simplify.cpp

Timestamp:

06/09/06 08:24:19 (18 years ago)

Author:

gumbau

Message:

 

File:

• GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/vmi/src/simplify.cpp (modified) (1 diff)

GTP/trunk/Lib/Geom/shared/GTGeometry/src/libs/vmi/src/simplify.cpp

@@ -424 +424 @@
 ///////////////////////////////////////////////////////////////////////////////
 
-void VMI::simplifyModel(Mesh *mesh, GLuint numDemandedTri) {
-    int e;
-    Change *c;
-    FILE* file = NULL;
-    char s[MAX_CHAR];
-    double cost = 0.0;
-    bheap_t *h = NULL;
-    float       percent =       (40.0) / (float)(mesh->numTriangles - numDemandedTri);
-
-    if (numDemandedTri > mesh->currentNumTriangles) return;
-
-    if (bSaveLog == TRUE) { // Save changes into a file
+void VMI::simplifyModel(Mesh *mesh, GLuint numDemandedTri)
+{
+        int                     e;
+        Change  *c;
+        FILE            *file                                   = NULL;
+        FILE            *file_simp_seq  =       NULL;
+        char            s[MAX_CHAR];
+        double  cost                                            = 0.0;
+        bheap_t *h                                                      = NULL;
+        float           percent =       (40.0) / (float)(mesh->numTriangles - numDemandedTri);
+
+        if (numDemandedTri > mesh->currentNumTriangles)
+        {
+                return;
+        }
+
+        // Save changes into a file
+        if (bSaveLog == TRUE)
+        {
 #ifdef SALIENCY       
 #ifdef KL // Kullback-Leibler
-        sprintf(s,"%s_VKL_S.log", filename);
+                sprintf(s,"%s_VKL_S.log", filename);
 #endif
 #ifdef MI // Mutual Information
-        sprintf(s,"%s_VMI_S.log", filename);
+                sprintf(s,"%s_VMI_S.log", filename);
 #endif
 #ifdef HE // Hellinger
-        sprintf(s,"%s_VHE_S.log", filename);
+                sprintf(s,"%s_VHE_S.log", filename);
 #endif
 #ifdef CS // Chi-Square
-        sprintf(s,"%s_VCS_S.log", filename);
+                sprintf(s,"%s_VCS_S.log", filename);
 #endif
 #else
 #ifdef KL // Kullback-Leibler
-        sprintf(s,"%s_VKL.log", filename);
+                sprintf(s,"%s_VKL.log", filename);
 #endif
 #ifdef MI // Mutual Information
-        sprintf(s,"%s_VMI.log", filename);
+                sprintf(s,"%s_VMI.log", filename);
 #endif
 #ifdef HE // Hellinger
-        sprintf(s,"%s_VHE.log", filename);
+                sprintf(s,"%s_VHE.log", filename);
 #endif
 #ifdef CS // Chi-Square
-        sprintf(s,"%s_VCS.log", filename);
-#endif
-#endif
-        
-        glmWriteOBJ(pmodel, s, GLM_NONE);
-        
-        /* open the file */
-        file = fopen(s, "a+");
-        if (!file) {
-            fprintf(stderr, "simplifyModel() failed: can't open file \"%s\" to write.\n", s);
-            exit(1);
-        }
-    }
-
-    h = initHeap(mesh);
-
-    printf("Starting simplification...\n");
-
-    while (mesh->currentNumTriangles > numDemandedTri /*&& cost == 0.0*/) {
-        
-        // Get the edge that has the minimum cost
-        e = bh_min(h);
-
-        c = createChange(mesh, e);
-
-        // Apply the edge collapse
-        computeChanges(mesh, c);
-        doChange(mesh, c);
-
-        if (bSaveLog == TRUE) writeChange(file, c);
-
-        cost = h->a[h->p[e]].key;
-        printf("Edge collapse e%d(%d,%d) %f MIN d %d m %d\n",c->e, c->u, c->v, cost, c->numDel, c->numMod);
-
-        mesh->edges[e].enable = FALSE;
-        bh_delete(h, e);
-        
-        // Get projected areas after the edge collapse
-        getProjectedAreas(histogram, numCameras);
-
-        computeCameraIs(histogram, numCameras, initialIs);
-
-        // Update the heap according to the edge collapse
-        h = updateHeap(h, mesh, c);
-
-                                mUPB((float)c->numDel * percent);
-
-        deleteChange(c);
-
-        printf("t %d\n", mesh->currentNumTriangles);
-    }
-
-    if (bSaveLog == TRUE) {
-        fclose(file);
-        printf("Log file written...Ok\n");
-    }
-
-    bh_free(h);
-}
+                sprintf(s,"%s_VCS.log", filename);
+#endif
+#endif
+
+                glmWriteOBJ(pmodel, s, GLM_NONE);
+
+                /* open the file */
+                file = fopen(s, "a+");
+                
+                if (!file)
+                {
+                        fprintf(stderr, "simplifyModel() failed: can't open file \"%s\" to write.\n", s);
+                        exit(1);
+                }
+        }
+
+        //      Open file of simplification sequence.
+        file_simp_seq   =       fopen("SimplifSequence.txt", "w");
+        
+        if (!file_simp_seq)
+        {
+                fprintf(stderr,
+                                                "simplifyModel() failed: ",
+                                                "can't open file \"SimplifSequence\" to write.\n");
+        }
+        
+        h = initHeap(mesh);
+
+        printf("Starting simplification...\n");
+
+        while (mesh->currentNumTriangles > numDemandedTri /*&& cost == 0.0*/)
+        {
+                // Get the edge that has the minimum cost
+                e = bh_min(h);
+
+                c = createChange(mesh, e);
+
+                // Apply the edge collapse
+                computeChanges(mesh, c);
+                doChange(mesh, c);
+
+                if (bSaveLog == TRUE) writeChange(file, c);
+
+                //      Write Simplification sequence.
+                saveSimplificationSequence(c);
+                
+                cost = h->a[h->p[e]].key;
+                
+                printf( "Edge collapse e%d(%d,%d) %f MIN d %d m %d\n",
+                                                c->e,
+                                                c->u,
+                                                c->v,
+                                                cost,
+                                                c->numDel,
+                                                c->numMod);
+
+                mesh->edges[e].enable = FALSE;
+                
+                bh_delete(h, e);
+
+                // Get projected areas after the edge collapse
+                getProjectedAreas(histogram, numCameras);
+
+                computeCameraIs(histogram, numCameras, initialIs);
+
+                // Update the heap according to the edge collapse
+                h = updateHeap(h, mesh, c);
+
+                mUPB((float)c->numDel * percent);
+
+                deleteChange(c);
+
+                printf("t %d\n", mesh->currentNumTriangles);
+        }
+
+        if (bSaveLog == TRUE)
+        {
+                fclose(file);
+                printf("Log file written...Ok\n");
+        }
+
+        //      Close simplification sequence file.
+        fclose(file_simp_seq);
+
+        bh_free(h);
+}
+