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

SCM Repository

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

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

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2143 - (view) (download)
Original Path: branches/vis12/test/mandel.diderot

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 : jhr 1971 input real escape = 4;
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 : jhr 2143 function real coord (real x, int idx) = lerp(x - fov, x + fov, 1, idx, reso);
21 :    
22 : glk 1893 // Strand definition
23 : jhr 1971 strand mandel (vec2 c) {
24 :     vec2 z = [0, 0];
25 : glk 1893 int iter = 0;
26 : jhr 1971 output vec3 rgb = [0, 0, 0];
27 : glk 1893 update {
28 :     z = [z[0]^2 - z[1]^2, 2.0*z[0]*z[1]] + c; // z = z^2 + c
29 :     iter += 1;
30 : glk 1948 if (|z| > escape) {
31 :     // point has escaped; set continuous color from iter and radius
32 : glk 1972 real t = 4*log(iter - log2(log(|z|)/log(escape)));
33 : glk 1948 rgb = [sin(-t), cos(t), sin(t)];
34 :     rgb = ([1.8, 1.8, 1.8] + rgb)/2.8;
35 : glk 1893 stabilize;
36 :     }
37 :     if (iter > maxiter) {
38 : jhr 1971 rgb = [0, 0, 0];
39 : glk 1893 stabilize;
40 :     }
41 :     }
42 :     }
43 :    
44 :     // Strand initialization
45 : jhr 2143 initially [ mandel([coord(center[0], crIdx), coord(center[1], ciIdx)])
46 : glk 1897 | ciIdx in 1..reso, crIdx in 1..reso ];

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