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

SCM Repository

[diderot] Annotation of /benchmarks/programs/lic2d/hact2.diderot
ViewVC logotype

Annotation of /benchmarks/programs/lic2d/hact2.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1615 - (view) (download)

1 : glk 1615 // lic-turb.diderot
2 :     //
3 :     // demo of line integral convolution (LIC) on 2D turbulent flow
4 :     //
5 :     // process output with:
6 :     // unu reshape -i lic-turb2d.txt -s 3 1020 561 | unu quantize -b 8 -o lic-turb2d.ppm
7 :    
8 :     int imgSizeX = 500;
9 :     int imgSizeY = 500;
10 :     real h = 0.005; // step size of integration
11 :     int stepNum = 10; // take this many steps both upstream and downstream
12 :    
13 :     // field#2(2)[2] V = load("imp-u.nrrd") ⊛ bspln3;
14 :     field#1(2)[2] V = load("imp-u.nrrd") ⊛ ctmr;
15 :     field#0(2)[] R = load("../../data/R.nrrd") ⊛ tent;
16 :    
17 :     strand LIC (int xi, int yi) {
18 :     real xx = lerp(0.0, 8.0, 0.0, real(xi), real(imgSizeX)-1.0);
19 :     real yy = lerp(0.0, 8.0, 0.0, real(yi), real(imgSizeY)-1.0);
20 :     vec2 pos0 = [xx,yy];
21 :     output vec3 RGB = [0.0, 0.0, 0.0];
22 :    
23 :     update {
24 :     tensor[2,2] J = ∇⊗V(pos0);
25 :     RGB = [J[0,0], J[0,1], J[0,0]];
26 :     stabilize;
27 :     }
28 :     }
29 :    
30 :     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