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 #ifndef __QSQRT_H__
00029 #define __QSQRT_H__
00030
00031 #if (!defined (CS_NO_QSQRT)) && defined (PROC_X86) && defined (COMP_GCC)
00032
00033
00034
00035
00036
00037
00038
00039
00040
00049 static inline float qsqrt (float x)
00050 {
00051 float ret;
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062 __asm__ (
00063 "flds %1\n"
00064 "movl $0xbe6f0000,%%eax\n"
00065 "subl %1,%%eax\n"
00066 "shrl $1,%%eax\n"
00067 "movl %%eax,%1\n"
00068 "flds %2\n"
00069 "fmul %%st(1)\n"
00070 "flds %3\n"
00071 "flds %1\n"
00072 "fld %%st\n"
00073 "fmul %%st\n"
00074 "fmul %%st(3)\n"
00075 "fsubr %%st(2)\n"
00076 "fmulp %%st(1)\n"
00077 "fld %%st\n"
00078 "fmul %%st\n"
00079 "fmulp %%st(3)\n"
00080 "fxch\n"
00081 "fsubp %%st,%%st(2)\n"
00082 "fmulp %%st(1)\n"
00083 "fmulp %%st(1)\n"
00084 : "=&t" (ret), "+m" (x) : "m" (0.5F), "m" (1.5F)
00085 : "eax", "st(1)", "st(2)", "st(3)", "st(4)", "st(5)", "st(6)", "st(7)"
00086 );
00087 return ret;
00088 }
00089
00097 static inline float qisqrt (float x)
00098 {
00099 float ret;
00100 __asm__ (
00101 "flds %1\n"
00102 "movl $0xbe6f0000,%%eax\n"
00103 "subl %1,%%eax\n"
00104 "shrl $1,%%eax\n"
00105 "movl %%eax,%1\n"
00106 "flds %2\n"
00107 "fmulp %%st(1)\n"
00108 "flds %3\n"
00109 "flds %1\n"
00110 "fld %%st\n"
00111 "fmul %%st\n"
00112 "fmul %%st(3)\n"
00113 "fsubr %%st(2)\n"
00114 "fmulp %%st(1)\n"
00115 "fld %%st\n"
00116 "fmul %%st\n"
00117 "fmulp %%st(3)\n"
00118 "fxch\n"
00119 "fsubp %%st,%%st(2)\n"
00120 "fmulp %%st(1)\n"
00121 : "=t" (ret), "+m" (x) : "m" (0.5F), "m" (1.5F)
00122 : "eax", "st(1)", "st(2)", "st(3)", "st(4)", "st(5)", "st(6)", "st(7)"
00123 );
00124 return ret;
00125 }
00126
00127 #elif (!defined (CS_NO_QSQRT)) && defined (PROC_POWERPC) && defined (COMP_GCC)
00128
00136 static inline float qsqrt(float x)
00137 {
00138 float y0 = 0.0;
00139
00140 if (x != 0.0)
00141 {
00142 float x0 = x * 0.5f;
00143
00144 asm ("frsqrte %0,%1" : "=f" (y0) : "f" (x));
00145
00146 y0 = y0 * (1.5f - x0 * y0 * y0);
00147 y0 = (y0 * (1.5f - x0 * y0 * y0)) * x;
00148 };
00149
00150 return y0;
00151 };
00152
00157 static inline float qisqrt(float x)
00158 {
00159 float x0 = x * 0.5f;
00160 float y0;
00161 asm ("frsqrte %0,%1" : "=f" (y0) : "f" (x));
00162
00163 y0 = y0 * (1.5f - x0 * y0 * y0);
00164 y0 = y0 * (1.5f - x0 * y0 * y0);
00165
00166 return y0;
00167 };
00168
00169 #elif (!defined (CS_NO_QSQRT)) && defined (PROC_X86) && defined (COMP_VC)
00170
00171 #include <math.h>
00172
00173 #define qsqrt(x) sqrtf(x)
00174 #define qisqrt(x) (1.0f / sqrtf(x))
00175
00176 #else
00177
00178 #include <math.h>
00179 #define qsqrt(x) sqrt(x)
00180 #define qisqrt(x) (1.0/sqrt(x))
00181
00182 #endif
00183
00184 #endif // __QSQRT_H__