Changeset 2288Tweet

Timestamp:

Apr 16, 2008, 12:11:38 AM (15 years ago)

Author:

Sam Hocevar

Message:

• Add mode 8, for the very first graphic in the paper.
• Add results for mode 8 in part0/lena-values.txt.
• Make the bytecode version use less memory.

Location:

research/2008-displacement

Files:

• README (modified) (2 diffs)
• main-cpushare (modified) (3 diffs)
• main.c (modified) (8 diffs)
• part0 (added)
• part0/lena-values.txt (added)

research/2008-displacement/README

@@ -16 +16 @@
 
 # Get phase 3 and phase 4 stuff
-ssh canard.zoy.org "cd test-20080329; tar cz *raster.txt *serp.txt" | tar xz
+# .. from canard
+ssh canard.zoy.org "cd test-20080329; tar cz *raster.txt" | tar xz
+ssh canard.zoy.org "cd test-20080329; tar cz *serp.txt" | tar xz
 for x in *-raster.txt; do y="$x"; y="${y%%-raster.txt}"; y="${y%%.tiff}"; y="${y##usc-sipi}"; \mv "$x" part3/"$y".txt; done
 for x in *-serp.txt; do y="$x"; y="${y%%-serp.txt}"; y="${y%%.tiff}"; y="${y##usc-sipi}"; \mv "$x" part4/"$y".txt; done
+# .. from poulet (cpushare)
+scp poulet.zoy.org:cpushare/'test*-*.txt' .
+for f in test6-*.txt; do grep '###' $f|cut -f2 -d'`'|cut -f1 -d"'"|tr / .|while read i ; do sed -ne '/`'$i'/,/limit/p' $f | cut -b61- >| part3/$(echo $i | cut -f5 -d.).txt ; done; done; rm -f test6-*.txt
+for f in test7-*.txt; do grep '###' $f|cut -f2 -d'`'|cut -f1 -d"'"|tr / .|while read i ; do sed -ne '/`'$i'/,/limit/p' $f | cut -b61- >| part4/$(echo $i | cut -f5 -d.).txt ; done; done; rm -f test7-*.txt
 
 # Condorcet voting for part 3 and 4
@@ -38 +44 @@
 cat part4/aerials2.1.01.txt | grep K | while read x k y ; do grep $k part4/* | awk '{ dx+=$10; dy+=$12; n+=1 } END { print dx/n, dy/n }' | read dx dy; echo "$k $(grep $k part4/* | awk 'BEGIN { dx='$dx'; dy='$dy' } { x=dx-$10; y=dy-$12; t+=2.^-5*(x*x+y*y); a+=t*$8; n+=t } END { print a/n, n }')"; done | sort -nk2 | head -20
 
+# Plot lena displacement
+set contour
+unset surface
+set view 0,0
+set cntrparam levels discrete 0.31, 0.35, 0.42, 0.52, 0.7, 1, 1.5
+set xrange [-1:1]
+set yrange [1:-1]
+set mxtics 2
+set mytics 2
+set size square
+set xzeroaxis
+set yzeroaxis
+set xlabel "dx" font "Italic,32"
+splot 'part0/lena-values.txt' with lines
+
+
+#set border 0
+set style line 1 lt 1 lw 1
+set style line 2 lt 2 lw 1 
+set style line 3 lt 3 lw 1 
+set style line 4 lt 6 lw 1 
+set style line 5 lt 1 lw 3 
+set style line 6 lt 2 lw 3 
+set style line 7 lt 3 lw 3 
+set style line 8 lt 6 lw 3
 ###
 ###

research/2008-displacement/main-cpushare

@@ -17 +17 @@
     print 'Usage: %s <buy_order.cpu> -<mode> [image list]' % (sys.argv[0],)
     sys.exit(-1)
-mode = sys.argv[2][1]
+mode = sys.argv[2]
 done = 0
 source = sys.argv[3:]
@@ -34 +34 @@
 
     def start(self):
-        global todo, source
+        global mode, todo, source
+        gd.gdMaxColors = 256 * 256 * 256
+        while todo < len(source):
+            try:
+                self.index = todo
+                im = gd.image(source[self.index])
+                break
+            except:
+                todo += 1
+                self.result = "Error\n"
         if todo >= len(source):
-            return # We're finished... FIXME: is the transaction kept stuck?
-        self.index = todo
+            return # We're finished... FIXME: isn't the transaction stuck?
         todo += 1
-        gd.gdMaxColors = 256 * 256 * 256
-        im = gd.image(source[self.index])
         # Send argument count
         self.protocol.sendString(PROTO_SECCOMP_FORWARD + chr(3))
@@ -46 +52 @@
         msg = chr(3)
         msg += "bytecode\0"
-        msg += "-1\0"
+        msg += mode + "\0"
         msg += source[self.index] + "\0"
         self.protocol.sendString(PROTO_SECCOMP_FORWARD + msg)

research/2008-displacement/main.c

@@ -25 +25 @@
 
 #ifdef BYTECODE
-#   define MAXIMAGES 128
+#   define MAXIMAGES 6
 static int slots[MAXIMAGES];
 static double slotbuf[MAXIMAGES * MAXWIDTH * MAXHEIGHT];
@@ -383 +383 @@
 {
     double *floats;
-    int x, y;
+    int w, h, x, y;
 
 #ifdef BYTECODE
@@ -390 +390 @@
     if(sys_read(0, &c, 1) != 1)
         sys_exit(-5);
-    WIDTH = ((int)(unsigned char)c) << 8;
+    w = ((int)(unsigned char)c) << 8;
     if(sys_read(0, &c, 1) != 1)
         sys_exit(-5);
-    WIDTH |= (int)(unsigned char)c;
+    w |= (int)(unsigned char)c;
 
     if(sys_read(0, &c, 1) != 1)
         sys_exit(-5);
-    HEIGHT = ((int)(unsigned char)c) << 8;
+    h = ((int)(unsigned char)c) << 8;
     if(sys_read(0, &c, 1) != 1)
         sys_exit(-5);
-    HEIGHT |= (int)(unsigned char)c;
+    h |= (int)(unsigned char)c;
+#else
+    SDL_Surface *tmp, *surface;
+    uint32_t *pixels;
+
+    tmp = IMG_Load(name);
+    if(!tmp)
+        return NULL;
+
+    w = tmp->w;
+    h = tmp->h;
+#endif
+
+    WIDTH = w > MAXWIDTH ? MAXWIDTH : w;
+    HEIGHT = h > MAXHEIGHT ? MAXHEIGHT : h;
 
     floats = new();
@@ -406 +420 @@
         return NULL;
 
-    for(y = 0; y < HEIGHT; y++)
-    for(x = 0; x < WIDTH; x++)
+#ifdef BYTECODE
+    for(y = 0; y < h; y++)
+    for(x = 0; x < w; x++)
     {
         if(sys_read(0, &c, 1) != 1)
             sys_exit(-5);
+        if(x >= WIDTH || y >= HEIGHT)
+            continue;
         put(floats, x, y, (double)(unsigned char)c / 0xff);
-    }   
+    }
 #else
-    SDL_Surface *tmp, *surface;
-    uint32_t *pixels;
-
-    tmp = IMG_Load(name);
-    if(!tmp)
-        return NULL;
-
-    WIDTH = tmp->w > MAXWIDTH ? MAXWIDTH : tmp->w;
-    HEIGHT = tmp->h > MAXHEIGHT ? MAXHEIGHT : tmp->h;
-    floats = new();
-    if(!floats)
-        return NULL;
-
     surface = SDL_CreateRGBSurface(SDL_SWSURFACE, WIDTH, HEIGHT, 32,
                                    0xff0000, 0xff00, 0xff, 0x0);
@@ -859 +863 @@
     msg(" -6 <src>           restrained (a,b,c,d) ED kernel analysis on src\n");
     msg(" -7 <src>           restrained displacement study on src\n");
+    msg(" -8 <src>           displacement values on src\n");
 }
 
@@ -891 +896 @@
         {
             double *dest = ed(src, false, 7, 0,
-                                0, 3, 5, 1, 0,
-                                0, 0, 0, 0, 0);
+                                 0, 3, 5, 1, 0,
+                                 0, 0, 0, 0, 0);
             study(src, dest, 1.2, 0.001, .16, .28);
             del(dest);
@@ -1068 +1073 @@
         break;
 
+        case 8:
+        {
+            const int STEP = 32;
+            double mat[NN][NN];
+            double *dest = ed(src, false, 7, 0,
+                                 0, 3, 5, 1, 0,
+                                 0, 0, 0, 0, 0);
+            double *tmp, *tmp2;
+            int dx, dy;
+
+            makegauss(mat, 1.2, 0., 0.);
+            tmp = gauss(src, mat);
+            for(dy = 0; dy <= STEP; dy++)
+            {
+                for(dx = 0; dx <= STEP; dx++)
+                {
+                    double fy = 2. / STEP * (dy - STEP / 2.);
+                    double fx = 2. / STEP * (dx - STEP / 2.);
+
+                    makegauss(mat, 1.2, fx, fy);
+                    tmp2 = gauss(dest, mat);
+                    msg("%g %g %g\n", fy, fx, 1000. * dist(tmp, tmp2, 1.));
+                    del(tmp2);
+                }
+                msg("\n");
+            }
+            del(tmp);
+            del(dest);
+            del(src);
+        }
+        break;
+
 #if 0
     tmp = ed(src, 7, 0, 0, 3, 5, 1, 0, 0, 0, 0, 0, 0);
@@ -1099 +1136 @@
 
 #if 0
-#   define STEP 32
-    dest = ed(src, 7, 0, 0, 3, 5, 1, 0, 0, 0, 0, 0, 0);
-    makegauss(mat, 1.2, 0., 0.);
-    tmp = gauss(src, mat);
-    for(dy = 0; dy < STEP; dy++)
-    {
-        for(dx = 0; dx < STEP; dx++)
-        {
-            double fy = 2. / STEP * (dy - STEP / 2.);
-            double fx = 2. / STEP * (dx - STEP / 2.);
-
-            makegauss(mat, 1.2, fx, fy);
-            tmp2 = gauss(dest, mat);
-            msg("%g %g %g\n", fy, fx, 1000. * dist(tmp, tmp2, 1.));
-            fflush(stdout);
-            del(tmp2);
-        }
-        msg("\n");
-    }
 
     save(dest, "output.bmp");