~ubuntu-branches/ubuntu/trusty/bmagic/trusty-proposed

<h1>bmutil.h</h1><a href="a00094.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 /*

<a name="l00003"></a>00003

<a name="l00004"></a>00004 Permission is hereby granted, free of charge, to any person

<a name="l00005"></a>00005 obtaining a copy of this software and associated documentation

<a name="l00006"></a>00006 files (the "Software"), to deal in the Software without restriction,

<a name="l00007"></a>00007 including without limitation the rights to use, copy, modify, merge,

<a name="l00008"></a>00008 publish, distribute, sublicense, and/or sell copies of the Software,

<a name="l00009"></a>00009 and to permit persons to whom the Software is furnished to do so,

<a name="l00010"></a>00010 subject to the following conditions:

<a name="l00011"></a>00011

<a name="l00012"></a>00012 The above copyright notice and this permission notice shall be included

<a name="l00013"></a>00013 in all copies or substantial portions of the Software.

<a name="l00014"></a>00014

<a name="l00015"></a>00015 THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,

<a name="l00016"></a>00016 EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES

<a name="l00017"></a>00017 OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.

<a name="l00018"></a>00018 IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM,

<a name="l00019"></a>00019 DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,

<a name="l00020"></a>00020 ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR

<a name="l00021"></a>00021 OTHER DEALINGS IN THE SOFTWARE.

<a name="l00022"></a>00022

<a name="l00023"></a>00023 For more information please visit: http://bmagic.sourceforge.net

<a name="l00024"></a>00024

<a name="l00025"></a>00025 */

<a name="l00026"></a>00026 #ifndef BMUTIL__H__INCLUDED__

<a name="l00027"></a>00027 #define BMUTIL__H__INCLUDED__

<a name="l00028"></a>00028

<a name="l00029"></a>00029 #include "<a class="code" href="a00087.html">bmdef.h</a>"

<a name="l00030"></a>00030 #include "<a class="code" href="a00086.html">bmconst.h</a>"

<a name="l00031"></a>00031

<a name="l00032"></a>00032 namespace bm

<a name="l00033"></a>00033 {

<a name="l00034"></a>00034

<a name="l00035"></a>00035

<a name="l00036"></a>00036 /*

<a name="l00037"></a>00037 // SSE2 vector equality comparison

<a name="l00038"></a>00038 int _mm_any_eq( vFloat a, vFloat b )

<a name="l00039"></a>00039 {

<a name="l00040"></a>00040

<a name="l00041"></a>00041 //test a==b for each float in a & b

<a name="l00042"></a>00042 vFloat mask = _mm_cmpeq_ps( a, b );

<a name="l00043"></a>00043

<a name="l00044"></a>00044 //copy top bit of each result to maskbits

<a name="l00045"></a>00045 int maskBits = _mm_movemask_ps( mask );

<a name="l00046"></a>00046 return maskBits != 0;

<a name="l00047"></a>00047 }

<a name="l00048"></a>00048 */

<a name="l00049"></a>00049

<a name="l00050"></a>00050 /**

<a name="l00051"></a>00051 Fast loop-less function to find LOG2

<a name="l00052"></a>00052 */

<a name="l00053"></a>00053 template<typename T>

<a name="l00054"></a><a class="code" href="a00106.html#46f258f92ab40e2a0c76aa775338faaa">00054</a> T <a class="code" href="a00106.html#46f258f92ab40e2a0c76aa775338faaa" title="Fast loop-less function to find LOG2.">ilog2</a>(T x)

<a name="l00055"></a>00055 {

<a name="l00056"></a>00056 unsigned int l = 0;

<a name="l00057"></a>00057 if(x >= 1<<16) { x>>=16; l |= 16; }

<a name="l00058"></a>00058 if(x >= 1<<8) { x>>=8; l |= 8; }

<a name="l00059"></a>00059 if(x >= 1<<4) { x>>=4; l |= 4; }

<a name="l00060"></a>00060 if(x >= 1<<2) { x>>=2; l |= 2; }

<a name="l00061"></a>00061 if(x >= 1<<1) l |=1;

<a name="l00062"></a>00062 return l;

<a name="l00063"></a>00063 }

<a name="l00064"></a>00064

<a name="l00065"></a>00065 template<>

<a name="l00066"></a><a class="code" href="a00106.html#5f818ac575622da8a9a75511e649ea90">00066</a> inline <a class="code" href="a00106.html#c654d6319039a86546d235a236fc7cf6">bm::gap_word_t</a> <a class="code" href="a00106.html#46f258f92ab40e2a0c76aa775338faaa" title="Fast loop-less function to find LOG2.">ilog2</a>(<a class="code" href="a00106.html#c654d6319039a86546d235a236fc7cf6">gap_word_t</a> x)

<a name="l00067"></a>00067 {

<a name="l00068"></a>00068 unsigned int l = 0;

<a name="l00069"></a>00069 if(x >= 1<<8) { x>>=8; l |= 8; }

<a name="l00070"></a>00070 if(x >= 1<<4) { x>>=4; l |= 4; }

<a name="l00071"></a>00071 if(x >= 1<<2) { x>>=2; l |= 2; }

<a name="l00072"></a>00072 if(x >= 1<<1) l |=1;

<a name="l00073"></a>00073 return l;

<a name="l00074"></a>00074 }

<a name="l00075"></a>00075

<a name="l00076"></a>00076

<a name="l00077"></a>00077 /**

<a name="l00078"></a>00078 Lookup table based integer LOG2

<a name="l00079"></a>00079 */

<a name="l00080"></a>00080 template<typename T>

<a name="l00081"></a><a class="code" href="a00106.html#3ab126a6fd7ffd872a4001638819a330">00081</a> T <a class="code" href="a00106.html#3ab126a6fd7ffd872a4001638819a330" title="Lookup table based integer LOG2.">ilog2_LUT</a>(T x)

100

<a name="l00082"></a>00082 {

101

<a name="l00083"></a>00083 unsigned l = 0;

102

<a name="l00084"></a>00084 if (x & 0xffff0000)

103

<a name="l00085"></a>00085 {

104

<a name="l00086"></a>00086 l += 16; x >>= 16;

105

<a name="l00087"></a>00087 }

106

<a name="l00088"></a>00088

107

<a name="l00089"></a>00089 if (x & 0xff00)

108

<a name="l00090"></a>00090 {

109

<a name="l00091"></a>00091 l += 8; x >>= 8;

110

<a name="l00092"></a>00092 }

111

<a name="l00093"></a>00093 return l + <a class="code" href="a00064.html" title="Structure keeps index of first ON bit for every byte.">first_bit_table<true>::_idx</a>[x];

112

<a name="l00094"></a>00094 }

113

<a name="l00095"></a>00095

114

<a name="l00096"></a>00096 /**

115

<a name="l00097"></a>00097 Lookup table based short integer LOG2

116

<a name="l00098"></a>00098 */

117

<a name="l00099"></a>00099 template<>

118

<a name="l00100"></a>00100 inline <a class="code" href="a00106.html#c654d6319039a86546d235a236fc7cf6">bm::gap_word_t</a> ilog2_LUT<bm::gap_word_t>(<a class="code" href="a00106.html#c654d6319039a86546d235a236fc7cf6">bm::gap_word_t</a> x)

119

<a name="l00101"></a>00101 {

120

<a name="l00102"></a>00102 unsigned l = 0;

121

<a name="l00103"></a>00103 if (x & 0xff00)

122

<a name="l00104"></a>00104 {

123

<a name="l00105"></a>00105 l += 8; x >>= 8;

124

<a name="l00106"></a>00106 }

125

<a name="l00107"></a>00107 return l + first_bit_table<true>::_idx[x];

126

<a name="l00108"></a>00108 }

127

<a name="l00109"></a>00109

128

<a name="l00110"></a>00110

129

<a name="l00111"></a>00111 // if we are running on x86 CPU we can use inline ASM

130

<a name="l00112"></a>00112

131

<a name="l00113"></a>00113 #ifdef BM_x86

132

<a name="l00114"></a>00114 #ifdef __GNUG__

133

<a name="l00115"></a>00115

134

<a name="l00116"></a>00116 inline

135

<a name="l00117"></a>00117 unsigned bsf_asm32(unsigned int v)

136

<a name="l00118"></a>00118 {

137

<a name="l00119"></a>00119 unsigned r;

138

<a name="l00120"></a>00120 asm volatile(" bsfl %1, %0": "=r"(r): "rm"(v) );

139

<a name="l00121"></a>00121 return r;

140

<a name="l00122"></a>00122 }

141

<a name="l00123"></a>00123

142

<a name="l00124"></a>00124 <a class="code" href="a00087.html#938734d014fb68dd8b2251fe8ec2b025">BMFORCEINLINE</a> unsigned bsr_asm32(register unsigned int v)

143

<a name="l00125"></a>00125 {

144

<a name="l00126"></a>00126 unsigned r;

145

<a name="l00127"></a>00127 asm volatile(" bsrl %1, %0": "=r"(r): "rm"(v) );

146

<a name="l00128"></a>00128 return r;

147

<a name="l00129"></a>00129 }

148

<a name="l00130"></a>00130

149

<a name="l00131"></a>00131 #endif // __GNUG__

150

<a name="l00132"></a>00132

151

<a name="l00133"></a>00133 #ifdef _MSC_VER

152

<a name="l00134"></a>00134 <a class="code" href="a00087.html#938734d014fb68dd8b2251fe8ec2b025">BMFORCEINLINE</a> unsigned int bsr_asm32(register unsigned int value)

153

<a name="l00135"></a>00135 {

154

<a name="l00136"></a>00136 __asm bsr eax, value

155

<a name="l00137"></a>00137 }

156

<a name="l00138"></a>00138 #endif

157

<a name="l00139"></a>00139

158

<a name="l00140"></a>00140 #endif // BM_x86

159

<a name="l00141"></a>00141

160

<a name="l00142"></a>00142

161

<a name="l00143"></a>00143

162

<a name="l00144"></a>00144 } // bm

163

<a name="l00145"></a>00145

164

<a name="l00146"></a>00146 #endif

165

</pre></div></div>

166

<hr size="1"><address style="text-align: right;">Generated on Mon Oct 12 12:16:39 2009 for BitMagic by

167

168

169

</body>

170

</html>

Older »