Fixed >scaled in example.

[forth.jl.git] / examples / mandelbrot.4th

diff --git a/examples/mandelbrot.4th b/examples/mandelbrot.4th
index 2767738..55f746e 100644 (file)
--- a/examples/mandelbrot.4th
+++ b/examples/mandelbrot.4th
@@ -1,7 +1,23 @@
-( --- Complex integer arithmetic --- )
+( --- Complex arithmetic --- )
+
+( Location of floating point. )
+: precision 10000 ;
+
+: sign dup abs / ;
+
+: >scaled
+    precision 10 / * over
+    ?dup 0<> if 
+        sign *
+    then
+    swap precision * +
+;
+
+( Redefine multiplication.  Yay forth! )
+: * precision */ ;
 
 : c* ( x1 y1 x2 y2 -- x3 y3 )
-        swap rot                ( x1 x2 y1 y2 )
+        swap -rot               ( x1 x2 y1 y2 )
         2dup * negate           ( x1 x2 y1 y2 -y1y2 )
         4 pick 4 pick * +       ( x1 x2 y1 y2 (x1x2-y1y2))
         4 roll 2 roll *         ( x2 y1 (x1x2-y1y2) x1y2 )
@@ -16,7 +32,7 @@
 
 : csq 2dup c* ;
 
-; cmagsq ( x1 y1 -- mag )
+: cmagsq ( x1 y1 -- mag )
         csq abs
 ;
 
@@ -28,15 +44,33 @@
 
 : inSet? ( cr ci -- res )
     
-    100 0 DO
+    100 0 do
 
-        2SWAP 2DUP 5 ROLL 5 ROLL
+        2swap 2dup 5 roll 5 roll
         iterate
-        2DUP cmagsq
-        100 > IF
-            LEAVE
-        THEN
+        2dup cmagsq
+        100 > if
+            leave
+        then
+
+    loop
+;
+
+: xsteps 100 ;
+: ysteps 50 ;
+
+( Draw the Mandelbrot Set!)
+: mandel ( x1 y1 x2 y2 -- )
 
-    LOOP
+    1 pick 4 pick -
+    2 pick 5 pick do
+        i 0 inSet? if
+            42 emit
+        else
+            '.' emit
+        then
+    dup +loop
 ;
 
+( Clean up - hide non-standard multiplication def. )
+hide *