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

SCM Repository

[diderot] Annotation of /trunk/test/bug-matrix.diderot
ViewVC logotype

Annotation of /trunk/test/bug-matrix.diderot

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1115 - (view) (download)

1 : jhr 1115
2 :     /*
3 :    
4 :     If 3x3 matrix subtraction works, mip.txt will be:
5 :    
6 :     0.928571 -0.142857 0.214286
7 :     -0.142857 0.714286 0.428571
8 :     0.214286 0.428571 0.357143
9 :     0.928571 -0.142857 0.214286
10 :     -0.142857 0.714286 0.428571
11 :     0.214286 0.428571 0.357143
12 :    
13 :     That is, the same three lines, repeated twice.
14 :    
15 :     */
16 :    
17 :     strand subtract (int bugvali) {
18 :     output vec3 out = [0.0, 0.0, 0.0];
19 :     int vali = bugvali;
20 :    
21 :     update {
22 :     vec3 norm = normalize([1.0, 2.0, -3.0]);
23 :     tensor[3,3] sub0 = identity[3] + (-1.0*norm⊗norm);
24 :     tensor[3,3] sub1 = identity[3] - norm⊗norm;
25 :    
26 :     if (0 == vali) {
27 :     out = sub0•[1.0, 0.0, 0.0];
28 :     } else if (1 == vali) {
29 :     out = sub0•[0.0, 1.0, 0.0];
30 :     } else if (2 == vali) {
31 :     out = sub0•[0.0, 0.0, 1.0];
32 :     }
33 :     else if (3 == vali) {
34 :     out = sub1•[1.0, 0.0, 0.0];
35 :     } else if (4 == vali) {
36 :     out = sub1•[0.0, 1.0, 0.0];
37 :     } else if (5 == vali) {
38 :     out = sub1•[0.0, 0.0, 1.0];
39 :     }
40 :    
41 :     stabilize;
42 :     }
43 :     }
44 :    
45 :     initially [ subtract(vali) | vali in 0..5 ];

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