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

SCM Repository

[diderot] Annotation of /branches/vis12-cl/src/include/Diderot/cl-inline-image.h
ViewVC logotype

Annotation of /branches/vis12-cl/src/include/Diderot/cl-inline-image.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3349 - (view) (download) (as text)

1 : jhr 1671 /*! \file cl-inline-image.h
2 :     *
3 :     * \author Lamont Samuels
4 :     */
5 :    
6 :     /*
7 : jhr 3349 * This code is part of the Diderot Project (http://diderot-language.cs.uchicago.edu)
8 :     *
9 :     * COPYRIGHT (c) 2015 The University of Chicago
10 : jhr 1671 * All rights reserved.
11 :     */
12 :    
13 :     #ifndef _DIDEROT_CL_INLINE_IMAGE_H_
14 :     #define _DIDEROT_CL_INLINE_IMAGE_H_
15 :    
16 :     #ifndef _DIDEROT_CL_TYPES_H_
17 :     # include "cl-types.h"
18 :     #endif
19 :     #ifndef _DIDEROT_CL_INLINE_MATRIX_H_
20 :     # include "cl-inline-matrix.h"
21 :     #endif
22 :    
23 :    
24 :     /***** 1D Image functions *****/
25 :    
26 :     // check if pos is inside the img, assuming that we have a border of width s.
27 :     //
28 : jhr 3118 inline bool Inside1 (Diderot_real_t pos, __global Diderot_image1D_t *img, int s)
29 : jhr 1671 {
30 :     return ((s-1 < pos) && (pos < (img->size[0] - s)));
31 :     }
32 :    
33 : jhr 3118 inline Diderot_real_t ToImageSpace1 ( __global Diderot_image1D_t *img, Diderot_real_t pos)
34 : jhr 1671 {
35 :     return img->s * pos + img->t;
36 :     }
37 :    
38 :     /***** 2D Image functions *****/
39 :    
40 :     // check if pos is inside the img, assuming that we have a border of width s.
41 :     //
42 : jhr 3118 inline bool Inside2 (Diderot_vec2_t posArg, __global Diderot_image2D_t *img, int s)
43 : jhr 1671 {
44 :     return ((s-1 < posArg.s0) && (posArg.s0 < (img->size[0] - s))
45 :     && (s-1 < posArg.s1) && (posArg.s1 < (img->size[1] - s)));
46 :     }
47 :    
48 : jhr 3118 inline Diderot_vec2_t ToImageSpace2 (__global Diderot_image2D_t *img, Diderot_vec2_t p)
49 : jhr 1671 {
50 : jhr 3118 return (Diderot_vec2_t)(dot(img->w2i[0], p), dot(img->w2i[1], p)) + img->tVec;
51 : jhr 1671 }
52 :    
53 : jhr 3118 inline Diderot_vec2_t ToWorldSpace2 (__global Diderot_image2D_t *img, Diderot_vec2_t g)
54 : jhr 1671 {
55 : jhr 3118 return (Diderot_vec2_t)(
56 : jhr 1671 dot(img->w2iT[0], g),
57 :     dot(img->w2iT[1], g));
58 :     }
59 :    
60 : jhr 3118 inline void ToWorldSpace2x2 ( Diderot_Mat2x2_t dst, __global Diderot_image2D_t *img, Diderot_Mat2x2_t src)
61 : jhr 1671 {
62 :     /* Diderot_Mat2x2_t tmp;
63 :     mulMat2x2Mat2x2f (tmp, img->w2iT, src);
64 :     mulMat2x2Mat2x2f (dst, tmp, img->w2i); */
65 :    
66 :     Diderot_Mat2x2_t tmp;
67 :    
68 : jhr 3118 tmp[0] = (Diderot_vec2_t)(
69 :     dot(img->w2iT[0], (Diderot_vec2_t)(src[0].s0, src[1].s0)),
70 :     dot(img->w2iT[0], (Diderot_vec2_t)(src[0].s1, src[1].s1)));
71 :     tmp[1] = (Diderot_vec2_t)(
72 :     dot(img->w2iT[1], (Diderot_vec2_t)(src[0].s0, src[1].s0)),
73 :     dot(img->w2iT[1], (Diderot_vec2_t)(src[0].s1, src[1].s1)));
74 : jhr 1671
75 : jhr 3118 dst[0] = (Diderot_vec2_t)(
76 :     dot(tmp[0], (Diderot_vec2_t)(img->w2i[0].s0, img->w2i[1].s0)),
77 :     dot(tmp[0], (Diderot_vec2_t)(img->w2i[0].s1, img->w2i[1].s1)));
78 :     dst[1] = (Diderot_vec2_t)(
79 :     dot(tmp[1], (Diderot_vec2_t)(img->w2i[0].s0, img->w2i[1].s0)),
80 :     dot(tmp[1], (Diderot_vec2_t)(img->w2i[0].s1, img->w2i[1].s1)));
81 : jhr 1671
82 :     }
83 :    
84 :     /***** 3D Image functions *****/
85 :    
86 :     // check if pos is inside the img, assuming that we have a border of width s.
87 :     //
88 : jhr 3134 inline bool Inside3 (Diderot_vec3_t posArg, __global Diderot_image3D_t *img, int s)
89 : jhr 1671 {
90 :     return ((s-1 < posArg.s0) && (posArg.s0 < (img->size[0] - s))
91 :     && (s-1 < posArg.s1) && (posArg.s1 < (img->size[1] - s))
92 :     && (s-1 < posArg.s2) && (posArg.s2 < (img->size[2] - s)));
93 :     }
94 :    
95 : jhr 3134 inline Diderot_vec3_t ToImageSpace3 (__global Diderot_image3D_t *img, Diderot_vec3_t p)
96 : jhr 1671 {
97 :     return VEC3(dot(img->w2i[0], p), dot(img->w2i[1], p), dot(img->w2i[2], p)) + img->tVec;
98 :     }
99 :    
100 : jhr 3134 inline Diderot_vec3_t ToWorldSpace3 (__global Diderot_image3D_t *img, Diderot_vec3_t g)
101 : jhr 1671 {
102 :     return VEC3(
103 :     dot(img->w2iT[0], g),
104 :     dot(img->w2iT[1], g),
105 :     dot(img->w2iT[2], g));
106 :     }
107 :    
108 : jhr 3134 inline void ToWorldSpace3x3 ( Diderot_Mat3x3_t dst, __global Diderot_image3D_t *img, Diderot_Mat3x3_t src)
109 : jhr 1671 {
110 :     Diderot_Mat3x3_t tmp;
111 :    
112 :     tmp[0] = VEC3(
113 :     dot(img->w2iT[0], VEC3(src[0].s0, src[1].s0, src[2].s0)),
114 :     dot(img->w2iT[0], VEC3(src[0].s1, src[1].s1, src[2].s1)),
115 :     dot(img->w2iT[0], VEC3(src[0].s2, src[1].s2, src[2].s2)));
116 :     tmp[1] = VEC3(
117 :     dot(img->w2iT[1], VEC3(src[0].s0, src[1].s0, src[2].s0)),
118 :     dot(img->w2iT[1], VEC3(src[0].s1, src[1].s1, src[2].s1)),
119 :     dot(img->w2iT[1], VEC3(src[0].s2, src[1].s2, src[2].s2)));
120 :     tmp[2] = VEC3(
121 :     dot(img->w2iT[2], VEC3(src[0].s0, src[1].s0, src[2].s0)),
122 :     dot(img->w2iT[2], VEC3(src[0].s1, src[1].s1, src[2].s1)),
123 :     dot(img->w2iT[2], VEC3(src[0].s2, src[1].s2, src[2].s2)));
124 :    
125 :     dst[0] = VEC3(
126 :     dot(tmp[0], VEC3(img->w2i[0].s0, img->w2i[1].s0, img->w2i[2].s0)),
127 :     dot(tmp[0], VEC3(img->w2i[0].s1, img->w2i[1].s1, img->w2i[2].s1)),
128 :     dot(tmp[0], VEC3(img->w2i[0].s2, img->w2i[1].s2, img->w2i[2].s2)));
129 :     dst[1] = VEC3(
130 :     dot(tmp[1], VEC3(img->w2i[0].s0, img->w2i[1].s0, img->w2i[2].s0)),
131 :     dot(tmp[1], VEC3(img->w2i[0].s1, img->w2i[1].s1, img->w2i[2].s1)),
132 :     dot(tmp[1], VEC3(img->w2i[0].s2, img->w2i[1].s2, img->w2i[2].s2)));
133 :     dst[2] = VEC3(
134 :     dot(tmp[2], VEC3(img->w2i[0].s0, img->w2i[1].s0, img->w2i[2].s0)),
135 :     dot(tmp[2], VEC3(img->w2i[0].s1, img->w2i[1].s1, img->w2i[2].s1)),
136 :     dot(tmp[2], VEC3(img->w2i[0].s2, img->w2i[1].s2, img->w2i[2].s2)));
137 :    
138 :     }
139 :    
140 :     #endif /* !_DIDEROT_CL_INLINE_IMAGE_H_ */

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