Updated.

[graphlib.git] / mandel / mandel.cpp

diff --git a/mandel/mandel.cpp b/mandel/mandel.cpp
index e57e2c8..1afd640 100644 (file)
--- a/mandel/mandel.cpp
+++ b/mandel/mandel.cpp
@@ -14,45 +14,45 @@ void do_mandel(DrawingWindow &w,
     double rouge, vert, bleu;
     int i;
 
-    pr = (Rmax - Rmin) / w.width;
-    pi = (Imax - Imin) / w.height;
+    pr = (Rmax - Rmin) / (w.width - 1);
+    pi = (Imax - Imin) / (w.height - 1);
     
     cr = Rmin;
     for (x = 0; x < w.width; x++) {
         ci = Imax;
         for (y = 0; y < w.height; y++) {
-            /* z_1 = c */
+            // z_1 = c
             zr = cr;
             zi = ci;
             for (i = 1; i <= maxiter; i++) {
                 zr2 = zr * zr;
                 zi2 = zi * zi;
                 if (zr2 + zi2 >= 4) {
-                    /* |z| >= 2 : on sort de la boucle */
+                    // |z| >= 2 : on sort de la boucle
                     break;
                 }
-                /* on calcule le z suivant */
-                zi = 2*zr*zi + ci;
+                // on calcule le z suivant
+                zi = 2 * zr * zi + ci;
                 zr = zr2 - zi2 + cr;
             }
             if (i > maxiter) {
                 rouge = vert = bleu = 0.0;
             } else {
-                /* on est sorti trop tôt du for(...): on affiche le pixel
-                   d'un couleur en fonction de i */
+                // on est sorti trop tôt du for(...): on affiche le
+                // pixel d'un couleur en fonction de i
                 int ii = (maxiter - i) % 96;
                 if (ii < 32) {
-                    /* vert -> bleu */
+                    // vert -> bleu
                     bleu = ii / 32.0;
                     vert = 1.0 - bleu;
                     rouge = 0.0;
                 } else if (ii < 64) {
-                    /* bleu -> rouge */
+                    // bleu -> rouge
                     rouge = (ii - 32) / 32.0;
                     bleu = 1.0 - rouge;
                     vert = 0.0;
                 } else {
-                    /* rouge -> vert */
+                    // rouge -> vert
                     vert = (ii - 64) / 32.0;
                     rouge = 1.0 - vert;
                     bleu = 0.0;
@@ -69,51 +69,48 @@ void do_mandel(DrawingWindow &w,
 
 void mandel(DrawingWindow &w)
 {
-    /* nombre max d'itérations */
+    // nombre max d'itérations
     int maxiter = 500;
-    /* zone d'intérêt par défaut */
+    // zone d'intérêt par défaut
     double Rmin = -2.05;
     double Rmax = 0.55;
     double Imin = -1.3;
     double Imax = 1.3;
-    /* cible du zoom */
+    // cible du zoom
     double Tr;
     double Ti;
 
     while (1) {
-        std::cout << "[ " << Rmin
-                  << (Imin >= 0.0 ? " + " : " - ")
-                  << (Imin >= 0.0 ? Imin : -Imin)
-                  << " i ; " << Rmax
-                  << (Imax >= 0.0 ? " + " : " - ")
-                  << (Imax >= 0.0 ? Imax : -Imax)
-                  << " i ] (" << maxiter << ")"
-                  << std::endl;
-
         do_mandel(w, Rmin, Rmax, Imin, Imax, maxiter);
 
+        w.setColor("white");
+        w.drawText(5, 5, "Cliquer sur l'image pour zoomer");
+
         int x, y;
         int button;
-        std::cout << "-=[ Cliquer sur l'image pour zoomer ]=-" << std::endl;
         w.waitMousePress(x, y, button);
-        Tr = Rmin + x * (Rmax - Rmin) / w.width;
-        Ti = Imax - y * (Imax - Imin) / w.height;
-        std::cout << "(" << x << ";" << y << ") -> ("        
-                  << Tr << ";" << Ti << ")" << std::endl;
 
-        double Rmin2 = Rmin / 2;
-        double Rmax2 = Rmax / 2;
-        double Imin2 = Imin / 2;
-        double Imax2 = Imax / 2;
-        double Rdecal = Tr - (Rmin2 + Rmax2) / 2;
-        double Idecal = Ti - (Imin2 + Imax2) / 2;
-        Rmin2 += Rdecal;
-        Rmax2 += Rdecal;
-        Imin2 += Idecal;
-        Imax2 += Idecal;
+        // calcul des coordonnées du point cliqué
+        Tr = Rmin + x * (Rmax - Rmin) / (w.width - 1);
+        Ti = Imax - y * (Imax - Imin) / (w.height - 1);
+
+        // calcul de la nouvelle zone d'intérêt :
+        // zoom ×2 en direction du point cliqué
+        const int zoom = 2;
+        double Rmin2 = Rmin / zoom;
+        double Rmax2 = Rmax / zoom;
+        double Imin2 = Imin / zoom;
+        double Imax2 = Imax / zoom;
+        double Rshift = Tr - (Rmin2 + Rmax2) / 2;
+        double Ishift = Ti - (Imin2 + Imax2) / 2;
+        Rmin2 += Rshift;
+        Rmax2 += Rshift;
+        Imin2 += Ishift;
+        Imax2 += Ishift;
 
-        double x_factor = w.width / (Rmax - Rmin);
-        double y_factor = w.height / (Imax - Imin);
+        // affichage d'un rectangle autour de la nouvelle zone d'intérêt
+        double x_factor = (w.width - 1) / (Rmax - Rmin);
+        double y_factor = (w.height - 1) / (Imax - Imin);
         w.setColor("white");
         w.drawRect((Rmin2 - Rmin) * x_factor, (Imax - Imin2) * y_factor,
                    (Rmax2 - Rmin) * x_factor, (Imax - Imax2) * y_factor);