a00064_source.html 6.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104
  1. <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
  2. <html xmlns="http://www.w3.org/1999/xhtml">
  3. <head>
  4. <meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
  5. <title>intrinsic_exponential.hpp Source File</title>
  6. <link href="tabs.css" rel="stylesheet" type="text/css"/>
  7. <link href="doxygen.css" rel="stylesheet" type="text/css"/>
  8. </head>
  9. <body>
  10. <!-- Generated by Doxygen 1.7.4 -->
  11. <div id="top">
  12. <div id="titlearea">
  13. <table cellspacing="0" cellpadding="0">
  14. <tbody>
  15. <tr style="height: 56px;">
  16. <td id="projectlogo"><img alt="Logo" src="logo-mini.png"/></td>
  17. </tr>
  18. </tbody>
  19. </table>
  20. </div>
  21. <div id="navrow1" class="tabs">
  22. <ul class="tablist">
  23. <li><a href="index.html"><span>Main&#160;Page</span></a></li>
  24. <li><a href="modules.html"><span>Modules</span></a></li>
  25. <li><a href="namespaces.html"><span>Namespaces</span></a></li>
  26. <li><a href="annotated.html"><span>Classes</span></a></li>
  27. <li class="current"><a href="files.html"><span>Files</span></a></li>
  28. </ul>
  29. </div>
  30. <div id="navrow2" class="tabs2">
  31. <ul class="tablist">
  32. <li><a href="files.html"><span>File&#160;List</span></a></li>
  33. </ul>
  34. </div>
  35. <div class="header">
  36. <div class="headertitle">
  37. <div class="title">intrinsic_exponential.hpp</div> </div>
  38. </div>
  39. <div class="contents">
  40. <div class="fragment"><pre class="fragment"><a name="l00001"></a>00001
  41. <a name="l00002"></a>00002 <span class="comment">// OpenGL Mathematics Copyright (c) 2005 - 2011 G-Truc Creation (www.g-truc.net)</span>
  42. <a name="l00004"></a>00004 <span class="comment"></span><span class="comment">// Created : 2009-05-11</span>
  43. <a name="l00005"></a>00005 <span class="comment">// Updated : 2009-05-11</span>
  44. <a name="l00006"></a>00006 <span class="comment">// Licence : This source is under MIT License</span>
  45. <a name="l00007"></a>00007 <span class="comment">// File : glm/core/intrinsic_exponential.hpp</span>
  46. <a name="l00009"></a>00009 <span class="comment"></span>
  47. <a name="l00010"></a>00010 <span class="preprocessor">#ifndef glm_detail_intrinsic_exponential</span>
  48. <a name="l00011"></a>00011 <span class="preprocessor"></span><span class="preprocessor">#define glm_detail_intrinsic_exponential</span>
  49. <a name="l00012"></a>00012 <span class="preprocessor"></span>
  50. <a name="l00013"></a>00013 <span class="preprocessor">#include &quot;setup.hpp&quot;</span>
  51. <a name="l00014"></a>00014
  52. <a name="l00015"></a>00015 <span class="preprocessor">#if((GLM_ARCH &amp; GLM_ARCH_SSE2) != GLM_ARCH_SSE2)</span>
  53. <a name="l00016"></a>00016 <span class="preprocessor"></span><span class="preprocessor"># error &quot;SSE2 instructions not supported or enabled&quot;</span>
  54. <a name="l00017"></a>00017 <span class="preprocessor"></span><span class="preprocessor">#else</span>
  55. <a name="l00018"></a>00018 <span class="preprocessor"></span>
  56. <a name="l00019"></a>00019 <span class="keyword">namespace </span>glm{
  57. <a name="l00020"></a>00020 <span class="keyword">namespace </span>detail
  58. <a name="l00021"></a>00021 {
  59. <a name="l00022"></a>00022 <span class="comment">/*</span>
  60. <a name="l00023"></a>00023 <span class="comment">GLM_FUNC_QUALIFIER __m128 sse_rsqrt_nr_ss(__m128 const x)</span>
  61. <a name="l00024"></a>00024 <span class="comment">{</span>
  62. <a name="l00025"></a>00025 <span class="comment"> __m128 recip = _mm_rsqrt_ss( x ); // &quot;estimate&quot; opcode</span>
  63. <a name="l00026"></a>00026 <span class="comment"> const static __m128 three = { 3, 3, 3, 3 }; // aligned consts for fast load</span>
  64. <a name="l00027"></a>00027 <span class="comment"> const static __m128 half = { 0.5,0.5,0.5,0.5 };</span>
  65. <a name="l00028"></a>00028 <span class="comment"> __m128 halfrecip = _mm_mul_ss( half, recip );</span>
  66. <a name="l00029"></a>00029 <span class="comment"> __m128 threeminus_xrr = _mm_sub_ss( three, _mm_mul_ss( x, _mm_mul_ss ( recip, recip ) ) );</span>
  67. <a name="l00030"></a>00030 <span class="comment"> return _mm_mul_ss( halfrecip, threeminus_xrr );</span>
  68. <a name="l00031"></a>00031 <span class="comment">}</span>
  69. <a name="l00032"></a>00032 <span class="comment"> </span>
  70. <a name="l00033"></a>00033 <span class="comment">GLM_FUNC_QUALIFIER __m128 sse_normalize_fast_ps( float * RESTRICT vOut, float * RESTRICT vIn )</span>
  71. <a name="l00034"></a>00034 <span class="comment">{</span>
  72. <a name="l00035"></a>00035 <span class="comment"> __m128 x = _mm_load_ss(&amp;vIn[0]);</span>
  73. <a name="l00036"></a>00036 <span class="comment"> __m128 y = _mm_load_ss(&amp;vIn[1]);</span>
  74. <a name="l00037"></a>00037 <span class="comment"> __m128 z = _mm_load_ss(&amp;vIn[2]);</span>
  75. <a name="l00038"></a>00038 <span class="comment"> </span>
  76. <a name="l00039"></a>00039 <span class="comment"> const __m128 l = // compute x*x + y*y + z*z</span>
  77. <a name="l00040"></a>00040 <span class="comment"> _mm_add_ss(</span>
  78. <a name="l00041"></a>00041 <span class="comment"> _mm_add_ss( _mm_mul_ss(x,x),</span>
  79. <a name="l00042"></a>00042 <span class="comment"> _mm_mul_ss(y,y)</span>
  80. <a name="l00043"></a>00043 <span class="comment"> ),</span>
  81. <a name="l00044"></a>00044 <span class="comment"> _mm_mul_ss( z, z )</span>
  82. <a name="l00045"></a>00045 <span class="comment"> );</span>
  83. <a name="l00046"></a>00046 <span class="comment"> </span>
  84. <a name="l00047"></a>00047 <span class="comment"> </span>
  85. <a name="l00048"></a>00048 <span class="comment"> const __m128 rsqt = _mm_rsqrt_nr_ss( l );</span>
  86. <a name="l00049"></a>00049 <span class="comment"> _mm_store_ss( &amp;vOut[0] , _mm_mul_ss( rsqt, x ) );</span>
  87. <a name="l00050"></a>00050 <span class="comment"> _mm_store_ss( &amp;vOut[1] , _mm_mul_ss( rsqt, y ) );</span>
  88. <a name="l00051"></a>00051 <span class="comment"> _mm_store_ss( &amp;vOut[2] , _mm_mul_ss( rsqt, z ) );</span>
  89. <a name="l00052"></a>00052 <span class="comment"> </span>
  90. <a name="l00053"></a>00053 <span class="comment"> return _mm_mul_ss( l , rsqt );</span>
  91. <a name="l00054"></a>00054 <span class="comment">}</span>
  92. <a name="l00055"></a>00055 <span class="comment">*/</span>
  93. <a name="l00056"></a>00056 }<span class="comment">//namespace detail</span>
  94. <a name="l00057"></a>00057 }<span class="comment">//namespace glm</span>
  95. <a name="l00058"></a>00058
  96. <a name="l00059"></a>00059 <span class="preprocessor">#endif//GLM_ARCH</span>
  97. <a name="l00060"></a>00060 <span class="preprocessor"></span><span class="preprocessor">#endif//glm_detail_intrinsic_exponential</span>
  98. </pre></div></div>
  99. </div>
  100. <hr class="footer"/><address class="footer"><small>Generated by&#160;
  101. <a href="http://www.doxygen.org/index.html">
  102. <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.7.4 </small></address>
  103. </body>
  104. </html>