svn merge ^/trunk/blender -r47023:HEAD
[blender-staging.git] / intern / cycles / util / util_boundbox.h
1 /*
2  * Copyright 2011, Blender Foundation.
3  *
4  * This program is free software; you can redistribute it and/or
5  * modify it under the terms of the GNU General Public License
6  * as published by the Free Software Foundation; either version 2
7  * of the License, or (at your option) any later version.
8  *
9  * This program is distributed in the hope that it will be useful,
10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
12  * GNU General Public License for more details.
13  *
14  * You should have received a copy of the GNU General Public License
15  * along with this program; if not, write to the Free Software Foundation,
16  * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
17  */
18
19 #ifndef __UTIL_BOUNDBOX_H__
20 #define __UTIL_BOUNDBOX_H__
21
22 #include <math.h>
23 #include <float.h>
24
25 #include "util_math.h"
26 #include "util_string.h"
27 #include "util_transform.h"
28 #include "util_types.h"
29
30 using namespace std;
31
32 CCL_NAMESPACE_BEGIN
33
34 class BoundBox
35 {
36 public:
37         float3 min, max;
38
39         __forceinline BoundBox()
40         {
41         }
42
43         __forceinline BoundBox(const float3& pt)
44         : min(pt), max(pt)
45         {
46         }
47
48         __forceinline BoundBox(const float3& min_, const float3& max_)
49         : min(min_), max(max_)
50         {
51         }
52
53         enum empty_t { empty = 0};
54
55         __forceinline BoundBox(empty_t)
56         : min(make_float3(FLT_MAX, FLT_MAX, FLT_MAX)), max(make_float3(-FLT_MAX, -FLT_MAX, -FLT_MAX))
57         {
58         }
59
60         __forceinline void grow(const float3& pt)  
61         {
62                 min = ccl::min(min, pt);
63                 max = ccl::max(max, pt);
64         }
65
66         __forceinline void grow(const BoundBox& bbox)
67         {
68                 grow(bbox.min);
69                 grow(bbox.max);
70         }
71
72         __forceinline void intersect(const BoundBox& bbox) 
73         {
74                 min = ccl::max(min, bbox.min);
75                 max = ccl::min(max, bbox.max);
76         }
77
78         /* todo: avoid using this */
79         __forceinline float safe_area() const
80         {
81                 if(!((min.x <= max.x) && (min.y <= max.y) && (min.z <= max.z)))
82                         return 0.0f;
83
84                 return area();
85         }
86
87         __forceinline float area() const
88         {
89                 return half_area()*2.0f;
90         }
91
92         __forceinline float half_area() const
93         {
94                 float3 d = max - min;
95                 return (d.x*d.z + d.y*d.z + d.x*d.y);
96         }
97
98         __forceinline float3 center() const
99         {
100                 return 0.5f*(min + max);
101         }
102
103         __forceinline float3 center2() const
104         {
105                 return min + max;
106         }
107
108         __forceinline float3 size() const
109         {
110                 return max - min;
111         }
112         
113         __forceinline bool valid() const
114         {
115                 return (min.x <= max.x) && (min.y <= max.y) && (min.z <= max.z) &&
116                        (isfinite(min.x) && isfinite(min.y) && isfinite(min.z)) &&
117                        (isfinite(max.x) && isfinite(max.y) && isfinite(max.z));
118         }
119
120         BoundBox transformed(const Transform *tfm)
121         {
122                 BoundBox result = BoundBox::empty;
123
124                 for(int i = 0; i < 8; i++) {
125                         float3 p;
126
127                         p.x = (i & 1)? min.x: max.x;
128                         p.y = (i & 2)? min.y: max.y;
129                         p.z = (i & 4)? min.z: max.z;
130
131                         result.grow(transform_point(tfm, p));
132                 }
133
134                 return result;
135         }
136 };
137
138 __forceinline BoundBox merge(const BoundBox& bbox, const float3& pt)
139 {
140         return BoundBox(min(bbox.min, pt), max(bbox.max, pt));
141 }
142
143 __forceinline BoundBox merge(const BoundBox& a, const BoundBox& b)
144 {
145         return BoundBox(min(a.min, b.min), max(a.max, b.max));
146 }
147
148 __forceinline BoundBox merge(const BoundBox& a, const BoundBox& b, const BoundBox& c, const BoundBox& d)
149 {
150         return merge(merge(a, b), merge(c, d));
151 }
152
153 __forceinline BoundBox intersect(const BoundBox& a, const BoundBox& b)
154 {
155         return BoundBox(max(a.min, b.min), min(a.max, b.max));
156 }
157
158 __forceinline BoundBox intersect(const BoundBox& a, const BoundBox& b, const BoundBox& c)
159 {
160         return intersect(a, intersect(b, c));
161 }
162
163 CCL_NAMESPACE_END
164
165 #endif /* __UTIL_BOUNDBOX_H__ */
166