Changeset 2278Tweet

Timestamp:

Apr 16, 2008, 12:09:44 AM (15 years ago)

Author:

Sam Hocevar

Message:

• Reactivate full ED analysis (raster FS, raster JaJuNi?, serpentine Ostromoukhov and serpentine optimum) as mode 5.
• Add a script to perform mode 5 on a directory's contents.

Location:

research/2008-displacement

Files:

• main.c (modified) (10 diffs)
• study-5.sh (added)

research/2008-displacement/main.c

@@ -512 +512 @@
 
 static void study(float const *src, float const *dest,
-                  float sigma, float precision)
+                  float sigma, float precision, float fdx, float fdy)
 {
 #   define Z 3
-    float mat[NN][NN];
+    float mat[NN][NN], mat2[NN][NN];
     float *tmp, *tmp2;
-    float e, e0, best = 1., fx = -1., fy = -1., step = 2., bfx = 0., bfy = 0.;
+    float e, e0, e1;
+    float best = 1., fx = -1., fy = -1., step = 2., bfx = 0., bfy = 0.;
     int dx, dy;
 
     makegauss(mat, sigma, 0., 0.);
     tmp = gauss(src, mat);
+
     tmp2 = gauss(dest, mat);
-
     e0 = dist(tmp, tmp2, 1.);
+    free(tmp2);
+
+    makegauss(mat2, sigma, fdx, fdy);
+    tmp2 = gauss(dest, mat2);
+    e1 = dist(tmp, tmp2, 1.);
     free(tmp2);
 
@@ -551 +557 @@
     free(tmp);
 
-    printf("E = %g E_min = %g dx = %g dy = %g\n",
-           1000 * e0, 1000 * best, fx, fy);
+    printf("E = %g E_fast = %g E_min = %g dx = %g dy = %g\n",
+           1000. * e0, 1000. * e1, 1000. * best, fx, fy);
     fflush(stdout);
 }
@@ -576 +582 @@
     fprintf(stderr, " -3 <src>           quick (a,b,c,d) ED kernel analysis on src\n");
     fprintf(stderr, " -4 <src>           exhaustive (a,b,c,d) ED kernel analysis on src\n");
+    fprintf(stderr, " -5 <src>           exhaustive displacement study on src\n");
 }
 
@@ -601 +608 @@
         return 2;
 
+    printf("### mode %i on `%s' ###\n", mode, argv[2]);
+
     switch(mode)
     {
         case 1:
         {
-            float *dest = ed(src, false, 7, 0, 1, 3, 5, 0, 0, 0, 0, 0, 0, 0);
-            study(src, dest, 1.2, 0.001);
+            float *dest = ed(src, false, 7, 0,
+                                0, 3, 5, 1, 0,
+                                0, 0, 0, 0, 0);
+            study(src, dest, 1.2, 0.001, .16, .28);
             free(dest);
             free(src);
@@ -625 +636 @@
             {
                 tmp = merge(src, src2, (float)i / 100.);
-                dest = ed(tmp, false, 7, 0, 1, 3, 5, 0, 0, 0, 0, 0, 0, 0);
-                study(tmp, dest, 1.2, 0.001);
+                dest = ed(tmp, false, 7, 0,
+                             0, 3, 5, 1, 0,
+                             0, 0, 0, 0, 0);
+                study(tmp, dest, 1.2, 0.001, .16, .28);
                 free(dest);
                 free(tmp);
@@ -638 +651 @@
         case 4:
         {
-            float mat0[NN][NN];
+            float mat[NN][NN];
             float *dest, *tmp, *tmp2;
             int a, b, c, d, e;
 
-            makegauss(mat0, 1.2, 0, 0);
-            for(a = 0; a <= 16; a++)
-            for(b = 0; b <= 16; b++)
-            for(c = 0; c <= 16; c++)
-            for(d = 0; d <= 16; d++)
-            for(e = 0; e <= 16; e++)
-            {
-                if(a + b + c + d + e != 16)
+#define TOTAL 16
+            makegauss(mat, 1.2, 0, 0);
+            for(a = 0; a <= TOTAL; a++)
+            for(b = 0; b <= TOTAL; b++)
+            for(c = 0; c <= TOTAL; c++)
+            for(d = 0; d <= TOTAL; d++)
+            for(e = 0; e <= TOTAL; e++)
+            {
+                if(a + b + c + d + e != TOTAL)
                     continue;
 
@@ -657 +671 @@
 #define SHUFFLE(p,q,n) \
     if(p+q) { int tmp = p+q; p = (p+n) % (tmp+1); q = tmp-p; }
-                SHUFFLE(c2, d2, 7); SHUFFLE(b2, d2, 5);
-                SHUFFLE(a2, d2, 3); SHUFFLE(b2, c2, 2);
-                SHUFFLE(a2, e2, 3); SHUFFLE(b2, e2, 2);
-                SHUFFLE(a2, c2, 4); SHUFFLE(a2, b2, 6);
-                SHUFFLE(c2, e2, 9); SHUFFLE(d2, e2, 7);
-                SHUFFLE(a2, d2, 9); SHUFFLE(a2, b2, 7);
+                SHUFFLE(c2, d2, 777); SHUFFLE(b2, d2, 555);
+                SHUFFLE(a2, d2, 333); SHUFFLE(b2, c2, 222);
+                SHUFFLE(a2, e2, 333); SHUFFLE(b2, e2, 222);
+                SHUFFLE(a2, c2, 444); SHUFFLE(a2, b2, 666);
+                SHUFFLE(c2, e2, 999); SHUFFLE(d2, e2, 777);
+                SHUFFLE(a2, d2, 999); SHUFFLE(a2, b2, 777);
 
 #if 0
@@ -674 +688 @@
                 printf("K = %d,%d,%d,%d ", a2, c2, d2, e2);
 
-                dest = ed(src, false, a2, 0, b2, c2, d2, e2, 0, 0, 0, 0, 0, 0);
-                if(mode == 3)
-                {
-                    tmp = gauss(src, mat0);
-                    tmp2 = gauss(dest, mat0);
+                dest = ed(src, false, a2, 0,
+                          b2, c2, d2, e2, 0,
+                           0,  0,  0,  0, 0);
+                if(mode == 4)
+                {
+                    study(src, dest, 1.2, 0.001, 0., 0.);
+                }
+                else
+                {
+                    tmp = gauss(src, mat);
+                    tmp2 = gauss(dest, mat);
                     printf("E = %.5g\n", 1000. * dist(tmp, tmp2, 1.));
                     free(tmp);
                     free(tmp2);
                 }
-                else
-                    study(src, dest, 1.2, 0.001);
                 fflush(stdout);
                 free(dest);
@@ -690 +708 @@
 
             free(src);
+        }
+        break;
+
+        case 5:
+        {
+            float *dest;
+
+            dest = ed(src, false, 7, 0,
+                         0, 3, 5, 1, 0,
+                         0, 0, 0, 0, 0);
+            study(src, dest, 1.2, 0.001, .16, .28);
+            free(dest);
+
+            dest = ed(src, false, 7, 5,
+                         3, 5, 7, 5, 3,
+                         1, 3, 5, 3, 1);
+            study(src, dest, 1.2, 0.001, .26, .76);
+            free(dest);
+
+            dest = ostromoukhov(src);
+            study(src, dest, 1.2, 0.001, .0, .19);
+            free(dest);
+
+            dest = ed(src, true, 2911, 0,
+                  0, 1373, 3457, 2258, 0,
+                  0,    0,    0,    0, 0);
+            study(src, dest, 1.2, 0.001, .0, .34);
         }
         break;
@@ -768 +813 @@
 
 #if 0
-    printf("%s: ", argv[1]);
-    makegauss(mat0, 1.2, 0, 0);
-    tmp = gauss(src, mat0);
-
-    dest = ed(src, false, 7, 0, 0, 3, 5, 1, 0, 0, 0, 0, 0, 0);
-    tmp2 = gauss(dest, mat0);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    makegauss(mat, 1.2, 0.16, 0.28);
-    tmp2 = gauss(dest, mat);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    free(dest);
-
-    dest = ed(src, false, 7, 5, 3, 5, 7, 5, 3, 1, 3, 5, 3, 1);
-    tmp2 = gauss(dest, mat0);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    makegauss(mat, 1.2, 0.26, 0.76);
-    tmp2 = gauss(dest, mat);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    free(dest);
-
-    dest = ostromoukhov(src);
-    tmp2 = gauss(dest, mat0);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    makegauss(mat, 1.2, 0.0, 0.19);
-    tmp2 = gauss(dest, mat);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    free(dest);
-
-    dest = ed(src, true, 2911, 0, 0, 1373, 3457, 2258, 0, 0, 0, 0, 0, 0);
-    tmp2 = gauss(dest, mat0);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    makegauss(mat, 1.2, 0.0, 0.34);
-    tmp2 = gauss(dest, mat);
-    printf("%.5g ", 1000. * dist(tmp, tmp2, 1.));
-    free(tmp2);
-    free(dest);
-
-    printf("\n");
-#endif
-
-#if 0
     printf("%s\n", argv[1]);