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

SCM Repository

[diderot] Annotation of /tests/DATm-bugs/tests/C8/observ.diderot
ViewVC logotype

Annotation of /tests/DATm-bugs/tests/C8/observ.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 5141 - (view) (download)

1 : cchiw 5141
2 :     tensor [2, 2]F0 = [[-1, -4], [3, -5]];
3 :     field#4(3)[]F1 = c4hexic⊛ image("rst/data/inputfile1.nrrd");
4 :     input real scale = 1;
5 :     field#4(3)[2, 2]G = inv((normalize(F0))*(F1));
6 :     int length =7;
7 :     strand f(int i) {
8 :    
9 :     output tensor [2, 2] out = [ [0.0,0.0], [0.0,0.0]];
10 :     update {
11 :    
12 :     tensor [3] pos=[0,0,0];
13 :     if(i==0){
14 :     pos = ([-0.04, -0.36, 0.03]);
15 :     }
16 :     else if(i==1){
17 :     pos = ([-0.04, -0.36, 0.03]);
18 :     }
19 :     else if(i==2){
20 :     pos = ([-0.2, 0.32, -0.35]);
21 :     }
22 :     else if(i==3){
23 :     pos = ([0.06, -0.3, -0.16]);
24 :     }
25 :     else if(i==4){
26 :     pos = ([0.0, 0.06, -0.39]);
27 :     }
28 :     else if(i==5){
29 :     pos = ([-0.39, -0.25, -0.32]);
30 :     }
31 :     else if(i==6){
32 :     pos = ([0.3, -0.24, -0.12]);
33 :     }
34 :     else if(i==7){
35 :     pos = ([-0.1, -0.17, 0.21]);
36 :     }
37 :    
38 :     if((inside(pos,F1)) && true){
39 :     if(((det(((normalize(F0))*(F1))(pos))>0.01) || (det(((normalize(F0))*(F1))(pos))< -0.01)))
40 :     { out = (G)(pos);}
41 :     else{out = [ [7.2,7.2], [7.2,7.2]];}
42 :     }
43 :     else{out = [ [7.2,7.2], [7.2,7.2]];}
44 :     stabilize;
45 :     }
46 :     }
47 :     initially [ f(i) | i in 0..length];

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