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