BLI_rect: Function to calculate a matrix from 2 rctf's
[blender.git] / source / blender / blenlib / intern / rct.c
1 /*
2  * ***** BEGIN GPL LICENSE BLOCK *****
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  * The Original Code is Copyright (C) 2001-2002 by NaN Holding BV.
19  * All rights reserved.
20  *
21  * The Original Code is: all of this file.
22  *
23  * Contributor(s): none yet.
24  *
25  * ***** END GPL LICENSE BLOCK *****
26  *
27  */
28
29 /** \file blender/blenlib/intern/rct.c
30  *  \ingroup bli
31  *
32  * A minimalist lib for functions doing stuff with rectangle structs.
33  */
34
35 #include <stdlib.h>
36 #include <stdio.h>
37 #include <math.h>
38
39 #include <limits.h>
40 #include <float.h>
41
42 #include "DNA_vec_types.h"
43 #include "BLI_rect.h"
44
45 /* avoid including BLI_math */
46 static void unit_m4(float m[4][4]);
47
48 /**
49  * Determine if a rect is empty. An empty
50  * rect is one with a zero (or negative)
51  * width or height.
52  *
53  * \return True if \a rect is empty.
54  */
55 bool BLI_rcti_is_empty(const rcti *rect)
56 {
57         return ((rect->xmax <= rect->xmin) || (rect->ymax <= rect->ymin));
58 }
59
60 bool BLI_rctf_is_empty(const rctf *rect)
61 {
62         return ((rect->xmax <= rect->xmin) || (rect->ymax <= rect->ymin));
63 }
64
65 bool BLI_rcti_isect_x(const rcti *rect, const int x)
66 {
67         if (x < rect->xmin) return false;
68         if (x > rect->xmax) return false;
69         return true;
70 }
71
72 bool BLI_rcti_isect_y(const rcti *rect, const int y)
73 {
74         if (y < rect->ymin) return false;
75         if (y > rect->ymax) return false;
76         return true;
77 }
78
79 bool BLI_rcti_isect_pt(const rcti *rect, const int x, const int y)
80 {
81         if (x < rect->xmin) return false;
82         if (x > rect->xmax) return false;
83         if (y < rect->ymin) return false;
84         if (y > rect->ymax) return false;
85         return true;
86 }
87
88 bool BLI_rcti_isect_pt_v(const rcti *rect, const int xy[2])
89 {
90         if (xy[0] < rect->xmin) return false;
91         if (xy[0] > rect->xmax) return false;
92         if (xy[1] < rect->ymin) return false;
93         if (xy[1] > rect->ymax) return false;
94         return true;
95 }
96
97 bool BLI_rctf_isect_x(const rctf *rect, const float x)
98 {
99         if (x < rect->xmin) return false;
100         if (x > rect->xmax) return false;
101         return true;
102 }
103
104 bool BLI_rctf_isect_y(const rctf *rect, const float y)
105 {
106         if (y < rect->ymin) return false;
107         if (y > rect->ymax) return false;
108         return true;
109 }
110
111 bool BLI_rctf_isect_pt(const rctf *rect, const float x, const float y)
112 {
113         if (x < rect->xmin) return false;
114         if (x > rect->xmax) return false;
115         if (y < rect->ymin) return false;
116         if (y > rect->ymax) return false;
117         return true;
118 }
119
120 bool BLI_rctf_isect_pt_v(const rctf *rect, const float xy[2])
121 {
122         if (xy[0] < rect->xmin) return false;
123         if (xy[0] > rect->xmax) return false;
124         if (xy[1] < rect->ymin) return false;
125         if (xy[1] > rect->ymax) return false;
126         return true;
127 }
128
129 /**
130  * \returns shortest distance from \a rect to x/y (0 if inside)
131 */
132
133 int BLI_rcti_length_x(const rcti *rect, const int x)
134 {
135         if (x < rect->xmin) return rect->xmin - x;
136         if (x > rect->xmax) return x - rect->xmax;
137         return 0;
138 }
139
140 int BLI_rcti_length_y(const rcti *rect, const int y)
141 {
142         if (y < rect->ymin) return rect->ymin - y;
143         if (y > rect->ymax) return y - rect->ymax;
144         return 0;
145 }
146
147 float BLI_rctf_length_x(const rctf *rect, const float x)
148 {
149         if (x < rect->xmin) return rect->xmin - x;
150         if (x > rect->xmax) return x - rect->xmax;
151         return 0.0f;
152 }
153
154 float BLI_rctf_length_y(const rctf *rect, const float y)
155 {
156         if (y < rect->ymin) return rect->ymin - y;
157         if (y > rect->ymax) return y - rect->ymax;
158         return 0.0f;
159 }
160
161 /**
162  * is \a rct_b inside \a rct_a
163  */
164 bool BLI_rctf_inside_rctf(rctf *rct_a, const rctf *rct_b)
165 {
166         return ((rct_a->xmin <= rct_b->xmin) &&
167                 (rct_a->xmax >= rct_b->xmax) &&
168                 (rct_a->ymin <= rct_b->ymin) &&
169                 (rct_a->ymax >= rct_b->ymax));
170 }
171 bool BLI_rcti_inside_rcti(rcti *rct_a, const rcti *rct_b)
172 {
173         return ((rct_a->xmin <= rct_b->xmin) &&
174                 (rct_a->xmax >= rct_b->xmax) &&
175                 (rct_a->ymin <= rct_b->ymin) &&
176                 (rct_a->ymax >= rct_b->ymax));
177 }
178
179
180 /* based closely on 'isect_seg_seg_v2_int', but in modified so corner cases are treated as intersections */
181 static int isect_segments_i(const int v1[2], const int v2[2], const int v3[2], const int v4[2])
182 {
183         const double div = (double)((v2[0] - v1[0]) * (v4[1] - v3[1]) - (v2[1] - v1[1]) * (v4[0] - v3[0]));
184         if (div == 0.0) {
185                 return 1; /* co-linear */
186         }
187         else {
188                 const double lambda = (double)((v1[1] - v3[1]) * (v4[0] - v3[0]) - (v1[0] - v3[0]) * (v4[1] - v3[1])) / div;
189                 const double mu    = (double)((v1[1] - v3[1]) * (v2[0] - v1[0]) - (v1[0] - v3[0]) * (v2[1] - v1[1])) / div;
190                 return (lambda >= 0.0 && lambda <= 1.0 && mu >= 0.0 && mu <= 1.0);
191         }
192 }
193 static int isect_segments_fl(const float v1[2], const float v2[2], const float v3[2], const float v4[2])
194 {
195         const double div = (double)((v2[0] - v1[0]) * (v4[1] - v3[1]) - (v2[1] - v1[1]) * (v4[0] - v3[0]));
196         if (div == 0.0) {
197                 return 1; /* co-linear */
198         }
199         else {
200                 const double lambda = (double)((v1[1] - v3[1]) * (v4[0] - v3[0]) - (v1[0] - v3[0]) * (v4[1] - v3[1])) / div;
201                 const double mu    = (double)((v1[1] - v3[1]) * (v2[0] - v1[0]) - (v1[0] - v3[0]) * (v2[1] - v1[1])) / div;
202                 return (lambda >= 0.0 && lambda <= 1.0 && mu >= 0.0 && mu <= 1.0);
203         }
204 }
205
206 bool BLI_rcti_isect_segment(const rcti *rect, const int s1[2], const int s2[2])
207 {
208         /* first do outside-bounds check for both points of the segment */
209         if (s1[0] < rect->xmin && s2[0] < rect->xmin) return false;
210         if (s1[0] > rect->xmax && s2[0] > rect->xmax) return false;
211         if (s1[1] < rect->ymin && s2[1] < rect->ymin) return false;
212         if (s1[1] > rect->ymax && s2[1] > rect->ymax) return false;
213
214         /* if either points intersect then we definetly intersect */
215         if (BLI_rcti_isect_pt_v(rect, s1) || BLI_rcti_isect_pt_v(rect, s2)) {
216                 return true;
217         }
218         else {
219                 /* both points are outside but may insersect the rect */
220                 int tvec1[2];
221                 int tvec2[2];
222                 /* diagonal: [/] */
223                 tvec1[0] = rect->xmin; tvec1[1] = rect->ymin;
224                 tvec2[0] = rect->xmin; tvec2[1] = rect->ymax;
225                 if (isect_segments_i(s1, s2, tvec1, tvec2)) {
226                         return true;
227                 }
228
229                 /* diagonal: [\] */
230                 tvec1[0] = rect->xmin; tvec1[1] = rect->ymax;
231                 tvec2[0] = rect->xmax; tvec2[1] = rect->ymin;
232                 if (isect_segments_i(s1, s2, tvec1, tvec2)) {
233                         return true;
234                 }
235
236                 /* no intersection */
237                 return false;
238         }
239 }
240
241 bool BLI_rctf_isect_segment(const rctf *rect, const float s1[2], const float s2[2])
242 {
243         /* first do outside-bounds check for both points of the segment */
244         if (s1[0] < rect->xmin && s2[0] < rect->xmin) return false;
245         if (s1[0] > rect->xmax && s2[0] > rect->xmax) return false;
246         if (s1[1] < rect->ymin && s2[1] < rect->ymin) return false;
247         if (s1[1] > rect->ymax && s2[1] > rect->ymax) return false;
248
249         /* if either points intersect then we definetly intersect */
250         if (BLI_rctf_isect_pt_v(rect, s1) || BLI_rctf_isect_pt_v(rect, s2)) {
251                 return true;
252         }
253         else {
254                 /* both points are outside but may insersect the rect */
255                 float tvec1[2];
256                 float tvec2[2];
257                 /* diagonal: [/] */
258                 tvec1[0] = rect->xmin; tvec1[1] = rect->ymin;
259                 tvec2[0] = rect->xmin; tvec2[1] = rect->ymax;
260                 if (isect_segments_fl(s1, s2, tvec1, tvec2)) {
261                         return true;
262                 }
263
264                 /* diagonal: [\] */
265                 tvec1[0] = rect->xmin; tvec1[1] = rect->ymax;
266                 tvec2[0] = rect->xmax; tvec2[1] = rect->ymin;
267                 if (isect_segments_fl(s1, s2, tvec1, tvec2)) {
268                         return true;
269                 }
270
271                 /* no intersection */
272                 return false;
273         }
274 }
275
276 bool BLI_rcti_isect_circle(const rcti *rect, const float xy[2], const float radius)
277 {
278         float dx, dy;
279
280         if (xy[0] >= rect->xmin && xy[0] <= rect->xmax) dx = 0;
281         else dx = (xy[0] < rect->xmin) ? (rect->xmin - xy[0]) : (xy[0] - rect->xmax);
282
283         if (xy[1] >= rect->ymin && xy[1] <= rect->ymax) dy = 0;
284         else dy = (xy[1] < rect->ymin) ? (rect->ymin - xy[1]) : (xy[1] - rect->ymax);
285
286         return dx * dx + dy * dy <= radius * radius;
287 }
288
289 bool BLI_rctf_isect_circle(const rctf *rect, const float xy[2], const float radius)
290 {
291         float dx, dy;
292
293         if (xy[0] >= rect->xmin && xy[0] <= rect->xmax) dx = 0;
294         else dx = (xy[0] < rect->xmin) ? (rect->xmin - xy[0]) : (xy[0] - rect->xmax);
295
296         if (xy[1] >= rect->ymin && xy[1] <= rect->ymax) dy = 0;
297         else dy = (xy[1] < rect->ymin) ? (rect->ymin - xy[1]) : (xy[1] - rect->ymax);
298
299         return dx * dx + dy * dy <= radius * radius;
300 }
301
302 void BLI_rctf_union(rctf *rct1, const rctf *rct2)
303 {
304         if (rct1->xmin > rct2->xmin) rct1->xmin = rct2->xmin;
305         if (rct1->xmax < rct2->xmax) rct1->xmax = rct2->xmax;
306         if (rct1->ymin > rct2->ymin) rct1->ymin = rct2->ymin;
307         if (rct1->ymax < rct2->ymax) rct1->ymax = rct2->ymax;
308 }
309
310 void BLI_rcti_union(rcti *rct1, const rcti *rct2)
311 {
312         if (rct1->xmin > rct2->xmin) rct1->xmin = rct2->xmin;
313         if (rct1->xmax < rct2->xmax) rct1->xmax = rct2->xmax;
314         if (rct1->ymin > rct2->ymin) rct1->ymin = rct2->ymin;
315         if (rct1->ymax < rct2->ymax) rct1->ymax = rct2->ymax;
316 }
317
318 void BLI_rctf_init(rctf *rect, float xmin, float xmax, float ymin, float ymax)
319 {
320         if (xmin <= xmax) {
321                 rect->xmin = xmin;
322                 rect->xmax = xmax;
323         }
324         else {
325                 rect->xmax = xmin;
326                 rect->xmin = xmax;
327         }
328         if (ymin <= ymax) {
329                 rect->ymin = ymin;
330                 rect->ymax = ymax;
331         }
332         else {
333                 rect->ymax = ymin;
334                 rect->ymin = ymax;
335         }
336 }
337
338 void BLI_rcti_init(rcti *rect, int xmin, int xmax, int ymin, int ymax)
339 {
340         if (xmin <= xmax) {
341                 rect->xmin = xmin;
342                 rect->xmax = xmax;
343         }
344         else {
345                 rect->xmax = xmin;
346                 rect->xmin = xmax;
347         }
348         if (ymin <= ymax) {
349                 rect->ymin = ymin;
350                 rect->ymax = ymax;
351         }
352         else {
353                 rect->ymax = ymin;
354                 rect->ymin = ymax;
355         }
356 }
357
358 void BLI_rctf_init_pt_radius(rctf *rect, const float xy[2], float size)
359 {
360         rect->xmin = xy[0] - size;
361         rect->xmax = xy[0] + size;
362         rect->ymin = xy[1] - size;
363         rect->ymax = xy[1] + size;
364 }
365
366 void BLI_rcti_init_pt_radius(rcti *rect, const int xy[2], int size)
367 {
368         rect->xmin = xy[0] - size;
369         rect->xmax = xy[0] + size;
370         rect->ymin = xy[1] - size;
371         rect->ymax = xy[1] + size;
372 }
373
374 void BLI_rcti_init_minmax(rcti *rect)
375 {
376         rect->xmin = rect->ymin = INT_MAX;
377         rect->xmax = rect->ymax = INT_MIN;
378 }
379
380 void BLI_rctf_init_minmax(rctf *rect)
381 {
382         rect->xmin = rect->ymin =  FLT_MAX;
383         rect->xmax = rect->ymax = -FLT_MAX;
384 }
385
386 void BLI_rcti_do_minmax_v(rcti *rect, const int xy[2])
387 {
388         if (xy[0] < rect->xmin) rect->xmin = xy[0];
389         if (xy[0] > rect->xmax) rect->xmax = xy[0];
390         if (xy[1] < rect->ymin) rect->ymin = xy[1];
391         if (xy[1] > rect->ymax) rect->ymax = xy[1];
392 }
393
394 void BLI_rctf_do_minmax_v(rctf *rect, const float xy[2])
395 {
396         if (xy[0] < rect->xmin) rect->xmin = xy[0];
397         if (xy[0] > rect->xmax) rect->xmax = xy[0];
398         if (xy[1] < rect->ymin) rect->ymin = xy[1];
399         if (xy[1] > rect->ymax) rect->ymax = xy[1];
400 }
401
402 /* given 2 rectangles - transform a point from one to another */
403 void BLI_rctf_transform_pt_v(const rctf *dst, const rctf *src, float xy_dst[2], const float xy_src[2])
404 {
405         xy_dst[0] = ((xy_src[0] - src->xmin) / (src->xmax - src->xmin));
406         xy_dst[0] =               dst->xmin + ((dst->xmax - dst->xmin) * xy_dst[0]);
407
408         xy_dst[1] = ((xy_src[1] - src->ymin) / (src->ymax - src->ymin));
409         xy_dst[1] =               dst->ymin + ((dst->ymax - dst->ymin) * xy_dst[1]);
410 }
411
412 /**
413  * Calculate a 4x4 matrix representing the transformation between two rectangles.
414  *
415  * \note Multiplying a vector by this matrix does *not* give the same value as #BLI_rctf_transform_pt_v.
416  */
417 void BLI_rctf_transform_calc_m4_pivot_min_ex(
418         const rctf *dst, const rctf *src, float matrix[4][4],
419         uint x, uint y)
420 {
421         BLI_assert(x < 3 && y < 3);
422
423         unit_m4(matrix);
424
425         matrix[x][x] = BLI_rctf_size_x(src) / BLI_rctf_size_x(dst);
426         matrix[y][y] = BLI_rctf_size_y(src) / BLI_rctf_size_y(dst);
427         matrix[3][x] = (src->xmin - dst->xmin) * matrix[x][x];
428         matrix[3][y] = (src->ymin - dst->ymin) * matrix[y][y];
429 }
430
431 void BLI_rctf_transform_calc_m4_pivot_min(
432         const rctf *dst, const rctf *src, float matrix[4][4])
433 {
434         BLI_rctf_transform_calc_m4_pivot_min_ex(dst, src, matrix, 0, 1);
435 }
436
437 void BLI_rcti_translate(rcti *rect, int x, int y)
438 {
439         rect->xmin += x;
440         rect->ymin += y;
441         rect->xmax += x;
442         rect->ymax += y;
443 }
444 void BLI_rctf_translate(rctf *rect, float x, float y)
445 {
446         rect->xmin += x;
447         rect->ymin += y;
448         rect->xmax += x;
449         rect->ymax += y;
450 }
451
452 void BLI_rcti_recenter(rcti *rect, int x, int y)
453 {
454         const int dx = x - BLI_rcti_cent_x(rect);
455         const int dy = y - BLI_rcti_cent_y(rect);
456         BLI_rcti_translate(rect, dx, dy);
457 }
458 void BLI_rctf_recenter(rctf *rect, float x, float y)
459 {
460         const float dx = x - BLI_rctf_cent_x(rect);
461         const float dy = y - BLI_rctf_cent_y(rect);
462         BLI_rctf_translate(rect, dx, dy);
463 }
464
465 /* change width & height around the central location */
466 void BLI_rcti_resize(rcti *rect, int x, int y)
467 {
468         rect->xmin = BLI_rcti_cent_x(rect) - (x / 2);
469         rect->ymin = BLI_rcti_cent_y(rect) - (y / 2);
470         rect->xmax = rect->xmin + x;
471         rect->ymax = rect->ymin + y;
472 }
473
474 void BLI_rctf_resize(rctf *rect, float x, float y)
475 {
476         rect->xmin = BLI_rctf_cent_x(rect) - (x * 0.5f);
477         rect->ymin = BLI_rctf_cent_y(rect) - (y * 0.5f);
478         rect->xmax = rect->xmin + x;
479         rect->ymax = rect->ymin + y;
480 }
481
482 void BLI_rcti_scale(rcti *rect, const float scale)
483 {
484         const int cent_x      = BLI_rcti_cent_x(rect);
485         const int cent_y      = BLI_rcti_cent_y(rect);
486         const int size_x_half = BLI_rcti_size_x(rect) * (scale * 0.5f);
487         const int size_y_half = BLI_rcti_size_y(rect) * (scale * 0.5f);
488         rect->xmin = cent_x - size_x_half;
489         rect->ymin = cent_y - size_y_half;
490         rect->xmax = cent_x + size_x_half;
491         rect->ymax = cent_y + size_y_half;
492 }
493
494 void BLI_rctf_scale(rctf *rect, const float scale)
495 {
496         const float cent_x      = BLI_rctf_cent_x(rect);
497         const float cent_y      = BLI_rctf_cent_y(rect);
498         const float size_x_half = BLI_rctf_size_x(rect) * (scale * 0.5f);
499         const float size_y_half = BLI_rctf_size_y(rect) * (scale * 0.5f);
500         rect->xmin = cent_x - size_x_half;
501         rect->ymin = cent_y - size_y_half;
502         rect->xmax = cent_x + size_x_half;
503         rect->ymax = cent_y + size_y_half;
504 }
505
506 void BLI_rctf_interp(rctf *rect, const rctf *rect_a, const rctf *rect_b, const float fac)
507 {
508         const float ifac = 1.0f - fac;
509         rect->xmin = (rect_a->xmin * ifac) + (rect_b->xmin * fac);
510         rect->xmax = (rect_a->xmax * ifac) + (rect_b->xmax * fac);
511         rect->ymin = (rect_a->ymin * ifac) + (rect_b->ymin * fac);
512         rect->ymax = (rect_a->ymax * ifac) + (rect_b->ymax * fac);
513 }
514
515 /* BLI_rcti_interp() not needed yet */
516
517
518 bool BLI_rctf_clamp_pt_v(const rctf *rect, float xy[2])
519 {
520         bool changed = false;
521         if (xy[0] < rect->xmin) { xy[0] = rect->xmin; changed = true; }
522         if (xy[0] > rect->xmax) { xy[0] = rect->xmax; changed = true; }
523         if (xy[1] < rect->ymin) { xy[1] = rect->ymin; changed = true; }
524         if (xy[1] > rect->ymax) { xy[1] = rect->ymax; changed = true; }
525         return changed;
526 }
527
528 bool BLI_rcti_clamp_pt_v(const rcti *rect, int xy[2])
529 {
530         bool changed = false;
531         if (xy[0] < rect->xmin) { xy[0] = rect->xmin; changed = true; }
532         if (xy[0] > rect->xmax) { xy[0] = rect->xmax; changed = true; }
533         if (xy[1] < rect->ymin) { xy[1] = rect->ymin; changed = true; }
534         if (xy[1] > rect->ymax) { xy[1] = rect->ymax; changed = true; }
535         return changed;
536 }
537
538 /**
539  * Clamp \a rect within \a rect_bounds, setting \a r_xy to the offset.
540  *
541  * \return true if a change is made.
542  */
543 bool BLI_rctf_clamp(rctf *rect, const rctf *rect_bounds, float r_xy[2])
544 {
545         bool changed = false;
546
547         r_xy[0] = 0.0f;
548         r_xy[1] = 0.0f;
549
550         if (rect->xmin < rect_bounds->xmin) {
551                 float ofs = rect_bounds->xmin - rect->xmin;
552                 rect->xmin += ofs;
553                 rect->xmax += ofs;
554                 r_xy[0] += ofs;
555                 changed = true;
556         }
557
558         if (rect->xmax > rect_bounds->xmax) {
559                 float ofs = rect_bounds->xmax - rect->xmax;
560                 rect->xmin += ofs;
561                 rect->xmax += ofs;
562                 r_xy[0] += ofs;
563                 changed = true;
564         }
565
566         if (rect->ymin < rect_bounds->ymin) {
567                 float ofs = rect_bounds->ymin - rect->ymin;
568                 rect->ymin += ofs;
569                 rect->ymax += ofs;
570                 r_xy[1] += ofs;
571                 changed = true;
572         }
573
574         if (rect->ymax > rect_bounds->ymax) {
575                 float ofs = rect_bounds->ymax - rect->ymax;
576                 rect->ymin += ofs;
577                 rect->ymax += ofs;
578                 r_xy[1] += ofs;
579                 changed = true;
580         }
581
582         return changed;
583 }
584
585 bool BLI_rcti_clamp(rcti *rect, const rcti *rect_bounds, int r_xy[2])
586 {
587         bool changed = false;
588
589         r_xy[0] = 0;
590         r_xy[1] = 0;
591
592         if (rect->xmin < rect_bounds->xmin) {
593                 int ofs = rect_bounds->xmin - rect->xmin;
594                 rect->xmin += ofs;
595                 rect->xmax += ofs;
596                 r_xy[0] += ofs;
597                 changed = true;
598         }
599
600         if (rect->xmax > rect_bounds->xmax) {
601                 int ofs = rect_bounds->xmax - rect->xmax;
602                 rect->xmin += ofs;
603                 rect->xmax += ofs;
604                 r_xy[0] += ofs;
605                 changed = true;
606         }
607
608         if (rect->ymin < rect_bounds->ymin) {
609                 int ofs = rect_bounds->ymin - rect->ymin;
610                 rect->ymin += ofs;
611                 rect->ymax += ofs;
612                 r_xy[1] += ofs;
613                 changed = true;
614         }
615
616         if (rect->ymax > rect_bounds->ymax) {
617                 int ofs = rect_bounds->ymax - rect->ymax;
618                 rect->ymin += ofs;
619                 rect->ymax += ofs;
620                 r_xy[1] += ofs;
621                 changed = true;
622         }
623
624         return changed;
625 }
626
627 bool BLI_rctf_compare(const rctf *rect_a, const rctf *rect_b, const float limit)
628 {
629         if (fabsf(rect_a->xmin - rect_b->xmin) < limit)
630                 if (fabsf(rect_a->xmax - rect_b->xmax) < limit)
631                         if (fabsf(rect_a->ymin - rect_b->ymin) < limit)
632                                 if (fabsf(rect_a->ymax - rect_b->ymax) < limit)
633                                         return true;
634
635         return false;
636 }
637
638 bool BLI_rcti_compare(const rcti *rect_a, const rcti *rect_b)
639 {
640         if (rect_a->xmin == rect_b->xmin)
641                 if (rect_a->xmax == rect_b->xmax)
642                         if (rect_a->ymin == rect_b->ymin)
643                                 if (rect_a->ymax == rect_b->ymax)
644                                         return true;
645
646         return false;
647 }
648
649 bool BLI_rctf_isect(const rctf *src1, const rctf *src2, rctf *dest)
650 {
651         float xmin, xmax;
652         float ymin, ymax;
653
654         xmin = (src1->xmin) > (src2->xmin) ? (src1->xmin) : (src2->xmin);
655         xmax = (src1->xmax) < (src2->xmax) ? (src1->xmax) : (src2->xmax);
656         ymin = (src1->ymin) > (src2->ymin) ? (src1->ymin) : (src2->ymin);
657         ymax = (src1->ymax) < (src2->ymax) ? (src1->ymax) : (src2->ymax);
658
659         if (xmax >= xmin && ymax >= ymin) {
660                 if (dest) {
661                         dest->xmin = xmin;
662                         dest->xmax = xmax;
663                         dest->ymin = ymin;
664                         dest->ymax = ymax;
665                 }
666                 return true;
667         }
668         else {
669                 if (dest) {
670                         dest->xmin = 0;
671                         dest->xmax = 0;
672                         dest->ymin = 0;
673                         dest->ymax = 0;
674                 }
675                 return false;
676         }
677 }
678
679 bool BLI_rcti_isect(const rcti *src1, const rcti *src2, rcti *dest)
680 {
681         int xmin, xmax;
682         int ymin, ymax;
683
684         xmin = (src1->xmin) > (src2->xmin) ? (src1->xmin) : (src2->xmin);
685         xmax = (src1->xmax) < (src2->xmax) ? (src1->xmax) : (src2->xmax);
686         ymin = (src1->ymin) > (src2->ymin) ? (src1->ymin) : (src2->ymin);
687         ymax = (src1->ymax) < (src2->ymax) ? (src1->ymax) : (src2->ymax);
688
689         if (xmax >= xmin && ymax >= ymin) {
690                 if (dest) {
691                         dest->xmin = xmin;
692                         dest->xmax = xmax;
693                         dest->ymin = ymin;
694                         dest->ymax = ymax;
695                 }
696                 return true;
697         }
698         else {
699                 if (dest) {
700                         dest->xmin = 0;
701                         dest->xmax = 0;
702                         dest->ymin = 0;
703                         dest->ymax = 0;
704                 }
705                 return false;
706         }
707 }
708
709 void BLI_rcti_rctf_copy(rcti *dst, const rctf *src)
710 {
711         dst->xmin = floorf(src->xmin + 0.5f);
712         dst->xmax = dst->xmin + floorf(BLI_rctf_size_x(src) + 0.5f);
713         dst->ymin = floorf(src->ymin + 0.5f);
714         dst->ymax = dst->ymin + floorf(BLI_rctf_size_y(src) + 0.5f);
715 }
716
717 void BLI_rcti_rctf_copy_floor(rcti *dst, const rctf *src)
718 {
719         dst->xmin = floorf(src->xmin);
720         dst->xmax = floorf(src->xmax);
721         dst->ymin = floorf(src->ymin);
722         dst->ymax = floorf(src->ymax);
723 }
724
725 void BLI_rcti_rctf_copy_round(rcti *dst, const rctf *src)
726 {
727         dst->xmin = floorf(src->xmin + 0.5f);
728         dst->xmax = floorf(src->xmax + 0.5f);
729         dst->ymin = floorf(src->ymin + 0.5f);
730         dst->ymax = floorf(src->ymax + 0.5f);
731 }
732
733 void BLI_rctf_rcti_copy(rctf *dst, const rcti *src)
734 {
735         dst->xmin = src->xmin;
736         dst->xmax = src->xmax;
737         dst->ymin = src->ymin;
738         dst->ymax = src->ymax;
739 }
740
741 void print_rctf(const char *str, const rctf *rect)
742 {
743         printf("%s: xmin %.8f, xmax %.8f, ymin %.8f, ymax %.8f (%.12fx%.12f)\n", str,
744                rect->xmin, rect->xmax, rect->ymin, rect->ymax, BLI_rctf_size_x(rect), BLI_rctf_size_y(rect));
745 }
746
747 void print_rcti(const char *str, const rcti *rect)
748 {
749         printf("%s: xmin %d, xmax %d, ymin %d, ymax %d (%dx%d)\n", str,
750                rect->xmin, rect->xmax, rect->ymin, rect->ymax, BLI_rcti_size_x(rect), BLI_rcti_size_y(rect));
751 }
752
753
754 /* -------------------------------------------------------------------- */
755 /* Comprehensive math (float only) */
756
757 /** \name Rect math functions
758  * \{ */
759
760 #define ROTATE_SINCOS(r_vec, mat2, vec) { \
761         (r_vec)[0] = (mat2)[1] * (vec)[0] + (+(mat2)[0]) * (vec)[1]; \
762         (r_vec)[1] = (mat2)[0] * (vec)[0] + (-(mat2)[1]) * (vec)[1]; \
763 } ((void)0)
764
765 /**
766  * Expand the rectangle to fit a rotated \a src.
767  */
768 void BLI_rctf_rotate_expand(rctf *dst, const rctf *src, const float angle)
769 {
770         const float mat2[2] = {sinf(angle), cosf(angle)};
771         const float cent[2] = {BLI_rctf_cent_x(src), BLI_rctf_cent_y(src)};
772         float corner[2], corner_rot[2], corder_max[2];
773
774         /* x is same for both corners */
775         corner[0] = src->xmax - cent[0];
776         corner[1] = src->ymax - cent[1];
777         ROTATE_SINCOS(corner_rot, mat2, corner);
778         corder_max[0] = fabsf(corner_rot[0]);
779         corder_max[1] = fabsf(corner_rot[1]);
780
781         corner[1] *= -1;
782         ROTATE_SINCOS(corner_rot, mat2, corner);
783         corder_max[0] = MAX2(corder_max[0], fabsf(corner_rot[0]));
784         corder_max[1] = MAX2(corder_max[1], fabsf(corner_rot[1]));
785
786         dst->xmin = cent[0] - corder_max[0];
787         dst->xmax = cent[0] + corder_max[0];
788         dst->ymin = cent[1] - corder_max[1];
789         dst->ymax = cent[1] + corder_max[1];
790 }
791
792 #undef ROTATE_SINCOS
793
794 /** \} */
795
796 static void unit_m4(float m[4][4])
797 {
798         m[0][0] = m[1][1] = m[2][2] = m[3][3] = 1.0f;
799         m[0][1] = m[0][2] = m[0][3] = 0.0f;
800         m[1][0] = m[1][2] = m[1][3] = 0.0f;
801         m[2][0] = m[2][1] = m[2][3] = 0.0f;
802         m[3][0] = m[3][1] = m[3][2] = 0.0f;
803 }