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-matrix.h
ViewVC logotype

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

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : jhr 1671 /*! \file inline-matrix.h
2 :     *
3 :     * \author John Reppy & Lamont Samuels
4 :     */
5 :    
6 :     /*
7 :     * COPYRIGHT (c) 2011 The Diderot Project (http://diderot-language.cs.uchicago.edu)
8 :     * All rights reserved.
9 :     */
10 :    
11 :     #ifndef _DIDEROT_CL_INLINE_MATRIX_H_
12 :     #define _DIDEROT_CL_INLINE_MATRIX_H_
13 :    
14 :     #ifndef _DIDEROT_CL_TYPES_H_
15 :     #include "cl-types.h"
16 :     #endif
17 :    
18 :     /********** 2x2 matrix functions **********/
19 :    
20 :    
21 : jhr 3164 inline void zero2x2 (Diderot_Mat2x2_t dst)
22 : jhr 1671 {
23 : jhr 3164 dst[0] = (Diderot_vec2_t)(0.0, 0.0);
24 :     dst[1] = (Diderot_vec2_t)(0.0, 0.0);
25 : jhr 1671 }
26 :    
27 : jhr 3164 inline void identity2x2 (Diderot_Mat2x2_t dst)
28 : jhr 1671 {
29 : jhr 3164 dst[0] = (Diderot_vec2_t)(1.0, 0.0);
30 :     dst[1] = (Diderot_vec2_t)(0.0, 1.0);
31 : jhr 1671 }
32 :    
33 : jhr 3164 inline void copy2x2 (Diderot_Mat2x2_t dst, Diderot_Mat2x2_t src)
34 : jhr 1671 {
35 :     dst[0] = src[0];
36 :     dst[1] = src[1];
37 :     }
38 :    
39 : jhr 3164 inline void scale2x2 (Diderot_Mat2x2_t dst, float s, Diderot_Mat2x2_t src)
40 : jhr 1671 {
41 : jhr 3164 Diderot_vec2_t scale = (Diderot_vec2_t)(s, s);
42 : jhr 1671 dst[0] = scale * src[0];
43 :     dst[1] = scale * src[1];
44 :     }
45 :    
46 : jhr 3164 inline void add2x2 (Diderot_Mat2x2_t dst, Diderot_Mat2x2_t a, Diderot_Mat2x2_t b)
47 : jhr 1671 {
48 :     dst[0] = a[0] + b[0];
49 :     dst[1] = a[1] + b[1];
50 :     }
51 :    
52 : jhr 3164 inline void sub2x2 (Diderot_Mat2x2_t dst, Diderot_Mat2x2_t a, Diderot_Mat2x2_t b)
53 : jhr 1671 {
54 :     dst[0] = a[0] - b[0];
55 :     dst[1] = a[1] - b[1];
56 :     }
57 :    
58 : jhr 3164 inline Diderot_vec2_t mulVec3Mat2x2 (Diderot_vec2_t v, Diderot_Mat2x2_t m)
59 : jhr 1671 {
60 : jhr 3164 return (Diderot_vec2_t)(
61 :     dot(v, (Diderot_vec2_t)(m[0].s0, m[1].s0)),
62 :     dot(v, (Diderot_vec2_t)(m[0].s1, m[1].s1)));
63 : jhr 1671 }
64 :    
65 : jhr 3164 inline Diderot_vec2_t mulMat2x2Vec2 (Diderot_Mat2x2_t m, Diderot_vec2_t v)
66 : jhr 1671 {
67 : jhr 3164 return (Diderot_vec2_t)(dot(m[0], v), dot(m[1], v));
68 : jhr 1671 }
69 :    
70 : jhr 3164 inline void mulMat2x2Mat2x2 (Diderot_Mat2x2_t dst, Diderot_Mat2x2_t m1, Diderot_Mat2x2_t m2)
71 : jhr 1671 {
72 : jhr 3164 dst[0] = (Diderot_vec2_t)(
73 :     dot(m1[0], (Diderot_vec2_t)(m2[0].s0, m2[1].s0)),
74 :     dot(m1[0], (Diderot_vec2_t)(m2[0].s1, m2[1].s1)));
75 :     dst[1] = (Diderot_vec2_t)(
76 :     dot(m1[1], (Diderot_vec2_t)(m2[0].s0, m2[1].s0)),
77 :     dot(m1[1], (Diderot_vec2_t)(m2[0].s1, m2[1].s1)));
78 : jhr 1671 }
79 :    
80 : jhr 3164 inline void transpose2x2 (Diderot_Mat2x2_t dst, Diderot_Mat2x2_t src)
81 : jhr 1671 {
82 : jhr 3164 dst[0] = (Diderot_vec2_t)(src[0].s0, src[1].s0);
83 :     dst[1] = (Diderot_vec2_t)(src[0].s1, src[1].s1);
84 : jhr 1671 }
85 :    
86 :     // The Frobenius norm of a matrix is the sqrt of the sum of the squares of the elements
87 : jhr 3164 inline float norm2x2 (Diderot_Mat2x2_t m)
88 : jhr 1671 {
89 :     return sqrt(dot(m[0],m[0]) + dot(m[1],m[1]));
90 :     }
91 :    
92 :    
93 :     /********** 3x3 matrix functions **********/
94 :    
95 : jhr 3164 inline void zero3x3 (Diderot_Mat3x3_t dst)
96 : jhr 1671 {
97 :     dst[0] = VEC3(0.0, 0.0, 0.0);
98 :     dst[1] = VEC3(0.0, 0.0, 0.0);
99 :     dst[2] = VEC3(0.0, 0.0, 0.0);
100 :     }
101 :    
102 : jhr 3164 inline void identity3x3 (Diderot_Mat3x3_t dst)
103 : jhr 1671 {
104 :     dst[0] = VEC3(1.0, 0.0, 0.0);
105 :     dst[1] = VEC3(0.0, 1.0, 0.0);
106 :     dst[2] = VEC3(0.0, 0.0, 1.0);
107 :     }
108 :    
109 : jhr 3164 inline void copy3x3 (Diderot_Mat3x3_t dst, Diderot_Mat3x3_t src)
110 : jhr 1671 {
111 :     dst[0] = src[0];
112 :     dst[1] = src[1];
113 :     dst[2] = src[2];
114 :     }
115 :    
116 : jhr 3164 inline void scale3x3 (Diderot_Mat3x3_t dst, float s, Diderot_Mat3x3_t src)
117 : jhr 1671 {
118 : jhr 3164 Diderot_vec3_t scale = VEC3(s, s, s);
119 : jhr 1671 dst[0] = scale * src[0];
120 :     dst[1] = scale * src[1];
121 :     dst[2] = scale * src[2];
122 :     }
123 :    
124 : jhr 3164 inline void add3x3 (Diderot_Mat3x3_t dst, Diderot_Mat3x3_t a, Diderot_Mat3x3_t b)
125 : jhr 1671 {
126 :     dst[0] = a[0] + b[0];
127 :     dst[1] = a[1] + b[1];
128 :     dst[2] = a[2] + b[2];
129 :     }
130 :    
131 : jhr 3164 inline void sub3x3 (Diderot_Mat3x3_t dst, Diderot_Mat3x3_t a, Diderot_Mat3x3_t b)
132 : jhr 1671 {
133 :     dst[0] = a[0] - b[0];
134 :     dst[1] = a[1] - b[1];
135 :     dst[2] = a[2] - b[2];
136 :     }
137 :    
138 : jhr 3164 inline Diderot_vec3_t mulVec3Mat3x3 (Diderot_vec3_t v, Diderot_Mat3x3_t m)
139 : jhr 1671 {
140 :     return VEC3(
141 :     dot(v, VEC3(m[0].s0, m[1].s0, m[2].s0)),
142 :     dot(v, VEC3(m[0].s1, m[1].s1, m[2].s1)),
143 :     dot(v, VEC3(m[0].s2, m[1].s2, m[2].s2)));
144 :     }
145 :    
146 : jhr 3164 inline Diderot_vec3_t mulMat3x3Vec3 (Diderot_Mat3x3_t m, Diderot_vec3_t v)
147 : jhr 1671 {
148 : jhr 1940 return VEC3(dot(m[0], v), dot(m[1], v), dot(m[2], v));
149 : jhr 1671 }
150 :    
151 : jhr 3164 inline void mulMat3x3Mat3x3 (Diderot_Mat3x3_t dst, Diderot_Mat3x3_t m1, Diderot_Mat3x3_t m2)
152 : jhr 1671 {
153 :     dst[0] = VEC3(
154 :     dot(m1[0], VEC3(m2[0].s0, m2[1].s0, m2[2].s0)),
155 :     dot(m1[0], VEC3(m2[0].s1, m2[1].s1, m2[2].s1)),
156 :     dot(m1[0], VEC3(m2[0].s2, m2[1].s2, m2[2].s2)));
157 :     dst[1] = VEC3(
158 :     dot(m1[1], VEC3(m2[0].s0, m2[1].s0, m2[2].s0)),
159 :     dot(m1[1], VEC3(m2[0].s1, m2[1].s1, m2[2].s1)),
160 :     dot(m1[1], VEC3(m2[0].s2, m2[1].s2, m2[2].s2)));
161 :     dst[2] = VEC3(
162 :     dot(m1[2], VEC3(m2[0].s0, m2[1].s0, m2[2].s0)),
163 :     dot(m1[2], VEC3(m2[0].s1, m2[1].s1, m2[2].s1)),
164 :     dot(m1[2], VEC3(m2[0].s2, m2[1].s2, m2[2].s2)));
165 :     }
166 : jhr 3164 inline void transpose3x3 (Diderot_Mat3x3_t dst, Diderot_Mat3x3_t src)
167 : jhr 1671 {
168 :     dst[0] = VEC3(src[0].s0, src[1].s0, src[2].s0);
169 :     dst[1] = VEC3(src[0].s1, src[1].s1, src[2].s1);
170 :     dst[2] = VEC3(src[0].s2, src[1].s2, src[2].s2);
171 :     }
172 :    
173 :     // The Frobenius norm of a matrix is the sqrt of the sum of the squares of the elements
174 : jhr 3164 inline float norm3x3 (Diderot_Mat3x3_t m)
175 : jhr 1671 {
176 :     return sqrt(dot(m[0],m[0]) + dot(m[1],m[1]) + dot(m[2],m[2]));
177 :     }
178 :    
179 :    
180 :     /********** 4x4 matrix functions **********/
181 :    
182 : jhr 3164 inline void copy4x4 (Diderot_Mat4x4_t dst, Diderot_Mat4x4_t src)
183 : jhr 1671 {
184 :     dst[0] = src[0];
185 :     dst[1] = src[1];
186 :     dst[2] = src[2];
187 :     dst[3] = src[3];
188 :     }
189 :    
190 : jhr 3164 inline void transpose4x4 (Diderot_Mat4x4_t dst, Diderot_Mat4x4_t src)
191 : jhr 1671 {
192 : jhr 3164 dst[0] = (Diderot_vec4_t)(src[0].s0, src[1].s0, src[2].s0, src[3].s0);
193 :     dst[1] = (Diderot_vec4_t)(src[0].s1, src[1].s1, src[2].s1, src[3].s1);
194 :     dst[2] = (Diderot_vec4_t)(src[0].s2, src[1].s2, src[2].s2, src[3].s2);
195 :     dst[3] = (Diderot_vec4_t)(src[0].s3, src[1].s3, src[2].s3, src[3].s3);
196 : jhr 1671 }
197 :    
198 :     // The Frobenius norm of a matrix is the sqrt of the sum of the squares of the elements
199 : jhr 3164 inline float norm4x4 (Diderot_Mat4x4_t m)
200 : jhr 1671 {
201 :     return sqrt(dot(m[0],m[0]) + dot(m[1],m[1]) + dot(m[2],m[2]) + dot(m[3],m[3]));
202 :     }
203 :    
204 :     #endif /* !_DIDEROT_CL_INLINE_MATRIX_H_ */

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