src/mesa/main/macros.h

   1 /**
   2  * \file macros.h
   3  * A collection of useful macros.
   4  */
   5
   6 /*
   7  * Mesa 3-D graphics library
   8  * Version:  4.0.3
   9  *
  10  * Copyright (C) 1999-2002  Brian Paul   All Rights Reserved.
  11  *
  12  * Permission is hereby granted, free of charge, to any person obtaining a
  13  * copy of this software and associated documentation files (the "Software"),
  14  * to deal in the Software without restriction, including without limitation
  15  * the rights to use, copy, modify, merge, publish, distribute, sublicense,
  16  * and/or sell copies of the Software, and to permit persons to whom the
  17  * Software is furnished to do so, subject to the following conditions:
  18  *
  19  * The above copyright notice and this permission notice shall be included
  20  * in all copies or substantial portions of the Software.
  21  *
  22  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
  23  * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  24  * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.  IN NO EVENT SHALL
  25  * BRIAN PAUL BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN
  26  * AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
  27  * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
  28  */
  29
  30
  31 #ifndef MACROS_H
  32 #define MACROS_H
  33
  34 #include "imports.h"
  35
  36
  37 /**
  38  * \name Integer / float conversion for colors, normals, etc.
  39  */
  40 /*@{*/
  41
  42 /** Convert GLubyte in [0,255] to GLfloat in [0.0,1.0] */
  43 extern GLfloat _mesa_ubyte_to_float_color_tab[256];
  44 #define UBYTE_TO_FLOAT(u) _mesa_ubyte_to_float_color_tab[(unsigned int)(u)]
  45
  46 /** Convert GLfloat in [0.0,1.0] to GLubyte in [0,255] */
  47 #define FLOAT_TO_UBYTE(X)       ((GLubyte) (GLint) ((X) * 255.0F))
  48
  49
  50 /** Convert GLbyte in [-128,127] to GLfloat in [-1.0,1.0] */
  51 #define BYTE_TO_FLOAT(B)        ((2.0F * (B) + 1.0F) * (1.0F/255.0F))
  52
  53 /** Convert GLfloat in [-1.0,1.0] to GLbyte in [-128,127] */
  54 #define FLOAT_TO_BYTE(X)        ( (((GLint) (255.0F * (X))) - 1) / 2 )
  55
  56
  57 /** Convert GLushort in [0,65536] to GLfloat in [0.0,1.0] */
  58 #define USHORT_TO_FLOAT(S)      ((GLfloat) (S) * (1.0F / 65535.0F))
  59
  60 /** Convert GLfloat in [0.0,1.0] to GLushort in [0,65536] */
  61 #define FLOAT_TO_USHORT(X)      ((GLushort) (GLint) ((X) * 65535.0F))
  62
  63 /** Convert GLshort in [-32768,32767] to GLfloat in [-1.0,1.0] */
  64 #define SHORT_TO_FLOAT(S)       ((2.0F * (S) + 1.0F) * (1.0F/65535.0F))
  65
  66 /** Convert GLfloat in [0.0,1.0] to GLshort in [-32768,32767] */
  67 #define FLOAT_TO_SHORT(X)       ( (((GLint) (65535.0F * (X))) - 1) / 2 )
  68
  69
  70 /** Convert GLuint in [0,4294967295] to GLfloat in [0.0,1.0] */
  71 #define UINT_TO_FLOAT(U)        ((GLfloat) (U) * (1.0F / 4294967295.0F))
  72
  73 /** Convert GLfloat in [0.0,1.0] to GLuint in [0,4294967295] */
  74 #define FLOAT_TO_UINT(X)        ((GLuint) ((X) * 4294967295.0))
  75
  76
  77 /** Convert GLint in [-2147483648,2147483647] to GLfloat in [-1.0,1.0] */
  78 #define INT_TO_FLOAT(I)         ((2.0F * (I) + 1.0F) * (1.0F/4294967294.0F))
  79
  80 /** Convert GLfloat in [-1.0,1.0] to GLint in [-2147483648,2147483647] */
  81 /* causes overflow:
  82 #define FLOAT_TO_INT(X)         ( (((GLint) (4294967294.0F * (X))) - 1) / 2 )
  83 */
  84 /* a close approximation: */
  85 #define FLOAT_TO_INT(X)         ( (GLint) (2147483647.0 * (X)) )
  86
  87
  88 #define BYTE_TO_UBYTE(b)   ((GLubyte) ((b) < 0 ? 0 : (GLubyte) (b)))
  89 #define SHORT_TO_UBYTE(s)  ((GLubyte) ((s) < 0 ? 0 : (GLubyte) ((s) >> 7)))
  90 #define USHORT_TO_UBYTE(s) ((GLubyte) ((s) >> 8))
  91 #define INT_TO_UBYTE(i)    ((GLubyte) ((i) < 0 ? 0 : (GLubyte) ((i) >> 23)))
  92 #define UINT_TO_UBYTE(i)   ((GLubyte) ((i) >> 24))
  93
  94
  95 #define BYTE_TO_USHORT(b)  ((b) < 0 ? 0 : ((GLushort) (((b) * 65535) / 255)))
  96 #define UBYTE_TO_USHORT(b) (((GLushort) (b) << 8) | (GLushort) (b))
  97 #define SHORT_TO_USHORT(s) ((s) < 0 ? 0 : ((GLushort) (((s) * 65535 / 32767))))
  98 #define INT_TO_USHORT(i)   ((i) < 0 ? 0 : ((GLushort) ((i) >> 15)))
  99 #define UINT_TO_USHORT(i)  ((i) < 0 ? 0 : ((GLushort) ((i) >> 16)))
 100 #define UNCLAMPED_FLOAT_TO_USHORT(us, f)  \
 101         us = ( (GLushort) IROUND( CLAMP((f), 0.0, 1.0) * 65535.0F) )
 102 #define CLAMPED_FLOAT_TO_USHORT(us, f)  \
 103         us = ( (GLushort) IROUND( (f) * 65535.0F) )
 104
 105
 106 /** Stepping a GLfloat pointer by a byte stride */
 107 #define STRIDE_F(p, i)  (p = (GLfloat *)((GLubyte *)p + i))
 108 /** Stepping a GLuint pointer by a byte stride */
 109 #define STRIDE_UI(p, i)  (p = (GLuint *)((GLubyte *)p + i))
 110 /** Stepping a GLubyte[4] pointer by a byte stride */
 111 #define STRIDE_4UB(p, i)  (p = (GLubyte (*)[4])((GLubyte *)p + i))
 112 /** Stepping a GLchan[4] pointer by a byte stride */
 113 #define STRIDE_4CHAN(p, i)  (p = (GLchan (*)[4])((GLubyte *)p + i))
 114 /** Stepping a GLchan pointer by a byte stride */
 115 #define STRIDE_CHAN(p, i)  (p = (GLchan *)((GLubyte *)p + i))
 116 /** Stepping a \p t pointer by a byte stride */
 117 #define STRIDE_T(p, t, i)  (p = (t)((GLubyte *)p + i))
 118
 119
 120 /**********************************************************************/
 121 /** \name 4-element vector operations */
 122 /*@{*/
 123
 124 /** Zero */
 125 #define ZERO_4V( DST )  (DST)[0] = (DST)[1] = (DST)[2] = (DST)[3] = 0
 126
 127 /** Test for equality */
 128 #define TEST_EQ_4V(a,b)  ((a)[0] == (b)[0] &&   \
 129                           (a)[1] == (b)[1] &&   \
 130                           (a)[2] == (b)[2] &&   \
 131                           (a)[3] == (b)[3])
 132
 133 /** Test for equality (unsigned bytes) */
 134 #if defined(__i386__)
 135 #define TEST_EQ_4UBV(DST, SRC) *((GLuint*)(DST)) == *((GLuint*)(SRC))
 136 #else
 137 #define TEST_EQ_4UBV(DST, SRC) TEST_EQ_4V(DST, SRC)
 138 #endif
 139
 140 /** Copy a 4-element vector */
 141 #define COPY_4V( DST, SRC )                     \
 142 do {                                            \
 143    (DST)[0] = (SRC)[0];                         \
 144    (DST)[1] = (SRC)[1];                         \
 145    (DST)[2] = (SRC)[2];                         \
 146    (DST)[3] = (SRC)[3];                         \
 147 } while (0)
 148
 149 /** Copy a 4-element vector with cast */
 150 #define COPY_4V_CAST( DST, SRC, CAST )          \
 151 do {                                            \
 152    (DST)[0] = (CAST)(SRC)[0];                   \
 153    (DST)[1] = (CAST)(SRC)[1];                   \
 154    (DST)[2] = (CAST)(SRC)[2];                   \
 155    (DST)[3] = (CAST)(SRC)[3];                   \
 156 } while (0)
 157
 158 /** Copy a 4-element unsigned byte vector */
 159 #if defined(__i386__)
 160 #define COPY_4UBV(DST, SRC)                     \
 161 do {                                            \
 162    *((GLuint*)(DST)) = *((GLuint*)(SRC));       \
 163 } while (0)
 164 #else
 165 /* The GLuint cast might fail if DST or SRC are not dword-aligned (RISC) */
 166 #define COPY_4UBV(DST, SRC)                     \
 167 do {                                            \
 168    (DST)[0] = (SRC)[0];                         \
 169    (DST)[1] = (SRC)[1];                         \
 170    (DST)[2] = (SRC)[2];                         \
 171    (DST)[3] = (SRC)[3];                         \
 172 } while (0)
 173 #endif
 174
 175 /** Copy a 4-element float vector */
 176 #define COPY_4FV( DST, SRC )                    \
 177 do {                                            \
 178    const GLfloat *_tmp = (SRC);                 \
 179    (DST)[0] = _tmp[0];                          \
 180    (DST)[1] = _tmp[1];                          \
 181    (DST)[2] = _tmp[2];                          \
 182    (DST)[3] = _tmp[3];                          \
 183 } while (0)
 184
 185
 186 /** Copy \p SZ elements into a 4-element vector */
 187 #define COPY_SZ_4V(DST, SZ, SRC)                \
 188 do {                                            \
 189    switch (SZ) {                                \
 190    case 4: (DST)[3] = (SRC)[3];                 \
 191    case 3: (DST)[2] = (SRC)[2];                 \
 192    case 2: (DST)[1] = (SRC)[1];                 \
 193    case 1: (DST)[0] = (SRC)[0];                 \
 194    }                                            \
 195 } while(0)
 196
 197 /** Copy \p SZ elements into a homegeneous (4-element) vector, giving
 198  * default values to the remaining */
 199 #define COPY_CLEAN_4V(DST, SZ, SRC)             \
 200 do {                                            \
 201       ASSIGN_4V( DST, 0, 0, 0, 1 );             \
 202       COPY_SZ_4V( DST, SZ, SRC );               \
 203 } while (0)
 204
 205 /** Subtraction */
 206 #define SUB_4V( DST, SRCA, SRCB )               \
 207 do {                                            \
 208       (DST)[0] = (SRCA)[0] - (SRCB)[0];         \
 209       (DST)[1] = (SRCA)[1] - (SRCB)[1];         \
 210       (DST)[2] = (SRCA)[2] - (SRCB)[2];         \
 211       (DST)[3] = (SRCA)[3] - (SRCB)[3];         \
 212 } while (0)
 213
 214 /** Addition */
 215 #define ADD_4V( DST, SRCA, SRCB )               \
 216 do {                                            \
 217       (DST)[0] = (SRCA)[0] + (SRCB)[0];         \
 218       (DST)[1] = (SRCA)[1] + (SRCB)[1];         \
 219       (DST)[2] = (SRCA)[2] + (SRCB)[2];         \
 220       (DST)[3] = (SRCA)[3] + (SRCB)[3];         \
 221 } while (0)
 222
 223 /** Element-wise multiplication */
 224 #define SCALE_4V( DST, SRCA, SRCB )             \
 225 do {                                            \
 226       (DST)[0] = (SRCA)[0] * (SRCB)[0];         \
 227       (DST)[1] = (SRCA)[1] * (SRCB)[1];         \
 228       (DST)[2] = (SRCA)[2] * (SRCB)[2];         \
 229       (DST)[3] = (SRCA)[3] * (SRCB)[3];         \
 230 } while (0)
 231
 232 /** In-place addition */
 233 #define ACC_4V( DST, SRC )                      \
 234 do {                                            \
 235       (DST)[0] += (SRC)[0];                     \
 236       (DST)[1] += (SRC)[1];                     \
 237       (DST)[2] += (SRC)[2];                     \
 238       (DST)[3] += (SRC)[3];                     \
 239 } while (0)
 240
 241 /** Element-wise multiplication and addition */
 242 #define ACC_SCALE_4V( DST, SRCA, SRCB )         \
 243 do {                                            \
 244       (DST)[0] += (SRCA)[0] * (SRCB)[0];        \
 245       (DST)[1] += (SRCA)[1] * (SRCB)[1];        \
 246       (DST)[2] += (SRCA)[2] * (SRCB)[2];        \
 247       (DST)[3] += (SRCA)[3] * (SRCB)[3];        \
 248 } while (0)
 249
 250 /** In-place scalar multiplication and addition */
 251 #define ACC_SCALE_SCALAR_4V( DST, S, SRCB )     \
 252 do {                                            \
 253       (DST)[0] += S * (SRCB)[0];                \
 254       (DST)[1] += S * (SRCB)[1];                \
 255       (DST)[2] += S * (SRCB)[2];                \
 256       (DST)[3] += S * (SRCB)[3];                \
 257 } while (0)
 258
 259 /** Scalar multiplication */
 260 #define SCALE_SCALAR_4V( DST, S, SRCB )         \
 261 do {                                            \
 262       (DST)[0] = S * (SRCB)[0];                 \
 263       (DST)[1] = S * (SRCB)[1];                 \
 264       (DST)[2] = S * (SRCB)[2];                 \
 265       (DST)[3] = S * (SRCB)[3];                 \
 266 } while (0)
 267
 268 /** In-place scalar multiplication */
 269 #define SELF_SCALE_SCALAR_4V( DST, S )          \
 270 do {                                            \
 271       (DST)[0] *= S;                            \
 272       (DST)[1] *= S;                            \
 273       (DST)[2] *= S;                            \
 274       (DST)[3] *= S;                            \
 275 } while (0)
 276
 277 /** Assignment */
 278 #define ASSIGN_4V( V, V0, V1, V2, V3 )          \
 279 do {                                            \
 280     V[0] = V0;                                  \
 281     V[1] = V1;                                  \
 282     V[2] = V2;                                  \
 283     V[3] = V3;                                  \
 284 } while(0)
 285
 286 /*@}*/
 287
 288
 289 /**********************************************************************/
 290 /** \name 3-element vector operations*/
 291 /*@{*/
 292
 293 /** Zero */
 294 #define ZERO_3V( DST )  (DST)[0] = (DST)[1] = (DST)[2] = 0
 295
 296 /** Test for equality */
 297 #define TEST_EQ_3V(a,b)  ((a)[0] == (b)[0] &&   \
 298                           (a)[1] == (b)[1] &&   \
 299                           (a)[2] == (b)[2])
 300
 301 /** Copy a 3-element vector */
 302 #define COPY_3V( DST, SRC )                     \
 303 do {                                            \
 304    (DST)[0] = (SRC)[0];                         \
 305    (DST)[1] = (SRC)[1];                         \
 306    (DST)[2] = (SRC)[2];                         \
 307 } while (0)
 308
 309 /** Copy a 3-element vector with cast */
 310 #define COPY_3V_CAST( DST, SRC, CAST )          \
 311 do {                                            \
 312    (DST)[0] = (CAST)(SRC)[0];                   \
 313    (DST)[1] = (CAST)(SRC)[1];                   \
 314    (DST)[2] = (CAST)(SRC)[2];                   \
 315 } while (0)
 316
 317 /** Copy a 3-element float vector */
 318 #define COPY_3FV( DST, SRC )                    \
 319 do {                                            \
 320    const GLfloat *_tmp = (SRC);                 \
 321    (DST)[0] = _tmp[0];                          \
 322    (DST)[1] = _tmp[1];                          \
 323    (DST)[2] = _tmp[2];                          \
 324 } while (0)
 325
 326 /** Subtraction */
 327 #define SUB_3V( DST, SRCA, SRCB )               \
 328 do {                                            \
 329       (DST)[0] = (SRCA)[0] - (SRCB)[0];         \
 330       (DST)[1] = (SRCA)[1] - (SRCB)[1];         \
 331       (DST)[2] = (SRCA)[2] - (SRCB)[2];         \
 332 } while (0)
 333
 334 /** Addition */
 335 #define ADD_3V( DST, SRCA, SRCB )               \
 336 do {                                            \
 337       (DST)[0] = (SRCA)[0] + (SRCB)[0];         \
 338       (DST)[1] = (SRCA)[1] + (SRCB)[1];         \
 339       (DST)[2] = (SRCA)[2] + (SRCB)[2];         \
 340 } while (0)
 341
 342 /** In-place scalar multiplication */
 343 #define SCALE_3V( DST, SRCA, SRCB )             \
 344 do {                                            \
 345       (DST)[0] = (SRCA)[0] * (SRCB)[0];         \
 346       (DST)[1] = (SRCA)[1] * (SRCB)[1];         \
 347       (DST)[2] = (SRCA)[2] * (SRCB)[2];         \
 348 } while (0)
 349
 350 /** In-place element-wise multiplication */
 351 #define SELF_SCALE_3V( DST, SRC )               \
 352 do {                                            \
 353       (DST)[0] *= (SRC)[0];                     \
 354       (DST)[1] *= (SRC)[1];                     \
 355       (DST)[2] *= (SRC)[2];                     \
 356 } while (0)
 357
 358 /** In-place addition */
 359 #define ACC_3V( DST, SRC )                      \
 360 do {                                            \
 361       (DST)[0] += (SRC)[0];                     \
 362       (DST)[1] += (SRC)[1];                     \
 363       (DST)[2] += (SRC)[2];                     \
 364 } while (0)
 365
 366 /** Element-wise multiplication and addition */
 367 #define ACC_SCALE_3V( DST, SRCA, SRCB )         \
 368 do {                                            \
 369       (DST)[0] += (SRCA)[0] * (SRCB)[0];        \
 370       (DST)[1] += (SRCA)[1] * (SRCB)[1];        \
 371       (DST)[2] += (SRCA)[2] * (SRCB)[2];        \
 372 } while (0)
 373
 374 /** Scalar multiplication */
 375 #define SCALE_SCALAR_3V( DST, S, SRCB )         \
 376 do {                                            \
 377       (DST)[0] = S * (SRCB)[0];                 \
 378       (DST)[1] = S * (SRCB)[1];                 \
 379       (DST)[2] = S * (SRCB)[2];                 \
 380 } while (0)
 381
 382 /** In-place scalar multiplication and addition */
 383 #define ACC_SCALE_SCALAR_3V( DST, S, SRCB )     \
 384 do {                                            \
 385       (DST)[0] += S * (SRCB)[0];                \
 386       (DST)[1] += S * (SRCB)[1];                \
 387       (DST)[2] += S * (SRCB)[2];                \
 388 } while (0)
 389
 390 /** In-place scalar multiplication */
 391 #define SELF_SCALE_SCALAR_3V( DST, S )          \
 392 do {                                            \
 393       (DST)[0] *= S;                            \
 394       (DST)[1] *= S;                            \
 395       (DST)[2] *= S;                            \
 396 } while (0)
 397
 398 /** In-place scalar addition */
 399 #define ACC_SCALAR_3V( DST, S )                 \
 400 do {                                            \
 401       (DST)[0] += S;                            \
 402       (DST)[1] += S;                            \
 403       (DST)[2] += S;                            \
 404 } while (0)
 405
 406 /** Assignment */
 407 #define ASSIGN_3V( V, V0, V1, V2 )      \
 408 do {                                    \
 409     V[0] = V0;                          \
 410     V[1] = V1;                          \
 411     V[2] = V2;                          \
 412 } while(0)
 413
 414 /*@}*/
 415
 416
 417 /**********************************************************************/
 418 /** \name 2-element vector operations*/
 419 /*@{*/
 420
 421 /** Zero */
 422 #define ZERO_2V( DST )  (DST)[0] = (DST)[1] = 0
 423
 424 /** Copy a 2-element vector */
 425 #define COPY_2V( DST, SRC )                     \
 426 do {                                            \
 427    (DST)[0] = (SRC)[0];                         \
 428    (DST)[1] = (SRC)[1];                         \
 429 } while (0)
 430
 431 /** Copy a 2-element vector with cast */
 432 #define COPY_2V_CAST( DST, SRC, CAST )          \
 433 do {                                            \
 434    (DST)[0] = (CAST)(SRC)[0];                   \
 435    (DST)[1] = (CAST)(SRC)[1];                   \
 436 } while (0)
 437
 438 /** Copy a 2-element float vector */
 439 #define COPY_2FV( DST, SRC )                    \
 440 do {                                            \
 441    const GLfloat *_tmp = (SRC);                 \
 442    (DST)[0] = _tmp[0];                          \
 443    (DST)[1] = _tmp[1];                          \
 444 } while (0)
 445
 446 /** Subtraction */
 447 #define SUB_2V( DST, SRCA, SRCB )               \
 448 do {                                            \
 449       (DST)[0] = (SRCA)[0] - (SRCB)[0];         \
 450       (DST)[1] = (SRCA)[1] - (SRCB)[1];         \
 451 } while (0)
 452
 453 /** Addition */
 454 #define ADD_2V( DST, SRCA, SRCB )               \
 455 do {                                            \
 456       (DST)[0] = (SRCA)[0] + (SRCB)[0];         \
 457       (DST)[1] = (SRCA)[1] + (SRCB)[1];         \
 458 } while (0)
 459
 460 /** In-place scalar multiplication */
 461 #define SCALE_2V( DST, SRCA, SRCB )             \
 462 do {                                            \
 463       (DST)[0] = (SRCA)[0] * (SRCB)[0];         \
 464       (DST)[1] = (SRCA)[1] * (SRCB)[1];         \
 465 } while (0)
 466
 467 /** In-place addition */
 468 #define ACC_2V( DST, SRC )                      \
 469 do {                                            \
 470       (DST)[0] += (SRC)[0];                     \
 471       (DST)[1] += (SRC)[1];                     \
 472 } while (0)
 473
 474 /** Element-wise multiplication and addition */
 475 #define ACC_SCALE_2V( DST, SRCA, SRCB )         \
 476 do {                                            \
 477       (DST)[0] += (SRCA)[0] * (SRCB)[0];        \
 478       (DST)[1] += (SRCA)[1] * (SRCB)[1];        \
 479 } while (0)
 480
 481 /** Scalar multiplication */
 482 #define SCALE_SCALAR_2V( DST, S, SRCB )         \
 483 do {                                            \
 484       (DST)[0] = S * (SRCB)[0];                 \
 485       (DST)[1] = S * (SRCB)[1];                 \
 486 } while (0)
 487
 488 /** In-place scalar multiplication and addition */
 489 #define ACC_SCALE_SCALAR_2V( DST, S, SRCB )     \
 490 do {                                            \
 491       (DST)[0] += S * (SRCB)[0];                \
 492       (DST)[1] += S * (SRCB)[1];                \
 493 } while (0)
 494
 495 /** In-place scalar multiplication */
 496 #define SELF_SCALE_SCALAR_2V( DST, S )          \
 497 do {                                            \
 498       (DST)[0] *= S;                            \
 499       (DST)[1] *= S;                            \
 500 } while (0)
 501
 502 /** In-place scalar addition */
 503 #define ACC_SCALAR_2V( DST, S )                 \
 504 do {                                            \
 505       (DST)[0] += S;                            \
 506       (DST)[1] += S;                            \
 507 } while (0)
 508
 509
 510
 511 /**
 512  * Linear interpolation
 513  *
 514  * \note \p OUT argument is evaluated twice!
 515  * \note Be wary of using *coord++ as an argument to any of these macros!
 516  */
 517 #define LINTERP(T, OUT, IN)     ((OUT) + (T) * ((IN) - (OUT)))
 518
 519 /* Can do better with integer math
 520  */
 521 #define INTERP_UB( t, dstub, outub, inub )      \
 522 do {                                            \
 523    GLfloat inf = UBYTE_TO_FLOAT( inub );        \
 524    GLfloat outf = UBYTE_TO_FLOAT( outub );      \
 525    GLfloat dstf = LINTERP( t, outf, inf );      \
 526    UNCLAMPED_FLOAT_TO_UBYTE( dstub, dstf );     \
 527 } while (0)
 528
 529 #define INTERP_CHAN( t, dstc, outc, inc )       \
 530 do {                                            \
 531    GLfloat inf = CHAN_TO_FLOAT( inc );          \
 532    GLfloat outf = CHAN_TO_FLOAT( outc );        \
 533    GLfloat dstf = LINTERP( t, outf, inf );      \
 534    UNCLAMPED_FLOAT_TO_CHAN( dstc, dstf );       \
 535 } while (0)
 536
 537 #define INTERP_UI( t, dstui, outui, inui )      \
 538    dstui = (GLuint) (GLint) LINTERP( (t), (GLfloat) (outui), (GLfloat) (inui) )
 539
 540 #define INTERP_F( t, dstf, outf, inf )          \
 541    dstf = LINTERP( t, outf, inf )
 542
 543 #define INTERP_4F( t, dst, out, in )            \
 544 do {                                            \
 545    dst[0] = LINTERP( (t), (out)[0], (in)[0] );  \
 546    dst[1] = LINTERP( (t), (out)[1], (in)[1] );  \
 547    dst[2] = LINTERP( (t), (out)[2], (in)[2] );  \
 548    dst[3] = LINTERP( (t), (out)[3], (in)[3] );  \
 549 } while (0)
 550
 551 #define INTERP_3F( t, dst, out, in )            \
 552 do {                                            \
 553    dst[0] = LINTERP( (t), (out)[0], (in)[0] );  \
 554    dst[1] = LINTERP( (t), (out)[1], (in)[1] );  \
 555    dst[2] = LINTERP( (t), (out)[2], (in)[2] );  \
 556 } while (0)
 557
 558 #define INTERP_4CHAN( t, dst, out, in )                 \
 559 do {                                                    \
 560    INTERP_CHAN( (t), (dst)[0], (out)[0], (in)[0] );     \
 561    INTERP_CHAN( (t), (dst)[1], (out)[1], (in)[1] );     \
 562    INTERP_CHAN( (t), (dst)[2], (out)[2], (in)[2] );     \
 563    INTERP_CHAN( (t), (dst)[3], (out)[3], (in)[3] );     \
 564 } while (0)
 565
 566 #define INTERP_3CHAN( t, dst, out, in )                 \
 567 do {                                                    \
 568    INTERP_CHAN( (t), (dst)[0], (out)[0], (in)[0] );     \
 569    INTERP_CHAN( (t), (dst)[1], (out)[1], (in)[1] );     \
 570    INTERP_CHAN( (t), (dst)[2], (out)[2], (in)[2] );     \
 571 } while (0)
 572
 573 #define INTERP_SZ( t, vec, to, out, in, sz )                            \
 574 do {                                                                    \
 575    switch (sz) {                                                        \
 576    case 4: vec[to][3] = LINTERP( (t), (vec)[out][3], (vec)[in][3] );    \
 577    case 3: vec[to][2] = LINTERP( (t), (vec)[out][2], (vec)[in][2] );    \
 578    case 2: vec[to][1] = LINTERP( (t), (vec)[out][1], (vec)[in][1] );    \
 579    case 1: vec[to][0] = LINTERP( (t), (vec)[out][0], (vec)[in][0] );    \
 580    }                                                                    \
 581 } while(0)
 582
 583
 584
 585 /** Assign scalers to short vectors */
 586 #define ASSIGN_2V( V, V0, V1 )  \
 587 do {                            \
 588     V[0] = V0;                  \
 589     V[1] = V1;                  \
 590 } while(0)
 591
 592 /*@}*/
 593
 594
 595
 596 /** Clamp X to [MIN,MAX] */
 597 #define CLAMP( X, MIN, MAX )  ( (X)<(MIN) ? (MIN) : ((X)>(MAX) ? (MAX) : (X)) )
 598
 599 /** Assign X to CLAMP(X, MIN, MAX) */
 600 #define CLAMP_SELF(x, mn, mx)  \
 601    ( (x)<(mn) ? ((x) = (mn)) : ((x)>(mx) ? ((x)=(mx)) : (x)) )
 602
 603
 604
 605 /** Minimum of two values: */
 606 #define MIN2( A, B )   ( (A)<(B) ? (A) : (B) )
 607
 608 /** Maximum of two values: */
 609 #define MAX2( A, B )   ( (A)>(B) ? (A) : (B) )
 610
 611 /** Dot product of two 2-element vectors */
 612 #define DOT2( a, b )  ( (a)[0]*(b)[0] + (a)[1]*(b)[1] )
 613
 614 /** Dot product of two 3-element vectors */
 615 #define DOT3( a, b )  ( (a)[0]*(b)[0] + (a)[1]*(b)[1] + (a)[2]*(b)[2] )
 616
 617 /** Dot product of two 4-element vectors */
 618 #define DOT4( a, b )  ( (a)[0]*(b)[0] + (a)[1]*(b)[1] + \
 619                         (a)[2]*(b)[2] + (a)[3]*(b)[3] )
 620
 621 /** Dot product of two 4-element vectors */
 622 #define DOT4V(v,a,b,c,d) (v[0]*(a) + v[1]*(b) + v[2]*(c) + v[3]*(d))
 623
 624
 625 /** Cross product of two 3-element vectors */
 626 #define CROSS3(n, u, v)                         \
 627 do {                                            \
 628    (n)[0] = (u)[1]*(v)[2] - (u)[2]*(v)[1];      \
 629    (n)[1] = (u)[2]*(v)[0] - (u)[0]*(v)[2];      \
 630    (n)[2] = (u)[0]*(v)[1] - (u)[1]*(v)[0];      \
 631 } while (0)
 632
 633
 634 /* Normalize a 3-element vector to unit length. */
 635 #define NORMALIZE_3FV( V )                      \
 636 do {                                            \
 637    GLfloat len = (GLfloat) LEN_SQUARED_3FV(V);  \
 638    if (len) {                                   \
 639       len = INV_SQRTF(len);                     \
 640       (V)[0] = (GLfloat) ((V)[0] * len);        \
 641       (V)[1] = (GLfloat) ((V)[1] * len);        \
 642       (V)[2] = (GLfloat) ((V)[2] * len);        \
 643    }                                            \
 644 } while(0)
 645
 646 #define LEN_3FV( V ) (SQRTF((V)[0]*(V)[0]+(V)[1]*(V)[1]+(V)[2]*(V)[2]))
 647 #define LEN_2FV( V ) (SQRTF((V)[0]*(V)[0]+(V)[1]*(V)[1]))
 648
 649 #define LEN_SQUARED_3FV( V ) ((V)[0]*(V)[0]+(V)[1]*(V)[1]+(V)[2]*(V)[2])
 650 #define LEN_SQUARED_2FV( V ) ((V)[0]*(V)[0]+(V)[1]*(V)[1])
 651
 652
 653 /*@}*/
 654
 655
 656 #endif