Index: /libpipi/trunk/examples/edd.c
===================================================================
--- /libpipi/trunk/examples/edd.c	(revision 2646)
+++ /libpipi/trunk/examples/edd.c	(revision 2647)
@@ -38,5 +38,5 @@
     double sigma = 1.2, precision = 0.001, step = 2.;
     double best = 1., fx = -1., fy = -1., bfx = 0., bfy = 0.;
-    double e, e0;
+    double e, e0, e1;
     pipi_image_t *img, *gauss, *dither, *tmp;
     int dx, dy;
@@ -60,4 +60,9 @@
     tmp = pipi_gaussian_blur(dither, sigma);
     e0 = pipi_measure_rmsd(gauss, tmp);
+    pipi_free(tmp);
+
+    /* Compute the fast error */
+    tmp = pipi_gaussian_blur_ext(dither, sigma, sigma, 0.16, 0.26);
+    e1 = pipi_measure_rmsd(gauss, tmp);
     pipi_free(tmp);
 
@@ -86,5 +91,5 @@
     }
 
-    printf("E: %g E_min: %g dx: %g dy: %g\n", e0, best, fx, fy);
+    printf("E: %g E_fast: %g E_min: %g dx: %g dy: %g\n", e0, e1, best, fx, fy);
 
     pipi_free(dither);