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 1685 - (view) (download)

1 : jhr 1671 // curl-turb2d.diderot
2 :     //
3 :     //
4 :     // process output with:
5 :     // unu reshape -i curl-turb2d.txt -s 1020 561 | unu quantize -b 8 -o curl-turb2d.png
6 :    
7 :    
8 :     int imgSizeX = 1020;
9 :     int imgSizeY = 561;
10 :    
11 :     field#1(2)[2] V = load("../data/turb2d.nrrd") ⊛ ctmr;
12 :    
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 :     output tensor[2,2] M = ∇⊗V(pos0);
19 :    
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