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

SCM Repository

[diderot] Annotation of /branches/charisee/test/EIN/curl3/x5.diderot
ViewVC logotype

Annotation of /branches/charisee/test/EIN/curl3/x5.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 2857 - (view) (download)

1 : cchiw 2856 field#2(3)[] s3 = load("../../../data/vfrhand-nohip.nhdr")⊛ bspln3;
2 :     field#2(3)[] t3 = bspln3 ⊛ load("../../../data/txs.nrrd");
3 :     field#2(3)[3] V3= load("../../../data/const-curl3d.nrrd") ⊛ bspln3;
4 :     field#2(3)[3] G3= load("../../../data/const-curl3d.nrrd") ⊛ bspln5;
5 :     int imgSizeX = 16;
6 :     int imgSizeY = 16;
7 :     int imgSizeZ = 16;
8 :    
9 :     strand tensorTest (int xi,int yi,int zi)
10 :     {
11 :    
12 :     real xx = lerp(0.3, 6.75, 0.0, real(xi), real(imgSizeX)-1.0);
13 :     real yy = lerp(0.6, 3.69, 0.0, real(yi), real(imgSizeY)-1.0);
14 :     real zz = lerp(2.00, 17.0, 0.0, real(zi), real(imgSizeZ)-1.0);
15 :     vec3 posv3 = [xx,yy,zz];
16 :     output tensor[3] out=[99,99,99];
17 :    
18 :     update
19 :     {
20 :     if (inside(posv3,V3))
21 : cchiw 2857 { if (inside(posv3,G3))
22 : cchiw 2856 {
23 :     out=∇×(G3×V3)(posv3);//lhs
24 : cchiw 2857 }
25 : cchiw 2856 }
26 :    
27 :     stabilize;
28 :     }
29 :     }
30 :    
31 :     initially [ tensorTest(xi, yi, zi) | zi in 0..(imgSizeZ-1), yi in 0..(imgSizeY-1), xi in 0..(imgSizeX-1) ];

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