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

SCM Repository

[diderot] Annotation of /trunk/test/probe/probe.cl
ViewVC logotype

Annotation of /trunk/test/probe/probe.cl

Parent Directory Parent Directory | Revision Log Revision Log


Revision 259 - (view) (download)

1 : lamonts 259 __kernel void probe (float * img,
2 :     float * h1,
3 :     float * h2,
4 :     float * posValues,
5 :     float * out,
6 :     float16 transformMatrix,
7 :     int ySize,
8 :     int zSize,
9 :     int xSize)
10 :     {
11 :     int location = get_global_id(0);
12 :     float4 pos = (float4)(posValues[location *3],
13 :     posValues[(location *3) + 1],
14 :     posValues[(location *3) + 2],1.0f);
15 :    
16 :     float probedVal;
17 :     float4 t_k, f, imgPt,nf,v, hx, hy, hz, h_tx, h_ty, hMult;
18 :     int4 n1;
19 :    
20 :     float4 d = (float4) (h2[0],h1[0],h1[0],h2[0]);
21 :     float4 c = (float4) (h2[1],h1[1],h1[1],h2[1]);
22 :     float4 b = (float4) (h2[2],h1[2],h1[2],h2[2]);
23 :     float4 a = (float4) (h2[3],h1[3],h1[3],h2[3]);
24 :    
25 :     // Transform the value to image space position.
26 :     imgPt = (float4) (dot(pos,transformMatrix.s0123),
27 :     dot(pos,transformMatrix.s4567),
28 :     dot(pos,transformMatrix.s89ab),
29 :     dot(pos,transformMatrix.scdef));
30 :    
31 :    
32 :     nf = floor(imgPt);
33 :     f = imgPt - nf;
34 :     n1 = convert_int4(nf);
35 :    
36 :     /* //Debugging
37 :     if(location == 1)
38 :     {
39 :     printf("\n============\nImage Point: (%f, %f, %f)\n",imgPt.x,imgPt.y,imgPt.z);
40 :     printf("N: (%d,%d,%d,%d)\n",n1.x,n1.y,n1.z,n1.w);
41 :     printf("f: (%f,%f,%f,%f)\n",f.x,f.y,f.z,f.w);
42 :     }
43 :     */
44 :    
45 :     // the t value for h(fx - i)
46 :     float t_i[4] = {-1.0 - f.x, -f.x, f.x - 1.0, f.x - 2.0};
47 :    
48 :     // the t value for h(fy - j)
49 :     float t_j[4] = {-1.0 - f.y, -f.y, f.y - 1.0, f.y - 2.0};
50 :    
51 :     t_k = (float4) (-1.0 - f.z, -f.z, f.z - 1.0, f.z - 2.0);
52 :    
53 :     float tx[4], ty[4];
54 :    
55 :     for (int i = -1; i <= 2; i++) {
56 :     for (int j = -1; j <= 2; j++) {
57 :     // compute z axis using vectors
58 :     int index = (n1.x+i) * ySize * zSize + (n1.y+j) * zSize + n1.z-1;
59 :     v = (float4)(
60 :     img[index],
61 :     img[index+1],
62 :     img[index+2],
63 :     img[index+3]);
64 :    
65 :    
66 :     h_tx = (float4) (t_i[i+1],t_i[i+1],t_i[i+1],t_i[i+1]);
67 :    
68 :     hx = (d + h_tx * (c + h_tx * (b + h_tx * a)));
69 :    
70 :     h_ty = (float4) (t_j[j+1],t_j[j+1],t_j[j+1],t_j[j+1]);
71 :    
72 :     hy = (d + h_ty * (c + h_ty * (b + h_ty * a)));
73 :    
74 :     hz = (d + t_k * (c + t_k * (b + t_k * a)));
75 :    
76 :     ty[j+1] = dot(v, (hx * hy * hz));
77 :    
78 :    
79 :     /* //Debugging
80 :     if(location==1)
81 :     {
82 :     printf("========J:%d======\nImage Data:(%f,%f,%f,%f)\n",j,v.x,v.y,v.z,v.w);
83 :     printf("hx:(%f,%f,%f,%f)\nhx:(%f,%f,%f,%f)\nhz:(%f,%f,%f,%f)\n",
84 :     hx.x,hx.y,hx.z,hx.w,
85 :     hy.x,hy.y,hy.z,hy.w,
86 :     hz.x,hz.y,hz.z,hz.w);
87 :     printf("hMult:(%f,%f,%f,%f)\n",hMult.x,hMult.y,hMult.z,hMult.w);
88 :     printf("ty[%d]: %f\n",j+1,ty[j+1]);
89 :    
90 :     }
91 :     */
92 :    
93 :     }
94 :     //tx[i+1] = dot((float4)(ty[0], ty[1], ty[2], ty[3]), d + t_j * (c + t_j * (b + t_j * a)));
95 :     tx[i+1] = ty[0] + ty[1] + ty[2] + ty[3];
96 :     }
97 :     //probedVal = dot((float4)(tx[0],tx[1],tx[2],tx[3]), d + t_i * (c + t_i * (b + t_i * a)));
98 :    
99 :     probedVal = tx[0] + tx[1] + tx[2] + tx[3];
100 :    
101 :     out[location] = probedVal;
102 :     }

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