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 #ifndef __QSQRT_H__
00028 #define __QSQRT_H__
00029
00030 #if (!defined (CS_NO_QSQRT)) && defined (PROC_X86) && defined (COMP_GCC)
00031
00032
00033
00034
00035
00036
00037
00038
00039
00048 static inline float qsqrt (float x)
00049 {
00050 float ret;
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061 __asm__ (
00062 "flds %1\n"
00063 "movl $0xbe6f0000,%%eax\n"
00064 "subl %1,%%eax\n"
00065 "shrl $1,%%eax\n"
00066 "movl %%eax,%1\n"
00067 "flds %2\n"
00068 "fmul %%st(1)\n"
00069 "flds %3\n"
00070 "flds %1\n"
00071 "fld %%st\n"
00072 "fmul %%st\n"
00073 "fmul %%st(3)\n"
00074 "fsubr %%st(2)\n"
00075 "fmulp %%st(1)\n"
00076 "fld %%st\n"
00077 "fmul %%st\n"
00078 "fmulp %%st(3)\n"
00079 "fxch\n"
00080 "fsubp %%st,%%st(2)\n"
00081 "fmulp %%st(1)\n"
00082 "fmulp %%st(1)\n"
00083 : "=&t" (ret), "+m" (x) : "m" (0.5F), "m" (1.5F)
00084 : "eax", "st", "st(1)", "st(2)", "st(3)", "st(4)", "st(5)", "st(6)", "st(7)"
00085 );
00086 return ret;
00087 }
00088
00096 static inline float qisqrt (float x)
00097 {
00098 float ret;
00099 __asm__ (
00100 "flds %1\n"
00101 "movl $0xbe6f0000,%%eax\n"
00102 "subl %1,%%eax\n"
00103 "shrl $1,%%eax\n"
00104 "movl %%eax,%1\n"
00105 "flds %2\n"
00106 "fmulp %%st(1)\n"
00107 "flds %3\n"
00108 "flds %1\n"
00109 "fld %%st\n"
00110 "fmul %%st\n"
00111 "fmul %%st(3)\n"
00112 "fsubr %%st(2)\n"
00113 "fmulp %%st(1)\n"
00114 "fld %%st\n"
00115 "fmul %%st\n"
00116 "fmulp %%st(3)\n"
00117 "fxch\n"
00118 "fsubp %%st,%%st(2)\n"
00119 "fmulp %%st(1)\n"
00120 : "=t" (ret), "+m" (x) : "m" (0.5F), "m" (1.5F)
00121 : "eax", "st", "st(1)", "st(2)", "st(3)", "st(4)", "st(5)", "st(6)", "st(7)"
00122 );
00123 return ret;
00124 }
00125
00126 #else
00127
00128 #include <math.h>
00129 #define qsqrt(x) sqrt(x)
00130 #define qisqrt(x) (1.0/sqrt(x))
00131
00132 #endif
00133
00134 #endif // __QSQRT_H__