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

Parent Directory | Revision Log

---

Revision 262 - (view) (download)

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

---

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