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

SCM Repository

[diderot] Annotation of /branches/lamont/src/include/Diderot/inline-vec3.h
ViewVC logotype

Annotation of /branches/lamont/src/include/Diderot/inline-vec3.h

Parent Directory Parent Directory | Revision Log Revision Log


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

1 : jhr 1115 /*! \file inline-vec3.h
2 :     *
3 :     * \author John Reppy
4 :     *
5 :     * \brief inline functions for 3D vectors
6 :     */
7 :    
8 :     /*
9 :     * COPYRIGHT (c) 2011 The Diderot Project (http://diderot-language.cs.uchicago.edu)
10 :     * All rights reserved.
11 :     */
12 :    
13 :     #ifndef _DIDEROT_INLINE_VEC3_H_
14 :     #define _DIDEROT_INLINE_VEC3_H_
15 :    
16 :     #ifndef _DIDEROT_TYPES_H_
17 :     #include "Diderot/types.h"
18 :     #endif
19 :    
20 : jhr 1640 STATIC_INLINE Diderot_vec3_t vec3 (Diderot_real_t a, Diderot_real_t b, Diderot_real_t c)
21 : jhr 1115 {
22 : jhr 1640 return __extension__ (Diderot_vec3_t){ a, b, c, 0.0f };
23 : jhr 1115 }
24 : lamonts 2301 STATIC_INLINE Diderot_ivec3_t vec3i (Diderot_int_t x, Diderot_int_t y, Diderot_int_t z)
25 :     {
26 :     return __extension__ (Diderot_ivec3_t){x,y,z, 0 };
27 :     }
28 : jhr 1640 STATIC_INLINE Diderot_ivec3_t vec3rtoi (Diderot_vec3_t v)
29 : jhr 1115 {
30 : jhr 1640 Diderot_union3_t u;
31 : jhr 1115 u.v = v;
32 : jhr 1640 return __extension__ (Diderot_ivec3_t){
33 :     (Diderot_int_t)u.r[0],
34 :     (Diderot_int_t)u.r[1],
35 :     (Diderot_int_t)u.r[2],
36 :     0
37 :     };
38 : jhr 1115 }
39 :    
40 : jhr 1640 STATIC_INLINE Diderot_vec3_t vec3itor (Diderot_ivec3_t v)
41 : jhr 1115 {
42 :     union3i_t u;
43 :     u.v = v;
44 : jhr 1640 return __extension__ (Diderot_vec3_t){ (Diderot_real_t)u.i[0], (Diderot_real_t)u.i[1], (Diderot_real_t)u.i[2], 0 };
45 : jhr 1115 }
46 :    
47 : jhr 1640 STATIC_INLINE Diderot_vec3_t scale3 (Diderot_real_t s, Diderot_vec3_t v)
48 : jhr 1115 {
49 : jhr 1640 return vec3(s, s, s) * v;
50 : jhr 1115 }
51 :    
52 : jhr 1640 STATIC_INLINE Diderot_vec3_t clamp3 (Diderot_vec3_t lo, Diderot_vec3_t hi, Diderot_vec3_t v)
53 : jhr 1295 {
54 :     // FIXME: there is probably a vectorized way to compute this
55 : jhr 1640 Diderot_union3_t a, b, c;
56 : jhr 1295 a.v = lo; b.v = hi; c.v = v;
57 : jhr 1640 return vec3(
58 :     clamp(a.r[0], b.r[0], c.r[0]),
59 :     clamp(a.r[1], b.r[1], c.r[1]),
60 :     clamp(a.r[2], b.r[2], c.r[2]));
61 : jhr 1295 }
62 :    
63 : jhr 1640 STATIC_INLINE Diderot_vec3_t lerp3 (Diderot_vec3_t a, Diderot_vec3_t b, Diderot_real_t t)
64 : jhr 1115 {
65 : jhr 1640 return a + scale3(t, b - a);
66 : jhr 1115 }
67 :    
68 : jhr 1640 STATIC_INLINE Diderot_vec3_t floor3 (Diderot_vec3_t v)
69 : jhr 1115 {
70 : jhr 1640 Diderot_union3_t u;
71 : jhr 1115 u.v = v;
72 : jhr 1640 return vec3(FLOOR(u.r[0]), FLOOR(u.r[1]), FLOOR(u.r[2]));
73 : jhr 1115 }
74 :    
75 : jhr 1640 STATIC_INLINE Diderot_ivec3_t truncToInt3 (Diderot_vec3_t v)
76 : jhr 1115 {
77 : jhr 1640 Diderot_union3_t t;
78 : jhr 1115 t.v = v;
79 : jhr 1640 return __extension__ (Diderot_ivec3_t){
80 :     (Diderot_int_t)TRUNC(t.r[0]),
81 :     (Diderot_int_t)TRUNC(t.r[1]),
82 :     (Diderot_int_t)TRUNC(t.r[2]),
83 : jhr 1115 0 };
84 :     }
85 :    
86 : jhr 1640 STATIC_INLINE Diderot_real_t dot3 (Diderot_vec3_t u, Diderot_vec3_t v)
87 : jhr 1115 {
88 : jhr 1640 Diderot_union3_t uv = __extension__ (Diderot_union3_t)(u*v);
89 : jhr 1115 return uv.r[0] + uv.r[1] + uv.r[2];
90 :     }
91 :    
92 : jhr 2418 STATIC_INLINE Diderot_real_t lengthSq3 (Diderot_vec3_t v)
93 :     {
94 :     return dot3(v, v);
95 :     }
96 :    
97 : jhr 1640 STATIC_INLINE Diderot_real_t length3 (Diderot_vec3_t v)
98 : jhr 1115 {
99 : jhr 2418 return SQRT(lengthSq3(v));
100 : jhr 1115 }
101 :    
102 : jhr 2418 STATIC_INLINE Diderot_real_t dist3 (Diderot_vec3_t u, Diderot_vec3_t v)
103 :     {
104 :     return length3(v-u);
105 :     }
106 :    
107 : jhr 1640 STATIC_INLINE Diderot_vec3_t normalize3 (Diderot_vec3_t v)
108 : jhr 1115 {
109 : jhr 1640 return scale3(1.0 / length3(v), v);
110 : jhr 1115 }
111 :    
112 : jhr 1640 STATIC_INLINE Diderot_vec3_t cross3 (Diderot_vec3_t u, Diderot_vec3_t v)
113 : jhr 1115 {
114 : jhr 1640 Diderot_union3_t uu = (Diderot_union3_t)u;
115 :     Diderot_union3_t uv = (Diderot_union3_t)v;
116 :     return __extension__ (Diderot_vec3_t){
117 : jhr 1115 uu.r[1]*uv.r[2] - uu.r[2]*uv.r[1],
118 :     uu.r[2]*uv.r[0] - uu.r[0]*uv.r[2],
119 :     uu.r[0]*uv.r[1] - uu.r[1]*uv.r[0],
120 :     0 };
121 :     }
122 :    
123 :     #endif /* !_DIDEROT_INLINE_VEC3_H_ */

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