Home My Page Projects Code Snippets Project Openings diderot
Summary Activity Tracker Tasks SCM

SCM Repository

[diderot] Annotation of /branches/vis12/test/mandel.diderot
ViewVC logotype

Annotation of /branches/vis12/test/mandel.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1970 - (view) (download)

1 : jhr 1938 /*! \file mandel.diderot
2 :     *
3 :     * \author Gordon Kindlmann
4 :     */
5 : glk 1904
6 : jhr 1938 /*
7 :     * COPYRIGHT (c) 2012 The Diderot Project (http://diderot-language.cs.uchicago.edu)
8 :     * All rights reserved.
9 :     */
10 : glk 1904
11 : glk 1894 // ../bin/diderotc --exec mandel.diderot && ./mandel && unu quantize -b 8 -i rgb.nrrd -o mandel.png
12 : glk 1893
13 :     // Global definitions
14 : glk 1904 input int reso = 2000;
15 : glk 1948 input real escape = 4.0;
16 : glk 1904 input int maxiter = 2000;
17 :     input vec2 center = [-1.2, -0.3]; // [0.0,0.0] for whole thing
18 : glk 1948 input real fov = 0.16; // 2.0 for whole thing
19 : glk 1893
20 :     // Strand definition
21 : glk 1897 strand mandel(vec2 c) {
22 : glk 1893 vec2 z = [0.0, 0.0];
23 :     int iter = 0;
24 :     output vec3 rgb = [0.0, 0.0, 0.0];
25 :     update {
26 :     z = [z[0]^2 - z[1]^2, 2.0*z[0]*z[1]] + c; // z = z^2 + c
27 :     iter += 1;
28 : glk 1948 if (|z| > escape) {
29 :     // point has escaped; set continuous color from iter and radius
30 : jhr 1970 real t = 4.0*log(real(iter) - log2(log(|z|)/log(escape)));
31 : glk 1948 rgb = [sin(-t), cos(t), sin(t)];
32 :     rgb = ([1.8, 1.8, 1.8] + rgb)/2.8;
33 : glk 1893 stabilize;
34 :     }
35 :     if (iter > maxiter) {
36 :     rgb = [0.0, 0.0, 0.0];
37 :     stabilize;
38 :     }
39 :     }
40 :     }
41 :    
42 :     // Strand initialization
43 : glk 1904 initially [ mandel([lerp(center[0]-fov, center[0]+fov,
44 :     1.0, real(crIdx), real(reso)),
45 :     lerp(center[1]-fov, center[1]+fov,
46 :     1.0, real(ciIdx), real(reso))])
47 : glk 1897 | ciIdx in 1..reso, crIdx in 1..reso ];

root@smlnj-gforge.cs.uchicago.edu
ViewVC Help
Powered by ViewVC 1.0.0