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

SCM Repository

[diderot] Annotation of /branches/vis12/test/curl-turb2d.diderot
ViewVC logotype

Annotation of /branches/vis12/test/curl-turb2d.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2164 - (view) (download)

1 : jhr 1671 // curl-turb2d.diderot
2 :     //
3 :     //
4 :     // process output with:
5 : jhr 2164 // unu quantize -i crl.nrrd -b 8 -o curl-turb2d.png
6 : jhr 1671
7 :    
8 :     int imgSizeX = 1020;
9 :     int imgSizeY = 561;
10 :    
11 : jhr 1927 field#1(2)[2] V = image("../data/turb2d.nrrd") ⊛ ctmr;
12 : jhr 1671
13 :     strand LIC (int xi, int yi) {
14 :     real xx = lerp(0.3, 6.75, 0.0, real(xi), real(imgSizeX)-1.0);
15 :     real yy = lerp(0.6, 3.69, 0.0, real(yi), real(imgSizeY)-1.0);
16 :     vec2 pos0 = [xx,yy];
17 :    
18 : glk 2163 output real crl = ∇×V(pos0);
19 : jhr 1671
20 :     update {
21 :    
22 :     stabilize;
23 :     }
24 :     }
25 :    
26 :     // strand LIC (int xi, int yi) {
27 :     // real xx = lerp(0.3, 6.75, 0.0, real(xi), real(imgSizeX)-1.0);
28 :     // real yy = lerp(0.6, 3.69, 0.0, real(yi), real(imgSizeY)-1.0);
29 :     // vec2 pos0 = [xx,yy];
30 :     //
31 :     // tensor[2,2] M = ∇V(pos0);
32 :     // output real curl = 0.0;
33 :     //
34 :     // update {
35 :     // // I'm not sure I remember the correct way to index into a tensor
36 :     // //curl = M[0,1] - M[1,0];
37 :     //
38 :     // // This should work, even if the above indexing was wrong
39 :     // vec2 v = [0.0, 1.0];
40 :     // vec2 u = [1.0, 0.0];
41 :     // curl = v•M•u - u•M•v;
42 :     //
43 :     // stabilize;
44 :     // }
45 :     // }
46 :    
47 :     initially [ LIC(xi, yi) | yi in 0..(imgSizeY-1), xi in 0..(imgSizeX-1) ];

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