xmmintrin.h 24 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008
  1. /*===---- xmmintrin.h - SSE intrinsics -------------------------------------===
  2. *
  3. * Permission is hereby granted, free of charge, to any person obtaining a copy
  4. * of this software and associated documentation files (the "Software"), to deal
  5. * in the Software without restriction, including without limitation the rights
  6. * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
  7. * copies of the Software, and to permit persons to whom the Software is
  8. * furnished to do so, subject to the following conditions:
  9. *
  10. * The above copyright notice and this permission notice shall be included in
  11. * all copies or substantial portions of the Software.
  12. *
  13. * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  14. * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
  15. * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
  16. * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
  17. * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
  18. * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
  19. * THE SOFTWARE.
  20. *
  21. *===-----------------------------------------------------------------------===
  22. */
  23. #ifndef __XMMINTRIN_H
  24. #define __XMMINTRIN_H
  25. #ifndef __SSE__
  26. #error "SSE instruction set not enabled"
  27. #else
  28. #include <mmintrin.h>
  29. typedef int __v4si __attribute__((__vector_size__(16)));
  30. typedef float __v4sf __attribute__((__vector_size__(16)));
  31. typedef float __m128 __attribute__((__vector_size__(16)));
  32. /* This header should only be included in a hosted environment as it depends on
  33. * a standard library to provide allocation routines. */
  34. #if __STDC_HOSTED__
  35. #include <mm_malloc.h>
  36. #endif
  37. /* Define the default attributes for the functions in this file. */
  38. #define __DEFAULT_FN_ATTRS __attribute__((__always_inline__, __nodebug__))
  39. static __inline__ __m128 __DEFAULT_FN_ATTRS
  40. _mm_add_ss(__m128 __a, __m128 __b)
  41. {
  42. __a[0] += __b[0];
  43. return __a;
  44. }
  45. static __inline__ __m128 __DEFAULT_FN_ATTRS
  46. _mm_add_ps(__m128 __a, __m128 __b)
  47. {
  48. return __a + __b;
  49. }
  50. static __inline__ __m128 __DEFAULT_FN_ATTRS
  51. _mm_sub_ss(__m128 __a, __m128 __b)
  52. {
  53. __a[0] -= __b[0];
  54. return __a;
  55. }
  56. static __inline__ __m128 __DEFAULT_FN_ATTRS
  57. _mm_sub_ps(__m128 __a, __m128 __b)
  58. {
  59. return __a - __b;
  60. }
  61. static __inline__ __m128 __DEFAULT_FN_ATTRS
  62. _mm_mul_ss(__m128 __a, __m128 __b)
  63. {
  64. __a[0] *= __b[0];
  65. return __a;
  66. }
  67. static __inline__ __m128 __DEFAULT_FN_ATTRS
  68. _mm_mul_ps(__m128 __a, __m128 __b)
  69. {
  70. return __a * __b;
  71. }
  72. static __inline__ __m128 __DEFAULT_FN_ATTRS
  73. _mm_div_ss(__m128 __a, __m128 __b)
  74. {
  75. __a[0] /= __b[0];
  76. return __a;
  77. }
  78. static __inline__ __m128 __DEFAULT_FN_ATTRS
  79. _mm_div_ps(__m128 __a, __m128 __b)
  80. {
  81. return __a / __b;
  82. }
  83. static __inline__ __m128 __DEFAULT_FN_ATTRS
  84. _mm_sqrt_ss(__m128 __a)
  85. {
  86. __m128 __c = __builtin_ia32_sqrtss(__a);
  87. return (__m128) { __c[0], __a[1], __a[2], __a[3] };
  88. }
  89. static __inline__ __m128 __DEFAULT_FN_ATTRS
  90. _mm_sqrt_ps(__m128 __a)
  91. {
  92. return __builtin_ia32_sqrtps(__a);
  93. }
  94. static __inline__ __m128 __DEFAULT_FN_ATTRS
  95. _mm_rcp_ss(__m128 __a)
  96. {
  97. __m128 __c = __builtin_ia32_rcpss(__a);
  98. return (__m128) { __c[0], __a[1], __a[2], __a[3] };
  99. }
  100. static __inline__ __m128 __DEFAULT_FN_ATTRS
  101. _mm_rcp_ps(__m128 __a)
  102. {
  103. return __builtin_ia32_rcpps(__a);
  104. }
  105. static __inline__ __m128 __DEFAULT_FN_ATTRS
  106. _mm_rsqrt_ss(__m128 __a)
  107. {
  108. __m128 __c = __builtin_ia32_rsqrtss(__a);
  109. return (__m128) { __c[0], __a[1], __a[2], __a[3] };
  110. }
  111. static __inline__ __m128 __DEFAULT_FN_ATTRS
  112. _mm_rsqrt_ps(__m128 __a)
  113. {
  114. return __builtin_ia32_rsqrtps(__a);
  115. }
  116. static __inline__ __m128 __DEFAULT_FN_ATTRS
  117. _mm_min_ss(__m128 __a, __m128 __b)
  118. {
  119. return __builtin_ia32_minss(__a, __b);
  120. }
  121. static __inline__ __m128 __DEFAULT_FN_ATTRS
  122. _mm_min_ps(__m128 __a, __m128 __b)
  123. {
  124. return __builtin_ia32_minps(__a, __b);
  125. }
  126. static __inline__ __m128 __DEFAULT_FN_ATTRS
  127. _mm_max_ss(__m128 __a, __m128 __b)
  128. {
  129. return __builtin_ia32_maxss(__a, __b);
  130. }
  131. static __inline__ __m128 __DEFAULT_FN_ATTRS
  132. _mm_max_ps(__m128 __a, __m128 __b)
  133. {
  134. return __builtin_ia32_maxps(__a, __b);
  135. }
  136. static __inline__ __m128 __DEFAULT_FN_ATTRS
  137. _mm_and_ps(__m128 __a, __m128 __b)
  138. {
  139. return (__m128)((__v4si)__a & (__v4si)__b);
  140. }
  141. static __inline__ __m128 __DEFAULT_FN_ATTRS
  142. _mm_andnot_ps(__m128 __a, __m128 __b)
  143. {
  144. return (__m128)(~(__v4si)__a & (__v4si)__b);
  145. }
  146. static __inline__ __m128 __DEFAULT_FN_ATTRS
  147. _mm_or_ps(__m128 __a, __m128 __b)
  148. {
  149. return (__m128)((__v4si)__a | (__v4si)__b);
  150. }
  151. static __inline__ __m128 __DEFAULT_FN_ATTRS
  152. _mm_xor_ps(__m128 __a, __m128 __b)
  153. {
  154. return (__m128)((__v4si)__a ^ (__v4si)__b);
  155. }
  156. static __inline__ __m128 __DEFAULT_FN_ATTRS
  157. _mm_cmpeq_ss(__m128 __a, __m128 __b)
  158. {
  159. return (__m128)__builtin_ia32_cmpeqss(__a, __b);
  160. }
  161. static __inline__ __m128 __DEFAULT_FN_ATTRS
  162. _mm_cmpeq_ps(__m128 __a, __m128 __b)
  163. {
  164. return (__m128)__builtin_ia32_cmpeqps(__a, __b);
  165. }
  166. static __inline__ __m128 __DEFAULT_FN_ATTRS
  167. _mm_cmplt_ss(__m128 __a, __m128 __b)
  168. {
  169. return (__m128)__builtin_ia32_cmpltss(__a, __b);
  170. }
  171. static __inline__ __m128 __DEFAULT_FN_ATTRS
  172. _mm_cmplt_ps(__m128 __a, __m128 __b)
  173. {
  174. return (__m128)__builtin_ia32_cmpltps(__a, __b);
  175. }
  176. static __inline__ __m128 __DEFAULT_FN_ATTRS
  177. _mm_cmple_ss(__m128 __a, __m128 __b)
  178. {
  179. return (__m128)__builtin_ia32_cmpless(__a, __b);
  180. }
  181. static __inline__ __m128 __DEFAULT_FN_ATTRS
  182. _mm_cmple_ps(__m128 __a, __m128 __b)
  183. {
  184. return (__m128)__builtin_ia32_cmpleps(__a, __b);
  185. }
  186. static __inline__ __m128 __DEFAULT_FN_ATTRS
  187. _mm_cmpgt_ss(__m128 __a, __m128 __b)
  188. {
  189. return (__m128)__builtin_shufflevector(__a,
  190. __builtin_ia32_cmpltss(__b, __a),
  191. 4, 1, 2, 3);
  192. }
  193. static __inline__ __m128 __DEFAULT_FN_ATTRS
  194. _mm_cmpgt_ps(__m128 __a, __m128 __b)
  195. {
  196. return (__m128)__builtin_ia32_cmpltps(__b, __a);
  197. }
  198. static __inline__ __m128 __DEFAULT_FN_ATTRS
  199. _mm_cmpge_ss(__m128 __a, __m128 __b)
  200. {
  201. return (__m128)__builtin_shufflevector(__a,
  202. __builtin_ia32_cmpless(__b, __a),
  203. 4, 1, 2, 3);
  204. }
  205. static __inline__ __m128 __DEFAULT_FN_ATTRS
  206. _mm_cmpge_ps(__m128 __a, __m128 __b)
  207. {
  208. return (__m128)__builtin_ia32_cmpleps(__b, __a);
  209. }
  210. static __inline__ __m128 __DEFAULT_FN_ATTRS
  211. _mm_cmpneq_ss(__m128 __a, __m128 __b)
  212. {
  213. return (__m128)__builtin_ia32_cmpneqss(__a, __b);
  214. }
  215. static __inline__ __m128 __DEFAULT_FN_ATTRS
  216. _mm_cmpneq_ps(__m128 __a, __m128 __b)
  217. {
  218. return (__m128)__builtin_ia32_cmpneqps(__a, __b);
  219. }
  220. static __inline__ __m128 __DEFAULT_FN_ATTRS
  221. _mm_cmpnlt_ss(__m128 __a, __m128 __b)
  222. {
  223. return (__m128)__builtin_ia32_cmpnltss(__a, __b);
  224. }
  225. static __inline__ __m128 __DEFAULT_FN_ATTRS
  226. _mm_cmpnlt_ps(__m128 __a, __m128 __b)
  227. {
  228. return (__m128)__builtin_ia32_cmpnltps(__a, __b);
  229. }
  230. static __inline__ __m128 __DEFAULT_FN_ATTRS
  231. _mm_cmpnle_ss(__m128 __a, __m128 __b)
  232. {
  233. return (__m128)__builtin_ia32_cmpnless(__a, __b);
  234. }
  235. static __inline__ __m128 __DEFAULT_FN_ATTRS
  236. _mm_cmpnle_ps(__m128 __a, __m128 __b)
  237. {
  238. return (__m128)__builtin_ia32_cmpnleps(__a, __b);
  239. }
  240. static __inline__ __m128 __DEFAULT_FN_ATTRS
  241. _mm_cmpngt_ss(__m128 __a, __m128 __b)
  242. {
  243. return (__m128)__builtin_shufflevector(__a,
  244. __builtin_ia32_cmpnltss(__b, __a),
  245. 4, 1, 2, 3);
  246. }
  247. static __inline__ __m128 __DEFAULT_FN_ATTRS
  248. _mm_cmpngt_ps(__m128 __a, __m128 __b)
  249. {
  250. return (__m128)__builtin_ia32_cmpnltps(__b, __a);
  251. }
  252. static __inline__ __m128 __DEFAULT_FN_ATTRS
  253. _mm_cmpnge_ss(__m128 __a, __m128 __b)
  254. {
  255. return (__m128)__builtin_shufflevector(__a,
  256. __builtin_ia32_cmpnless(__b, __a),
  257. 4, 1, 2, 3);
  258. }
  259. static __inline__ __m128 __DEFAULT_FN_ATTRS
  260. _mm_cmpnge_ps(__m128 __a, __m128 __b)
  261. {
  262. return (__m128)__builtin_ia32_cmpnleps(__b, __a);
  263. }
  264. static __inline__ __m128 __DEFAULT_FN_ATTRS
  265. _mm_cmpord_ss(__m128 __a, __m128 __b)
  266. {
  267. return (__m128)__builtin_ia32_cmpordss(__a, __b);
  268. }
  269. static __inline__ __m128 __DEFAULT_FN_ATTRS
  270. _mm_cmpord_ps(__m128 __a, __m128 __b)
  271. {
  272. return (__m128)__builtin_ia32_cmpordps(__a, __b);
  273. }
  274. static __inline__ __m128 __DEFAULT_FN_ATTRS
  275. _mm_cmpunord_ss(__m128 __a, __m128 __b)
  276. {
  277. return (__m128)__builtin_ia32_cmpunordss(__a, __b);
  278. }
  279. static __inline__ __m128 __DEFAULT_FN_ATTRS
  280. _mm_cmpunord_ps(__m128 __a, __m128 __b)
  281. {
  282. return (__m128)__builtin_ia32_cmpunordps(__a, __b);
  283. }
  284. static __inline__ int __DEFAULT_FN_ATTRS
  285. _mm_comieq_ss(__m128 __a, __m128 __b)
  286. {
  287. return __builtin_ia32_comieq(__a, __b);
  288. }
  289. static __inline__ int __DEFAULT_FN_ATTRS
  290. _mm_comilt_ss(__m128 __a, __m128 __b)
  291. {
  292. return __builtin_ia32_comilt(__a, __b);
  293. }
  294. static __inline__ int __DEFAULT_FN_ATTRS
  295. _mm_comile_ss(__m128 __a, __m128 __b)
  296. {
  297. return __builtin_ia32_comile(__a, __b);
  298. }
  299. static __inline__ int __DEFAULT_FN_ATTRS
  300. _mm_comigt_ss(__m128 __a, __m128 __b)
  301. {
  302. return __builtin_ia32_comigt(__a, __b);
  303. }
  304. static __inline__ int __DEFAULT_FN_ATTRS
  305. _mm_comige_ss(__m128 __a, __m128 __b)
  306. {
  307. return __builtin_ia32_comige(__a, __b);
  308. }
  309. static __inline__ int __DEFAULT_FN_ATTRS
  310. _mm_comineq_ss(__m128 __a, __m128 __b)
  311. {
  312. return __builtin_ia32_comineq(__a, __b);
  313. }
  314. static __inline__ int __DEFAULT_FN_ATTRS
  315. _mm_ucomieq_ss(__m128 __a, __m128 __b)
  316. {
  317. return __builtin_ia32_ucomieq(__a, __b);
  318. }
  319. static __inline__ int __DEFAULT_FN_ATTRS
  320. _mm_ucomilt_ss(__m128 __a, __m128 __b)
  321. {
  322. return __builtin_ia32_ucomilt(__a, __b);
  323. }
  324. static __inline__ int __DEFAULT_FN_ATTRS
  325. _mm_ucomile_ss(__m128 __a, __m128 __b)
  326. {
  327. return __builtin_ia32_ucomile(__a, __b);
  328. }
  329. static __inline__ int __DEFAULT_FN_ATTRS
  330. _mm_ucomigt_ss(__m128 __a, __m128 __b)
  331. {
  332. return __builtin_ia32_ucomigt(__a, __b);
  333. }
  334. static __inline__ int __DEFAULT_FN_ATTRS
  335. _mm_ucomige_ss(__m128 __a, __m128 __b)
  336. {
  337. return __builtin_ia32_ucomige(__a, __b);
  338. }
  339. static __inline__ int __DEFAULT_FN_ATTRS
  340. _mm_ucomineq_ss(__m128 __a, __m128 __b)
  341. {
  342. return __builtin_ia32_ucomineq(__a, __b);
  343. }
  344. static __inline__ int __DEFAULT_FN_ATTRS
  345. _mm_cvtss_si32(__m128 __a)
  346. {
  347. return __builtin_ia32_cvtss2si(__a);
  348. }
  349. static __inline__ int __DEFAULT_FN_ATTRS
  350. _mm_cvt_ss2si(__m128 __a)
  351. {
  352. return _mm_cvtss_si32(__a);
  353. }
  354. #ifdef __x86_64__
  355. static __inline__ long long __DEFAULT_FN_ATTRS
  356. _mm_cvtss_si64(__m128 __a)
  357. {
  358. return __builtin_ia32_cvtss2si64(__a);
  359. }
  360. #endif
  361. static __inline__ __m64 __DEFAULT_FN_ATTRS
  362. _mm_cvtps_pi32(__m128 __a)
  363. {
  364. return (__m64)__builtin_ia32_cvtps2pi(__a);
  365. }
  366. static __inline__ __m64 __DEFAULT_FN_ATTRS
  367. _mm_cvt_ps2pi(__m128 __a)
  368. {
  369. return _mm_cvtps_pi32(__a);
  370. }
  371. static __inline__ int __DEFAULT_FN_ATTRS
  372. _mm_cvttss_si32(__m128 __a)
  373. {
  374. return __a[0];
  375. }
  376. static __inline__ int __DEFAULT_FN_ATTRS
  377. _mm_cvtt_ss2si(__m128 __a)
  378. {
  379. return _mm_cvttss_si32(__a);
  380. }
  381. static __inline__ long long __DEFAULT_FN_ATTRS
  382. _mm_cvttss_si64(__m128 __a)
  383. {
  384. return __a[0];
  385. }
  386. static __inline__ __m64 __DEFAULT_FN_ATTRS
  387. _mm_cvttps_pi32(__m128 __a)
  388. {
  389. return (__m64)__builtin_ia32_cvttps2pi(__a);
  390. }
  391. static __inline__ __m64 __DEFAULT_FN_ATTRS
  392. _mm_cvtt_ps2pi(__m128 __a)
  393. {
  394. return _mm_cvttps_pi32(__a);
  395. }
  396. static __inline__ __m128 __DEFAULT_FN_ATTRS
  397. _mm_cvtsi32_ss(__m128 __a, int __b)
  398. {
  399. __a[0] = __b;
  400. return __a;
  401. }
  402. static __inline__ __m128 __DEFAULT_FN_ATTRS
  403. _mm_cvt_si2ss(__m128 __a, int __b)
  404. {
  405. return _mm_cvtsi32_ss(__a, __b);
  406. }
  407. #ifdef __x86_64__
  408. static __inline__ __m128 __DEFAULT_FN_ATTRS
  409. _mm_cvtsi64_ss(__m128 __a, long long __b)
  410. {
  411. __a[0] = __b;
  412. return __a;
  413. }
  414. #endif
  415. static __inline__ __m128 __DEFAULT_FN_ATTRS
  416. _mm_cvtpi32_ps(__m128 __a, __m64 __b)
  417. {
  418. return __builtin_ia32_cvtpi2ps(__a, (__v2si)__b);
  419. }
  420. static __inline__ __m128 __DEFAULT_FN_ATTRS
  421. _mm_cvt_pi2ps(__m128 __a, __m64 __b)
  422. {
  423. return _mm_cvtpi32_ps(__a, __b);
  424. }
  425. static __inline__ float __DEFAULT_FN_ATTRS
  426. _mm_cvtss_f32(__m128 __a)
  427. {
  428. return __a[0];
  429. }
  430. static __inline__ __m128 __DEFAULT_FN_ATTRS
  431. _mm_loadh_pi(__m128 __a, const __m64 *__p)
  432. {
  433. typedef float __mm_loadh_pi_v2f32 __attribute__((__vector_size__(8)));
  434. struct __mm_loadh_pi_struct {
  435. __mm_loadh_pi_v2f32 __u;
  436. } __attribute__((__packed__, __may_alias__));
  437. __mm_loadh_pi_v2f32 __b = ((struct __mm_loadh_pi_struct*)__p)->__u;
  438. __m128 __bb = __builtin_shufflevector(__b, __b, 0, 1, 0, 1);
  439. return __builtin_shufflevector(__a, __bb, 0, 1, 4, 5);
  440. }
  441. static __inline__ __m128 __DEFAULT_FN_ATTRS
  442. _mm_loadl_pi(__m128 __a, const __m64 *__p)
  443. {
  444. typedef float __mm_loadl_pi_v2f32 __attribute__((__vector_size__(8)));
  445. struct __mm_loadl_pi_struct {
  446. __mm_loadl_pi_v2f32 __u;
  447. } __attribute__((__packed__, __may_alias__));
  448. __mm_loadl_pi_v2f32 __b = ((struct __mm_loadl_pi_struct*)__p)->__u;
  449. __m128 __bb = __builtin_shufflevector(__b, __b, 0, 1, 0, 1);
  450. return __builtin_shufflevector(__a, __bb, 4, 5, 2, 3);
  451. }
  452. static __inline__ __m128 __DEFAULT_FN_ATTRS
  453. _mm_load_ss(const float *__p)
  454. {
  455. struct __mm_load_ss_struct {
  456. float __u;
  457. } __attribute__((__packed__, __may_alias__));
  458. float __u = ((struct __mm_load_ss_struct*)__p)->__u;
  459. return (__m128){ __u, 0, 0, 0 };
  460. }
  461. static __inline__ __m128 __DEFAULT_FN_ATTRS
  462. _mm_load1_ps(const float *__p)
  463. {
  464. struct __mm_load1_ps_struct {
  465. float __u;
  466. } __attribute__((__packed__, __may_alias__));
  467. float __u = ((struct __mm_load1_ps_struct*)__p)->__u;
  468. return (__m128){ __u, __u, __u, __u };
  469. }
  470. #define _mm_load_ps1(p) _mm_load1_ps(p)
  471. static __inline__ __m128 __DEFAULT_FN_ATTRS
  472. _mm_load_ps(const float *__p)
  473. {
  474. return *(__m128*)__p;
  475. }
  476. static __inline__ __m128 __DEFAULT_FN_ATTRS
  477. _mm_loadu_ps(const float *__p)
  478. {
  479. struct __loadu_ps {
  480. __m128 __v;
  481. } __attribute__((__packed__, __may_alias__));
  482. return ((struct __loadu_ps*)__p)->__v;
  483. }
  484. static __inline__ __m128 __DEFAULT_FN_ATTRS
  485. _mm_loadr_ps(const float *__p)
  486. {
  487. __m128 __a = _mm_load_ps(__p);
  488. return __builtin_shufflevector(__a, __a, 3, 2, 1, 0);
  489. }
  490. static __inline__ __m128 __DEFAULT_FN_ATTRS
  491. _mm_set_ss(float __w)
  492. {
  493. return (__m128){ __w, 0, 0, 0 };
  494. }
  495. static __inline__ __m128 __DEFAULT_FN_ATTRS
  496. _mm_set1_ps(float __w)
  497. {
  498. return (__m128){ __w, __w, __w, __w };
  499. }
  500. /* Microsoft specific. */
  501. static __inline__ __m128 __DEFAULT_FN_ATTRS
  502. _mm_set_ps1(float __w)
  503. {
  504. return _mm_set1_ps(__w);
  505. }
  506. static __inline__ __m128 __DEFAULT_FN_ATTRS
  507. _mm_set_ps(float __z, float __y, float __x, float __w)
  508. {
  509. return (__m128){ __w, __x, __y, __z };
  510. }
  511. static __inline__ __m128 __DEFAULT_FN_ATTRS
  512. _mm_setr_ps(float __z, float __y, float __x, float __w)
  513. {
  514. return (__m128){ __z, __y, __x, __w };
  515. }
  516. static __inline__ __m128 __DEFAULT_FN_ATTRS
  517. _mm_setzero_ps(void)
  518. {
  519. return (__m128){ 0, 0, 0, 0 };
  520. }
  521. static __inline__ void __DEFAULT_FN_ATTRS
  522. _mm_storeh_pi(__m64 *__p, __m128 __a)
  523. {
  524. __builtin_ia32_storehps((__v2si *)__p, __a);
  525. }
  526. static __inline__ void __DEFAULT_FN_ATTRS
  527. _mm_storel_pi(__m64 *__p, __m128 __a)
  528. {
  529. __builtin_ia32_storelps((__v2si *)__p, __a);
  530. }
  531. static __inline__ void __DEFAULT_FN_ATTRS
  532. _mm_store_ss(float *__p, __m128 __a)
  533. {
  534. struct __mm_store_ss_struct {
  535. float __u;
  536. } __attribute__((__packed__, __may_alias__));
  537. ((struct __mm_store_ss_struct*)__p)->__u = __a[0];
  538. }
  539. static __inline__ void __DEFAULT_FN_ATTRS
  540. _mm_storeu_ps(float *__p, __m128 __a)
  541. {
  542. __builtin_ia32_storeups(__p, __a);
  543. }
  544. static __inline__ void __DEFAULT_FN_ATTRS
  545. _mm_store1_ps(float *__p, __m128 __a)
  546. {
  547. __a = __builtin_shufflevector(__a, __a, 0, 0, 0, 0);
  548. _mm_storeu_ps(__p, __a);
  549. }
  550. static __inline__ void __DEFAULT_FN_ATTRS
  551. _mm_store_ps1(float *__p, __m128 __a)
  552. {
  553. return _mm_store1_ps(__p, __a);
  554. }
  555. static __inline__ void __DEFAULT_FN_ATTRS
  556. _mm_store_ps(float *__p, __m128 __a)
  557. {
  558. *(__m128 *)__p = __a;
  559. }
  560. static __inline__ void __DEFAULT_FN_ATTRS
  561. _mm_storer_ps(float *__p, __m128 __a)
  562. {
  563. __a = __builtin_shufflevector(__a, __a, 3, 2, 1, 0);
  564. _mm_store_ps(__p, __a);
  565. }
  566. #define _MM_HINT_T0 3
  567. #define _MM_HINT_T1 2
  568. #define _MM_HINT_T2 1
  569. #define _MM_HINT_NTA 0
  570. #ifndef _MSC_VER
  571. /* FIXME: We have to #define this because "sel" must be a constant integer, and
  572. Sema doesn't do any form of constant propagation yet. */
  573. #define _mm_prefetch(a, sel) (__builtin_prefetch((void *)(a), 0, (sel)))
  574. #endif
  575. static __inline__ void __DEFAULT_FN_ATTRS
  576. _mm_stream_pi(__m64 *__p, __m64 __a)
  577. {
  578. __builtin_ia32_movntq(__p, __a);
  579. }
  580. static __inline__ void __DEFAULT_FN_ATTRS
  581. _mm_stream_ps(float *__p, __m128 __a)
  582. {
  583. __builtin_ia32_movntps(__p, __a);
  584. }
  585. static __inline__ void __DEFAULT_FN_ATTRS
  586. _mm_sfence(void)
  587. {
  588. __builtin_ia32_sfence();
  589. }
  590. static __inline__ int __DEFAULT_FN_ATTRS
  591. _mm_extract_pi16(__m64 __a, int __n)
  592. {
  593. __v4hi __b = (__v4hi)__a;
  594. return (unsigned short)__b[__n & 3];
  595. }
  596. static __inline__ __m64 __DEFAULT_FN_ATTRS
  597. _mm_insert_pi16(__m64 __a, int __d, int __n)
  598. {
  599. __v4hi __b = (__v4hi)__a;
  600. __b[__n & 3] = __d;
  601. return (__m64)__b;
  602. }
  603. static __inline__ __m64 __DEFAULT_FN_ATTRS
  604. _mm_max_pi16(__m64 __a, __m64 __b)
  605. {
  606. return (__m64)__builtin_ia32_pmaxsw((__v4hi)__a, (__v4hi)__b);
  607. }
  608. static __inline__ __m64 __DEFAULT_FN_ATTRS
  609. _mm_max_pu8(__m64 __a, __m64 __b)
  610. {
  611. return (__m64)__builtin_ia32_pmaxub((__v8qi)__a, (__v8qi)__b);
  612. }
  613. static __inline__ __m64 __DEFAULT_FN_ATTRS
  614. _mm_min_pi16(__m64 __a, __m64 __b)
  615. {
  616. return (__m64)__builtin_ia32_pminsw((__v4hi)__a, (__v4hi)__b);
  617. }
  618. static __inline__ __m64 __DEFAULT_FN_ATTRS
  619. _mm_min_pu8(__m64 __a, __m64 __b)
  620. {
  621. return (__m64)__builtin_ia32_pminub((__v8qi)__a, (__v8qi)__b);
  622. }
  623. static __inline__ int __DEFAULT_FN_ATTRS
  624. _mm_movemask_pi8(__m64 __a)
  625. {
  626. return __builtin_ia32_pmovmskb((__v8qi)__a);
  627. }
  628. static __inline__ __m64 __DEFAULT_FN_ATTRS
  629. _mm_mulhi_pu16(__m64 __a, __m64 __b)
  630. {
  631. return (__m64)__builtin_ia32_pmulhuw((__v4hi)__a, (__v4hi)__b);
  632. }
  633. #define _mm_shuffle_pi16(a, n) __extension__ ({ \
  634. __m64 __a = (a); \
  635. (__m64)__builtin_ia32_pshufw((__v4hi)__a, (n)); })
  636. static __inline__ void __DEFAULT_FN_ATTRS
  637. _mm_maskmove_si64(__m64 __d, __m64 __n, char *__p)
  638. {
  639. __builtin_ia32_maskmovq((__v8qi)__d, (__v8qi)__n, __p);
  640. }
  641. static __inline__ __m64 __DEFAULT_FN_ATTRS
  642. _mm_avg_pu8(__m64 __a, __m64 __b)
  643. {
  644. return (__m64)__builtin_ia32_pavgb((__v8qi)__a, (__v8qi)__b);
  645. }
  646. static __inline__ __m64 __DEFAULT_FN_ATTRS
  647. _mm_avg_pu16(__m64 __a, __m64 __b)
  648. {
  649. return (__m64)__builtin_ia32_pavgw((__v4hi)__a, (__v4hi)__b);
  650. }
  651. static __inline__ __m64 __DEFAULT_FN_ATTRS
  652. _mm_sad_pu8(__m64 __a, __m64 __b)
  653. {
  654. return (__m64)__builtin_ia32_psadbw((__v8qi)__a, (__v8qi)__b);
  655. }
  656. static __inline__ unsigned int __DEFAULT_FN_ATTRS
  657. _mm_getcsr(void)
  658. {
  659. return __builtin_ia32_stmxcsr();
  660. }
  661. static __inline__ void __DEFAULT_FN_ATTRS
  662. _mm_setcsr(unsigned int __i)
  663. {
  664. __builtin_ia32_ldmxcsr(__i);
  665. }
  666. #define _mm_shuffle_ps(a, b, mask) __extension__ ({ \
  667. __m128 __a = (a); \
  668. __m128 __b = (b); \
  669. (__m128)__builtin_shufflevector((__v4sf)__a, (__v4sf)__b, \
  670. (mask) & 0x3, ((mask) & 0xc) >> 2, \
  671. (((mask) & 0x30) >> 4) + 4, \
  672. (((mask) & 0xc0) >> 6) + 4); })
  673. static __inline__ __m128 __DEFAULT_FN_ATTRS
  674. _mm_unpackhi_ps(__m128 __a, __m128 __b)
  675. {
  676. return __builtin_shufflevector(__a, __b, 2, 6, 3, 7);
  677. }
  678. static __inline__ __m128 __DEFAULT_FN_ATTRS
  679. _mm_unpacklo_ps(__m128 __a, __m128 __b)
  680. {
  681. return __builtin_shufflevector(__a, __b, 0, 4, 1, 5);
  682. }
  683. static __inline__ __m128 __DEFAULT_FN_ATTRS
  684. _mm_move_ss(__m128 __a, __m128 __b)
  685. {
  686. return __builtin_shufflevector(__a, __b, 4, 1, 2, 3);
  687. }
  688. static __inline__ __m128 __DEFAULT_FN_ATTRS
  689. _mm_movehl_ps(__m128 __a, __m128 __b)
  690. {
  691. return __builtin_shufflevector(__a, __b, 6, 7, 2, 3);
  692. }
  693. static __inline__ __m128 __DEFAULT_FN_ATTRS
  694. _mm_movelh_ps(__m128 __a, __m128 __b)
  695. {
  696. return __builtin_shufflevector(__a, __b, 0, 1, 4, 5);
  697. }
  698. static __inline__ __m128 __DEFAULT_FN_ATTRS
  699. _mm_cvtpi16_ps(__m64 __a)
  700. {
  701. __m64 __b, __c;
  702. __m128 __r;
  703. __b = _mm_setzero_si64();
  704. __b = _mm_cmpgt_pi16(__b, __a);
  705. __c = _mm_unpackhi_pi16(__a, __b);
  706. __r = _mm_setzero_ps();
  707. __r = _mm_cvtpi32_ps(__r, __c);
  708. __r = _mm_movelh_ps(__r, __r);
  709. __c = _mm_unpacklo_pi16(__a, __b);
  710. __r = _mm_cvtpi32_ps(__r, __c);
  711. return __r;
  712. }
  713. static __inline__ __m128 __DEFAULT_FN_ATTRS
  714. _mm_cvtpu16_ps(__m64 __a)
  715. {
  716. __m64 __b, __c;
  717. __m128 __r;
  718. __b = _mm_setzero_si64();
  719. __c = _mm_unpackhi_pi16(__a, __b);
  720. __r = _mm_setzero_ps();
  721. __r = _mm_cvtpi32_ps(__r, __c);
  722. __r = _mm_movelh_ps(__r, __r);
  723. __c = _mm_unpacklo_pi16(__a, __b);
  724. __r = _mm_cvtpi32_ps(__r, __c);
  725. return __r;
  726. }
  727. static __inline__ __m128 __DEFAULT_FN_ATTRS
  728. _mm_cvtpi8_ps(__m64 __a)
  729. {
  730. __m64 __b;
  731. __b = _mm_setzero_si64();
  732. __b = _mm_cmpgt_pi8(__b, __a);
  733. __b = _mm_unpacklo_pi8(__a, __b);
  734. return _mm_cvtpi16_ps(__b);
  735. }
  736. static __inline__ __m128 __DEFAULT_FN_ATTRS
  737. _mm_cvtpu8_ps(__m64 __a)
  738. {
  739. __m64 __b;
  740. __b = _mm_setzero_si64();
  741. __b = _mm_unpacklo_pi8(__a, __b);
  742. return _mm_cvtpi16_ps(__b);
  743. }
  744. static __inline__ __m128 __DEFAULT_FN_ATTRS
  745. _mm_cvtpi32x2_ps(__m64 __a, __m64 __b)
  746. {
  747. __m128 __c;
  748. __c = _mm_setzero_ps();
  749. __c = _mm_cvtpi32_ps(__c, __b);
  750. __c = _mm_movelh_ps(__c, __c);
  751. return _mm_cvtpi32_ps(__c, __a);
  752. }
  753. static __inline__ __m64 __DEFAULT_FN_ATTRS
  754. _mm_cvtps_pi16(__m128 __a)
  755. {
  756. __m64 __b, __c;
  757. __b = _mm_cvtps_pi32(__a);
  758. __a = _mm_movehl_ps(__a, __a);
  759. __c = _mm_cvtps_pi32(__a);
  760. return _mm_packs_pi32(__b, __c);
  761. }
  762. static __inline__ __m64 __DEFAULT_FN_ATTRS
  763. _mm_cvtps_pi8(__m128 __a)
  764. {
  765. __m64 __b, __c;
  766. __b = _mm_cvtps_pi16(__a);
  767. __c = _mm_setzero_si64();
  768. return _mm_packs_pi16(__b, __c);
  769. }
  770. static __inline__ int __DEFAULT_FN_ATTRS
  771. _mm_movemask_ps(__m128 __a)
  772. {
  773. return __builtin_ia32_movmskps(__a);
  774. }
  775. #define _MM_SHUFFLE(z, y, x, w) (((z) << 6) | ((y) << 4) | ((x) << 2) | (w))
  776. #define _MM_EXCEPT_INVALID (0x0001)
  777. #define _MM_EXCEPT_DENORM (0x0002)
  778. #define _MM_EXCEPT_DIV_ZERO (0x0004)
  779. #define _MM_EXCEPT_OVERFLOW (0x0008)
  780. #define _MM_EXCEPT_UNDERFLOW (0x0010)
  781. #define _MM_EXCEPT_INEXACT (0x0020)
  782. #define _MM_EXCEPT_MASK (0x003f)
  783. #define _MM_MASK_INVALID (0x0080)
  784. #define _MM_MASK_DENORM (0x0100)
  785. #define _MM_MASK_DIV_ZERO (0x0200)
  786. #define _MM_MASK_OVERFLOW (0x0400)
  787. #define _MM_MASK_UNDERFLOW (0x0800)
  788. #define _MM_MASK_INEXACT (0x1000)
  789. #define _MM_MASK_MASK (0x1f80)
  790. #define _MM_ROUND_NEAREST (0x0000)
  791. #define _MM_ROUND_DOWN (0x2000)
  792. #define _MM_ROUND_UP (0x4000)
  793. #define _MM_ROUND_TOWARD_ZERO (0x6000)
  794. #define _MM_ROUND_MASK (0x6000)
  795. #define _MM_FLUSH_ZERO_MASK (0x8000)
  796. #define _MM_FLUSH_ZERO_ON (0x8000)
  797. #define _MM_FLUSH_ZERO_OFF (0x0000)
  798. #define _MM_GET_EXCEPTION_MASK() (_mm_getcsr() & _MM_MASK_MASK)
  799. #define _MM_GET_EXCEPTION_STATE() (_mm_getcsr() & _MM_EXCEPT_MASK)
  800. #define _MM_GET_FLUSH_ZERO_MODE() (_mm_getcsr() & _MM_FLUSH_ZERO_MASK)
  801. #define _MM_GET_ROUNDING_MODE() (_mm_getcsr() & _MM_ROUND_MASK)
  802. #define _MM_SET_EXCEPTION_MASK(x) (_mm_setcsr((_mm_getcsr() & ~_MM_MASK_MASK) | (x)))
  803. #define _MM_SET_EXCEPTION_STATE(x) (_mm_setcsr((_mm_getcsr() & ~_MM_EXCEPT_MASK) | (x)))
  804. #define _MM_SET_FLUSH_ZERO_MODE(x) (_mm_setcsr((_mm_getcsr() & ~_MM_FLUSH_ZERO_MASK) | (x)))
  805. #define _MM_SET_ROUNDING_MODE(x) (_mm_setcsr((_mm_getcsr() & ~_MM_ROUND_MASK) | (x)))
  806. #define _MM_TRANSPOSE4_PS(row0, row1, row2, row3) \
  807. do { \
  808. __m128 tmp3, tmp2, tmp1, tmp0; \
  809. tmp0 = _mm_unpacklo_ps((row0), (row1)); \
  810. tmp2 = _mm_unpacklo_ps((row2), (row3)); \
  811. tmp1 = _mm_unpackhi_ps((row0), (row1)); \
  812. tmp3 = _mm_unpackhi_ps((row2), (row3)); \
  813. (row0) = _mm_movelh_ps(tmp0, tmp2); \
  814. (row1) = _mm_movehl_ps(tmp2, tmp0); \
  815. (row2) = _mm_movelh_ps(tmp1, tmp3); \
  816. (row3) = _mm_movehl_ps(tmp3, tmp1); \
  817. } while (0)
  818. /* Aliases for compatibility. */
  819. #define _m_pextrw _mm_extract_pi16
  820. #define _m_pinsrw _mm_insert_pi16
  821. #define _m_pmaxsw _mm_max_pi16
  822. #define _m_pmaxub _mm_max_pu8
  823. #define _m_pminsw _mm_min_pi16
  824. #define _m_pminub _mm_min_pu8
  825. #define _m_pmovmskb _mm_movemask_pi8
  826. #define _m_pmulhuw _mm_mulhi_pu16
  827. #define _m_pshufw _mm_shuffle_pi16
  828. #define _m_maskmovq _mm_maskmove_si64
  829. #define _m_pavgb _mm_avg_pu8
  830. #define _m_pavgw _mm_avg_pu16
  831. #define _m_psadbw _mm_sad_pu8
  832. #define _m_ _mm_
  833. #define _m_ _mm_
  834. #undef __DEFAULT_FN_ATTRS
  835. /* Ugly hack for backwards-compatibility (compatible with gcc) */
  836. #if defined(__SSE2__) && !__has_feature(modules)
  837. #include <emmintrin.h>
  838. #endif
  839. #endif /* __SSE__ */
  840. #endif /* __XMMINTRIN_H */