workbench/libs/mesa/src/gallium/auxiliary/util/u_math.c

   1 /**************************************************************************
   2  *
   3  * Copyright 2008 Tungsten Graphics, Inc., Cedar Park, Texas.
   4  * All Rights Reserved.
   5  *
   6  * Permission is hereby granted, free of charge, to any person obtaining a
   7  * copy of this software and associated documentation files (the
   8  * "Software"), to deal in the Software without restriction, including
   9  * without limitation the rights to use, copy, modify, merge, publish,
  10  * distribute, sub license, and/or sell copies of the Software, and to
  11  * permit persons to whom the Software is furnished to do so, subject to
  12  * the following conditions:
  13  *
  14  * The above copyright notice and this permission notice (including the
  15  * next paragraph) shall be included in all copies or substantial portions
  16  * of the Software.
  17  *
  18  * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
  19  * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
  20  * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NON-INFRINGEMENT.
  21  * IN NO EVENT SHALL TUNGSTEN GRAPHICS AND/OR ITS SUPPLIERS BE LIABLE FOR
  22  * ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
  23  * TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE
  24  * SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
  25  *
  26  **************************************************************************/
  27
  28
  29
  30 #include "util/u_math.h"
  31
  32
  33 /** 2^x, for x in [-1.0, 1.0) */
  34 float pow2_table[POW2_TABLE_SIZE];
  35
  36
  37 static void
  38 init_pow2_table(void)
  39 {
  40    int i;
  41    for (i = 0; i < POW2_TABLE_SIZE; i++)
  42       pow2_table[i] = (float) pow(2.0, (i - POW2_TABLE_OFFSET) / POW2_TABLE_SCALE);
  43 }
  44
  45 #if defined(PIPE_OS_AROS)
  46 /*
  47  * NOTE: log_base_2(x) = log(x) / log(2)
  48  * NOTE: 1.442695 = 1/log(2).
  49  */
  50 #define log2(x)  ((float) (log(x) * 1.442695f))
  51 #endif
  52
  53 /** log2(x), for x in [1.0, 2.0) */
  54 float log2_table[LOG2_TABLE_SIZE];
  55
  56
  57 static void
  58 init_log2_table(void)
  59 {
  60    unsigned i;
  61    for (i = 0; i < LOG2_TABLE_SIZE; i++)
  62       log2_table[i] = (float) log2(1.0 + i * (1.0 / LOG2_TABLE_SCALE));
  63 }
  64
  65
  66 /**
  67  * One time init for math utilities.
  68  */
  69 void
  70 util_init_math(void)
  71 {
  72    static boolean initialized = FALSE;
  73    if (!initialized) {
  74       init_pow2_table();
  75       init_log2_table();
  76       initialized = TRUE;
  77    }
  78 }
  79
  80