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

SCM Repository

[diderot] Annotation of /tests/vis15-bugs/pbfs1.diderot
ViewVC logotype

Annotation of /tests/vis15-bugs/pbfs1.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 5149 - (view) (download)

1 : jhr 5149 #version 1
2 : glk 5146
3 : jhr 5149 image(3)[] img = image("pbfs-vol.nrrd");
4 :    
5 : glk 5146 field#2(3)[] F = bspln3 ⊛ clamp(img);
6 :    
7 : jhr 5149 vec3 cent = [0,0,0];
8 :     int sz = 10;
9 :     real width = 2;
10 : glk 5146
11 :     function bool featureLost(vec3 x) {
12 :     if (!inside(x, F)) { return true; }
13 : jhr 5149 return F(x) < 10;
14 : glk 5146 }
15 :    
16 : jhr 5149 strand point (vec3 pos0)
17 :     {
18 :     output real out = -1.0;
19 :     vec3 pos = pos0;
20 : glk 5146
21 :     update {
22 :     // can't proceed if outside field or have zero gradient
23 :     if (featureLost(pos)) { die; }
24 : jhr 5149 out = F (pos);
25 :     stabilize;
26 : glk 5146 }
27 :    
28 :     }
29 :    
30 :     initially { point([cent[0] + width*lerp(-1, 1, 0, ii, sz-1)/2,
31 :     cent[1] + width*lerp(-1, 1, 0, jj, sz-1)/2,
32 : jhr 5149 cent[2] + width*lerp(-1, 1, 0, kk, sz-1)/2])
33 : glk 5146 | ii in 0 .. sz-1,
34 :     jj in 0 .. sz-1,
35 :     kk in 0 .. sz-1 };

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