00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044 #ifdef HAVE_CONFIG_H
00045 #include <config.h>
00046 #endif
00047
00048 #include "prim_type.h"
00049 #include "fixpoint.h"
00050 #include "fe.h"
00051 #include "fe_internal.h"
00052
00053
00054
00055 static uint32 logtable[] = {
00056 (uint32) (0 * (1 << DEFAULT_RADIX)),
00057 (uint32) (0.0223678130284545 * (1 << DEFAULT_RADIX)),
00058 (uint32) (0.0443941193584534 * (1 << DEFAULT_RADIX)),
00059 (uint32) (0.0660891904577724 * (1 << DEFAULT_RADIX)),
00060 (uint32) (0.0874628412503394 * (1 << DEFAULT_RADIX)),
00061 (uint32) (0.108524456778169 * (1 << DEFAULT_RADIX)),
00062 (uint32) (0.129283016944966 * (1 << DEFAULT_RADIX)),
00063 (uint32) (0.149747119504682 * (1 << DEFAULT_RADIX)),
00064 (uint32) (0.169925001442312 * (1 << DEFAULT_RADIX)),
00065 (uint32) (0.189824558880017 * (1 << DEFAULT_RADIX)),
00066 (uint32) (0.20945336562895 * (1 << DEFAULT_RADIX)),
00067 (uint32) (0.228818690495881 * (1 << DEFAULT_RADIX)),
00068 (uint32) (0.247927513443586 * (1 << DEFAULT_RADIX)),
00069 (uint32) (0.266786540694901 * (1 << DEFAULT_RADIX)),
00070 (uint32) (0.285402218862248 * (1 << DEFAULT_RADIX)),
00071 (uint32) (0.303780748177103 * (1 << DEFAULT_RADIX)),
00072 (uint32) (0.321928094887362 * (1 << DEFAULT_RADIX)),
00073 (uint32) (0.339850002884625 * (1 << DEFAULT_RADIX)),
00074 (uint32) (0.357552004618084 * (1 << DEFAULT_RADIX)),
00075 (uint32) (0.375039431346925 * (1 << DEFAULT_RADIX)),
00076 (uint32) (0.39231742277876 * (1 << DEFAULT_RADIX)),
00077 (uint32) (0.409390936137702 * (1 << DEFAULT_RADIX)),
00078 (uint32) (0.426264754702098 * (1 << DEFAULT_RADIX)),
00079 (uint32) (0.442943495848728 * (1 << DEFAULT_RADIX)),
00080 (uint32) (0.459431618637297 * (1 << DEFAULT_RADIX)),
00081 (uint32) (0.475733430966398 * (1 << DEFAULT_RADIX)),
00082 (uint32) (0.491853096329675 * (1 << DEFAULT_RADIX)),
00083 (uint32) (0.507794640198696 * (1 << DEFAULT_RADIX)),
00084 (uint32) (0.523561956057013 * (1 << DEFAULT_RADIX)),
00085 (uint32) (0.539158811108031 * (1 << DEFAULT_RADIX)),
00086 (uint32) (0.554588851677637 * (1 << DEFAULT_RADIX)),
00087 (uint32) (0.569855608330948 * (1 << DEFAULT_RADIX)),
00088 (uint32) (0.584962500721156 * (1 << DEFAULT_RADIX)),
00089 (uint32) (0.599912842187128 * (1 << DEFAULT_RADIX)),
00090 (uint32) (0.614709844115208 * (1 << DEFAULT_RADIX)),
00091 (uint32) (0.62935662007961 * (1 << DEFAULT_RADIX)),
00092 (uint32) (0.643856189774725 * (1 << DEFAULT_RADIX)),
00093 (uint32) (0.658211482751795 * (1 << DEFAULT_RADIX)),
00094 (uint32) (0.672425341971496 * (1 << DEFAULT_RADIX)),
00095 (uint32) (0.686500527183218 * (1 << DEFAULT_RADIX)),
00096 (uint32) (0.700439718141092 * (1 << DEFAULT_RADIX)),
00097 (uint32) (0.714245517666123 * (1 << DEFAULT_RADIX)),
00098 (uint32) (0.727920454563199 * (1 << DEFAULT_RADIX)),
00099 (uint32) (0.741466986401147 * (1 << DEFAULT_RADIX)),
00100 (uint32) (0.754887502163469 * (1 << DEFAULT_RADIX)),
00101 (uint32) (0.768184324776926 * (1 << DEFAULT_RADIX)),
00102 (uint32) (0.78135971352466 * (1 << DEFAULT_RADIX)),
00103 (uint32) (0.794415866350106 * (1 << DEFAULT_RADIX)),
00104 (uint32) (0.807354922057604 * (1 << DEFAULT_RADIX)),
00105 (uint32) (0.820178962415188 * (1 << DEFAULT_RADIX)),
00106 (uint32) (0.832890014164742 * (1 << DEFAULT_RADIX)),
00107 (uint32) (0.845490050944375 * (1 << DEFAULT_RADIX)),
00108 (uint32) (0.857980995127572 * (1 << DEFAULT_RADIX)),
00109 (uint32) (0.870364719583405 * (1 << DEFAULT_RADIX)),
00110 (uint32) (0.882643049361841 * (1 << DEFAULT_RADIX)),
00111 (uint32) (0.894817763307944 * (1 << DEFAULT_RADIX)),
00112 (uint32) (0.906890595608518 * (1 << DEFAULT_RADIX)),
00113 (uint32) (0.918863237274595 * (1 << DEFAULT_RADIX)),
00114 (uint32) (0.930737337562886 * (1 << DEFAULT_RADIX)),
00115 (uint32) (0.94251450533924 * (1 << DEFAULT_RADIX)),
00116 (uint32) (0.954196310386875 * (1 << DEFAULT_RADIX)),
00117 (uint32) (0.965784284662087 * (1 << DEFAULT_RADIX)),
00118 (uint32) (0.977279923499917 * (1 << DEFAULT_RADIX)),
00119 (uint32) (0.988684686772166 * (1 << DEFAULT_RADIX)),
00120 };
00121
00122 int32
00123 fixlog2(uint32 x)
00124 {
00125 uint32 y;
00126
00127 if (x == 0)
00128 return MIN_FIXLOG2;
00129
00130
00131 #ifdef __i386__
00132 asm("bsrl %1, %0\n": "=r"(y): "g"(x):"cc");
00133 x <<= (31 - y);
00134 #elif defined(__ppc__)
00135 asm("cntlzw %0, %1\n": "=r"(y):"r"(x));
00136 x <<= y;
00137 y = 31 - y;
00138 #elif defined(XSCALE)
00139 asm("clz %0, %1\n": "=r"(y):"r"(x));
00140 x <<= y;
00141 y = 31 - y;
00142 #else
00143 for (y = 31; y >= 0; --y) {
00144 if (x & 0x80000000)
00145 break;
00146 x <<= 1;
00147 }
00148 #endif
00149 y <<= DEFAULT_RADIX;
00150
00151 x = (x >> 25) & 0x3f;
00152 return y + logtable[x];
00153 }
00154
00155 int
00156 fixlog(uint32 x)
00157 {
00158 int32 y;
00159
00160 if (x == 0)
00161 return MIN_FIXLOG;
00162
00163 y = fixlog2(x);
00164 return FIXMUL(y, FIXLN_2);
00165 }