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

SCM Repository

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

Annotation of /trunk/src/include/Diderot/inline-vec3.h

Parent Directory Parent Directory | Revision Log Revision Log


Revision 3349 - (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 : jhr 3349 * This code is part of the Diderot Project (http://diderot-language.cs.uchicago.edu)
10 :     *
11 :     * COPYRIGHT (c) 2015 The University of Chicago
12 : jhr 1115 * All rights reserved.
13 :     */
14 :    
15 :     #ifndef _DIDEROT_INLINE_VEC3_H_
16 :     #define _DIDEROT_INLINE_VEC3_H_
17 :    
18 :     #ifndef _DIDEROT_TYPES_H_
19 :     #include "Diderot/types.h"
20 :     #endif
21 :    
22 : jhr 1640 STATIC_INLINE Diderot_vec3_t vec3 (Diderot_real_t a, Diderot_real_t b, Diderot_real_t c)
23 : jhr 1115 {
24 : jhr 1640 return __extension__ (Diderot_vec3_t){ a, b, c, 0.0f };
25 : jhr 1115 }
26 :    
27 : jhr 1640 STATIC_INLINE Diderot_ivec3_t vec3rtoi (Diderot_vec3_t v)
28 : jhr 1115 {
29 : jhr 1640 Diderot_union3_t u;
30 : jhr 1115 u.v = v;
31 : jhr 1640 return __extension__ (Diderot_ivec3_t){
32 :     (Diderot_int_t)u.r[0],
33 :     (Diderot_int_t)u.r[1],
34 :     (Diderot_int_t)u.r[2],
35 :     0
36 :     };
37 : jhr 1115 }
38 :    
39 : jhr 1640 STATIC_INLINE Diderot_vec3_t vec3itor (Diderot_ivec3_t v)
40 : jhr 1115 {
41 :     union3i_t u;
42 :     u.v = v;
43 : 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 };
44 : jhr 1115 }
45 :    
46 : jhr 1640 STATIC_INLINE Diderot_vec3_t scale3 (Diderot_real_t s, Diderot_vec3_t v)
47 : jhr 1115 {
48 : jhr 1640 return vec3(s, s, s) * v;
49 : jhr 1115 }
50 :    
51 : jhr 1640 STATIC_INLINE Diderot_vec3_t clamp3 (Diderot_vec3_t lo, Diderot_vec3_t hi, Diderot_vec3_t v)
52 : jhr 1295 {
53 :     // FIXME: there is probably a vectorized way to compute this
54 : jhr 1640 Diderot_union3_t a, b, c;
55 : jhr 1295 a.v = lo; b.v = hi; c.v = v;
56 : jhr 1640 return vec3(
57 :     clamp(a.r[0], b.r[0], c.r[0]),
58 :     clamp(a.r[1], b.r[1], c.r[1]),
59 :     clamp(a.r[2], b.r[2], c.r[2]));
60 : jhr 1295 }
61 :    
62 : jhr 1640 STATIC_INLINE Diderot_vec3_t lerp3 (Diderot_vec3_t a, Diderot_vec3_t b, Diderot_real_t t)
63 : jhr 1115 {
64 : jhr 1640 return a + scale3(t, b - a);
65 : jhr 1115 }
66 :    
67 : jhr 1640 STATIC_INLINE Diderot_vec3_t floor3 (Diderot_vec3_t v)
68 : jhr 1115 {
69 : jhr 1640 Diderot_union3_t u;
70 : jhr 1115 u.v = v;
71 : jhr 1640 return vec3(FLOOR(u.r[0]), FLOOR(u.r[1]), FLOOR(u.r[2]));
72 : jhr 1115 }
73 :    
74 : jhr 1640 STATIC_INLINE Diderot_ivec3_t truncToInt3 (Diderot_vec3_t v)
75 : jhr 1115 {
76 : jhr 1640 Diderot_union3_t t;
77 : jhr 1115 t.v = v;
78 : jhr 1640 return __extension__ (Diderot_ivec3_t){
79 :     (Diderot_int_t)TRUNC(t.r[0]),
80 :     (Diderot_int_t)TRUNC(t.r[1]),
81 :     (Diderot_int_t)TRUNC(t.r[2]),
82 : jhr 1115 0 };
83 :     }
84 :    
85 : jhr 1640 STATIC_INLINE Diderot_real_t dot3 (Diderot_vec3_t u, Diderot_vec3_t v)
86 : jhr 1115 {
87 : jhr 1640 Diderot_union3_t uv = __extension__ (Diderot_union3_t)(u*v);
88 : jhr 1115 return uv.r[0] + uv.r[1] + uv.r[2];
89 :     }
90 :    
91 : jhr 1640 STATIC_INLINE Diderot_real_t length3 (Diderot_vec3_t v)
92 : jhr 1115 {
93 : jhr 1640 return SQRT(dot3(v, v));
94 : jhr 1115 }
95 :    
96 : jhr 1640 STATIC_INLINE Diderot_vec3_t normalize3 (Diderot_vec3_t v)
97 : jhr 1115 {
98 : jhr 1640 return scale3(1.0 / length3(v), v);
99 : jhr 1115 }
100 :    
101 : jhr 1640 STATIC_INLINE Diderot_vec3_t cross3 (Diderot_vec3_t u, Diderot_vec3_t v)
102 : jhr 1115 {
103 : jhr 1640 Diderot_union3_t uu = (Diderot_union3_t)u;
104 :     Diderot_union3_t uv = (Diderot_union3_t)v;
105 :     return __extension__ (Diderot_vec3_t){
106 : jhr 1115 uu.r[1]*uv.r[2] - uu.r[2]*uv.r[1],
107 :     uu.r[2]*uv.r[0] - uu.r[0]*uv.r[2],
108 :     uu.r[0]*uv.r[1] - uu.r[1]*uv.r[0],
109 :     0 };
110 :     }
111 :    
112 :     #endif /* !_DIDEROT_INLINE_VEC3_H_ */

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