summaryrefslogtreecommitdiff
path: root/Master/texmf-dist/doc/latex/lisp-on-tex/tug2013/bench/luatex/mandel.tex
blob: 2efa5d0ec7e29798f0775286e9f686b46ee7dfcb (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
\documentclass{article}
\usepackage[pdftex,a3paper,margin=1pt,landscape]{geometry}
%% draw macros
\newcount\mlength
\mlength0
\newcount\cstate
\newdimen\mandelunit
\mandelunit=0.5pt


\def\wp{%
  \ifnum\cstate=1 \global\advance\mlength1
  \else\ifnum\cstate=-1 \vrule width \mlength\mandelunit height \mandelunit depth 0pt \global\mlength1 \fi\fi
  \global\cstate1}
\def\bp{%
  \ifnum\cstate=-1 \global\advance\mlength1
  \else\ifnum\cstate=1 \hspace*{\mlength\mandelunit}\global\mlength1 \fi\fi
  \global\cstate-1}
\def\r{%
  \ifnum\cstate=1\hspace*{\mlength\mandelunit}\else\vrule width \mlength\mandelunit height \mandelunit depth 0pt \fi
  \global\cstate0 \global\mlength0}


\begin{document}
  \thispagestyle{empty}\noindent
  \leavevmode\baselineskip=\mandelunit
  \directlua{
    maxloop=20
    scale=0.002
    function isMandel(a,b)
      x=0
      y=0
      xx=0
      yy=0
      for k=0,maxloop do
        xx=x*x
        yy=y*y
        if (xx + yy) > 4.0 then 
          return false
        else
          y=2*x*y+b
          x=xx-yy+a
        end 
      end
      return true
    end
    for b=1.0,-1.0,-scale do
      for a=-1.5,0.5,scale do
        if isMandel(a,b) then
          tex.print("\noexpand\\bp")
        else
          tex.print("\noexpand\\wp")
        end
      end
      tex.print("\noexpand\\r\noexpand\\\noexpand\\")
    end
  }
\end{document}