sse2.h 97 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233123412351236123712381239124012411242124312441245124612471248124912501251125212531254125512561257125812591260126112621263126412651266126712681269127012711272127312741275127612771278127912801281128212831284128512861287128812891290129112921293129412951296129712981299130013011302130313041305130613071308130913101311131213131314131513161317131813191320132113221323132413251326132713281329133013311332133313341335133613371338133913401341134213431344134513461347134813491350135113521353135413551356135713581359136013611362136313641365136613671368136913701371137213731374137513761377137813791380138113821383138413851386138713881389139013911392139313941395139613971398139914001401140214031404140514061407140814091410141114121413141414151416141714181419142014211422142314241425142614271428142914301431143214331434143514361437143814391440144114421443144414451446144714481449145014511452145314541455145614571458145914601461146214631464146514661467146814691470147114721473147414751476147714781479148014811482148314841485148614871488148914901491149214931494149514961497149814991500150115021503150415051506150715081509151015111512151315141515151615171518151915201521152215231524152515261527152815291530153115321533153415351536153715381539154015411542154315441545154615471548154915501551155215531554155515561557155815591560156115621563156415651566156715681569157015711572157315741575157615771578157915801581158215831584158515861587158815891590159115921593159415951596159715981599160016011602160316041605160616071608160916101611161216131614161516161617161816191620162116221623162416251626162716281629163016311632163316341635163616371638163916401641164216431644164516461647164816491650165116521653165416551656165716581659166016611662166316641665166616671668166916701671167216731674167516761677167816791680168116821683168416851686168716881689169016911692169316941695169616971698169917001701170217031704170517061707170817091710171117121713171417151716171717181719172017211722172317241725172617271728172917301731173217331734173517361737173817391740174117421743174417451746174717481749175017511752175317541755175617571758175917601761176217631764176517661767176817691770177117721773177417751776177717781779178017811782178317841785178617871788178917901791179217931794179517961797179817991800180118021803180418051806180718081809181018111812181318141815181618171818181918201821182218231824182518261827182818291830183118321833183418351836183718381839184018411842184318441845184618471848184918501851185218531854185518561857185818591860186118621863186418651866186718681869187018711872187318741875187618771878187918801881188218831884188518861887188818891890189118921893189418951896189718981899190019011902190319041905190619071908190919101911191219131914191519161917191819191920192119221923192419251926192719281929193019311932193319341935193619371938193919401941194219431944194519461947194819491950195119521953195419551956195719581959196019611962196319641965196619671968196919701971197219731974197519761977197819791980198119821983198419851986198719881989199019911992199319941995199619971998199920002001200220032004200520062007200820092010201120122013201420152016201720182019202020212022202320242025202620272028202920302031203220332034203520362037203820392040204120422043204420452046204720482049205020512052205320542055205620572058205920602061206220632064206520662067206820692070207120722073207420752076207720782079208020812082208320842085208620872088208920902091209220932094209520962097209820992100210121022103210421052106210721082109211021112112211321142115211621172118211921202121212221232124212521262127212821292130213121322133213421352136213721382139214021412142214321442145214621472148214921502151215221532154215521562157215821592160216121622163216421652166216721682169217021712172217321742175217621772178217921802181218221832184218521862187218821892190219121922193219421952196219721982199220022012202220322042205220622072208220922102211221222132214221522162217221822192220222122222223222422252226222722282229223022312232223322342235223622372238223922402241224222432244224522462247224822492250225122522253225422552256225722582259226022612262226322642265226622672268226922702271227222732274227522762277227822792280228122822283228422852286228722882289229022912292229322942295229622972298229923002301230223032304230523062307230823092310231123122313231423152316231723182319232023212322232323242325232623272328232923302331233223332334233523362337233823392340234123422343234423452346234723482349235023512352235323542355235623572358235923602361236223632364236523662367236823692370237123722373237423752376237723782379238023812382238323842385238623872388238923902391239223932394239523962397239823992400240124022403240424052406240724082409241024112412241324142415241624172418241924202421242224232424242524262427242824292430243124322433243424352436243724382439244024412442244324442445244624472448244924502451245224532454245524562457245824592460246124622463246424652466246724682469247024712472247324742475247624772478247924802481248224832484248524862487248824892490249124922493249424952496249724982499250025012502250325042505250625072508250925102511251225132514251525162517251825192520252125222523252425252526252725282529253025312532253325342535253625372538253925402541254225432544254525462547254825492550255125522553255425552556255725582559256025612562256325642565256625672568256925702571257225732574257525762577257825792580258125822583258425852586258725882589259025912592259325942595259625972598259926002601260226032604260526062607260826092610261126122613261426152616261726182619262026212622262326242625262626272628262926302631263226332634263526362637263826392640264126422643264426452646264726482649265026512652265326542655265626572658265926602661266226632664266526662667266826692670267126722673267426752676267726782679268026812682268326842685268626872688268926902691269226932694269526962697269826992700270127022703270427052706270727082709271027112712271327142715271627172718271927202721272227232724272527262727272827292730273127322733273427352736273727382739274027412742274327442745274627472748274927502751275227532754275527562757275827592760276127622763276427652766276727682769277027712772277327742775277627772778277927802781278227832784278527862787278827892790279127922793279427952796279727982799280028012802280328042805280628072808280928102811281228132814281528162817281828192820282128222823282428252826282728282829283028312832283328342835283628372838283928402841284228432844284528462847284828492850285128522853285428552856285728582859286028612862286328642865286628672868286928702871287228732874287528762877287828792880288128822883288428852886288728882889289028912892289328942895289628972898289929002901290229032904290529062907290829092910291129122913291429152916291729182919292029212922292329242925292629272928292929302931293229332934293529362937293829392940294129422943294429452946294729482949295029512952295329542955295629572958295929602961296229632964296529662967296829692970297129722973297429752976297729782979298029812982298329842985298629872988298929902991299229932994299529962997299829993000300130023003300430053006300730083009301030113012301330143015301630173018301930203021302230233024302530263027302830293030303130323033303430353036303730383039304030413042304330443045304630473048304930503051305230533054305530563057305830593060306130623063306430653066306730683069307030713072307330743075307630773078307930803081308230833084308530863087308830893090309130923093309430953096309730983099310031013102310331043105310631073108310931103111311231133114311531163117311831193120312131223123312431253126312731283129313031313132313331343135313631373138313931403141314231433144314531463147314831493150315131523153315431553156315731583159316031613162316331643165316631673168316931703171317231733174317531763177317831793180318131823183318431853186318731883189319031913192319331943195319631973198319932003201320232033204320532063207320832093210321132123213321432153216321732183219322032213222322332243225322632273228322932303231323232333234323532363237323832393240324132423243324432453246324732483249325032513252325332543255325632573258325932603261326232633264326532663267326832693270327132723273327432753276327732783279328032813282328332843285328632873288328932903291329232933294329532963297329832993300330133023303330433053306330733083309331033113312331333143315331633173318331933203321332233233324332533263327332833293330333133323333333433353336333733383339334033413342334333443345334633473348334933503351335233533354335533563357335833593360336133623363336433653366336733683369337033713372337333743375337633773378337933803381338233833384338533863387338833893390339133923393339433953396339733983399340034013402340334043405340634073408340934103411341234133414341534163417341834193420342134223423342434253426342734283429343034313432343334343435343634373438343934403441344234433444344534463447344834493450345134523453345434553456345734583459346034613462346334643465346634673468346934703471347234733474347534763477347834793480348134823483348434853486348734883489349034913492349334943495349634973498349935003501350235033504350535063507350835093510351135123513351435153516351735183519352035213522352335243525352635273528352935303531353235333534353535363537353835393540354135423543354435453546354735483549355035513552355335543555355635573558355935603561356235633564356535663567356835693570357135723573357435753576357735783579358035813582358335843585358635873588358935903591359235933594359535963597359835993600360136023603360436053606360736083609361036113612361336143615361636173618361936203621362236233624362536263627362836293630363136323633363436353636363736383639364036413642364336443645364636473648364936503651365236533654365536563657365836593660366136623663366436653666366736683669367036713672367336743675367636773678367936803681368236833684368536863687368836893690369136923693369436953696369736983699370037013702370337043705370637073708370937103711371237133714371537163717371837193720372137223723372437253726372737283729373037313732373337343735373637373738373937403741374237433744374537463747374837493750375137523753375437553756375737583759376037613762376337643765376637673768376937703771377237733774377537763777377837793780378137823783378437853786378737883789379037913792379337943795379637973798379938003801380238033804380538063807380838093810381138123813381438153816381738183819382038213822382338243825382638273828382938303831383238333834383538363837383838393840384138423843384438453846384738483849385038513852385338543855385638573858385938603861386238633864386538663867386838693870387138723873387438753876387738783879388038813882388338843885388638873888388938903891389238933894389538963897389838993900390139023903390439053906390739083909391039113912391339143915391639173918391939203921392239233924392539263927392839293930393139323933393439353936393739383939394039413942394339443945394639473948394939503951395239533954395539563957395839593960396139623963396439653966396739683969397039713972397339743975397639773978397939803981398239833984398539863987398839893990399139923993399439953996399739983999400040014002400340044005400640074008400940104011401240134014401540164017401840194020402140224023402440254026402740284029403040314032403340344035403640374038403940404041404240434044404540464047404840494050405140524053405440554056405740584059406040614062406340644065406640674068406940704071407240734074407540764077407840794080408140824083408440854086408740884089409040914092409340944095409640974098409941004101410241034104410541064107410841094110411141124113411441154116411741184119412041214122412341244125412641274128412941304131413241334134413541364137413841394140414141424143414441454146414741484149415041514152415341544155415641574158415941604161416241634164416541664167416841694170417141724173417441754176417741784179418041814182418341844185418641874188418941904191419241934194419541964197
  1. /* Permission is hereby granted, free of charge, to any person
  2. * obtaining a copy of this software and associated documentation
  3. * files (the "Software"), to deal in the Software without
  4. * restriction, including without limitation the rights to use, copy,
  5. * modify, merge, publish, distribute, sublicense, and/or sell copies
  6. * of the Software, and to permit persons to whom the Software is
  7. * furnished to do so, subject to the following conditions:
  8. *
  9. * The above copyright notice and this permission notice shall be
  10. * included in all copies or substantial portions of the Software.
  11. *
  12. * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
  13. * EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
  14. * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
  15. * NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS
  16. * BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN
  17. * ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
  18. * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
  19. * SOFTWARE.
  20. *
  21. * Copyright:
  22. * 2017 Evan Nemerson <[email protected]>
  23. * 2015-2017 John W. Ratcliff <[email protected]>
  24. * 2015 Brandon Rowlett <[email protected]>
  25. * 2015 Ken Fast <[email protected]>
  26. * 2017 Hasindu Gamaarachchi <[email protected]>
  27. * 2018 Jeff Daily <[email protected]>
  28. */
  29. #if !defined(SIMDE__SSE2_H)
  30. #if !defined(SIMDE__SSE2_H)
  31. #define SIMDE__SSE2_H
  32. #endif
  33. #include "sse.h"
  34. #if defined(SIMDE_SSE2_NATIVE)
  35. #undef SIMDE_SSE2_NATIVE
  36. #endif
  37. #if defined(SIMDE_SSE2_FORCE_NATIVE)
  38. #define SIMDE_SSE2_NATIVE
  39. #elif defined(__SSE2__) && !defined(SIMDE_SSE2_NO_NATIVE) && \
  40. !defined(SIMDE_NO_NATIVE)
  41. #define SIMDE_SSE2_NATIVE
  42. #elif defined(__ARM_NEON) && !defined(SIMDE_SSE2_NO_NEON) && \
  43. !defined(SIMDE_NO_NEON)
  44. #define SIMDE_SSE2_NEON
  45. #endif
  46. #if defined(SIMDE_SSE2_NATIVE) && !defined(SIMDE_SSE_NATIVE)
  47. #if defined(SIMDE_SSE2_FORCE_NATIVE)
  48. #error Native SSE2 support requires native SSE support
  49. #else
  50. #warning Native SSE2 support requires native SSE support, disabling
  51. #undef SIMDE_SSE2_NATIVE
  52. #endif
  53. #elif defined(SIMDE_SSE2_NEON) && !defined(SIMDE_SSE_NEON)
  54. #warning SSE2 NEON support requires SSE NEON support, disabling
  55. #undef SIMDE_SSE_NEON
  56. #endif
  57. #if defined(SIMDE_SSE2_NATIVE)
  58. #include <emmintrin.h>
  59. #else
  60. #if defined(SIMDE_SSE2_NEON)
  61. #include <arm_neon.h>
  62. #endif
  63. #endif
  64. #include <stdint.h>
  65. #include <limits.h>
  66. #include <string.h>
  67. #define vreinterpretq_m128i_s32(v) \
  68. (simde__m128i) { .neon_i32 = v }
  69. #define vreinterpretq_m128i_u64(v) \
  70. (simde__m128i) { .neon_u64 = v }
  71. #define vreinterpretq_s32_m128i(a) a.neon_i32
  72. #define vreinterpretq_u64_m128i(a) a.neon_u64
  73. SIMDE__BEGIN_DECLS
  74. typedef SIMDE_ALIGN(16) union {
  75. #if defined(SIMDE__ENABLE_GCC_VEC_EXT)
  76. int8_t i8 __attribute__((__vector_size__(16), __may_alias__));
  77. int16_t i16 __attribute__((__vector_size__(16), __may_alias__));
  78. int32_t i32 __attribute__((__vector_size__(16), __may_alias__));
  79. int64_t i64 __attribute__((__vector_size__(16), __may_alias__));
  80. uint8_t u8 __attribute__((__vector_size__(16), __may_alias__));
  81. uint16_t u16 __attribute__((__vector_size__(16), __may_alias__));
  82. uint32_t u32 __attribute__((__vector_size__(16), __may_alias__));
  83. uint64_t u64 __attribute__((__vector_size__(16), __may_alias__));
  84. #if defined(SIMDE__HAVE_INT128)
  85. simde_int128 i128 __attribute__((__vector_size__(16), __may_alias__));
  86. simde_uint128 u128 __attribute__((__vector_size__(16), __may_alias__));
  87. #endif
  88. simde_float32 f32 __attribute__((__vector_size__(16), __may_alias__));
  89. simde_float64 f64 __attribute__((__vector_size__(16), __may_alias__));
  90. #else
  91. int8_t i8[16];
  92. int16_t i16[8];
  93. int32_t i32[4];
  94. int64_t i64[2];
  95. uint8_t u8[16];
  96. uint16_t u16[8];
  97. uint32_t u32[4];
  98. uint64_t u64[2];
  99. #if defined(SIMDE__HAVE_INT128)
  100. simde_int128 i128[1];
  101. simde_uint128 u128[1];
  102. #endif
  103. simde_float32 f32[4];
  104. simde_float64 f64[2];
  105. #endif
  106. #if defined(SIMDE_SSE2_NATIVE)
  107. __m128i n;
  108. #elif defined(SIMDE_SSE2_NEON)
  109. int8x16_t neon_i8;
  110. int16x8_t neon_i16;
  111. int32x4_t neon_i32;
  112. int64x2_t neon_i64;
  113. uint8x16_t neon_u8;
  114. uint16x8_t neon_u16;
  115. uint32x4_t neon_u32;
  116. uint64x2_t neon_u64;
  117. float32x4_t neon_f32;
  118. #if defined(SIMDE_ARCH_AMD64)
  119. float64x2_t neon_f64;
  120. #endif
  121. #endif
  122. } simde__m128i;
  123. typedef SIMDE_ALIGN(16) union {
  124. #if defined(SIMDE__ENABLE_GCC_VEC_EXT)
  125. int8_t i8 __attribute__((__vector_size__(16), __may_alias__));
  126. int16_t i16 __attribute__((__vector_size__(16), __may_alias__));
  127. int32_t i32 __attribute__((__vector_size__(16), __may_alias__));
  128. int64_t i64 __attribute__((__vector_size__(16), __may_alias__));
  129. uint8_t u8 __attribute__((__vector_size__(16), __may_alias__));
  130. uint16_t u16 __attribute__((__vector_size__(16), __may_alias__));
  131. uint32_t u32 __attribute__((__vector_size__(16), __may_alias__));
  132. uint64_t u64 __attribute__((__vector_size__(16), __may_alias__));
  133. simde_float32 f32 __attribute__((__vector_size__(16), __may_alias__));
  134. simde_float64 f64 __attribute__((__vector_size__(16), __may_alias__));
  135. #else
  136. int8_t i8[16];
  137. int16_t i16[8];
  138. int32_t i32[4];
  139. int64_t i64[2];
  140. uint8_t u8[16];
  141. uint16_t u16[8];
  142. uint32_t u32[4];
  143. uint64_t u64[2];
  144. simde_float32 f32[4];
  145. simde_float64 f64[2];
  146. #endif
  147. #if defined(SIMDE_SSE2_NATIVE)
  148. __m128d n;
  149. #elif defined(SIMDE_SSE2_NEON)
  150. int8x16_t neon_i8;
  151. int16x8_t neon_i16;
  152. int32x4_t neon_i32;
  153. int64x2_t neon_i64;
  154. uint8x16_t neon_u8;
  155. uint16x8_t neon_u16;
  156. uint32x4_t neon_u32;
  157. uint64x2_t neon_u64;
  158. float32x4_t neon_f32;
  159. #if defined(SIMDE_ARCH_AMD64)
  160. float64x2_t neon_f64;
  161. #endif
  162. #endif
  163. } simde__m128d;
  164. #if defined(SIMDE_SSE2_NATIVE)
  165. HEDLEY_STATIC_ASSERT(sizeof(__m128i) == sizeof(simde__m128i),
  166. "__m128i size doesn't match simde__m128i size");
  167. HEDLEY_STATIC_ASSERT(sizeof(__m128d) == sizeof(simde__m128d),
  168. "__m128d size doesn't match simde__m128d size");
  169. SIMDE__FUNCTION_ATTRIBUTES simde__m128i SIMDE__M128I_C(__m128i v)
  170. {
  171. simde__m128i r;
  172. r.n = v;
  173. return r;
  174. }
  175. SIMDE__FUNCTION_ATTRIBUTES simde__m128d SIMDE__M128D_C(__m128d v)
  176. {
  177. simde__m128d r;
  178. r.n = v;
  179. return r;
  180. }
  181. #elif defined(SIMDE_SSE_NEON)
  182. #define SIMDE__M128I_NEON_C(T, expr) \
  183. (simde__m128i) { .neon_##T = expr }
  184. #define SIMDE__M128D_NEON_C(T, expr) \
  185. (simde__m128d) { .neon_##T = expr }
  186. #endif
  187. HEDLEY_STATIC_ASSERT(16 == sizeof(simde__m128i), "simde__m128i size incorrect");
  188. HEDLEY_STATIC_ASSERT(16 == sizeof(simde__m128d), "simde__m128d size incorrect");
  189. SIMDE__FUNCTION_ATTRIBUTES
  190. simde__m128i simde_mm_add_epi8(simde__m128i a, simde__m128i b)
  191. {
  192. #if defined(SIMDE_SSE2_NATIVE)
  193. return SIMDE__M128I_C(_mm_add_epi8(a.n, b.n));
  194. #elif defined(SIMDE_SSE2_NEON)
  195. return SIMDE__M128I_NEON_C(i8, vaddq_s8(a.neon_i8, b.neon_i8));
  196. #else
  197. simde__m128i r;
  198. SIMDE__VECTORIZE
  199. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  200. r.i8[i] = a.i8[i] + b.i8[i];
  201. }
  202. return r;
  203. #endif
  204. }
  205. SIMDE__FUNCTION_ATTRIBUTES
  206. simde__m128i simde_mm_add_epi16(simde__m128i a, simde__m128i b)
  207. {
  208. #if defined(SIMDE_SSE2_NATIVE)
  209. return SIMDE__M128I_C(_mm_add_epi16(a.n, b.n));
  210. #elif defined(SIMDE_SSE2_NEON)
  211. return SIMDE__M128I_NEON_C(i16, vaddq_s16(a.neon_i16, b.neon_i16));
  212. #else
  213. simde__m128i r;
  214. SIMDE__VECTORIZE
  215. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  216. r.i16[i] = a.i16[i] + b.i16[i];
  217. }
  218. return r;
  219. #endif
  220. }
  221. SIMDE__FUNCTION_ATTRIBUTES
  222. simde__m128i simde_mm_add_epi32(simde__m128i a, simde__m128i b)
  223. {
  224. #if defined(SIMDE_SSE2_NATIVE)
  225. return SIMDE__M128I_C(_mm_add_epi32(a.n, b.n));
  226. #elif defined(SIMDE_SSE2_NEON)
  227. return SIMDE__M128I_NEON_C(i32, vaddq_s32(a.neon_i32, b.neon_i32));
  228. #else
  229. simde__m128i r;
  230. SIMDE__VECTORIZE
  231. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  232. r.i32[i] = a.i32[i] + b.i32[i];
  233. }
  234. return r;
  235. #endif
  236. }
  237. SIMDE__FUNCTION_ATTRIBUTES
  238. simde__m128i simde_mm_add_epi64(simde__m128i a, simde__m128i b)
  239. {
  240. #if defined(SIMDE_SSE2_NATIVE)
  241. return SIMDE__M128I_C(_mm_add_epi64(a.n, b.n));
  242. #elif defined(SIMDE_SSE2_NEON)
  243. return SIMDE__M128I_NEON_C(i64, vaddq_s64(a.neon_i64, b.neon_i64));
  244. #else
  245. simde__m128i r;
  246. SIMDE__VECTORIZE
  247. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  248. r.i64[i] = a.i64[i] + b.i64[i];
  249. }
  250. return r;
  251. #endif
  252. }
  253. SIMDE__FUNCTION_ATTRIBUTES
  254. simde__m128d simde_mm_add_pd(simde__m128d a, simde__m128d b)
  255. {
  256. #if defined(SIMDE_SSE2_NATIVE)
  257. return SIMDE__M128D_C(_mm_add_pd(a.n, b.n));
  258. #elif defined(SIMDE_SSE2_NEON) && defined(SIMDE_ARCH_AMD64)
  259. return SIMDE__M128I_NEON_C(f64, vaddq_f64(a.neon_f64, b.neon_f64));
  260. #else
  261. simde__m128d r;
  262. SIMDE__VECTORIZE
  263. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  264. r.f64[i] = a.f64[i] + b.f64[i];
  265. }
  266. return r;
  267. #endif
  268. }
  269. SIMDE__FUNCTION_ATTRIBUTES
  270. simde__m128d simde_mm_add_sd(simde__m128d a, simde__m128d b)
  271. {
  272. #if defined(SIMDE_SSE2_NATIVE)
  273. return SIMDE__M128D_C(_mm_add_sd(a.n, b.n));
  274. #else
  275. simde__m128d r;
  276. r.f64[0] = a.f64[0] + b.f64[0];
  277. r.f64[1] = a.f64[1];
  278. return r;
  279. #endif
  280. }
  281. SIMDE__FUNCTION_ATTRIBUTES
  282. simde__m64 simde_mm_add_si64(simde__m64 a, simde__m64 b)
  283. {
  284. #if defined(SIMDE_SSE2_NATIVE)
  285. return SIMDE__M64_C(_mm_add_si64(a.n, b.n));
  286. #elif defined(SIMDE_SSE2_NEON)
  287. return SIMDE__M64_NEON_C(i64, vadd_s64(a.neon_i64, b.neon_i64));
  288. #else
  289. simde__m64 r;
  290. r.i64[0] = a.i64[0] + b.i64[0];
  291. return r;
  292. #endif
  293. }
  294. SIMDE__FUNCTION_ATTRIBUTES
  295. simde__m128i simde_mm_adds_epi8(simde__m128i a, simde__m128i b)
  296. {
  297. #if defined(SIMDE_SSE2_NATIVE)
  298. return SIMDE__M128I_C(_mm_adds_epi8(a.n, b.n));
  299. #elif defined(SIMDE_SSE2_NEON)
  300. return SIMDE__M128I_NEON_C(i8, vqaddq_s8(a.neon_i8, b.neon_i8));
  301. #else
  302. simde__m128i r;
  303. SIMDE__VECTORIZE
  304. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  305. if ((((b.i8[i]) > 0) && ((a.i8[i]) > (INT8_MAX - (b.i8[i]))))) {
  306. r.i8[i] = INT8_MAX;
  307. } else if ((((b.i8[i]) < 0) &&
  308. ((a.i8[i]) < (INT8_MIN - (b.i8[i]))))) {
  309. r.i8[i] = INT8_MIN;
  310. } else {
  311. r.i8[i] = (a.i8[i]) + (b.i8[i]);
  312. }
  313. }
  314. return r;
  315. #endif
  316. }
  317. SIMDE__FUNCTION_ATTRIBUTES
  318. simde__m128i simde_mm_adds_epi16(simde__m128i a, simde__m128i b)
  319. {
  320. #if defined(SIMDE_SSE2_NATIVE)
  321. return SIMDE__M128I_C(_mm_adds_epi16(a.n, b.n));
  322. #elif defined(SIMDE_SSE2_NEON)
  323. return SIMDE__M128I_NEON_C(i16, vqaddq_s16(a.neon_i16, b.neon_i16));
  324. #else
  325. simde__m128i r;
  326. SIMDE__VECTORIZE
  327. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  328. if ((((b.i16[i]) > 0) &&
  329. ((a.i16[i]) > (INT16_MAX - (b.i16[i]))))) {
  330. r.i16[i] = INT16_MAX;
  331. } else if ((((b.i16[i]) < 0) &&
  332. ((a.i16[i]) < (INT16_MIN - (b.i16[i]))))) {
  333. r.i16[i] = INT16_MIN;
  334. } else {
  335. r.i16[i] = (a.i16[i]) + (b.i16[i]);
  336. }
  337. }
  338. return r;
  339. #endif
  340. }
  341. SIMDE__FUNCTION_ATTRIBUTES
  342. simde__m128i simde_mm_adds_epu8(simde__m128i a, simde__m128i b)
  343. {
  344. #if defined(SIMDE_SSE2_NATIVE)
  345. return SIMDE__M128I_C(_mm_adds_epu8(a.n, b.n));
  346. #elif defined(SIMDE_SSE2_NEON)
  347. return SIMDE__M128I_NEON_C(u8, vqaddq_u8(a.neon_u8, b.neon_u8));
  348. #else
  349. simde__m128i r;
  350. SIMDE__VECTORIZE
  351. for (size_t i = 0; i < (sizeof(r.u8) / sizeof(r.u8[0])); i++) {
  352. r.u8[i] = ((UINT8_MAX - a.u8[i]) > b.u8[i])
  353. ? (a.u8[i] + b.u8[i])
  354. : UINT8_MAX;
  355. }
  356. return r;
  357. #endif
  358. }
  359. SIMDE__FUNCTION_ATTRIBUTES
  360. simde__m128i simde_mm_adds_epu16(simde__m128i a, simde__m128i b)
  361. {
  362. #if defined(SIMDE_SSE2_NATIVE)
  363. return SIMDE__M128I_C(_mm_adds_epu16(a.n, b.n));
  364. #elif defined(SIMDE_SSE2_NEON)
  365. return SIMDE__M128I_NEON_C(u16, vqaddq_u16(a.neon_u16, b.neon_u16));
  366. #else
  367. simde__m128i r;
  368. SIMDE__VECTORIZE
  369. for (size_t i = 0; i < (sizeof(r.u16) / sizeof(r.u16[0])); i++) {
  370. r.u16[i] = ((UINT16_MAX - a.u16[i]) > b.u16[i])
  371. ? (a.u16[i] + b.u16[i])
  372. : UINT16_MAX;
  373. }
  374. return r;
  375. #endif
  376. }
  377. SIMDE__FUNCTION_ATTRIBUTES
  378. simde__m128d simde_mm_and_pd(simde__m128d a, simde__m128d b)
  379. {
  380. #if defined(SIMDE_SSE2_NATIVE)
  381. return SIMDE__M128D_C(_mm_and_pd(a.n, b.n));
  382. #elif defined(SIMDE_SSE2_NEON)
  383. return SIMDE__M128D_NEON_C(i32, vandq_s32(a.neon_i32, b.neon_i32));
  384. #else
  385. simde__m128d r;
  386. SIMDE__VECTORIZE
  387. for (size_t i = 0; i < (sizeof(r.u64) / sizeof(r.u64[0])); i++) {
  388. r.u64[i] = a.u64[i] & b.u64[i];
  389. }
  390. return r;
  391. #endif
  392. }
  393. SIMDE__FUNCTION_ATTRIBUTES
  394. simde__m128i simde_mm_and_si128(simde__m128i a, simde__m128i b)
  395. {
  396. #if defined(SIMDE_SSE2_NATIVE)
  397. return SIMDE__M128I_C(_mm_and_si128(a.n, b.n));
  398. #elif defined(SIMDE_SSE_NEON)
  399. return SIMDE__M128I_NEON_C(i32, vandq_s32(b.neon_i32, a.neon_i32));
  400. #else
  401. simde__m128i r;
  402. SIMDE__VECTORIZE
  403. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  404. r.i64[i] = a.i64[i] & b.i64[i];
  405. }
  406. return r;
  407. #endif
  408. }
  409. SIMDE__FUNCTION_ATTRIBUTES
  410. simde__m128d simde_mm_andnot_pd(simde__m128d a, simde__m128d b)
  411. {
  412. #if defined(SIMDE_SSE2_NATIVE)
  413. return SIMDE__M128D_C(_mm_andnot_pd(a.n, b.n));
  414. #elif defined(SIMDE_SSE2_NEON)
  415. return SIMDE__M128D_NEON_C(i32, vbicq_s32(a.neon_i32, b.neon_i32));
  416. #else
  417. simde__m128d r;
  418. SIMDE__VECTORIZE
  419. for (size_t i = 0; i < (sizeof(r.u64) / sizeof(r.u64[0])); i++) {
  420. r.u64[i] = ~a.u64[i] & b.u64[i];
  421. }
  422. return r;
  423. #endif
  424. }
  425. SIMDE__FUNCTION_ATTRIBUTES
  426. simde__m128i simde_mm_andnot_si128(simde__m128i a, simde__m128i b)
  427. {
  428. #if defined(SIMDE_SSE2_NATIVE)
  429. return SIMDE__M128I_C(_mm_andnot_si128(a.n, b.n));
  430. #elif defined(SIMDE_SSE2_NEON)
  431. return SIMDE__M128I_NEON_C(i32, vbicq_s32(b.neon_i32, a.neon_i32));
  432. #else
  433. simde__m128i r;
  434. SIMDE__VECTORIZE
  435. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  436. r.i64[i] = ~(a.i64[i]) & b.i64[i];
  437. }
  438. return r;
  439. #endif
  440. }
  441. SIMDE__FUNCTION_ATTRIBUTES
  442. simde__m128i simde_mm_avg_epu8(simde__m128i a, simde__m128i b)
  443. {
  444. #if defined(SIMDE_SSE2_NATIVE)
  445. return SIMDE__M128I_C(_mm_avg_epu8(a.n, b.n));
  446. #elif defined(SIMDE_SSE2_NEON)
  447. return SIMDE__M128I_NEON_C(u8, vrhaddq_u8(b.neon_u8, a.neon_u8));
  448. #else
  449. simde__m128i r;
  450. SIMDE__VECTORIZE
  451. for (size_t i = 0; i < (sizeof(r.u8) / sizeof(r.u8[0])); i++) {
  452. r.u8[i] = (a.u8[i] + b.u8[i] + 1) >> 1;
  453. }
  454. return r;
  455. #endif
  456. }
  457. SIMDE__FUNCTION_ATTRIBUTES
  458. simde__m128i simde_mm_avg_epu16(simde__m128i a, simde__m128i b)
  459. {
  460. #if defined(SIMDE_SSE2_NATIVE)
  461. return SIMDE__M128I_C(_mm_avg_epu16(a.n, b.n));
  462. #elif defined(SIMDE_SSE2_NEON)
  463. return SIMDE__M128I_NEON_C(u16, vrhaddq_u16(b.neon_u16, a.neon_u16));
  464. #else
  465. simde__m128i r;
  466. SIMDE__VECTORIZE
  467. for (size_t i = 0; i < (sizeof(r.u16) / sizeof(r.u16[0])); i++) {
  468. r.u16[i] = (a.u16[i] + b.u16[i] + 1) >> 1;
  469. }
  470. return r;
  471. #endif
  472. }
  473. SIMDE__FUNCTION_ATTRIBUTES
  474. simde__m128i simde_mm_bslli_si128(simde__m128i a, const int imm8)
  475. {
  476. simde__m128i r;
  477. if (HEDLEY_UNLIKELY(imm8 > 15)) {
  478. r.u64[0] = 0;
  479. r.u64[1] = 0;
  480. return r;
  481. }
  482. const int s = imm8 * 8;
  483. #if defined(SIMDE__HAVE_INT128)
  484. r.u128[0] = a.u128[0] << s;
  485. #else
  486. if (s < 64) {
  487. r.u64[0] = (a.u64[0] << s);
  488. r.u64[1] = (a.u64[1] << s) | (a.u64[0] >> (64 - s));
  489. } else {
  490. r.u64[0] = 0;
  491. r.u64[1] = a.u64[0] << (s - 64);
  492. }
  493. #endif
  494. return r;
  495. }
  496. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  497. #define simde_mm_bslli_si128(a, imm8) SIMDE__M128I_C(_mm_slli_si128(a.n, imm8))
  498. #elif defined(SIMDE_SSE2_NEON)
  499. #define simde_mm_bslli_si128(a, imm8) \
  500. SIMDE__M128I_NEON_C( \
  501. i8, \
  502. (((imm8) <= 0) ? ((a).neon_i8) \
  503. : (((imm8) > 15) ? (vdupq_n_s8(0)) \
  504. : (vextq_s8(vdupq_n_s8(0), \
  505. (a).neon_i8, \
  506. 16 - (imm8))))))
  507. #endif
  508. #define simde_mm_slli_si128(a, imm8) simde_mm_bslli_si128(a, imm8)
  509. SIMDE__FUNCTION_ATTRIBUTES
  510. simde__m128i simde_mm_bsrli_si128(simde__m128i a, const int imm8)
  511. {
  512. simde__m128i r;
  513. if (HEDLEY_UNLIKELY(imm8 > 15)) {
  514. r.u64[0] = 0;
  515. r.u64[1] = 0;
  516. return r;
  517. }
  518. const int s = imm8 * 8;
  519. #if defined(SIMDE__HAVE_INT128)
  520. r.u128[0] = a.u128[0] >> s;
  521. #else
  522. if (s < 64) {
  523. r.u64[0] = (a.u64[0] >> s) | (a.u64[1] << (64 - s));
  524. r.u64[1] = (a.u64[1] >> s);
  525. } else {
  526. r.u64[0] = a.u64[1] >> (s - 64);
  527. r.u64[1] = 0;
  528. }
  529. #endif
  530. return r;
  531. }
  532. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  533. #define simde_mm_bsrli_si128(a, imm8) SIMDE__M128I_C(_mm_srli_si128(a.n, imm8))
  534. #elif defined(SIMDE_SSE2_NEON)
  535. #define simde_mm_bsrli_si128(a, imm8) \
  536. SIMDE__M128I_NEON_C( \
  537. i8, \
  538. ((imm8) <= 0) \
  539. ? ((a).neon_i8) \
  540. : (((imm8) > 15) ? (vdupq_n_s8(0)) \
  541. : (vextq_s8((a).neon_i8, \
  542. vdupq_n_s8(0), (imm8)))))
  543. #endif
  544. #define simde_mm_srli_si128(a, imm8) simde_mm_bsrli_si128(a, imm8)
  545. SIMDE__FUNCTION_ATTRIBUTES
  546. void simde_mm_clflush(void const *p)
  547. {
  548. #if defined(SIMDE_SSE2_NATIVE)
  549. _mm_clflush(p);
  550. #else
  551. (void)p;
  552. #endif
  553. }
  554. SIMDE__FUNCTION_ATTRIBUTES
  555. int simde_mm_comieq_sd(simde__m128d a, simde__m128d b)
  556. {
  557. #if defined(SIMDE_SSE2_NATIVE)
  558. return _mm_comieq_sd(a.n, b.n);
  559. #else
  560. return a.f64[0] == b.f64[0];
  561. #endif
  562. }
  563. SIMDE__FUNCTION_ATTRIBUTES
  564. int simde_mm_comige_sd(simde__m128d a, simde__m128d b)
  565. {
  566. #if defined(SIMDE_SSE2_NATIVE)
  567. return _mm_comige_sd(a.n, b.n);
  568. #else
  569. return a.f64[0] >= b.f64[0];
  570. #endif
  571. }
  572. SIMDE__FUNCTION_ATTRIBUTES
  573. int simde_mm_comigt_sd(simde__m128d a, simde__m128d b)
  574. {
  575. #if defined(SIMDE_SSE2_NATIVE)
  576. return _mm_comigt_sd(a.n, b.n);
  577. #else
  578. return a.f64[0] > b.f64[0];
  579. #endif
  580. }
  581. SIMDE__FUNCTION_ATTRIBUTES
  582. int simde_mm_comile_sd(simde__m128d a, simde__m128d b)
  583. {
  584. #if defined(SIMDE_SSE2_NATIVE)
  585. return _mm_comile_sd(a.n, b.n);
  586. #else
  587. return a.f64[0] <= b.f64[0];
  588. #endif
  589. }
  590. SIMDE__FUNCTION_ATTRIBUTES
  591. int simde_mm_comilt_sd(simde__m128d a, simde__m128d b)
  592. {
  593. #if defined(SIMDE_SSE2_NATIVE)
  594. return _mm_comilt_sd(a.n, b.n);
  595. #else
  596. return a.f64[0] < b.f64[0];
  597. #endif
  598. }
  599. SIMDE__FUNCTION_ATTRIBUTES
  600. int simde_mm_comineq_sd(simde__m128d a, simde__m128d b)
  601. {
  602. #if defined(SIMDE_SSE2_NATIVE)
  603. return _mm_comineq_sd(a.n, b.n);
  604. #else
  605. return a.f64[0] != b.f64[0];
  606. #endif
  607. }
  608. SIMDE__FUNCTION_ATTRIBUTES
  609. simde__m128 simde_mm_castpd_ps(simde__m128d a)
  610. {
  611. #if defined(SIMDE_SSE2_NATIVE)
  612. return SIMDE__M128_C(_mm_castpd_ps(a.n));
  613. #else
  614. union {
  615. simde__m128d pd;
  616. simde__m128 ps;
  617. } r;
  618. r.pd = a;
  619. return r.ps;
  620. #endif
  621. }
  622. SIMDE__FUNCTION_ATTRIBUTES
  623. simde__m128i simde_mm_castpd_si128(simde__m128d a)
  624. {
  625. #if defined(SIMDE_SSE2_NATIVE)
  626. return SIMDE__M128I_C(_mm_castpd_si128(a.n));
  627. #else
  628. union {
  629. simde__m128d pd;
  630. simde__m128i si128;
  631. } r;
  632. r.pd = a;
  633. return r.si128;
  634. #endif
  635. }
  636. SIMDE__FUNCTION_ATTRIBUTES
  637. simde__m128d simde_mm_castps_pd(simde__m128 a)
  638. {
  639. #if defined(SIMDE_SSE2_NATIVE)
  640. return SIMDE__M128D_C(_mm_castps_pd(a.n));
  641. #else
  642. union {
  643. simde__m128 ps;
  644. simde__m128d pd;
  645. } r;
  646. r.ps = a;
  647. return r.pd;
  648. #endif
  649. }
  650. SIMDE__FUNCTION_ATTRIBUTES
  651. simde__m128i simde_mm_castps_si128(simde__m128 a)
  652. {
  653. #if defined(SIMDE_SSE2_NATIVE)
  654. return SIMDE__M128I_C(_mm_castps_si128(a.n));
  655. #elif defined(SIMDE_SSE2_NEON)
  656. return SIMDE__M128I_NEON_C(i32, a.neon_i32);
  657. #else
  658. union {
  659. simde__m128 ps;
  660. simde__m128i si128;
  661. } r;
  662. r.ps = a;
  663. return r.si128;
  664. #endif
  665. }
  666. SIMDE__FUNCTION_ATTRIBUTES
  667. simde__m128d simde_mm_castsi128_pd(simde__m128i a)
  668. {
  669. #if defined(SIMDE_SSE2_NATIVE)
  670. return SIMDE__M128D_C(_mm_castsi128_pd(a.n));
  671. #else
  672. union {
  673. simde__m128i si128;
  674. simde__m128d pd;
  675. } r;
  676. r.si128 = a;
  677. return r.pd;
  678. #endif
  679. }
  680. SIMDE__FUNCTION_ATTRIBUTES
  681. simde__m128 simde_mm_castsi128_ps(simde__m128i a)
  682. {
  683. #if defined(SIMDE_SSE2_NATIVE)
  684. return SIMDE__M128_C(_mm_castsi128_ps(a.n));
  685. #elif defined(SIMDE_SSE2_NEON)
  686. return SIMDE__M128_NEON_C(f32, a.neon_f32);
  687. #else
  688. union {
  689. simde__m128i si128;
  690. simde__m128 ps;
  691. } r;
  692. r.si128 = a;
  693. return r.ps;
  694. #endif
  695. }
  696. SIMDE__FUNCTION_ATTRIBUTES
  697. simde__m128i simde_mm_cmpeq_epi8(simde__m128i a, simde__m128i b)
  698. {
  699. #if defined(SIMDE_SSE2_NATIVE)
  700. return SIMDE__M128I_C(_mm_cmpeq_epi8(a.n, b.n));
  701. #elif defined(SIMDE_SSE2_NEON)
  702. return SIMDE__M128I_NEON_C(
  703. i8, vreinterpretq_s8_u8(vceqq_s8(b.neon_i8, a.neon_i8)));
  704. #else
  705. simde__m128i r;
  706. SIMDE__VECTORIZE
  707. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  708. r.i8[i] = (a.i8[i] == b.i8[i]) ? 0xff : 0x00;
  709. }
  710. return r;
  711. #endif
  712. }
  713. SIMDE__FUNCTION_ATTRIBUTES
  714. simde__m128i simde_mm_cmpeq_epi16(simde__m128i a, simde__m128i b)
  715. {
  716. #if defined(SIMDE_SSE2_NATIVE)
  717. return SIMDE__M128I_C(_mm_cmpeq_epi16(a.n, b.n));
  718. #elif defined(SIMDE_SSE2_NEON)
  719. return SIMDE__M128I_NEON_C(
  720. i16, vreinterpretq_s16_u16(vceqq_s16(b.neon_i16, a.neon_i16)));
  721. #else
  722. simde__m128i r;
  723. SIMDE__VECTORIZE
  724. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  725. r.i16[i] = (a.i16[i] == b.i16[i]) ? 0xffff : 0x0000;
  726. }
  727. return r;
  728. #endif
  729. }
  730. SIMDE__FUNCTION_ATTRIBUTES
  731. simde__m128i simde_mm_cmpeq_epi32(simde__m128i a, simde__m128i b)
  732. {
  733. #if defined(SIMDE_SSE2_NATIVE)
  734. return SIMDE__M128I_C(_mm_cmpeq_epi32(a.n, b.n));
  735. #elif defined(SIMDE_SSE2_NEON)
  736. return SIMDE__M128I_NEON_C(
  737. i32, vreinterpretq_s32_u32(vceqq_s32(b.neon_i32, a.neon_i32)));
  738. #else
  739. simde__m128i r;
  740. SIMDE__VECTORIZE
  741. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  742. r.i32[i] = (a.i32[i] == b.i32[i]) ? 0xffffffff : 0x00000000;
  743. }
  744. return r;
  745. #endif
  746. }
  747. SIMDE__FUNCTION_ATTRIBUTES
  748. simde__m128d simde_mm_cmpeq_pd(simde__m128d a, simde__m128d b)
  749. {
  750. #if defined(SIMDE_SSE2_NATIVE)
  751. return SIMDE__M128D_C(_mm_cmpeq_pd(a.n, b.n));
  752. #elif defined(SIMDE_SSE2_NEON)
  753. return SIMDE__M128D_NEON_C(
  754. i32, vreinterpretq_s32_u32(
  755. vceqq_s32(vreinterpretq_s32_f32(b.neon_f32),
  756. vreinterpretq_s32_f32(a.neon_f32))));
  757. #else
  758. simde__m128d r;
  759. SIMDE__VECTORIZE
  760. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  761. r.u64[i] = (a.f64[i] == b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  762. }
  763. return r;
  764. #endif
  765. }
  766. SIMDE__FUNCTION_ATTRIBUTES
  767. simde__m128d simde_mm_cmpeq_sd(simde__m128d a, simde__m128d b)
  768. {
  769. #if defined(SIMDE_SSE2_NATIVE)
  770. return SIMDE__M128D_C(_mm_cmpeq_sd(a.n, b.n));
  771. #else
  772. simde__m128d r;
  773. r.u64[0] = (a.f64[0] == b.f64[0]) ? ~UINT64_C(0) : 0;
  774. r.u64[1] = a.u64[1];
  775. return r;
  776. #endif
  777. }
  778. SIMDE__FUNCTION_ATTRIBUTES
  779. simde__m128d simde_mm_cmpneq_pd(simde__m128d a, simde__m128d b)
  780. {
  781. #if defined(SIMDE_SSE2_NATIVE)
  782. return SIMDE__M128D_C(_mm_cmpneq_pd(a.n, b.n));
  783. #elif defined(SIMDE_SSE2_NEON)
  784. return SIMDE__M128D_NEON_C(f32,
  785. vreinterpretq_f32_u16(vmvnq_u16(
  786. vceqq_s16(b.neon_i16, a.neon_i16))));
  787. #else
  788. simde__m128d r;
  789. SIMDE__VECTORIZE
  790. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  791. r.u64[i] = (a.f64[i] != b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  792. }
  793. return r;
  794. #endif
  795. }
  796. SIMDE__FUNCTION_ATTRIBUTES
  797. simde__m128d simde_mm_cmpneq_sd(simde__m128d a, simde__m128d b)
  798. {
  799. #if defined(SIMDE_SSE2_NATIVE)
  800. return SIMDE__M128D_C(_mm_cmpneq_sd(a.n, b.n));
  801. #else
  802. simde__m128d r;
  803. r.u64[0] = (a.f64[0] != b.f64[0]) ? ~UINT64_C(0) : UINT64_C(0);
  804. r.u64[1] = a.u64[1];
  805. return r;
  806. #endif
  807. }
  808. SIMDE__FUNCTION_ATTRIBUTES
  809. simde__m128i simde_mm_cmplt_epi8(simde__m128i a, simde__m128i b)
  810. {
  811. #if defined(SIMDE_SSE2_NATIVE)
  812. return SIMDE__M128I_C(_mm_cmplt_epi8(a.n, b.n));
  813. #elif defined(SIMDE_SSE2_NEON)
  814. return SIMDE__M128I_NEON_C(
  815. i8, vreinterpretq_s8_u8(vcltq_s8(a.neon_i8, b.neon_i8)));
  816. #else
  817. simde__m128i r;
  818. SIMDE__VECTORIZE
  819. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  820. r.i8[i] = (a.i8[i] < b.i8[i]) ? 0xff : 0x00;
  821. }
  822. return r;
  823. #endif
  824. }
  825. SIMDE__FUNCTION_ATTRIBUTES
  826. simde__m128i simde_mm_cmplt_epi16(simde__m128i a, simde__m128i b)
  827. {
  828. #if defined(SIMDE_SSE2_NATIVE)
  829. return SIMDE__M128I_C(_mm_cmplt_epi16(a.n, b.n));
  830. #elif defined(SIMDE_SSE2_NEON)
  831. return SIMDE__M128I_NEON_C(
  832. i16, vreinterpretq_s16_u16(vcltq_s16(a.neon_i16, b.neon_i16)));
  833. #else
  834. simde__m128i r;
  835. SIMDE__VECTORIZE
  836. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  837. r.i16[i] = (a.i16[i] < b.i16[i]) ? 0xffff : 0x0000;
  838. }
  839. return r;
  840. #endif
  841. }
  842. SIMDE__FUNCTION_ATTRIBUTES
  843. simde__m128i simde_mm_cmplt_epi32(simde__m128i a, simde__m128i b)
  844. {
  845. #if defined(SIMDE_SSE2_NATIVE)
  846. return SIMDE__M128I_C(_mm_cmplt_epi32(a.n, b.n));
  847. #elif defined(SIMDE_SSE2_NEON)
  848. return SIMDE__M128I_NEON_C(
  849. i32, vreinterpretq_s32_u32(vcltq_s32(a.neon_i32, b.neon_i32)));
  850. #else
  851. simde__m128i r;
  852. SIMDE__VECTORIZE
  853. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  854. r.i32[i] = (a.i32[i] < b.i32[i]) ? 0xffffffff : 0x00000000;
  855. }
  856. return r;
  857. #endif
  858. }
  859. SIMDE__FUNCTION_ATTRIBUTES
  860. simde__m128d simde_mm_cmplt_pd(simde__m128d a, simde__m128d b)
  861. {
  862. #if defined(SIMDE_SSE2_NATIVE)
  863. return SIMDE__M128D_C(_mm_cmplt_pd(a.n, b.n));
  864. #else
  865. simde__m128d r;
  866. SIMDE__VECTORIZE
  867. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  868. r.u64[i] = (a.f64[i] < b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  869. }
  870. return r;
  871. #endif
  872. }
  873. SIMDE__FUNCTION_ATTRIBUTES
  874. simde__m128d simde_mm_cmplt_sd(simde__m128d a, simde__m128d b)
  875. {
  876. #if defined(SIMDE_SSE2_NATIVE)
  877. return SIMDE__M128D_C(_mm_cmplt_sd(a.n, b.n));
  878. #else
  879. simde__m128d r;
  880. r.u64[0] = (a.f64[0] < b.f64[0]) ? ~UINT64_C(0) : UINT64_C(0);
  881. r.u64[1] = a.u64[1];
  882. return r;
  883. #endif
  884. }
  885. SIMDE__FUNCTION_ATTRIBUTES
  886. simde__m128d simde_mm_cmple_pd(simde__m128d a, simde__m128d b)
  887. {
  888. #if defined(SIMDE_SSE2_NATIVE)
  889. return SIMDE__M128D_C(_mm_cmple_pd(a.n, b.n));
  890. #else
  891. simde__m128d r;
  892. SIMDE__VECTORIZE
  893. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  894. r.u64[i] = (a.f64[i] <= b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  895. }
  896. return r;
  897. #endif
  898. }
  899. SIMDE__FUNCTION_ATTRIBUTES
  900. simde__m128d simde_mm_cmple_sd(simde__m128d a, simde__m128d b)
  901. {
  902. #if defined(SIMDE_SSE2_NATIVE)
  903. return SIMDE__M128D_C(_mm_cmple_sd(a.n, b.n));
  904. #else
  905. simde__m128d r;
  906. r.u64[0] = (a.f64[0] <= b.f64[0]) ? ~UINT64_C(0) : UINT64_C(0);
  907. r.u64[1] = a.u64[1];
  908. return r;
  909. #endif
  910. }
  911. SIMDE__FUNCTION_ATTRIBUTES
  912. simde__m128i simde_mm_cmpgt_epi8(simde__m128i a, simde__m128i b)
  913. {
  914. #if defined(SIMDE_SSE2_NATIVE)
  915. return SIMDE__M128I_C(_mm_cmpgt_epi8(a.n, b.n));
  916. #elif defined(SIMDE_SSE2_NEON)
  917. return SIMDE__M128I_NEON_C(
  918. i8, vreinterpretq_s8_u8(vcgtq_s8(a.neon_i8, b.neon_i8)));
  919. #else
  920. simde__m128i r;
  921. SIMDE__VECTORIZE
  922. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  923. r.i8[i] = (a.i8[i] > b.i8[i]) ? 0xff : 0x00;
  924. }
  925. return r;
  926. #endif
  927. }
  928. SIMDE__FUNCTION_ATTRIBUTES
  929. simde__m128i simde_mm_cmpgt_epi16(simde__m128i a, simde__m128i b)
  930. {
  931. #if defined(SIMDE_SSE2_NATIVE)
  932. return SIMDE__M128I_C(_mm_cmpgt_epi16(a.n, b.n));
  933. #elif defined(SIMDE_SSE2_NEON)
  934. return SIMDE__M128I_NEON_C(
  935. i16, vreinterpretq_s16_u16(vcgtq_s16(a.neon_i16, b.neon_i16)));
  936. #else
  937. simde__m128i r;
  938. SIMDE__VECTORIZE
  939. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  940. r.i16[i] = (a.i16[i] > b.i16[i]) ? 0xffff : 0x0000;
  941. }
  942. return r;
  943. #endif
  944. }
  945. SIMDE__FUNCTION_ATTRIBUTES
  946. simde__m128i simde_mm_cmpgt_epi32(simde__m128i a, simde__m128i b)
  947. {
  948. #if defined(SIMDE_SSE2_NATIVE)
  949. return SIMDE__M128I_C(_mm_cmpgt_epi32(a.n, b.n));
  950. #elif defined(SIMDE_SSE2_NEON)
  951. return SIMDE__M128I_NEON_C(
  952. i32, vreinterpretq_s32_u32(vcgtq_s32(a.neon_i32, b.neon_i32)));
  953. #else
  954. simde__m128i r;
  955. SIMDE__VECTORIZE
  956. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  957. r.i32[i] = (a.i32[i] > b.i32[i]) ? 0xffffffff : 0x00000000;
  958. }
  959. return r;
  960. #endif
  961. }
  962. SIMDE__FUNCTION_ATTRIBUTES
  963. simde__m128d simde_mm_cmpgt_pd(simde__m128d a, simde__m128d b)
  964. {
  965. #if defined(SIMDE_SSE2_NATIVE)
  966. return SIMDE__M128D_C(_mm_cmpgt_pd(a.n, b.n));
  967. #else
  968. simde__m128d r;
  969. SIMDE__VECTORIZE
  970. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  971. r.u64[i] = (a.f64[i] > b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  972. }
  973. return r;
  974. #endif
  975. }
  976. SIMDE__FUNCTION_ATTRIBUTES
  977. simde__m128d simde_mm_cmpgt_sd(simde__m128d a, simde__m128d b)
  978. {
  979. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  980. return SIMDE__M128D_C(_mm_cmpgt_sd(a.n, b.n));
  981. #else
  982. simde__m128d r;
  983. r.u64[0] = (a.f64[0] > b.f64[0]) ? ~UINT64_C(0) : UINT64_C(0);
  984. r.u64[1] = a.u64[1];
  985. return r;
  986. #endif
  987. }
  988. SIMDE__FUNCTION_ATTRIBUTES
  989. simde__m128d simde_mm_cmpge_pd(simde__m128d a, simde__m128d b)
  990. {
  991. #if defined(SIMDE_SSE2_NATIVE)
  992. return SIMDE__M128D_C(_mm_cmpge_pd(a.n, b.n));
  993. #else
  994. simde__m128d r;
  995. SIMDE__VECTORIZE
  996. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  997. r.u64[i] = (a.f64[i] >= b.f64[i]) ? ~UINT64_C(0) : UINT64_C(0);
  998. }
  999. return r;
  1000. #endif
  1001. }
  1002. SIMDE__FUNCTION_ATTRIBUTES
  1003. simde__m128d simde_mm_cmpge_sd(simde__m128d a, simde__m128d b)
  1004. {
  1005. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  1006. return SIMDE__M128D_C(_mm_cmpge_sd(a.n, b.n));
  1007. #else
  1008. simde__m128d r;
  1009. r.u64[0] = (a.f64[0] >= b.f64[0]) ? ~UINT64_C(0) : UINT64_C(0);
  1010. r.u64[1] = a.u64[1];
  1011. return r;
  1012. #endif
  1013. }
  1014. SIMDE__FUNCTION_ATTRIBUTES
  1015. simde__m128d simde_mm_cmpnge_pd(simde__m128d a, simde__m128d b)
  1016. {
  1017. #if defined(SIMDE_SSE2_NATIVE)
  1018. return SIMDE__M128D_C(_mm_cmpnge_pd(a.n, b.n));
  1019. #else
  1020. return simde_mm_cmplt_pd(a, b);
  1021. #endif
  1022. }
  1023. SIMDE__FUNCTION_ATTRIBUTES
  1024. simde__m128d simde_mm_cmpnge_sd(simde__m128d a, simde__m128d b)
  1025. {
  1026. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  1027. return SIMDE__M128D_C(_mm_cmpnge_sd(a.n, b.n));
  1028. #else
  1029. return simde_mm_cmplt_sd(a, b);
  1030. #endif
  1031. }
  1032. SIMDE__FUNCTION_ATTRIBUTES
  1033. simde__m128d simde_mm_cmpnlt_pd(simde__m128d a, simde__m128d b)
  1034. {
  1035. #if defined(SIMDE_SSE2_NATIVE)
  1036. return SIMDE__M128D_C(_mm_cmpnlt_pd(a.n, b.n));
  1037. #else
  1038. return simde_mm_cmpge_pd(a, b);
  1039. #endif
  1040. }
  1041. SIMDE__FUNCTION_ATTRIBUTES
  1042. simde__m128d simde_mm_cmpnlt_sd(simde__m128d a, simde__m128d b)
  1043. {
  1044. #if defined(SIMDE_SSE2_NATIVE)
  1045. return SIMDE__M128D_C(_mm_cmpnlt_sd(a.n, b.n));
  1046. #else
  1047. return simde_mm_cmpge_sd(a, b);
  1048. #endif
  1049. }
  1050. SIMDE__FUNCTION_ATTRIBUTES
  1051. simde__m128d simde_mm_cmpnle_pd(simde__m128d a, simde__m128d b)
  1052. {
  1053. #if defined(SIMDE_SSE2_NATIVE)
  1054. return SIMDE__M128D_C(_mm_cmpnle_pd(a.n, b.n));
  1055. #else
  1056. return simde_mm_cmpgt_pd(a, b);
  1057. #endif
  1058. }
  1059. SIMDE__FUNCTION_ATTRIBUTES
  1060. simde__m128d simde_mm_cmpnle_sd(simde__m128d a, simde__m128d b)
  1061. {
  1062. #if defined(SIMDE_SSE2_NATIVE)
  1063. return SIMDE__M128D_C(_mm_cmpnle_sd(a.n, b.n));
  1064. #else
  1065. return simde_mm_cmpgt_sd(a, b);
  1066. #endif
  1067. }
  1068. SIMDE__FUNCTION_ATTRIBUTES
  1069. simde__m128d simde_mm_cmpord_pd(simde__m128d a, simde__m128d b)
  1070. {
  1071. #if defined(SIMDE_SSE2_NATIVE)
  1072. return SIMDE__M128D_C(_mm_cmpord_pd(a.n, b.n));
  1073. #else
  1074. simde__m128d r;
  1075. SIMDE__VECTORIZE
  1076. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1077. r.u64[i] = (!isnan(a.f64[i]) && !isnan(b.f64[i])) ? ~UINT64_C(0)
  1078. : UINT64_C(0);
  1079. }
  1080. return r;
  1081. #endif
  1082. }
  1083. SIMDE__FUNCTION_ATTRIBUTES
  1084. simde__m128d simde_mm_cmpord_sd(simde__m128d a, simde__m128d b)
  1085. {
  1086. #if defined(SIMDE_SSE2_NATIVE)
  1087. return SIMDE__M128D_C(_mm_cmpord_sd(a.n, b.n));
  1088. #else
  1089. simde__m128d r;
  1090. r.u64[0] = (!isnan(a.f64[0]) && !isnan(b.f64[0])) ? ~UINT64_C(0)
  1091. : UINT64_C(0);
  1092. r.u64[1] = a.u64[1];
  1093. return r;
  1094. #endif
  1095. }
  1096. SIMDE__FUNCTION_ATTRIBUTES
  1097. simde__m128d simde_mm_cmpunord_pd(simde__m128d a, simde__m128d b)
  1098. {
  1099. #if defined(SIMDE_SSE2_NATIVE)
  1100. return SIMDE__M128D_C(_mm_cmpunord_pd(a.n, b.n));
  1101. #else
  1102. simde__m128d r;
  1103. SIMDE__VECTORIZE
  1104. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1105. r.u64[i] = (isnan(a.f64[i]) || isnan(b.f64[i])) ? ~UINT64_C(0)
  1106. : UINT64_C(0);
  1107. }
  1108. return r;
  1109. #endif
  1110. }
  1111. SIMDE__FUNCTION_ATTRIBUTES
  1112. simde__m128d simde_mm_cmpunord_sd(simde__m128d a, simde__m128d b)
  1113. {
  1114. #if defined(SIMDE_SSE2_NATIVE)
  1115. return SIMDE__M128D_C(_mm_cmpunord_sd(a.n, b.n));
  1116. #else
  1117. simde__m128d r;
  1118. r.u64[0] = (isnan(a.f64[0]) || isnan(b.f64[0])) ? ~UINT64_C(0)
  1119. : UINT64_C(0);
  1120. r.u64[1] = a.u64[1];
  1121. return r;
  1122. #endif
  1123. }
  1124. SIMDE__FUNCTION_ATTRIBUTES
  1125. simde__m128d simde_mm_cvtepi32_pd(simde__m128i a)
  1126. {
  1127. #if defined(SIMDE_SSE2_NATIVE)
  1128. return SIMDE__M128D_C(_mm_cvtepi32_pd(a.n));
  1129. #else
  1130. simde__m128d r;
  1131. SIMDE__VECTORIZE
  1132. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1133. r.f64[i] = (simde_float64)a.i32[i];
  1134. }
  1135. return r;
  1136. #endif
  1137. }
  1138. SIMDE__FUNCTION_ATTRIBUTES
  1139. simde__m128 simde_mm_cvtepi32_ps(simde__m128i a)
  1140. {
  1141. #if defined(SIMDE_SSE2_NATIVE)
  1142. return SIMDE__M128_C(_mm_cvtepi32_ps(a.n));
  1143. #elif defined(SIMDE_SSE2_NEON)
  1144. return SIMDE__M128_NEON_C(f32, vcvtq_f32_s32(a.neon_i32));
  1145. #else
  1146. simde__m128 r;
  1147. SIMDE__VECTORIZE
  1148. for (size_t i = 0; i < (sizeof(r.f32) / sizeof(r.f32[0])); i++) {
  1149. r.f32[i] = (simde_float32)a.i32[i];
  1150. }
  1151. return r;
  1152. #endif
  1153. }
  1154. SIMDE__FUNCTION_ATTRIBUTES
  1155. simde__m128i simde_mm_cvtpd_epi32(simde__m128d a)
  1156. {
  1157. #if defined(SIMDE_SSE2_NATIVE)
  1158. return SIMDE__M128I_C(_mm_cvtpd_epi32(a.n));
  1159. #else
  1160. simde__m128i r;
  1161. SIMDE__VECTORIZE
  1162. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1163. r.i32[i] = (int32_t)a.f64[i];
  1164. }
  1165. return r;
  1166. #endif
  1167. }
  1168. SIMDE__FUNCTION_ATTRIBUTES
  1169. simde__m64 simde_mm_cvtpd_pi32(simde__m128d a)
  1170. {
  1171. #if defined(SIMDE_SSE2_NATIVE)
  1172. return SIMDE__M64_C(_mm_cvtpd_pi32(a.n));
  1173. #else
  1174. simde__m64 r;
  1175. SIMDE__VECTORIZE
  1176. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  1177. r.i32[i] = (int32_t)a.f64[i];
  1178. }
  1179. return r;
  1180. #endif
  1181. }
  1182. SIMDE__FUNCTION_ATTRIBUTES
  1183. simde__m128 simde_mm_cvtpd_ps(simde__m128d a)
  1184. {
  1185. #if defined(SIMDE_SSE2_NATIVE)
  1186. return SIMDE__M128_C(_mm_cvtpd_ps(a.n));
  1187. #else
  1188. simde__m128 r;
  1189. SIMDE__VECTORIZE
  1190. for (size_t i = 0; i < (sizeof(a.f64) / sizeof(a.f64[0])); i++) {
  1191. r.f32[i] = (simde_float32)a.f64[i];
  1192. }
  1193. return r;
  1194. #endif
  1195. }
  1196. SIMDE__FUNCTION_ATTRIBUTES
  1197. simde__m128d simde_mm_cvtpi32_pd(simde__m64 a)
  1198. {
  1199. #if defined(SIMDE_SSE2_NATIVE)
  1200. return SIMDE__M128D_C(_mm_cvtpi32_pd(a.n));
  1201. #else
  1202. simde__m128d r;
  1203. SIMDE__VECTORIZE
  1204. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1205. r.f64[i] = (simde_float64)a.i32[i];
  1206. }
  1207. return r;
  1208. #endif
  1209. }
  1210. SIMDE__FUNCTION_ATTRIBUTES
  1211. simde__m128i simde_mm_cvtps_epi32(simde__m128 a)
  1212. {
  1213. #if defined(SIMDE_SSE2_NATIVE)
  1214. return SIMDE__M128I_C(_mm_cvtps_epi32(a.n));
  1215. #elif defined(SIMDE_SSE2_NEON)
  1216. /* The default rounding mode on SSE is 'round to even', which ArmV7
  1217. does not support! It is supported on ARMv8 however. */
  1218. #if defined(SIMDE_ARCH_AARCH64)
  1219. return SIMDE__M128I_NEON_C(i32, vcvtnq_s32_f32(a.neon_f32));
  1220. #else
  1221. uint32x4_t signmask = vdupq_n_u32(0x80000000);
  1222. float32x4_t half = vbslq_f32(signmask, a.neon_f32,
  1223. vdupq_n_f32(0.5f)); /* +/- 0.5 */
  1224. int32x4_t r_normal = vcvtq_s32_f32(
  1225. vaddq_f32(a.neon_f32, half)); /* round to integer: [a + 0.5]*/
  1226. int32x4_t r_trunc =
  1227. vcvtq_s32_f32(a.neon_f32); /* truncate to integer: [a] */
  1228. int32x4_t plusone = vshrq_n_s32(vnegq_s32(r_trunc), 31); /* 1 or 0 */
  1229. int32x4_t r_even = vbicq_s32(vaddq_s32(r_trunc, plusone),
  1230. vdupq_n_s32(1)); /* ([a] + {0,1}) & ~1 */
  1231. float32x4_t delta = vsubq_f32(
  1232. a.neon_f32,
  1233. vcvtq_f32_s32(r_trunc)); /* compute delta: delta = (a - [a]) */
  1234. uint32x4_t is_delta_half =
  1235. vceqq_f32(delta, half); /* delta == +/- 0.5 */
  1236. return SIMDE__M128I_NEON_C(i32,
  1237. vbslq_s32(is_delta_half, r_even, r_normal));
  1238. #endif
  1239. #else
  1240. simde__m128i r;
  1241. SIMDE__VECTORIZE
  1242. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  1243. r.i32[i] = (int32_t)a.f32[i];
  1244. }
  1245. return r;
  1246. #endif
  1247. }
  1248. SIMDE__FUNCTION_ATTRIBUTES
  1249. simde__m128d simde_mm_cvtps_pd(simde__m128 a)
  1250. {
  1251. #if defined(SIMDE_SSE2_NATIVE)
  1252. return SIMDE__M128D_C(_mm_cvtps_pd(a.n));
  1253. #else
  1254. simde__m128d r;
  1255. SIMDE__VECTORIZE
  1256. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1257. r.f64[i] = a.f32[i];
  1258. }
  1259. return r;
  1260. #endif
  1261. }
  1262. SIMDE__FUNCTION_ATTRIBUTES
  1263. double simde_mm_cvtsd_f64(simde__m128d a)
  1264. {
  1265. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  1266. return _mm_cvtsd_f64(a.n);
  1267. #else
  1268. return a.f64[0];
  1269. #endif
  1270. }
  1271. SIMDE__FUNCTION_ATTRIBUTES
  1272. int32_t simde_mm_cvtsd_si32(simde__m128d a)
  1273. {
  1274. #if defined(SIMDE_SSE2_NATIVE)
  1275. return _mm_cvtsd_si32(a.n);
  1276. #else
  1277. return (int32_t)a.f64[0];
  1278. #endif
  1279. }
  1280. SIMDE__FUNCTION_ATTRIBUTES
  1281. int32_t simde_mm_cvtsd_si64(simde__m128d a)
  1282. {
  1283. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  1284. #if defined(__PGI)
  1285. return _mm_cvtsd_si64x(a.n);
  1286. #else
  1287. return _mm_cvtsd_si64(a.n);
  1288. #endif
  1289. #else
  1290. return (int32_t)a.f64[0];
  1291. #endif
  1292. }
  1293. #define simde_mm_cvtsd_si64x(a) simde_mm_cvtsd_si64(a)
  1294. SIMDE__FUNCTION_ATTRIBUTES
  1295. simde__m128 simde_mm_cvtsd_ss(simde__m128 a, simde__m128d b)
  1296. {
  1297. #if defined(SIMDE_SSE2_NATIVE)
  1298. return SIMDE__M128_C(_mm_cvtsd_ss(a.n, b.n));
  1299. #else
  1300. simde__m128 r;
  1301. r.f32[0] = (simde_float32)b.f64[0];
  1302. SIMDE__VECTORIZE
  1303. for (size_t i = 1; i < (sizeof(r) / sizeof(r.i32[0])); i++) {
  1304. r.i32[i] = a.i32[i];
  1305. }
  1306. return r;
  1307. #endif
  1308. }
  1309. SIMDE__FUNCTION_ATTRIBUTES
  1310. int32_t simde_mm_cvtsi128_si32(simde__m128i a)
  1311. {
  1312. #if defined(SIMDE_SSE2_NATIVE)
  1313. return _mm_cvtsi128_si32(a.n);
  1314. #elif defined(SIMDE_SSE2_NEON)
  1315. return vgetq_lane_s32(a.neon_i32, 0);
  1316. #else
  1317. return a.i32[0];
  1318. #endif
  1319. }
  1320. SIMDE__FUNCTION_ATTRIBUTES
  1321. int64_t simde_mm_cvtsi128_si64(simde__m128i a)
  1322. {
  1323. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  1324. #if defined(__PGI)
  1325. return _mm_cvtsi128_si64x(a.n);
  1326. #else
  1327. return _mm_cvtsi128_si64(a.n);
  1328. #endif
  1329. #else
  1330. return a.i64[0];
  1331. #endif
  1332. }
  1333. #define simde_mm_cvtsi128_si64x(a) simde_mm_cvtsi128_si64(a)
  1334. SIMDE__FUNCTION_ATTRIBUTES
  1335. simde__m128d simde_mm_cvtsi32_sd(simde__m128d a, int32_t b)
  1336. {
  1337. #if defined(SIMDE_SSE2_NATIVE)
  1338. return SIMDE__M128D_C(_mm_cvtsi32_sd(a.n, b));
  1339. #else
  1340. simde__m128d r;
  1341. r.f64[0] = (simde_float64)b;
  1342. r.i64[1] = a.i64[1];
  1343. return r;
  1344. #endif
  1345. }
  1346. SIMDE__FUNCTION_ATTRIBUTES
  1347. simde__m128i simde_mm_cvtsi32_si128(int32_t a)
  1348. {
  1349. simde__m128i r;
  1350. #if defined(SIMDE_SSE2_NATIVE)
  1351. r.n = _mm_cvtsi32_si128(a);
  1352. #elif defined(SIMDE_SSE2_NEON)
  1353. r.neon_i32 = vsetq_lane_s32(a, vdupq_n_s32(0), 0);
  1354. #else
  1355. r.i32[0] = a;
  1356. r.i32[1] = 0;
  1357. r.i32[2] = 0;
  1358. r.i32[3] = 0;
  1359. #endif
  1360. return r;
  1361. }
  1362. SIMDE__FUNCTION_ATTRIBUTES
  1363. simde__m128d simde_mm_cvtsi64_sd(simde__m128d a, int32_t b)
  1364. {
  1365. simde__m128d r;
  1366. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  1367. #if !defined(__PGI)
  1368. r.n = _mm_cvtsi64_sd(a.n, b);
  1369. #else
  1370. r.n = _mm_cvtsi64x_sd(a.n, b);
  1371. #endif
  1372. #else
  1373. r.f64[0] = (simde_float64)b;
  1374. r.f64[1] = a.f64[1];
  1375. #endif
  1376. return r;
  1377. }
  1378. #define simde_mm_cvtsi64x_sd(a, b) simde_mm_cvtsi64(a, b)
  1379. SIMDE__FUNCTION_ATTRIBUTES
  1380. simde__m128i simde_mm_cvtsi64_si128(int64_t a)
  1381. {
  1382. simde__m128i r;
  1383. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  1384. #if !defined(__PGI)
  1385. r.n = _mm_cvtsi64_si128(a);
  1386. #else
  1387. r.n = _mm_cvtsi64x_si128(a);
  1388. #endif
  1389. #else
  1390. r.i64[0] = a;
  1391. r.i64[1] = 0;
  1392. #endif
  1393. return r;
  1394. }
  1395. #define simde_mm_cvtsi64x_si128(a) simde_mm_cvtsi64_si128(a)
  1396. SIMDE__FUNCTION_ATTRIBUTES
  1397. simde__m128d simde_mm_cvtss_sd(simde__m128d a, simde__m128 b)
  1398. {
  1399. simde__m128d r;
  1400. #if defined(SIMDE_SSE2_NATIVE)
  1401. r.n = _mm_cvtss_sd(a.n, b.n);
  1402. #else
  1403. r.f64[0] = b.f32[0];
  1404. r.i64[1] = a.i64[1];
  1405. #endif
  1406. return r;
  1407. }
  1408. SIMDE__FUNCTION_ATTRIBUTES
  1409. simde__m128i simde_mm_cvttpd_epi32(simde__m128d a)
  1410. {
  1411. simde__m128i r;
  1412. #if defined(SIMDE_SSE2_NATIVE)
  1413. r.n = _mm_cvttpd_epi32(a.n);
  1414. #else
  1415. for (size_t i = 0; i < (sizeof(a.f64) / sizeof(a.f64[0])); i++) {
  1416. r.i32[i] = (int32_t)trunc(a.f64[i]);
  1417. }
  1418. #endif
  1419. return r;
  1420. }
  1421. SIMDE__FUNCTION_ATTRIBUTES
  1422. simde__m64 simde_mm_cvttpd_pi32(simde__m128d a)
  1423. {
  1424. simde__m64 r;
  1425. #if defined(SIMDE_SSE2_NATIVE)
  1426. r.n = _mm_cvttpd_pi32(a.n);
  1427. #else
  1428. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  1429. r.i32[i] = (int32_t)trunc(a.f64[i]);
  1430. }
  1431. #endif
  1432. return r;
  1433. }
  1434. SIMDE__FUNCTION_ATTRIBUTES
  1435. simde__m128i simde_mm_cvttps_epi32(simde__m128 a)
  1436. {
  1437. simde__m128i r;
  1438. #if defined(SIMDE_SSE2_NATIVE)
  1439. r.n = _mm_cvttps_epi32(a.n);
  1440. #elif defined(SIMDE_SSE2_NEON)
  1441. r.neon_i32 = vcvtq_s32_f32(a.neon_f32);
  1442. #else
  1443. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  1444. r.i32[i] = (int32_t)truncf(a.f32[i]);
  1445. }
  1446. #endif
  1447. return r;
  1448. }
  1449. SIMDE__FUNCTION_ATTRIBUTES
  1450. int32_t simde_mm_cvttsd_si32(simde__m128d a)
  1451. {
  1452. #if defined(SIMDE_SSE2_NATIVE)
  1453. return _mm_cvttsd_si32(a.n);
  1454. #else
  1455. return (int32_t)trunc(a.f64[0]);
  1456. #endif
  1457. }
  1458. SIMDE__FUNCTION_ATTRIBUTES
  1459. int64_t simde_mm_cvttsd_si64(simde__m128d a)
  1460. {
  1461. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  1462. #if !defined(__PGI)
  1463. return _mm_cvttsd_si64(a.n);
  1464. #else
  1465. return _mm_cvttsd_si64x(a.n);
  1466. #endif
  1467. #else
  1468. return (int64_t)trunc(a.f64[0]);
  1469. #endif
  1470. }
  1471. #define simde_mm_cvttsd_si64x(a) simde_mm_cvttsd_si64(a)
  1472. SIMDE__FUNCTION_ATTRIBUTES
  1473. simde__m128d simde_mm_div_pd(simde__m128d a, simde__m128d b)
  1474. {
  1475. simde__m128d r;
  1476. #if defined(SIMDE_SSE2_NATIVE)
  1477. r.n = _mm_div_pd(a.n, b.n);
  1478. #else
  1479. SIMDE__VECTORIZE
  1480. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1481. r.f64[i] = a.f64[i] / b.f64[i];
  1482. }
  1483. #endif
  1484. return r;
  1485. }
  1486. SIMDE__FUNCTION_ATTRIBUTES
  1487. simde__m128d simde_mm_div_sd(simde__m128d a, simde__m128d b)
  1488. {
  1489. simde__m128d r;
  1490. #if defined(SIMDE_SSE2_NATIVE)
  1491. r.n = _mm_div_sd(a.n, b.n);
  1492. #else
  1493. r.f64[0] = a.f64[0] / b.f64[0];
  1494. r.f64[1] = a.f64[1];
  1495. #endif
  1496. return r;
  1497. }
  1498. SIMDE__FUNCTION_ATTRIBUTES
  1499. int32_t simde_mm_extract_epi16(simde__m128i a, const int imm8)
  1500. {
  1501. return a.u16[imm8 & 7];
  1502. }
  1503. #if defined(SIMDE_SSE2_NATIVE) && \
  1504. (!defined(SIMDE__REALLY_GCC) || HEDLEY_GCC_VERSION_CHECK(4, 6, 0))
  1505. #define simde_mm_extract_epi16(a, imm8) _mm_extract_epi16(a.n, imm8)
  1506. #elif defined(SIMDE_SSE2_NEON)
  1507. #define simde_mm_extract_epi16(a, imm8) \
  1508. (vgetq_lane_s16((a).neon_i16, (imm8)) & ((int32_t)UINT32_C(0x0000ffff)))
  1509. #endif
  1510. SIMDE__FUNCTION_ATTRIBUTES
  1511. simde__m128i simde_mm_insert_epi16(simde__m128i a, int32_t i, const int imm8)
  1512. {
  1513. a.u16[imm8 & 7] = (int16_t)i;
  1514. return a;
  1515. }
  1516. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  1517. #define simde_mm_insert_epi16(a, i, imm8) \
  1518. SIMDE__M128I_C(_mm_insert_epi16((a).n, (i), (imm8)))
  1519. #elif defined(SIMDE_SSE2_NEON)
  1520. #define simde_mm_insert_epi16(a, i, imm8) \
  1521. SIMDE__M128I_NEON_C(i16, vsetq_lane_s16((i), a.neon_i16, (imm8)))
  1522. #endif
  1523. SIMDE__FUNCTION_ATTRIBUTES
  1524. simde__m128d
  1525. simde_mm_load_pd(simde_float64 const mem_addr[HEDLEY_ARRAY_PARAM(2)])
  1526. {
  1527. simde__m128d r;
  1528. simde_assert_aligned(16, mem_addr);
  1529. #if defined(SIMDE_SSE2_NATIVE)
  1530. r.n = _mm_load_pd(mem_addr);
  1531. #elif defined(SIMDE_SSE2_NEON)
  1532. r.neon_u32 = vld1q_u32((uint32_t const *)mem_addr);
  1533. #else
  1534. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  1535. memcpy(&r, mem_addr, sizeof(r));
  1536. #endif
  1537. return r;
  1538. }
  1539. SIMDE__FUNCTION_ATTRIBUTES
  1540. simde__m128d simde_mm_load_pd1(simde_float64 const *mem_addr)
  1541. {
  1542. simde__m128d r;
  1543. #if defined(SIMDE_SSE2_NATIVE)
  1544. r.n = _mm_load_pd1(mem_addr);
  1545. #else
  1546. r.f64[0] = *mem_addr;
  1547. r.f64[1] = *mem_addr;
  1548. #endif
  1549. return r;
  1550. }
  1551. #define simde_mm_load1_pd(mem_addr) simde_mm_load_pd1(mem_addr)
  1552. SIMDE__FUNCTION_ATTRIBUTES
  1553. simde__m128d simde_mm_load_sd(simde_float64 const *mem_addr)
  1554. {
  1555. simde__m128d r;
  1556. #if defined(SIMDE_SSE2_NATIVE)
  1557. r.n = _mm_load_sd(mem_addr);
  1558. #else
  1559. memcpy(&r, mem_addr, sizeof(simde_float64));
  1560. r.u64[1] = 0;
  1561. #endif
  1562. return r;
  1563. }
  1564. SIMDE__FUNCTION_ATTRIBUTES
  1565. simde__m128i simde_mm_load_si128(simde__m128i const *mem_addr)
  1566. {
  1567. simde__m128i r;
  1568. simde_assert_aligned(16, mem_addr);
  1569. #if defined(SIMDE_SSE2_NATIVE)
  1570. r.n = _mm_load_si128(&(mem_addr->n));
  1571. #elif defined(SIMDE_SSE2_NEON)
  1572. r.neon_i32 = vld1q_s32((int32_t const *)mem_addr);
  1573. #else
  1574. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  1575. memcpy(&r, mem_addr, sizeof(r));
  1576. #endif
  1577. return r;
  1578. }
  1579. SIMDE__FUNCTION_ATTRIBUTES
  1580. simde__m128d simde_mm_loadh_pd(simde__m128d a, simde_float64 const *mem_addr)
  1581. {
  1582. simde__m128d r;
  1583. #if defined(SIMDE_SSE2_NATIVE)
  1584. r.n = _mm_loadh_pd(a.n, mem_addr);
  1585. #else
  1586. simde_float64 t;
  1587. memcpy(&t, mem_addr, sizeof(t));
  1588. r.f64[0] = a.f64[0];
  1589. r.f64[1] = t;
  1590. #endif
  1591. return r;
  1592. }
  1593. SIMDE__FUNCTION_ATTRIBUTES
  1594. simde__m128i simde_mm_loadl_epi64(simde__m128i const *mem_addr)
  1595. {
  1596. simde__m128i r;
  1597. #if defined(SIMDE_SSE2_NATIVE)
  1598. r.n = _mm_loadl_epi64(&mem_addr->n);
  1599. #elif defined(SIMDE_SSE2_NEON)
  1600. r.neon_i32 = vcombine_s32(vld1_s32((int32_t const *)mem_addr),
  1601. vcreate_s32(0));
  1602. #else
  1603. r.u64[0] = mem_addr->u64[0];
  1604. r.u64[1] = 0;
  1605. #endif
  1606. return r;
  1607. }
  1608. SIMDE__FUNCTION_ATTRIBUTES
  1609. simde__m128d simde_mm_loadl_pd(simde__m128d a, simde_float64 const *mem_addr)
  1610. {
  1611. simde__m128d r;
  1612. #if defined(SIMDE_SSE2_NATIVE)
  1613. r.n = _mm_loadl_pd(a.n, mem_addr);
  1614. #else
  1615. memcpy(&r, mem_addr, sizeof(simde_float64));
  1616. r.u64[1] = a.u64[1];
  1617. #endif
  1618. return r;
  1619. }
  1620. SIMDE__FUNCTION_ATTRIBUTES
  1621. simde__m128d
  1622. simde_mm_loadr_pd(simde_float64 const mem_addr[HEDLEY_ARRAY_PARAM(2)])
  1623. {
  1624. simde__m128d r;
  1625. simde_assert_aligned(16, mem_addr);
  1626. #if defined(SIMDE_SSE2_NATIVE)
  1627. r.n = _mm_loadr_pd(mem_addr);
  1628. #else
  1629. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  1630. r.f64[0] = mem_addr[1];
  1631. r.f64[1] = mem_addr[0];
  1632. #endif
  1633. return r;
  1634. }
  1635. SIMDE__FUNCTION_ATTRIBUTES
  1636. simde__m128d
  1637. simde_mm_loadu_pd(simde_float64 const mem_addr[HEDLEY_ARRAY_PARAM(2)])
  1638. {
  1639. simde__m128d r;
  1640. #if defined(SIMDE_SSE2_NATIVE)
  1641. r.n = _mm_loadu_pd(mem_addr);
  1642. #else
  1643. simde_float64 l, h;
  1644. memcpy(&l, &mem_addr[0], sizeof(l));
  1645. memcpy(&h, &mem_addr[1], sizeof(h));
  1646. r.f64[0] = l;
  1647. r.f64[1] = h;
  1648. #endif
  1649. return r;
  1650. }
  1651. SIMDE__FUNCTION_ATTRIBUTES
  1652. simde__m128i simde_mm_loadu_si128(simde__m128i const *mem_addr)
  1653. {
  1654. simde__m128i r;
  1655. #if defined(SIMDE_SSE2_NATIVE)
  1656. r.n = _mm_loadu_si128(&((*mem_addr).n));
  1657. #elif defined(SIMDE_SSE2_NEON)
  1658. r.neon_i32 = vld1q_s32((int32_t const *)mem_addr);
  1659. #else
  1660. memcpy(&r, mem_addr, sizeof(r));
  1661. #endif
  1662. return r;
  1663. }
  1664. SIMDE__FUNCTION_ATTRIBUTES
  1665. simde__m128i simde_mm_madd_epi16(simde__m128i a, simde__m128i b)
  1666. {
  1667. simde__m128i r;
  1668. #if defined(SIMDE_SSE2_NATIVE)
  1669. r.n = _mm_madd_epi16(a.n, b.n);
  1670. #elif defined(SIMDE_SSE2_NEON)
  1671. int32x4_t pl =
  1672. vmull_s16(vget_low_s16(a.neon_i16), vget_low_s16(b.neon_i16));
  1673. int32x4_t ph =
  1674. vmull_s16(vget_high_s16(a.neon_i16), vget_high_s16(b.neon_i16));
  1675. int32x2_t rl = vpadd_s32(vget_low_s32(pl), vget_high_s32(pl));
  1676. int32x2_t rh = vpadd_s32(vget_low_s32(ph), vget_high_s32(ph));
  1677. r.neon_i32 = vcombine_s32(rl, rh);
  1678. #else
  1679. SIMDE__VECTORIZE
  1680. for (size_t i = 0; i < (sizeof(r) / sizeof(r.i16[0])); i += 2) {
  1681. r.i32[i / 2] =
  1682. (a.i16[i] * b.i16[i]) + (a.i16[i + 1] * b.i16[i + 1]);
  1683. }
  1684. #endif
  1685. return r;
  1686. }
  1687. SIMDE__FUNCTION_ATTRIBUTES
  1688. void simde_mm_maskmoveu_si128(simde__m128i a, simde__m128i mask,
  1689. int8_t mem_addr[HEDLEY_ARRAY_PARAM(16)])
  1690. {
  1691. #if defined(SIMDE_SSE2_NATIVE)
  1692. _mm_maskmoveu_si128(a.n, mask.n, (char *)mem_addr);
  1693. #else
  1694. for (size_t i = 0; i < 16; i++) {
  1695. if (mask.u8[i] & 0x80) {
  1696. mem_addr[i] = a.i8[i];
  1697. }
  1698. }
  1699. #endif
  1700. }
  1701. SIMDE__FUNCTION_ATTRIBUTES
  1702. int32_t simde_mm_movemask_epi8(simde__m128i a)
  1703. {
  1704. #if defined(SIMDE_SSE2_NATIVE)
  1705. return _mm_movemask_epi8(a.n);
  1706. #elif defined(SIMDE_SSE2_NEON)
  1707. uint8x16_t input = a.neon_u8;
  1708. SIMDE_ALIGN(16)
  1709. static const int8_t xr[8] = {-7, -6, -5, -4, -3, -2, -1, 0};
  1710. uint8x8_t mask_and = vdup_n_u8(0x80);
  1711. int8x8_t mask_shift = vld1_s8(xr);
  1712. uint8x8_t lo = vget_low_u8(input);
  1713. uint8x8_t hi = vget_high_u8(input);
  1714. lo = vand_u8(lo, mask_and);
  1715. lo = vshl_u8(lo, mask_shift);
  1716. hi = vand_u8(hi, mask_and);
  1717. hi = vshl_u8(hi, mask_shift);
  1718. lo = vpadd_u8(lo, lo);
  1719. lo = vpadd_u8(lo, lo);
  1720. lo = vpadd_u8(lo, lo);
  1721. hi = vpadd_u8(hi, hi);
  1722. hi = vpadd_u8(hi, hi);
  1723. hi = vpadd_u8(hi, hi);
  1724. return ((hi[0] << 8) | (lo[0] & 0xFF));
  1725. #else
  1726. int32_t r = 0;
  1727. SIMDE__VECTORIZE_REDUCTION(| : r)
  1728. for (size_t i = 0; i < 16; i++) {
  1729. r |= (a.u8[15 - i] >> 7) << (15 - i);
  1730. }
  1731. return r;
  1732. #endif
  1733. }
  1734. SIMDE__FUNCTION_ATTRIBUTES
  1735. int32_t simde_mm_movemask_pd(simde__m128d a)
  1736. {
  1737. #if defined(SIMDE_SSE2_NATIVE)
  1738. return _mm_movemask_pd(a.n);
  1739. #else
  1740. int32_t r = 0;
  1741. SIMDE__VECTORIZE
  1742. for (size_t i = 0; i < (sizeof(a.u64) / sizeof(a.u64[0])); i++) {
  1743. r |= (a.u64[i] >> 63) << i;
  1744. }
  1745. return r;
  1746. #endif
  1747. }
  1748. SIMDE__FUNCTION_ATTRIBUTES
  1749. simde__m64 simde_mm_movepi64_pi64(simde__m128i a)
  1750. {
  1751. simde__m64 r;
  1752. #if defined(SIMDE_SSE2_NATIVE)
  1753. r.n = _mm_movepi64_pi64(a.n);
  1754. #else
  1755. r.i64[0] = a.i64[0];
  1756. #endif
  1757. return r;
  1758. }
  1759. SIMDE__FUNCTION_ATTRIBUTES
  1760. simde__m128i simde_mm_movpi64_epi64(simde__m64 a)
  1761. {
  1762. simde__m128i r;
  1763. #if defined(SIMDE_SSE2_NATIVE)
  1764. r.n = _mm_movpi64_epi64(a.n);
  1765. #else
  1766. r.i64[0] = a.i64[0];
  1767. r.i64[1] = 0;
  1768. #endif
  1769. return r;
  1770. }
  1771. SIMDE__FUNCTION_ATTRIBUTES
  1772. simde__m128i simde_mm_min_epi16(simde__m128i a, simde__m128i b)
  1773. {
  1774. simde__m128i r;
  1775. #if defined(SIMDE_SSE2_NATIVE)
  1776. r.n = _mm_min_epi16(a.n, b.n);
  1777. #elif defined(SIMDE_SSE2_NEON)
  1778. r.neon_i16 = vminq_s16(a.neon_i16, b.neon_i16);
  1779. #else
  1780. SIMDE__VECTORIZE
  1781. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  1782. r.i16[i] = (a.i16[i] < b.i16[i]) ? a.i16[i] : b.i16[i];
  1783. }
  1784. #endif
  1785. return r;
  1786. }
  1787. SIMDE__FUNCTION_ATTRIBUTES
  1788. simde__m128i simde_mm_min_epu8(simde__m128i a, simde__m128i b)
  1789. {
  1790. simde__m128i r;
  1791. #if defined(SIMDE_SSE2_NATIVE)
  1792. r.n = _mm_min_epu8(a.n, b.n);
  1793. #elif defined(SIMDE_SSE2_NEON)
  1794. r.neon_u8 = vminq_u8(a.neon_u8, b.neon_u8);
  1795. #else
  1796. SIMDE__VECTORIZE
  1797. for (size_t i = 0; i < (sizeof(r.u8) / sizeof(r.u8[0])); i++) {
  1798. r.u8[i] = (a.u8[i] < b.u8[i]) ? a.u8[i] : b.u8[i];
  1799. }
  1800. #endif
  1801. return r;
  1802. }
  1803. SIMDE__FUNCTION_ATTRIBUTES
  1804. simde__m128d simde_mm_min_pd(simde__m128d a, simde__m128d b)
  1805. {
  1806. simde__m128d r;
  1807. #if defined(SIMDE_SSE2_NATIVE)
  1808. r.n = _mm_min_pd(a.n, b.n);
  1809. #else
  1810. SIMDE__VECTORIZE
  1811. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1812. r.f64[i] = (a.f64[i] < b.f64[i]) ? a.f64[i] : b.f64[i];
  1813. }
  1814. #endif
  1815. return r;
  1816. }
  1817. SIMDE__FUNCTION_ATTRIBUTES
  1818. simde__m128d simde_mm_min_sd(simde__m128d a, simde__m128d b)
  1819. {
  1820. simde__m128d r;
  1821. #if defined(SIMDE_SSE2_NATIVE)
  1822. r.n = _mm_min_sd(a.n, b.n);
  1823. #else
  1824. r.f64[0] = (a.f64[0] < b.f64[0]) ? a.f64[0] : b.f64[0];
  1825. r.f64[1] = a.f64[1];
  1826. #endif
  1827. return r;
  1828. }
  1829. SIMDE__FUNCTION_ATTRIBUTES
  1830. simde__m128i simde_mm_max_epi16(simde__m128i a, simde__m128i b)
  1831. {
  1832. simde__m128i r;
  1833. #if defined(SIMDE_SSE2_NATIVE)
  1834. r.n = _mm_max_epi16(a.n, b.n);
  1835. #elif defined(SIMDE_SSE2_NEON)
  1836. r.neon_i16 = vmaxq_s16(a.neon_i16, b.neon_i16);
  1837. #else
  1838. SIMDE__VECTORIZE
  1839. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  1840. r.i16[i] = (a.i16[i] > b.i16[i]) ? a.i16[i] : b.i16[i];
  1841. }
  1842. #endif
  1843. return r;
  1844. }
  1845. SIMDE__FUNCTION_ATTRIBUTES
  1846. simde__m128i simde_mm_max_epu8(simde__m128i a, simde__m128i b)
  1847. {
  1848. simde__m128i r;
  1849. #if defined(SIMDE_SSE2_NATIVE)
  1850. r.n = _mm_max_epu8(a.n, b.n);
  1851. #elif defined(SIMDE_SSE2_NEON)
  1852. r.neon_u8 = vmaxq_u8(a.neon_u8, b.neon_u8);
  1853. #else
  1854. SIMDE__VECTORIZE
  1855. for (size_t i = 0; i < (sizeof(r.u8) / sizeof(r.u8[0])); i++) {
  1856. r.u8[i] = (a.u8[i] > b.u8[i]) ? a.u8[i] : b.u8[i];
  1857. }
  1858. #endif
  1859. return r;
  1860. }
  1861. SIMDE__FUNCTION_ATTRIBUTES
  1862. simde__m128d simde_mm_max_pd(simde__m128d a, simde__m128d b)
  1863. {
  1864. simde__m128d r;
  1865. #if defined(SIMDE_SSE2_NATIVE)
  1866. r.n = _mm_max_pd(a.n, b.n);
  1867. #else
  1868. SIMDE__VECTORIZE
  1869. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1870. r.f64[i] = (a.f64[i] > b.f64[i]) ? a.f64[i] : b.f64[i];
  1871. }
  1872. #endif
  1873. return r;
  1874. }
  1875. SIMDE__FUNCTION_ATTRIBUTES
  1876. simde__m128d simde_mm_max_sd(simde__m128d a, simde__m128d b)
  1877. {
  1878. simde__m128d r;
  1879. #if defined(SIMDE_SSE2_NATIVE)
  1880. r.n = _mm_max_sd(a.n, b.n);
  1881. #else
  1882. r.f64[0] = (a.f64[0] > b.f64[0]) ? a.f64[0] : b.f64[0];
  1883. r.f64[1] = a.f64[1];
  1884. #endif
  1885. return r;
  1886. }
  1887. SIMDE__FUNCTION_ATTRIBUTES
  1888. simde__m128i simde_mm_move_epi64(simde__m128i a)
  1889. {
  1890. simde__m128i r;
  1891. #if defined(SIMDE_SSE2_NATIVE)
  1892. r.n = _mm_move_epi64(a.n);
  1893. #elif defined(SIMDE_SSE2_NEON)
  1894. r.neon_i64 = vsetq_lane_s64(0, a.neon_i64, 1);
  1895. #else
  1896. r.i64[0] = a.i64[0];
  1897. r.i64[1] = 0;
  1898. #endif
  1899. return r;
  1900. }
  1901. SIMDE__FUNCTION_ATTRIBUTES
  1902. simde__m128d simde_mm_move_sd(simde__m128d a, simde__m128d b)
  1903. {
  1904. simde__m128d r;
  1905. #if defined(SIMDE_SSE2_NATIVE)
  1906. r.n = _mm_move_sd(a.n, b.n);
  1907. #else
  1908. r.f64[0] = b.f64[0];
  1909. r.f64[1] = a.f64[1];
  1910. #endif
  1911. return r;
  1912. }
  1913. SIMDE__FUNCTION_ATTRIBUTES
  1914. simde__m128i simde_mm_mul_epu32(simde__m128i a, simde__m128i b)
  1915. {
  1916. simde__m128i r;
  1917. #if defined(SIMDE_SSE2_NATIVE)
  1918. r.n = _mm_mul_epu32(a.n, b.n);
  1919. #else
  1920. SIMDE__VECTORIZE
  1921. for (size_t i = 0; i < (sizeof(r.u64) / sizeof(r.u64[0])); i++) {
  1922. r.u64[i] = ((uint64_t)a.u32[i * 2]) * ((uint64_t)b.u32[i * 2]);
  1923. }
  1924. #endif
  1925. return r;
  1926. }
  1927. SIMDE__FUNCTION_ATTRIBUTES
  1928. simde__m128i simde_x_mm_mul_epi64(simde__m128i a, simde__m128i b)
  1929. {
  1930. simde__m128i r;
  1931. SIMDE__VECTORIZE
  1932. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  1933. r.i64[i] = a.i64[i] * b.i64[i];
  1934. }
  1935. return r;
  1936. }
  1937. SIMDE__FUNCTION_ATTRIBUTES
  1938. simde__m128i simde_x_mm_mod_epi64(simde__m128i a, simde__m128i b)
  1939. {
  1940. simde__m128i r;
  1941. SIMDE__VECTORIZE
  1942. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  1943. r.i64[i] = a.i64[i] % b.i64[i];
  1944. }
  1945. return r;
  1946. }
  1947. SIMDE__FUNCTION_ATTRIBUTES
  1948. simde__m128d simde_mm_mul_pd(simde__m128d a, simde__m128d b)
  1949. {
  1950. simde__m128d r;
  1951. #if defined(SIMDE_SSE2_NATIVE)
  1952. r.n = _mm_mul_pd(a.n, b.n);
  1953. #else
  1954. SIMDE__VECTORIZE
  1955. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  1956. r.f64[i] = a.f64[i] * b.f64[i];
  1957. }
  1958. #endif
  1959. return r;
  1960. }
  1961. SIMDE__FUNCTION_ATTRIBUTES
  1962. simde__m128d simde_mm_mul_sd(simde__m128d a, simde__m128d b)
  1963. {
  1964. simde__m128d r;
  1965. #if defined(SIMDE_SSE2_NATIVE)
  1966. r.n = _mm_mul_sd(a.n, b.n);
  1967. #else
  1968. r.f64[0] = a.f64[0] * b.f64[0];
  1969. r.f64[1] = a.f64[1];
  1970. #endif
  1971. return r;
  1972. }
  1973. SIMDE__FUNCTION_ATTRIBUTES
  1974. simde__m64 simde_mm_mul_su32(simde__m64 a, simde__m64 b)
  1975. {
  1976. simde__m64 r;
  1977. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  1978. r.n = _mm_mul_su32(a.n, b.n);
  1979. #else
  1980. r.u64[0] = ((uint64_t)a.u32[0]) * ((uint64_t)b.u32[0]);
  1981. #endif
  1982. return r;
  1983. }
  1984. SIMDE__FUNCTION_ATTRIBUTES
  1985. simde__m128i simde_mm_mulhi_epi16(simde__m128i a, simde__m128i b)
  1986. {
  1987. simde__m128i r;
  1988. #if defined(SIMDE_SSE2_NATIVE)
  1989. r.n = _mm_mulhi_epi16(a.n, b.n);
  1990. #elif defined(SIMDE_SSE2_NEON)
  1991. int16x4_t a3210 = vget_low_s16(a.neon_i16);
  1992. int16x4_t b3210 = vget_low_s16(b.neon_i16);
  1993. int32x4_t ab3210 = vmull_s16(a3210, b3210); /* 3333222211110000 */
  1994. int16x4_t a7654 = vget_high_s16(a.neon_i16);
  1995. int16x4_t b7654 = vget_high_s16(b.neon_i16);
  1996. int32x4_t ab7654 = vmull_s16(a7654, b7654); /* 7777666655554444 */
  1997. uint16x8x2_t rv = vuzpq_u16(vreinterpretq_u16_s32(ab3210),
  1998. vreinterpretq_u16_s32(ab7654));
  1999. r.neon_u16 = rv.val[1];
  2000. #else
  2001. SIMDE__VECTORIZE
  2002. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2003. r.u16[i] = (uint16_t)(((uint32_t)(((int32_t)a.i16[i]) *
  2004. ((int32_t)b.i16[i]))) >>
  2005. 16);
  2006. }
  2007. #endif
  2008. return r;
  2009. }
  2010. SIMDE__FUNCTION_ATTRIBUTES
  2011. simde__m128i simde_mm_mulhi_epu16(simde__m128i a, simde__m128i b)
  2012. {
  2013. simde__m128i r;
  2014. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  2015. r.n = _mm_mulhi_epu16(a.n, b.n);
  2016. #else
  2017. SIMDE__VECTORIZE
  2018. for (size_t i = 0; i < (sizeof(r.u16) / sizeof(r.u16[0])); i++) {
  2019. r.u16[i] = (uint16_t)(
  2020. (((uint32_t)a.u16[i]) * ((uint32_t)b.u16[i])) >> 16);
  2021. }
  2022. #endif
  2023. return r;
  2024. }
  2025. SIMDE__FUNCTION_ATTRIBUTES
  2026. simde__m128i simde_mm_mullo_epi16(simde__m128i a, simde__m128i b)
  2027. {
  2028. simde__m128i r;
  2029. #if defined(SIMDE_SSE2_NATIVE)
  2030. r.n = _mm_mullo_epi16(a.n, b.n);
  2031. #elif defined(SIMDE_SSE2_NEON)
  2032. r.neon_i16 = vmulq_s16(a.neon_i16, b.neon_i16);
  2033. #else
  2034. SIMDE__VECTORIZE
  2035. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2036. r.u16[i] = (uint16_t)(((uint32_t)(((int32_t)a.i16[i]) *
  2037. ((int32_t)b.i16[i]))) &
  2038. 0xffff);
  2039. }
  2040. #endif
  2041. return r;
  2042. }
  2043. SIMDE__FUNCTION_ATTRIBUTES
  2044. simde__m128d simde_mm_or_pd(simde__m128d a, simde__m128d b)
  2045. {
  2046. simde__m128d r;
  2047. #if defined(SIMDE_SSE2_NATIVE)
  2048. r.n = _mm_or_pd(a.n, b.n);
  2049. #else
  2050. SIMDE__VECTORIZE
  2051. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2052. r.i64[i] = a.i64[i] | b.i64[i];
  2053. }
  2054. #endif
  2055. return r;
  2056. }
  2057. SIMDE__FUNCTION_ATTRIBUTES
  2058. simde__m128i simde_mm_or_si128(simde__m128i a, simde__m128i b)
  2059. {
  2060. simde__m128i r;
  2061. #if defined(SIMDE_SSE2_NATIVE)
  2062. r.n = _mm_or_si128(a.n, b.n);
  2063. #elif defined(SIMDE_SSE2_NEON)
  2064. r.neon_i32 = vorrq_s32(a.neon_i32, b.neon_i32);
  2065. #else
  2066. SIMDE__VECTORIZE
  2067. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2068. r.i64[i] = a.i64[i] | b.i64[i];
  2069. }
  2070. #endif
  2071. return r;
  2072. }
  2073. SIMDE__FUNCTION_ATTRIBUTES
  2074. simde__m128i simde_mm_packs_epi16(simde__m128i a, simde__m128i b)
  2075. {
  2076. simde__m128i r;
  2077. #if defined(SIMDE_SSE2_NATIVE)
  2078. r.n = _mm_packs_epi16(a.n, b.n);
  2079. #elif defined(SIMDE_SSE2_NEON)
  2080. r.neon_i8 = vcombine_s8(vqmovn_s16(a.neon_i16), vqmovn_s16(b.neon_i16));
  2081. #else
  2082. SIMDE__VECTORIZE
  2083. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2084. r.i8[i] = (a.i16[i] > INT8_MAX)
  2085. ? INT8_MAX
  2086. : ((a.i16[i] < INT8_MIN)
  2087. ? INT8_MIN
  2088. : ((int8_t)a.i16[i]));
  2089. r.i8[i + 8] = (b.i16[i] > INT8_MAX)
  2090. ? INT8_MAX
  2091. : ((b.i16[i] < INT8_MIN)
  2092. ? INT8_MIN
  2093. : ((int8_t)b.i16[i]));
  2094. }
  2095. #endif
  2096. return r;
  2097. }
  2098. SIMDE__FUNCTION_ATTRIBUTES
  2099. simde__m128i simde_mm_packs_epi32(simde__m128i a, simde__m128i b)
  2100. {
  2101. simde__m128i r;
  2102. #if defined(SIMDE_SSE2_NATIVE)
  2103. r.n = _mm_packs_epi32(a.n, b.n);
  2104. #elif defined(SIMDE_SSE2_NEON)
  2105. r.neon_i16 =
  2106. vcombine_s16(vqmovn_s32(a.neon_i32), vqmovn_s32(b.neon_i32));
  2107. #else
  2108. SIMDE__VECTORIZE
  2109. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2110. r.i16[i] = (a.i32[i] > INT16_MAX)
  2111. ? INT16_MAX
  2112. : ((a.i32[i] < INT16_MIN)
  2113. ? INT16_MIN
  2114. : ((int16_t)a.i32[i]));
  2115. r.i16[i + 4] = (b.i32[i] > INT16_MAX)
  2116. ? INT16_MAX
  2117. : ((b.i32[i] < INT16_MIN)
  2118. ? INT16_MIN
  2119. : ((int16_t)b.i32[i]));
  2120. }
  2121. #endif
  2122. return r;
  2123. }
  2124. SIMDE__FUNCTION_ATTRIBUTES
  2125. simde__m128i simde_mm_packus_epi16(simde__m128i a, simde__m128i b)
  2126. {
  2127. simde__m128i r;
  2128. #if defined(SIMDE_SSE2_NATIVE)
  2129. r.n = _mm_packus_epi16(a.n, b.n);
  2130. #elif defined(SIMDE_SSE2_NEON)
  2131. r.neon_u8 =
  2132. vcombine_u8(vqmovun_s16(a.neon_i16), vqmovun_s16(b.neon_i16));
  2133. #else
  2134. SIMDE__VECTORIZE
  2135. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2136. r.u8[i] = (a.i16[i] > UINT8_MAX)
  2137. ? UINT8_MAX
  2138. : ((a.i16[i] < 0) ? 0 : ((int8_t)a.i16[i]));
  2139. r.u8[i + 8] =
  2140. (b.i16[i] > UINT8_MAX)
  2141. ? UINT8_MAX
  2142. : ((b.i16[i] < 0) ? 0 : ((int8_t)b.i16[i]));
  2143. }
  2144. #endif
  2145. return r;
  2146. }
  2147. SIMDE__FUNCTION_ATTRIBUTES
  2148. void simde_mm_pause(void)
  2149. {
  2150. #if defined(SIMDE_SSE2_NATIVE)
  2151. _mm_pause();
  2152. #endif
  2153. }
  2154. SIMDE__FUNCTION_ATTRIBUTES
  2155. simde__m128i simde_mm_sad_epu8(simde__m128i a, simde__m128i b)
  2156. {
  2157. simde__m128i r;
  2158. #if defined(SIMDE_SSE2_NATIVE)
  2159. r.n = _mm_sad_epu8(a.n, b.n);
  2160. #else
  2161. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2162. uint16_t tmp = 0;
  2163. SIMDE__VECTORIZE_REDUCTION(+ : tmp)
  2164. for (size_t j = 0; j < ((sizeof(r.u8) / sizeof(r.u8[0])) / 2);
  2165. j++) {
  2166. const size_t e = j + (i * 8);
  2167. tmp += (a.u8[e] > b.u8[e]) ? (a.u8[e] - b.u8[e])
  2168. : (b.u8[e] - a.u8[e]);
  2169. }
  2170. r.i64[i] = tmp;
  2171. }
  2172. #endif
  2173. return r;
  2174. }
  2175. SIMDE__FUNCTION_ATTRIBUTES
  2176. simde__m128i simde_mm_set_epi8(int8_t e15, int8_t e14, int8_t e13, int8_t e12,
  2177. int8_t e11, int8_t e10, int8_t e9, int8_t e8,
  2178. int8_t e7, int8_t e6, int8_t e5, int8_t e4,
  2179. int8_t e3, int8_t e2, int8_t e1, int8_t e0)
  2180. {
  2181. simde__m128i r;
  2182. #if defined(SIMDE_SSE2_NATIVE)
  2183. r.n = _mm_set_epi8(e15, e14, e13, e12, e11, e10, e9, e8, e7, e6, e5, e4,
  2184. e3, e2, e1, e0);
  2185. #else
  2186. r.i8[0] = e0;
  2187. r.i8[1] = e1;
  2188. r.i8[2] = e2;
  2189. r.i8[3] = e3;
  2190. r.i8[4] = e4;
  2191. r.i8[5] = e5;
  2192. r.i8[6] = e6;
  2193. r.i8[7] = e7;
  2194. r.i8[8] = e8;
  2195. r.i8[9] = e9;
  2196. r.i8[10] = e10;
  2197. r.i8[11] = e11;
  2198. r.i8[12] = e12;
  2199. r.i8[13] = e13;
  2200. r.i8[14] = e14;
  2201. r.i8[15] = e15;
  2202. #endif
  2203. return r;
  2204. }
  2205. SIMDE__FUNCTION_ATTRIBUTES
  2206. simde__m128i simde_mm_set_epi16(int16_t e7, int16_t e6, int16_t e5, int16_t e4,
  2207. int16_t e3, int16_t e2, int16_t e1, int16_t e0)
  2208. {
  2209. simde__m128i r;
  2210. #if defined(SIMDE_SSE2_NATIVE)
  2211. r.n = _mm_set_epi16(e7, e6, e5, e4, e3, e2, e1, e0);
  2212. #elif defined(SIMDE_SSE2_NEON)
  2213. SIMDE_ALIGN(16) int16_t data[8] = {e0, e1, e2, e3, e4, e5, e6, e7};
  2214. r.neon_i16 = vld1q_s16(data);
  2215. #else
  2216. r.i16[0] = e0;
  2217. r.i16[1] = e1;
  2218. r.i16[2] = e2;
  2219. r.i16[3] = e3;
  2220. r.i16[4] = e4;
  2221. r.i16[5] = e5;
  2222. r.i16[6] = e6;
  2223. r.i16[7] = e7;
  2224. #endif
  2225. return r;
  2226. }
  2227. SIMDE__FUNCTION_ATTRIBUTES
  2228. simde__m128i simde_mm_set_epi32(int32_t e3, int32_t e2, int32_t e1, int32_t e0)
  2229. {
  2230. simde__m128i r;
  2231. #if defined(SIMDE_SSE2_NATIVE)
  2232. r.n = _mm_set_epi32(e3, e2, e1, e0);
  2233. #elif defined(SIMDE_SSE2_NEON)
  2234. SIMDE_ALIGN(16) int32_t data[4] = {e0, e1, e2, e3};
  2235. r.neon_i32 = vld1q_s32(data);
  2236. #else
  2237. r.i32[0] = e0;
  2238. r.i32[1] = e1;
  2239. r.i32[2] = e2;
  2240. r.i32[3] = e3;
  2241. #endif
  2242. return r;
  2243. }
  2244. SIMDE__FUNCTION_ATTRIBUTES
  2245. simde__m128i simde_mm_set_epi64(simde__m64 e1, simde__m64 e0)
  2246. {
  2247. simde__m128i r;
  2248. #if defined(SIMDE_SSE2_NATIVE)
  2249. r.n = _mm_set_epi64(e1.n, e0.n);
  2250. #else
  2251. r.i64[0] = e0.i64[0];
  2252. r.i64[1] = e1.i64[0];
  2253. #endif
  2254. return r;
  2255. }
  2256. SIMDE__FUNCTION_ATTRIBUTES
  2257. simde__m128i simde_mm_set_epi64x(int64_t e1, int64_t e0)
  2258. {
  2259. simde__m128i r;
  2260. #if defined(SIMDE_SSE2_NATIVE)
  2261. r.n = _mm_set_epi64x(e1, e0);
  2262. #elif defined(SIMDE_SSE2_NEON)
  2263. r = SIMDE__M128I_NEON_C(i64,
  2264. vcombine_s64(vdup_n_s64(e0), vdup_n_s64(e1)));
  2265. #else
  2266. r.i64[0] = e0;
  2267. r.i64[1] = e1;
  2268. #endif
  2269. return r;
  2270. }
  2271. SIMDE__FUNCTION_ATTRIBUTES
  2272. simde__m128i simde_x_mm_set_epu8(uint8_t e15, uint8_t e14, uint8_t e13,
  2273. uint8_t e12, uint8_t e11, uint8_t e10,
  2274. uint8_t e9, uint8_t e8, uint8_t e7, uint8_t e6,
  2275. uint8_t e5, uint8_t e4, uint8_t e3, uint8_t e2,
  2276. uint8_t e1, uint8_t e0)
  2277. {
  2278. simde__m128i r;
  2279. r.u8[0] = e0;
  2280. r.u8[1] = e1;
  2281. r.u8[2] = e2;
  2282. r.u8[3] = e3;
  2283. r.u8[4] = e4;
  2284. r.u8[5] = e5;
  2285. r.u8[6] = e6;
  2286. r.u8[7] = e7;
  2287. r.u8[8] = e8;
  2288. r.u8[9] = e9;
  2289. r.u8[10] = e10;
  2290. r.u8[11] = e11;
  2291. r.u8[12] = e12;
  2292. r.u8[13] = e13;
  2293. r.u8[14] = e14;
  2294. r.u8[15] = e15;
  2295. return r;
  2296. }
  2297. SIMDE__FUNCTION_ATTRIBUTES
  2298. simde__m128i simde_x_mm_set_epu16(uint16_t e7, uint16_t e6, uint16_t e5,
  2299. uint16_t e4, uint16_t e3, uint16_t e2,
  2300. uint16_t e1, uint16_t e0)
  2301. {
  2302. simde__m128i r;
  2303. r.u16[0] = e0;
  2304. r.u16[1] = e1;
  2305. r.u16[2] = e2;
  2306. r.u16[3] = e3;
  2307. r.u16[4] = e4;
  2308. r.u16[5] = e5;
  2309. r.u16[6] = e6;
  2310. r.u16[7] = e7;
  2311. return r;
  2312. }
  2313. SIMDE__FUNCTION_ATTRIBUTES
  2314. simde__m128i simde_x_mm_set_epu32(uint32_t e3, uint32_t e2, uint32_t e1,
  2315. uint32_t e0)
  2316. {
  2317. simde__m128i r;
  2318. r.u32[0] = e0;
  2319. r.u32[1] = e1;
  2320. r.u32[2] = e2;
  2321. r.u32[3] = e3;
  2322. return r;
  2323. }
  2324. SIMDE__FUNCTION_ATTRIBUTES
  2325. simde__m128i simde_x_mm_set_epu64x(uint64_t e1, uint64_t e0)
  2326. {
  2327. simde__m128i r;
  2328. r.u64[0] = e0;
  2329. r.u64[1] = e1;
  2330. return r;
  2331. }
  2332. SIMDE__FUNCTION_ATTRIBUTES
  2333. simde__m128d simde_mm_set_pd(simde_float64 e1, simde_float64 e0)
  2334. {
  2335. simde__m128d r;
  2336. #if defined(SIMDE_SSE2_NATIVE)
  2337. r.n = _mm_set_pd(e1, e0);
  2338. #else
  2339. r.f64[0] = e0;
  2340. r.f64[1] = e1;
  2341. #endif
  2342. return r;
  2343. }
  2344. SIMDE__FUNCTION_ATTRIBUTES
  2345. simde__m128d simde_mm_set_pd1(simde_float64 a)
  2346. {
  2347. simde__m128d r;
  2348. #if defined(SIMDE_SSE2_NATIVE)
  2349. r.n = _mm_set1_pd(a);
  2350. #else
  2351. r.f64[0] = a;
  2352. r.f64[1] = a;
  2353. #endif
  2354. return r;
  2355. }
  2356. SIMDE__FUNCTION_ATTRIBUTES
  2357. simde__m128d simde_mm_set_sd(simde_float64 a)
  2358. {
  2359. simde__m128d r;
  2360. #if defined(SIMDE_SSE2_NATIVE)
  2361. r.n = _mm_set_sd(a);
  2362. #else
  2363. r.f64[0] = a;
  2364. r.u64[1] = 0;
  2365. #endif
  2366. return r;
  2367. }
  2368. SIMDE__FUNCTION_ATTRIBUTES
  2369. simde__m128i simde_mm_set1_epi8(int8_t a)
  2370. {
  2371. simde__m128i r;
  2372. #if defined(SIMDE_SSE2_NATIVE)
  2373. r.n = _mm_set1_epi8(a);
  2374. #elif defined(SIMDE_SSE2_NEON)
  2375. r.neon_i8 = vdupq_n_s8(a);
  2376. #else
  2377. SIMDE__VECTORIZE
  2378. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  2379. r.i8[i] = a;
  2380. }
  2381. #endif
  2382. return r;
  2383. }
  2384. SIMDE__FUNCTION_ATTRIBUTES
  2385. simde__m128i simde_mm_set1_epi16(int16_t a)
  2386. {
  2387. simde__m128i r;
  2388. #if defined(SIMDE_SSE2_NATIVE)
  2389. r.n = _mm_set1_epi16(a);
  2390. #elif defined(SIMDE_SSE2_NEON)
  2391. r.neon_i16 = vdupq_n_s16(a);
  2392. #else
  2393. SIMDE__VECTORIZE
  2394. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2395. r.i16[i] = a;
  2396. }
  2397. #endif
  2398. return r;
  2399. }
  2400. SIMDE__FUNCTION_ATTRIBUTES
  2401. simde__m128i simde_mm_set1_epi32(int32_t a)
  2402. {
  2403. simde__m128i r;
  2404. #if defined(SIMDE_SSE2_NATIVE)
  2405. r.n = _mm_set1_epi32(a);
  2406. #elif defined(SIMDE_SSE2_NEON)
  2407. r.neon_i32 = vdupq_n_s32(a);
  2408. #else
  2409. SIMDE__VECTORIZE
  2410. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2411. r.i32[i] = a;
  2412. }
  2413. #endif
  2414. return r;
  2415. }
  2416. SIMDE__FUNCTION_ATTRIBUTES
  2417. simde__m128i simde_mm_set1_epi64x(int64_t a)
  2418. {
  2419. simde__m128i r;
  2420. #if defined(SIMDE_SSE2_NATIVE)
  2421. r.n = _mm_set1_epi64x(a);
  2422. #elif defined(SIMDE_SSE2_NEON)
  2423. r.neon_i64 = vmovq_n_s64(a);
  2424. #else
  2425. SIMDE__VECTORIZE
  2426. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2427. r.i64[i] = a;
  2428. }
  2429. #endif
  2430. return r;
  2431. }
  2432. SIMDE__FUNCTION_ATTRIBUTES
  2433. simde__m128i simde_mm_set1_epi64(simde__m64 a)
  2434. {
  2435. simde__m128i r;
  2436. #if defined(SIMDE_SSE2_NATIVE)
  2437. r.n = _mm_set1_epi64(a.n);
  2438. #else
  2439. SIMDE__VECTORIZE
  2440. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2441. r.i64[i] = a.i64[0];
  2442. }
  2443. #endif
  2444. return r;
  2445. }
  2446. SIMDE__FUNCTION_ATTRIBUTES
  2447. simde__m128d simde_mm_set1_pd(simde_float64 a)
  2448. {
  2449. simde__m128d r;
  2450. #if defined(SIMDE_SSE2_NATIVE)
  2451. r.n = _mm_set1_pd(a);
  2452. #else
  2453. SIMDE__VECTORIZE
  2454. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2455. r.f64[i] = a;
  2456. }
  2457. #endif
  2458. return r;
  2459. }
  2460. SIMDE__FUNCTION_ATTRIBUTES
  2461. simde__m128i simde_mm_setr_epi8(int8_t e15, int8_t e14, int8_t e13, int8_t e12,
  2462. int8_t e11, int8_t e10, int8_t e9, int8_t e8,
  2463. int8_t e7, int8_t e6, int8_t e5, int8_t e4,
  2464. int8_t e3, int8_t e2, int8_t e1, int8_t e0)
  2465. {
  2466. simde__m128i r;
  2467. #if defined(SIMDE_SSE2_NATIVE)
  2468. r.n = _mm_setr_epi8(e15, e14, e13, e12, e11, e10, e9, e8, e7, e6, e5,
  2469. e4, e3, e2, e1, e0);
  2470. #elif defined(SIMDE_SSE2_NEON)
  2471. int8_t t[] = {e15, e14, e13, e12, e11, e10, e9, e8,
  2472. e7, e6, e5, e4, e3, e2, e1, e0};
  2473. r.neon_i8 = vld1q_s8(t);
  2474. #else
  2475. r.i8[0] = e15;
  2476. r.i8[1] = e14;
  2477. r.i8[2] = e13;
  2478. r.i8[3] = e12;
  2479. r.i8[4] = e11;
  2480. r.i8[5] = e10;
  2481. r.i8[6] = e9;
  2482. r.i8[7] = e8;
  2483. r.i8[8] = e7;
  2484. r.i8[9] = e6;
  2485. r.i8[10] = e5;
  2486. r.i8[11] = e4;
  2487. r.i8[12] = e3;
  2488. r.i8[13] = e2;
  2489. r.i8[14] = e1;
  2490. r.i8[15] = e0;
  2491. #endif
  2492. return r;
  2493. }
  2494. SIMDE__FUNCTION_ATTRIBUTES
  2495. simde__m128i simde_mm_setr_epi16(int16_t e7, int16_t e6, int16_t e5, int16_t e4,
  2496. int16_t e3, int16_t e2, int16_t e1, int16_t e0)
  2497. {
  2498. simde__m128i r;
  2499. #if defined(SIMDE_SSE2_NATIVE)
  2500. r.n = _mm_setr_epi16(e7, e6, e5, e4, e3, e2, e1, e0);
  2501. #elif defined(SIMDE_SSE2_NEON)
  2502. int16_t t[] = {e7, e6, e5, e4, e3, e2, e1, e0};
  2503. r.neon_i16 = vld1q_s16(t);
  2504. #else
  2505. r.i16[0] = e7;
  2506. r.i16[1] = e6;
  2507. r.i16[2] = e5;
  2508. r.i16[3] = e4;
  2509. r.i16[4] = e3;
  2510. r.i16[5] = e2;
  2511. r.i16[6] = e1;
  2512. r.i16[7] = e0;
  2513. #endif
  2514. return r;
  2515. }
  2516. SIMDE__FUNCTION_ATTRIBUTES
  2517. simde__m128i simde_mm_setr_epi32(int32_t e3, int32_t e2, int32_t e1, int32_t e0)
  2518. {
  2519. simde__m128i r;
  2520. #if defined(SIMDE_SSE2_NATIVE)
  2521. r.n = _mm_setr_epi32(e3, e2, e1, e0);
  2522. #elif defined(SIMDE_SSE2_NEON)
  2523. int32_t t[] = {e3, e2, e1, e0};
  2524. r.neon_i32 = vld1q_s32(t);
  2525. #else
  2526. r.i32[0] = e3;
  2527. r.i32[1] = e2;
  2528. r.i32[2] = e1;
  2529. r.i32[3] = e0;
  2530. #endif
  2531. return r;
  2532. }
  2533. SIMDE__FUNCTION_ATTRIBUTES
  2534. simde__m128i simde_mm_setr_epi64(simde__m64 e1, simde__m64 e0)
  2535. {
  2536. simde__m128i r;
  2537. #if defined(SIMDE_SSE2_NATIVE)
  2538. r.n = _mm_setr_epi64(e1.n, e0.n);
  2539. #elif defined(SIMDE_SSE2_NEON)
  2540. r.neon_i64 = vcombine_s64(e1.neon_i64, e0.neon_i64);
  2541. #else
  2542. r.i64[0] = e1.i64[0];
  2543. r.i64[1] = e0.i64[0];
  2544. #endif
  2545. return r;
  2546. }
  2547. SIMDE__FUNCTION_ATTRIBUTES
  2548. simde__m128d simde_mm_setr_pd(simde_float64 e1, simde_float64 e0)
  2549. {
  2550. simde__m128d r;
  2551. #if defined(SIMDE_SSE2_NATIVE)
  2552. r.n = _mm_setr_pd(e1, e0);
  2553. #else
  2554. r.f64[0] = e1;
  2555. r.f64[1] = e0;
  2556. #endif
  2557. return r;
  2558. }
  2559. SIMDE__FUNCTION_ATTRIBUTES
  2560. simde__m128d simde_mm_setzero_pd(void)
  2561. {
  2562. simde__m128d r;
  2563. #if defined(SIMDE_SSE2_NATIVE)
  2564. r.n = _mm_setzero_pd();
  2565. #else
  2566. r.u64[0] = 0;
  2567. r.u64[1] = 0;
  2568. #endif
  2569. return r;
  2570. }
  2571. SIMDE__FUNCTION_ATTRIBUTES
  2572. simde__m128i simde_mm_setzero_si128(void)
  2573. {
  2574. simde__m128i r;
  2575. #if defined(SIMDE_SSE2_NATIVE)
  2576. r.n = _mm_setzero_si128();
  2577. #elif defined(SIMDE_SSE2_NEON)
  2578. r.neon_i32 = vdupq_n_s32(0);
  2579. #else
  2580. r.u64[0] = 0;
  2581. r.u64[1] = 0;
  2582. #endif
  2583. return r;
  2584. }
  2585. SIMDE__FUNCTION_ATTRIBUTES
  2586. simde__m128i simde_mm_shuffle_epi32(simde__m128i a, const int imm8)
  2587. {
  2588. simde__m128i r;
  2589. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2590. r.i32[i] = a.i32[(imm8 >> (i * 2)) & 3];
  2591. }
  2592. return r;
  2593. }
  2594. #if defined(SIMDE_SSE2_NATIVE)
  2595. #define simde_mm_shuffle_epi32(a, imm8) \
  2596. SIMDE__M128I_C(_mm_shuffle_epi32((a).n, (imm8)))
  2597. #elif defined(SIMDE__SHUFFLE_VECTOR)
  2598. #define simde_mm_shuffle_epi32(a, imm8) \
  2599. ({ \
  2600. const simde__m128i simde__tmp_a_ = a; \
  2601. (simde__m128i){.i32 = SIMDE__SHUFFLE_VECTOR( \
  2602. 32, 16, (simde__tmp_a_).i32, \
  2603. (simde__tmp_a_).i32, ((imm8)) & 3, \
  2604. ((imm8) >> 2) & 3, ((imm8) >> 4) & 3, \
  2605. ((imm8) >> 6) & 3)}; \
  2606. })
  2607. #endif
  2608. SIMDE__FUNCTION_ATTRIBUTES
  2609. simde__m128d simde_mm_shuffle_pd(simde__m128d a, simde__m128d b, const int imm8)
  2610. {
  2611. simde__m128d r;
  2612. r.f64[0] = ((imm8 & 1) == 0) ? a.f64[0] : a.f64[1];
  2613. r.f64[1] = ((imm8 & 2) == 0) ? b.f64[0] : b.f64[1];
  2614. return r;
  2615. }
  2616. #if defined(SIMDE_SSE2_NATIVE) && !defined(__PGI)
  2617. #define simde_mm_shuffle_pd(a, b, imm8) \
  2618. SIMDE__M128D_C(_mm_shuffle_pd((a).n, (b).n, (imm8)))
  2619. #elif defined(SIMDE__SHUFFLE_VECTOR)
  2620. #define simde_mm_shuffle_pd(a, b, imm8) \
  2621. ({ \
  2622. (simde__m128d){.f64 = SIMDE__SHUFFLE_VECTOR( \
  2623. 64, 16, (a).f64, (b).f64, \
  2624. (((imm8)) & 1), \
  2625. (((imm8) >> 1) & 1) + 2)}; \
  2626. })
  2627. #endif
  2628. SIMDE__FUNCTION_ATTRIBUTES
  2629. simde__m128i simde_mm_shufflehi_epi16(simde__m128i a, const int imm8)
  2630. {
  2631. simde__m128i r;
  2632. r.i64[0] = a.i64[0];
  2633. for (size_t i = 4; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2634. r.i16[i] = a.i16[((imm8 >> ((i - 4) * 2)) & 3) + 4];
  2635. }
  2636. return r;
  2637. }
  2638. #if defined(SIMDE_SSE2_NATIVE)
  2639. #define simde_mm_shufflehi_epi16(a, imm8) \
  2640. SIMDE__M128I_C(_mm_shufflehi_epi16((a).n, (imm8)))
  2641. #elif defined(SIMDE__SHUFFLE_VECTOR)
  2642. #define simde_mm_shufflehi_epi16(a, imm8) \
  2643. ({ \
  2644. const simde__m128i simde__tmp_a_ = a; \
  2645. (simde__m128i){.i16 = SIMDE__SHUFFLE_VECTOR( \
  2646. 16, 16, (simde__tmp_a_).i16, \
  2647. (simde__tmp_a_).i16, 0, 1, 2, 3, \
  2648. (((imm8)) & 3) + 4, \
  2649. (((imm8) >> 2) & 3) + 4, \
  2650. (((imm8) >> 4) & 3) + 4, \
  2651. (((imm8) >> 6) & 3) + 4)}; \
  2652. })
  2653. #endif
  2654. SIMDE__FUNCTION_ATTRIBUTES
  2655. simde__m128i simde_mm_shufflelo_epi16(simde__m128i a, const int imm8)
  2656. {
  2657. simde__m128i r;
  2658. for (size_t i = 0; i < ((sizeof(r.i16) / sizeof(r.i16[0])) / 2); i++) {
  2659. r.i16[i] = a.i16[((imm8 >> (i * 2)) & 3)];
  2660. }
  2661. r.i64[1] = a.i64[1];
  2662. return r;
  2663. }
  2664. #if defined(SIMDE_SSE2_NATIVE)
  2665. #define simde_mm_shufflelo_epi16(a, imm8) \
  2666. SIMDE__M128I_C(_mm_shufflelo_epi16((a).n, (imm8)))
  2667. #elif defined(SIMDE__SHUFFLE_VECTOR)
  2668. #define simde_mm_shufflelo_epi16(a, imm8) \
  2669. ({ \
  2670. const simde__m128i simde__tmp_a_ = a; \
  2671. (simde__m128i){.i16 = SIMDE__SHUFFLE_VECTOR( \
  2672. 16, 16, (simde__tmp_a_).i16, \
  2673. (simde__tmp_a_).i16, (((imm8)) & 3), \
  2674. (((imm8) >> 2) & 3), \
  2675. (((imm8) >> 4) & 3), \
  2676. (((imm8) >> 6) & 3), 4, 5, 6, 7)}; \
  2677. })
  2678. #endif
  2679. SIMDE__FUNCTION_ATTRIBUTES
  2680. simde__m128i simde_mm_sll_epi16(simde__m128i a, simde__m128i count)
  2681. {
  2682. #if defined(SIMDE_SSE2_NATIVE)
  2683. return SIMDE__M128I_C(_mm_sll_epi16(a.n, count.n));
  2684. #else
  2685. simde__m128i r;
  2686. if (count.u64[0] > 15)
  2687. return simde_mm_setzero_si128();
  2688. const int s = (int)(count.u64[0]);
  2689. SIMDE__VECTORIZE
  2690. for (size_t i = 0; i < (sizeof(r.u16) / sizeof(r.u16[0])); i++) {
  2691. r.u16[i] = a.u16[i] << s;
  2692. }
  2693. return r;
  2694. #endif
  2695. }
  2696. SIMDE__FUNCTION_ATTRIBUTES
  2697. simde__m128i simde_mm_sll_epi32(simde__m128i a, simde__m128i count)
  2698. {
  2699. #if defined(SIMDE_SSE2_NATIVE)
  2700. return SIMDE__M128I_C(_mm_sll_epi32(a.n, count.n));
  2701. #else
  2702. simde__m128i r;
  2703. if (count.u64[0] > 31)
  2704. return simde_mm_setzero_si128();
  2705. const int s = (int)(count.u64[0]);
  2706. SIMDE__VECTORIZE
  2707. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2708. r.i32[i] = a.i32[i] << s;
  2709. }
  2710. return r;
  2711. #endif
  2712. }
  2713. SIMDE__FUNCTION_ATTRIBUTES
  2714. simde__m128i simde_mm_sll_epi64(simde__m128i a, simde__m128i count)
  2715. {
  2716. #if defined(SIMDE_SSE2_NATIVE)
  2717. return SIMDE__M128I_C(_mm_sll_epi64(a.n, count.n));
  2718. #else
  2719. simde__m128i r;
  2720. if (count.u64[0] > 63)
  2721. return simde_mm_setzero_si128();
  2722. const int s = (int)(count.u64[0]);
  2723. SIMDE__VECTORIZE
  2724. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2725. r.i64[i] = a.i64[i] << s;
  2726. }
  2727. return r;
  2728. #endif
  2729. }
  2730. SIMDE__FUNCTION_ATTRIBUTES
  2731. simde__m128d simde_mm_sqrt_pd(simde__m128d a)
  2732. {
  2733. #if defined(SIMDE_SSE2_NATIVE)
  2734. return SIMDE__M128D_C(_mm_sqrt_pd(a.n));
  2735. #else
  2736. simde__m128d r;
  2737. SIMDE__VECTORIZE
  2738. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  2739. r.f64[i] = sqrt(a.f64[i]);
  2740. }
  2741. return r;
  2742. #endif
  2743. }
  2744. SIMDE__FUNCTION_ATTRIBUTES
  2745. simde__m128d simde_mm_sqrt_sd(simde__m128d a, simde__m128d b)
  2746. {
  2747. #if defined(SIMDE_SSE2_NATIVE)
  2748. return SIMDE__M128D_C(_mm_sqrt_sd(a.n, b.n));
  2749. #else
  2750. simde__m128d r;
  2751. r.f64[0] = sqrt(b.f64[0]);
  2752. r.f64[1] = a.f64[1];
  2753. return r;
  2754. #endif
  2755. }
  2756. SIMDE__FUNCTION_ATTRIBUTES
  2757. simde__m128i simde_mm_srl_epi16(simde__m128i a, simde__m128i count)
  2758. {
  2759. #if defined(SIMDE_SSE2_NATIVE)
  2760. return SIMDE__M128I_C(_mm_srl_epi16(a.n, count.n));
  2761. #else
  2762. simde__m128i r;
  2763. if (count.u64[0] > 15)
  2764. return simde_mm_setzero_si128();
  2765. const int s = (int)(count.u64[0]);
  2766. SIMDE__VECTORIZE
  2767. for (size_t i = 0; i < (sizeof(r.u16) / sizeof(r.u16[0])); i++) {
  2768. r.u16[i] = a.u16[i] >> s;
  2769. }
  2770. return r;
  2771. #endif
  2772. }
  2773. SIMDE__FUNCTION_ATTRIBUTES
  2774. simde__m128i simde_mm_srl_epi32(simde__m128i a, simde__m128i count)
  2775. {
  2776. #if defined(SIMDE_SSE2_NATIVE)
  2777. return SIMDE__M128I_C(_mm_srl_epi32(a.n, count.n));
  2778. #else
  2779. simde__m128i r;
  2780. if (count.u64[0] > 31)
  2781. return simde_mm_setzero_si128();
  2782. const int s = (int)(count.u64[0]);
  2783. SIMDE__VECTORIZE
  2784. for (size_t i = 0; i < (sizeof(r.u32) / sizeof(r.u32[0])); i++) {
  2785. r.u32[i] = a.u32[i] >> s;
  2786. }
  2787. return r;
  2788. #endif
  2789. }
  2790. SIMDE__FUNCTION_ATTRIBUTES
  2791. simde__m128i simde_mm_srl_epi64(simde__m128i a, simde__m128i count)
  2792. {
  2793. #if defined(SIMDE_SSE2_NATIVE)
  2794. return SIMDE__M128I_C(_mm_srl_epi64(a.n, count.n));
  2795. #else
  2796. simde__m128i r;
  2797. if (count.u64[0] > 31)
  2798. return simde_mm_setzero_si128();
  2799. const int s = (int)(count.u64[0]);
  2800. SIMDE__VECTORIZE
  2801. for (size_t i = 0; i < (sizeof(r.u64) / sizeof(r.u64[0])); i++) {
  2802. r.u64[i] = a.u64[i] >> s;
  2803. }
  2804. return r;
  2805. #endif
  2806. }
  2807. SIMDE__FUNCTION_ATTRIBUTES
  2808. simde__m128i simde_mm_srai_epi16(simde__m128i a, int imm8)
  2809. {
  2810. simde__m128i r;
  2811. const uint16_t m =
  2812. (uint16_t)((~0U) << ((sizeof(int16_t) * CHAR_BIT) - imm8));
  2813. SIMDE__VECTORIZE
  2814. for (size_t i = 0; i < (sizeof(r) / sizeof(r.u16[0])); i++) {
  2815. const uint16_t is_neg = ((uint16_t)(
  2816. ((a.u16[i]) >> ((sizeof(int16_t) * CHAR_BIT) - 1))));
  2817. r.u16[i] = (a.u16[i] >> imm8) | (m * is_neg);
  2818. }
  2819. return r;
  2820. }
  2821. #if defined(SIMDE_SSE2_NATIVE)
  2822. #define simde_mm_srai_epi16(a, imm8) \
  2823. SIMDE__M128I_C(_mm_srai_epi16((a).n, (imm8)));
  2824. #endif
  2825. SIMDE__FUNCTION_ATTRIBUTES
  2826. simde__m128i simde_mm_srai_epi32(simde__m128i a, int imm8)
  2827. {
  2828. simde__m128i r;
  2829. const uint32_t m =
  2830. (uint32_t)((~0U) << ((sizeof(int) * CHAR_BIT) - imm8));
  2831. SIMDE__VECTORIZE
  2832. for (size_t i = 0; i < (sizeof(r) / sizeof(r.u32[0])); i++) {
  2833. uint32_t is_neg = ((uint32_t)(
  2834. ((a.u32[i]) >> ((sizeof(int32_t) * CHAR_BIT) - 1))));
  2835. r.u32[i] = (a.u32[i] >> imm8) | (m * is_neg);
  2836. }
  2837. return r;
  2838. }
  2839. #if defined(SIMDE_SSE2_NATIVE)
  2840. #define simde_mm_srai_epi32(a, imm8) \
  2841. SIMDE__M128I_C(_mm_srai_epi32((a).n, (imm8)))
  2842. #elif defined(SIMDE_SSE2_NEON)
  2843. #define simde_mm_srai_epi32(a, imm8) \
  2844. SIMDE__M128I_NEON_C( \
  2845. i32, \
  2846. ((imm8) <= 0) \
  2847. ? (a.neon_i32) \
  2848. : (((imm8) > 31) \
  2849. ? (vshrq_n_s32(vshrq_n_s32(a.neon_i32, 16), \
  2850. 16)) \
  2851. : (vshrq_n_s32(a.neon_i32, (imm8)))))
  2852. #endif
  2853. SIMDE__FUNCTION_ATTRIBUTES
  2854. simde__m128i simde_mm_sra_epi16(simde__m128i a, simde__m128i count)
  2855. {
  2856. #if defined(SIMDE_SSE2_NATIVE)
  2857. return SIMDE__M128I_C(_mm_sra_epi16(a.n, count.n));
  2858. #else
  2859. simde__m128i r;
  2860. int cnt = (int)count.i64[0];
  2861. if (cnt > 15 || cnt < 0) {
  2862. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0]));
  2863. i++) {
  2864. r.u16[i] = (a.i16[i] < 0) ? 0xffff : 0x0000;
  2865. }
  2866. } else {
  2867. const uint16_t m = (uint16_t)(
  2868. (~0U) << ((sizeof(int16_t) * CHAR_BIT) - cnt));
  2869. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0]));
  2870. i++) {
  2871. const uint16_t is_neg = a.i16[i] < 0;
  2872. r.u16[i] = (a.u16[i] >> cnt) | (m * is_neg);
  2873. }
  2874. }
  2875. return r;
  2876. #endif
  2877. }
  2878. SIMDE__FUNCTION_ATTRIBUTES
  2879. simde__m128i simde_mm_sra_epi32(simde__m128i a, simde__m128i count)
  2880. {
  2881. #if defined(SIMDE_SSE2_NATIVE) && !defined(SIMDE_BUG_GCC_BAD_MM_SRA_EPI32)
  2882. return SIMDE__M128I_C(_mm_sra_epi32(a.n, count.n));
  2883. #else
  2884. simde__m128i r;
  2885. const uint64_t cnt = count.u64[0];
  2886. if (cnt > 31) {
  2887. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0]));
  2888. i++) {
  2889. r.u32[i] = (a.i32[i] < 0) ? UINT32_MAX : 0;
  2890. }
  2891. } else if (cnt == 0) {
  2892. memcpy(&r, &a, sizeof(r));
  2893. } else {
  2894. const uint32_t m = (uint32_t)(
  2895. (~0U) << ((sizeof(int32_t) * CHAR_BIT) - cnt));
  2896. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0]));
  2897. i++) {
  2898. const uint32_t is_neg = a.i32[i] < 0;
  2899. r.u32[i] = (a.u32[i] >> cnt) | (m * is_neg);
  2900. }
  2901. }
  2902. return r;
  2903. #endif
  2904. }
  2905. SIMDE__FUNCTION_ATTRIBUTES
  2906. simde__m128i simde_mm_slli_epi16(simde__m128i a, const int imm8)
  2907. {
  2908. simde__m128i r;
  2909. const int s = (imm8 > ((int)sizeof(r.i16[0]) * CHAR_BIT) - 1) ? 0
  2910. : imm8;
  2911. SIMDE__VECTORIZE
  2912. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2913. r.i16[i] = a.i16[i] << s;
  2914. }
  2915. return r;
  2916. }
  2917. #if defined(SIMDE_SSE2_NATIVE)
  2918. #define simde_mm_slli_epi16(a, imm8) SIMDE__M128I_C(_mm_slli_epi16(a.n, imm8));
  2919. #elif defined(SIMDE_SSE2_NEON)
  2920. #define simde_mm_slli_epi16(a, imm8) \
  2921. SIMDE__M128I_NEON_C( \
  2922. i16, ((imm8) <= 0) \
  2923. ? ((a).neon_i16) \
  2924. : (((imm8) > 31) ? (vdupq_n_s16(0)) \
  2925. : (vshlq_n_s16((a).neon_i16, \
  2926. (imm8)))))
  2927. #endif
  2928. SIMDE__FUNCTION_ATTRIBUTES
  2929. simde__m128i simde_mm_slli_epi32(simde__m128i a, const int imm8)
  2930. {
  2931. simde__m128i r;
  2932. const int s = (imm8 > ((int)sizeof(r.i32[0]) * CHAR_BIT) - 1) ? 0
  2933. : imm8;
  2934. SIMDE__VECTORIZE
  2935. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2936. r.i32[i] = a.i32[i] << s;
  2937. }
  2938. return r;
  2939. }
  2940. #if defined(SIMDE_SSE2_NATIVE)
  2941. #define simde_mm_slli_epi32(a, imm8) SIMDE__M128I_C(_mm_slli_epi32(a.n, imm8));
  2942. #elif defined(SIMDE_SSE2_NEON)
  2943. #define simde_mm_slli_epi32(a, imm8) \
  2944. SIMDE__M128I_NEON_C( \
  2945. i32, ((imm8) <= 0) \
  2946. ? ((a).neon_i32) \
  2947. : (((imm8) > 31) ? (vdupq_n_s32(0)) \
  2948. : (vshlq_n_s32((a).neon_i32, \
  2949. (imm8)))))
  2950. #endif
  2951. SIMDE__FUNCTION_ATTRIBUTES
  2952. simde__m128i simde_mm_slli_epi64(simde__m128i a, const int imm8)
  2953. {
  2954. simde__m128i r;
  2955. const int s = (imm8 > ((int)sizeof(r.i64[0]) * CHAR_BIT) - 1) ? 0
  2956. : imm8;
  2957. SIMDE__VECTORIZE
  2958. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  2959. r.i64[i] = a.i64[i] << s;
  2960. }
  2961. return r;
  2962. }
  2963. #if defined(SIMDE_SSE2_NATIVE)
  2964. #define simde_mm_slli_epi64(a, imm8) SIMDE__M128I_C(_mm_slli_epi64(a.n, imm8));
  2965. #endif
  2966. SIMDE__FUNCTION_ATTRIBUTES
  2967. simde__m128i simde_mm_srli_epi16(simde__m128i a, const int imm8)
  2968. {
  2969. simde__m128i r;
  2970. const int s = (imm8 > ((int)sizeof(r.i16[0]) * CHAR_BIT) - 1) ? 0
  2971. : imm8;
  2972. SIMDE__VECTORIZE
  2973. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  2974. r.u16[i] = a.u16[i] >> s;
  2975. }
  2976. return r;
  2977. }
  2978. #if defined(SIMDE_SSE2_NATIVE)
  2979. #define simde_mm_srli_epi16(a, imm8) SIMDE__M128I_C(_mm_srli_epi16(a.n, imm8));
  2980. #elif defined(SIMDE_SSE2_NEON)
  2981. #define simde_mm_srli_epi16(a, imm8) \
  2982. SIMDE__M128I_NEON_C( \
  2983. u16, ((imm8) <= 0) \
  2984. ? ((a).neon_u16) \
  2985. : (((imm8) > 31) ? (vdupq_n_u16(0)) \
  2986. : (vshrq_n_u16((a).neon_u16, \
  2987. (imm8)))))
  2988. #endif
  2989. SIMDE__FUNCTION_ATTRIBUTES
  2990. simde__m128i simde_mm_srli_epi32(simde__m128i a, const int imm8)
  2991. {
  2992. simde__m128i r;
  2993. const int s = (imm8 > ((int)sizeof(r.i32[0]) * CHAR_BIT) - 1) ? 0
  2994. : imm8;
  2995. SIMDE__VECTORIZE
  2996. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  2997. r.u32[i] = a.u32[i] >> s;
  2998. }
  2999. return r;
  3000. }
  3001. #if defined(SIMDE_SSE2_NATIVE)
  3002. #define simde_mm_srli_epi32(a, imm8) SIMDE__M128I_C(_mm_srli_epi32(a.n, imm8))
  3003. #elif defined(SIMDE_SSE2_NEON)
  3004. #define simde_mm_srli_epi32(a, imm8) \
  3005. SIMDE__M128I_NEON_C( \
  3006. u32, ((imm8) <= 0) \
  3007. ? ((a).neon_u32) \
  3008. : (((imm8) > 31) ? (vdupq_n_u32(0)) \
  3009. : (vshrq_n_u32((a).neon_u32, \
  3010. (imm8)))))
  3011. #endif
  3012. SIMDE__FUNCTION_ATTRIBUTES
  3013. simde__m128i simde_mm_srli_epi64(simde__m128i a, const int imm8)
  3014. {
  3015. simde__m128i r;
  3016. const unsigned char s = imm8 & 255;
  3017. SIMDE__VECTORIZE
  3018. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  3019. if (s > 63) {
  3020. r.u64[i] = 0;
  3021. } else {
  3022. r.u64[i] = a.u64[i] >> s;
  3023. }
  3024. }
  3025. return r;
  3026. }
  3027. #if defined(SIMDE_SSE2_NATIVE)
  3028. #define simde_mm_srli_epi64(a, imm8) SIMDE__M128I_C(_mm_srli_epi64(a.n, imm8))
  3029. #elif defined(SIMDE_SSE2_NEON)
  3030. #define simde_mm_srli_epi64(a, imm8) \
  3031. SIMDE__M128I_NEON_C( \
  3032. u64, \
  3033. (((imm8)&255) < 0 || ((imm8)&255) > 63) \
  3034. ? (vdupq_n_u64(0)) \
  3035. : ((((imm8)&255) == 0) \
  3036. ? (a.neon_u64) \
  3037. : (vshrq_n_u64((a).neon_u64, (imm8)&255))))
  3038. #endif
  3039. SIMDE__FUNCTION_ATTRIBUTES
  3040. void simde_mm_store_pd(simde_float64 mem_addr[HEDLEY_ARRAY_PARAM(2)],
  3041. simde__m128d a)
  3042. {
  3043. simde_assert_aligned(16, mem_addr);
  3044. #if defined(SIMDE_SSE2_NATIVE)
  3045. _mm_store_pd(mem_addr, a.n);
  3046. #else
  3047. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3048. memcpy(mem_addr, &a, sizeof(a));
  3049. #endif
  3050. }
  3051. SIMDE__FUNCTION_ATTRIBUTES
  3052. void simde_mm_store1_pd(simde_float64 mem_addr[HEDLEY_ARRAY_PARAM(2)],
  3053. simde__m128d a)
  3054. {
  3055. simde_assert_aligned(16, mem_addr);
  3056. #if defined(SIMDE_SSE2_NATIVE)
  3057. _mm_store1_pd(mem_addr, a.n);
  3058. #else
  3059. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3060. mem_addr[0] = a.f64[0];
  3061. mem_addr[1] = a.f64[0];
  3062. #endif
  3063. }
  3064. #define simde_mm_store_pd1(mem_addr, a) simde_mm_store1_pd(mem_addr, a)
  3065. SIMDE__FUNCTION_ATTRIBUTES
  3066. void simde_mm_store_sd(simde_float64 *mem_addr, simde__m128d a)
  3067. {
  3068. #if defined(SIMDE_SSE2_NATIVE)
  3069. _mm_store_sd(mem_addr, a.n);
  3070. #else
  3071. memcpy(mem_addr, &a, sizeof(a.f64[0]));
  3072. #endif
  3073. }
  3074. SIMDE__FUNCTION_ATTRIBUTES
  3075. void simde_mm_store_si128(simde__m128i *mem_addr, simde__m128i a)
  3076. {
  3077. #if defined(SIMDE_SSE2_NATIVE)
  3078. _mm_store_si128(&mem_addr->n, a.n);
  3079. #elif defined(SIMDE_SSE2_NEON)
  3080. vst1q_s32((int32_t *)mem_addr, a.neon_i32);
  3081. #else
  3082. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3083. memcpy(mem_addr, &a, sizeof(a));
  3084. #endif
  3085. }
  3086. SIMDE__FUNCTION_ATTRIBUTES
  3087. void simde_mm_storeh_pd(simde_float64 *mem_addr, simde__m128d a)
  3088. {
  3089. #if defined(SIMDE_SSE2_NATIVE)
  3090. _mm_storeh_pd(mem_addr, a.n);
  3091. #else
  3092. *mem_addr = a.f64[1];
  3093. #endif
  3094. }
  3095. SIMDE__FUNCTION_ATTRIBUTES
  3096. void simde_mm_storel_epi64(simde__m128i *mem_addr, simde__m128i a)
  3097. {
  3098. #if defined(SIMDE_SSE2_NATIVE)
  3099. _mm_storel_epi64(&(mem_addr->n), a.n);
  3100. #elif defined(SIMDE_SSE2_NEON)
  3101. mem_addr->i64[0] = vgetq_lane_s64(a.neon_i64, 0);
  3102. #else
  3103. mem_addr->i64[0] = a.i64[0];
  3104. #endif
  3105. }
  3106. SIMDE__FUNCTION_ATTRIBUTES
  3107. void simde_mm_storel_pd(simde_float64 *mem_addr, simde__m128d a)
  3108. {
  3109. #if defined(SIMDE_SSE2_NATIVE)
  3110. _mm_storel_pd(mem_addr, a.n);
  3111. #else
  3112. *mem_addr = a.f64[0];
  3113. #endif
  3114. }
  3115. SIMDE__FUNCTION_ATTRIBUTES
  3116. void simde_mm_storer_pd(simde_float64 mem_addr[2], simde__m128d a)
  3117. {
  3118. simde_assert_aligned(16, mem_addr);
  3119. #if defined(SIMDE_SSE2_NATIVE)
  3120. _mm_storer_pd(mem_addr, a.n);
  3121. #else
  3122. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3123. mem_addr[0] = a.f64[1];
  3124. mem_addr[1] = a.f64[0];
  3125. #endif
  3126. }
  3127. SIMDE__FUNCTION_ATTRIBUTES
  3128. void simde_mm_storeu_pd(simde_float64 *mem_addr, simde__m128d a)
  3129. {
  3130. #if defined(SIMDE_SSE2_NATIVE)
  3131. _mm_storeu_pd(mem_addr, a.n);
  3132. #else
  3133. memcpy(mem_addr, &a, sizeof(a));
  3134. #endif
  3135. }
  3136. SIMDE__FUNCTION_ATTRIBUTES
  3137. void simde_mm_storeu_si128(simde__m128i *mem_addr, simde__m128i a)
  3138. {
  3139. #if defined(SIMDE_SSE2_NATIVE)
  3140. _mm_storeu_si128(&mem_addr->n, a.n);
  3141. #elif defined(SIMDE_SSE2_NEON)
  3142. int32_t v[4];
  3143. vst1q_s32(v, a.neon_i32);
  3144. memcpy(mem_addr, v, sizeof(v));
  3145. #else
  3146. memcpy(mem_addr, &a, sizeof(a));
  3147. #endif
  3148. }
  3149. SIMDE__FUNCTION_ATTRIBUTES
  3150. void simde_mm_stream_pd(simde_float64 mem_addr[HEDLEY_ARRAY_PARAM(2)],
  3151. simde__m128d a)
  3152. {
  3153. #if defined(SIMDE_SSE2_NATIVE)
  3154. _mm_stream_pd(mem_addr, a.n);
  3155. #else
  3156. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3157. memcpy(mem_addr, &a, sizeof(a));
  3158. #endif
  3159. }
  3160. SIMDE__FUNCTION_ATTRIBUTES
  3161. void simde_mm_stream_si128(simde__m128i *mem_addr, simde__m128i a)
  3162. {
  3163. #if defined(SIMDE_SSE2_NATIVE)
  3164. _mm_stream_si128(&mem_addr->n, a.n);
  3165. #else
  3166. SIMDE__ASSUME_ALIGNED(mem_addr, 16);
  3167. memcpy(mem_addr, &a, sizeof(a));
  3168. #endif
  3169. }
  3170. SIMDE__FUNCTION_ATTRIBUTES
  3171. void simde_mm_stream_si32(int32_t *mem_addr, int32_t a)
  3172. {
  3173. #if defined(SIMDE_SSE2_NATIVE)
  3174. _mm_stream_si32(mem_addr, a);
  3175. #else
  3176. *mem_addr = a;
  3177. #endif
  3178. }
  3179. SIMDE__FUNCTION_ATTRIBUTES
  3180. void simde_mm_stream_si64(int64_t *mem_addr, int64_t a)
  3181. {
  3182. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE_ARCH_AMD64)
  3183. #if defined(SIMDE__REALLY_GCC) && !HEDLEY_GCC_VERSION_CHECK(4, 8, 0)
  3184. *mem_addr = a;
  3185. #elif defined(__GNUC__)
  3186. _mm_stream_si64((long long *)mem_addr, a);
  3187. #else
  3188. _mm_stream_si64(mem_addr, a);
  3189. #endif
  3190. #else
  3191. *mem_addr = a;
  3192. #endif
  3193. }
  3194. SIMDE__FUNCTION_ATTRIBUTES
  3195. simde__m128i simde_mm_sub_epi8(simde__m128i a, simde__m128i b)
  3196. {
  3197. #if defined(SIMDE_SSE2_NATIVE)
  3198. return SIMDE__M128I_C(_mm_sub_epi8(a.n, b.n));
  3199. #elif defined(SIMDE_SSE2_NEON)
  3200. return SIMDE__M128I_NEON_C(i8, vsubq_s8(a.neon_i8, b.neon_i8));
  3201. #else
  3202. simde__m128i r;
  3203. SIMDE__VECTORIZE
  3204. for (size_t i = 0; i < (sizeof(r.i8) / sizeof(r.i8[0])); i++) {
  3205. r.i8[i] = a.i8[i] - b.i8[i];
  3206. }
  3207. return r;
  3208. #endif
  3209. }
  3210. SIMDE__FUNCTION_ATTRIBUTES
  3211. simde__m128i simde_mm_sub_epi16(simde__m128i a, simde__m128i b)
  3212. {
  3213. #if defined(SIMDE_SSE2_NATIVE)
  3214. return SIMDE__M128I_C(_mm_sub_epi16(a.n, b.n));
  3215. #elif defined(SIMDE_SSE2_NEON)
  3216. return SIMDE__M128I_NEON_C(i16, vsubq_s16(a.neon_i16, b.neon_i16));
  3217. #else
  3218. simde__m128i r;
  3219. SIMDE__VECTORIZE
  3220. for (size_t i = 0; i < (sizeof(r.i16) / sizeof(r.i16[0])); i++) {
  3221. r.i16[i] = a.i16[i] - b.i16[i];
  3222. }
  3223. return r;
  3224. #endif
  3225. }
  3226. SIMDE__FUNCTION_ATTRIBUTES
  3227. simde__m128i simde_mm_sub_epi32(simde__m128i a, simde__m128i b)
  3228. {
  3229. #if defined(SIMDE_SSE2_NATIVE)
  3230. return SIMDE__M128I_C(_mm_sub_epi32(a.n, b.n));
  3231. #elif defined(SIMDE_SSE2_NEON)
  3232. return SIMDE__M128I_NEON_C(i32, vsubq_s32(a.neon_i32, b.neon_i32));
  3233. #else
  3234. simde__m128i r;
  3235. SIMDE__VECTORIZE
  3236. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  3237. r.i32[i] = a.i32[i] - b.i32[i];
  3238. }
  3239. return r;
  3240. #endif
  3241. }
  3242. SIMDE__FUNCTION_ATTRIBUTES
  3243. simde__m128i simde_mm_sub_epi64(simde__m128i a, simde__m128i b)
  3244. {
  3245. #if defined(SIMDE_SSE2_NATIVE)
  3246. return SIMDE__M128I_C(_mm_sub_epi64(a.n, b.n));
  3247. #elif defined(SIMDE_SSE2_NEON)
  3248. return SIMDE__M128I_NEON_C(i64, vsubq_s64(a.neon_i64, b.neon_i64));
  3249. #else
  3250. simde__m128i r;
  3251. SIMDE__VECTORIZE
  3252. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  3253. r.i64[i] = a.i64[i] - b.i64[i];
  3254. }
  3255. return r;
  3256. #endif
  3257. }
  3258. SIMDE__FUNCTION_ATTRIBUTES
  3259. simde__m128d simde_mm_sub_pd(simde__m128d a, simde__m128d b)
  3260. {
  3261. #if defined(SIMDE_SSE2_NATIVE)
  3262. return SIMDE__M128D_C(_mm_sub_pd(a.n, b.n));
  3263. #else
  3264. simde__m128d r;
  3265. SIMDE__VECTORIZE
  3266. for (size_t i = 0; i < (sizeof(r.f64) / sizeof(r.f64[0])); i++) {
  3267. r.f64[i] = a.f64[i] - b.f64[i];
  3268. }
  3269. return r;
  3270. #endif
  3271. }
  3272. SIMDE__FUNCTION_ATTRIBUTES
  3273. simde__m128d simde_mm_sub_sd(simde__m128d a, simde__m128d b)
  3274. {
  3275. #if defined(SIMDE_SSE2_NATIVE)
  3276. return SIMDE__M128D_C(_mm_sub_sd(a.n, b.n));
  3277. #else
  3278. simde__m128d r;
  3279. r.f64[0] = a.f64[0] - b.f64[0];
  3280. r.f64[1] = a.f64[1];
  3281. return r;
  3282. #endif
  3283. }
  3284. SIMDE__FUNCTION_ATTRIBUTES
  3285. simde__m64 simde_mm_sub_si64(simde__m64 a, simde__m64 b)
  3286. {
  3287. #if defined(SIMDE_SSE2_NATIVE)
  3288. return SIMDE__M64_C(_mm_sub_si64(a.n, b.n));
  3289. #else
  3290. simde__m64 r;
  3291. r.i64[0] = a.i64[0] - b.i64[0];
  3292. return r;
  3293. #endif
  3294. }
  3295. SIMDE__FUNCTION_ATTRIBUTES
  3296. simde__m128i simde_mm_subs_epi8(simde__m128i a, simde__m128i b)
  3297. {
  3298. #if defined(SIMDE_SSE2_NATIVE)
  3299. return SIMDE__M128I_C(_mm_subs_epi8(a.n, b.n));
  3300. #elif defined(SIMDE_SSE2_NEON)
  3301. return SIMDE__M128I_NEON_C(i8, vqsubq_s8(a.neon_i8, b.neon_i8));
  3302. #else
  3303. simde__m128i r;
  3304. SIMDE__VECTORIZE
  3305. for (size_t i = 0; i < (sizeof(r) / sizeof(r.i8[0])); i++) {
  3306. if (((b.i8[i]) > 0 && (a.i8[i]) < INT8_MIN + (b.i8[i]))) {
  3307. r.i8[i] = INT8_MIN;
  3308. } else if ((b.i8[i]) < 0 && (a.i8[i]) > INT8_MAX + (b.i8[i])) {
  3309. r.i8[i] = INT8_MAX;
  3310. } else {
  3311. r.i8[i] = (a.i8[i]) - (b.i8[i]);
  3312. }
  3313. }
  3314. return r;
  3315. #endif
  3316. }
  3317. SIMDE__FUNCTION_ATTRIBUTES
  3318. simde__m128i simde_mm_subs_epi16(simde__m128i a, simde__m128i b)
  3319. {
  3320. #if defined(SIMDE_SSE2_NATIVE)
  3321. return SIMDE__M128I_C(_mm_subs_epi16(a.n, b.n));
  3322. #elif defined(SIMDE_SSE2_NEON)
  3323. return SIMDE__M128I_NEON_C(i16, vqsubq_s16(a.neon_i16, b.neon_i16));
  3324. #else
  3325. simde__m128i r;
  3326. SIMDE__VECTORIZE
  3327. for (size_t i = 0; i < (sizeof(r) / sizeof(r.i16[0])); i++) {
  3328. if (((b.i16[i]) > 0 && (a.i16[i]) < INT16_MIN + (b.i16[i]))) {
  3329. r.i16[i] = INT16_MIN;
  3330. } else if ((b.i16[i]) < 0 &&
  3331. (a.i16[i]) > INT16_MAX + (b.i16[i])) {
  3332. r.i16[i] = INT16_MAX;
  3333. } else {
  3334. r.i16[i] = (a.i16[i]) - (b.i16[i]);
  3335. }
  3336. }
  3337. return r;
  3338. #endif
  3339. }
  3340. SIMDE__FUNCTION_ATTRIBUTES
  3341. simde__m128i simde_mm_subs_epu8(simde__m128i a, simde__m128i b)
  3342. {
  3343. #if defined(SIMDE_SSE2_NATIVE)
  3344. return SIMDE__M128I_C(_mm_subs_epu8(a.n, b.n));
  3345. #elif defined(SIMDE_SSE2_NEON)
  3346. return SIMDE__M128I_NEON_C(u8, vqsubq_u8(a.neon_u8, b.neon_u8));
  3347. #else
  3348. simde__m128i r;
  3349. SIMDE__VECTORIZE
  3350. for (size_t i = 0; i < (sizeof(r) / sizeof(r.i8[0])); i++) {
  3351. const int32_t x = a.u8[i] - b.u8[i];
  3352. if (x < 0) {
  3353. r.u8[i] = 0;
  3354. } else if (x > UINT8_MAX) {
  3355. r.u8[i] = UINT8_MAX;
  3356. } else {
  3357. r.u8[i] = (uint8_t)x;
  3358. }
  3359. }
  3360. return r;
  3361. #endif
  3362. }
  3363. SIMDE__FUNCTION_ATTRIBUTES
  3364. simde__m128i simde_mm_subs_epu16(simde__m128i a, simde__m128i b)
  3365. {
  3366. #if defined(SIMDE_SSE2_NATIVE)
  3367. return SIMDE__M128I_C(_mm_subs_epu16(a.n, b.n));
  3368. #elif defined(SIMDE_SSE2_NEON)
  3369. return SIMDE__M128I_NEON_C(u16, vqsubq_u16(a.neon_u16, b.neon_u16));
  3370. #else
  3371. simde__m128i r;
  3372. SIMDE__VECTORIZE
  3373. for (size_t i = 0; i < (sizeof(r) / sizeof(r.i16[0])); i++) {
  3374. const int32_t x = a.u16[i] - b.u16[i];
  3375. if (x < 0) {
  3376. r.u16[i] = 0;
  3377. } else if (x > UINT16_MAX) {
  3378. r.u16[i] = UINT16_MAX;
  3379. } else {
  3380. r.u16[i] = (uint16_t)x;
  3381. }
  3382. }
  3383. return r;
  3384. #endif
  3385. }
  3386. SIMDE__FUNCTION_ATTRIBUTES
  3387. int simde_mm_ucomieq_sd(simde__m128d a, simde__m128d b)
  3388. {
  3389. #if defined(SIMDE_SSE2_NATIVE)
  3390. return _mm_ucomieq_sd(a.n, b.n);
  3391. #else
  3392. fenv_t envp;
  3393. int x = feholdexcept(&envp);
  3394. int r = a.f64[0] == b.f64[0];
  3395. if (HEDLEY_LIKELY(x == 0))
  3396. fesetenv(&envp);
  3397. return r;
  3398. #endif
  3399. }
  3400. SIMDE__FUNCTION_ATTRIBUTES
  3401. int simde_mm_ucomige_sd(simde__m128d a, simde__m128d b)
  3402. {
  3403. #if defined(SIMDE_SSE2_NATIVE)
  3404. return _mm_ucomige_sd(a.n, b.n);
  3405. #else
  3406. fenv_t envp;
  3407. int x = feholdexcept(&envp);
  3408. int r = a.f64[0] >= b.f64[0];
  3409. if (HEDLEY_LIKELY(x == 0))
  3410. fesetenv(&envp);
  3411. return r;
  3412. #endif
  3413. }
  3414. SIMDE__FUNCTION_ATTRIBUTES
  3415. int simde_mm_ucomigt_sd(simde__m128d a, simde__m128d b)
  3416. {
  3417. #if defined(SIMDE_SSE2_NATIVE)
  3418. return _mm_ucomigt_sd(a.n, b.n);
  3419. #else
  3420. fenv_t envp;
  3421. int x = feholdexcept(&envp);
  3422. int r = a.f64[0] > b.f64[0];
  3423. if (HEDLEY_LIKELY(x == 0))
  3424. fesetenv(&envp);
  3425. return r;
  3426. #endif
  3427. }
  3428. SIMDE__FUNCTION_ATTRIBUTES
  3429. int simde_mm_ucomile_sd(simde__m128d a, simde__m128d b)
  3430. {
  3431. #if defined(SIMDE_SSE2_NATIVE)
  3432. return _mm_ucomile_sd(a.n, b.n);
  3433. #else
  3434. fenv_t envp;
  3435. int x = feholdexcept(&envp);
  3436. int r = a.f64[0] <= b.f64[0];
  3437. if (HEDLEY_LIKELY(x == 0))
  3438. fesetenv(&envp);
  3439. return r;
  3440. #endif
  3441. }
  3442. SIMDE__FUNCTION_ATTRIBUTES
  3443. int simde_mm_ucomilt_sd(simde__m128d a, simde__m128d b)
  3444. {
  3445. #if defined(SIMDE_SSE2_NATIVE)
  3446. return _mm_ucomilt_sd(a.n, b.n);
  3447. #else
  3448. fenv_t envp;
  3449. int x = feholdexcept(&envp);
  3450. int r = a.f64[0] < b.f64[0];
  3451. if (HEDLEY_LIKELY(x == 0))
  3452. fesetenv(&envp);
  3453. return r;
  3454. #endif
  3455. }
  3456. SIMDE__FUNCTION_ATTRIBUTES
  3457. int simde_mm_ucomineq_sd(simde__m128d a, simde__m128d b)
  3458. {
  3459. #if defined(SIMDE_SSE2_NATIVE)
  3460. return _mm_ucomineq_sd(a.n, b.n);
  3461. #else
  3462. fenv_t envp;
  3463. int x = feholdexcept(&envp);
  3464. int r = a.f64[0] != b.f64[0];
  3465. if (HEDLEY_LIKELY(x == 0))
  3466. fesetenv(&envp);
  3467. return r;
  3468. #endif
  3469. }
  3470. SIMDE__FUNCTION_ATTRIBUTES
  3471. simde__m128d simde_mm_undefined_pd(void)
  3472. {
  3473. simde__m128d r;
  3474. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE__HAVE_UNDEFINED128)
  3475. r.n = _mm_undefined_pd();
  3476. #else
  3477. r = simde_mm_setzero_pd();
  3478. #endif
  3479. return r;
  3480. }
  3481. SIMDE__FUNCTION_ATTRIBUTES
  3482. simde__m128i simde_mm_undefined_si128(void)
  3483. {
  3484. simde__m128i r;
  3485. #if defined(SIMDE_SSE2_NATIVE) && defined(SIMDE__HAVE_UNDEFINED128)
  3486. r.n = _mm_undefined_si128();
  3487. #else
  3488. r = simde_mm_setzero_si128();
  3489. #endif
  3490. return r;
  3491. }
  3492. SIMDE__FUNCTION_ATTRIBUTES
  3493. void simde_mm_lfence(void)
  3494. {
  3495. #if defined(SIMDE_SSE2_NATIVE)
  3496. _mm_lfence();
  3497. #else
  3498. simde_mm_sfence();
  3499. #endif
  3500. }
  3501. SIMDE__FUNCTION_ATTRIBUTES
  3502. void simde_mm_mfence(void)
  3503. {
  3504. #if defined(SIMDE_SSE2_NATIVE)
  3505. _mm_mfence();
  3506. #else
  3507. simde_mm_sfence();
  3508. #endif
  3509. }
  3510. SIMDE__FUNCTION_ATTRIBUTES
  3511. simde__m128i simde_mm_unpackhi_epi8(simde__m128i a, simde__m128i b)
  3512. {
  3513. #if defined(SIMDE_SSE2_NATIVE)
  3514. return SIMDE__M128I_C(_mm_unpackhi_epi8(a.n, b.n));
  3515. #elif defined(SIMDE_SSE2_NEON)
  3516. int8x8_t a1 = vreinterpret_s8_s16(vget_high_s16(a.neon_i16));
  3517. int8x8_t b1 = vreinterpret_s8_s16(vget_high_s16(b.neon_i16));
  3518. int8x8x2_t result = vzip_s8(a1, b1);
  3519. return SIMDE__M128I_NEON_C(i8,
  3520. vcombine_s8(result.val[0], result.val[1]));
  3521. #else
  3522. simde__m128i r;
  3523. SIMDE__VECTORIZE
  3524. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i8[0])) / 2); i++) {
  3525. r.i8[(i * 2)] = a.i8[i + ((sizeof(r) / sizeof(r.i8[0])) / 2)];
  3526. r.i8[(i * 2) + 1] =
  3527. b.i8[i + ((sizeof(r) / sizeof(r.i8[0])) / 2)];
  3528. }
  3529. return r;
  3530. #endif
  3531. }
  3532. SIMDE__FUNCTION_ATTRIBUTES
  3533. simde__m128i simde_mm_unpackhi_epi16(simde__m128i a, simde__m128i b)
  3534. {
  3535. #if defined(SIMDE_SSE2_NATIVE)
  3536. return SIMDE__M128I_C(_mm_unpackhi_epi16(a.n, b.n));
  3537. #elif defined(SIMDE_SSE2_NEON)
  3538. int16x4_t a1 = vget_high_s16(a.neon_i16);
  3539. int16x4_t b1 = vget_high_s16(b.neon_i16);
  3540. int16x4x2_t result = vzip_s16(a1, b1);
  3541. return SIMDE__M128I_NEON_C(i16,
  3542. vcombine_s16(result.val[0], result.val[1]));
  3543. #else
  3544. simde__m128i r;
  3545. SIMDE__VECTORIZE
  3546. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i16[0])) / 2); i++) {
  3547. r.i16[(i * 2)] =
  3548. a.i16[i + ((sizeof(r) / sizeof(r.i16[0])) / 2)];
  3549. r.i16[(i * 2) + 1] =
  3550. b.i16[i + ((sizeof(r) / sizeof(r.i16[0])) / 2)];
  3551. }
  3552. return r;
  3553. #endif
  3554. }
  3555. SIMDE__FUNCTION_ATTRIBUTES
  3556. simde__m128i simde_mm_unpackhi_epi32(simde__m128i a, simde__m128i b)
  3557. {
  3558. #if defined(SIMDE_SSE2_NATIVE)
  3559. return SIMDE__M128I_C(_mm_unpackhi_epi32(a.n, b.n));
  3560. #elif defined(SIMDE_SSE2_NEON)
  3561. int32x2_t a1 = vget_high_s32(a.neon_i32);
  3562. int32x2_t b1 = vget_high_s32(b.neon_i32);
  3563. int32x2x2_t result = vzip_s32(a1, b1);
  3564. return SIMDE__M128I_NEON_C(i32,
  3565. vcombine_s32(result.val[0], result.val[1]));
  3566. #else
  3567. simde__m128i r;
  3568. SIMDE__VECTORIZE
  3569. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i32[0])) / 2); i++) {
  3570. r.i32[(i * 2)] =
  3571. a.i32[i + ((sizeof(r) / sizeof(r.i32[0])) / 2)];
  3572. r.i32[(i * 2) + 1] =
  3573. b.i32[i + ((sizeof(r) / sizeof(r.i32[0])) / 2)];
  3574. }
  3575. return r;
  3576. #endif
  3577. }
  3578. SIMDE__FUNCTION_ATTRIBUTES
  3579. simde__m128i simde_mm_unpackhi_epi64(simde__m128i a, simde__m128i b)
  3580. {
  3581. #if defined(SIMDE_SSE2_NATIVE)
  3582. return SIMDE__M128I_C(_mm_unpackhi_epi64(a.n, b.n));
  3583. #else
  3584. simde__m128i r;
  3585. SIMDE__VECTORIZE
  3586. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i64[0])) / 2); i++) {
  3587. r.i64[(i * 2)] =
  3588. a.i64[i + ((sizeof(r) / sizeof(r.i64[0])) / 2)];
  3589. r.i64[(i * 2) + 1] =
  3590. b.i64[i + ((sizeof(r) / sizeof(r.i64[0])) / 2)];
  3591. }
  3592. return r;
  3593. #endif
  3594. }
  3595. SIMDE__FUNCTION_ATTRIBUTES
  3596. simde__m128d simde_mm_unpackhi_pd(simde__m128d a, simde__m128d b)
  3597. {
  3598. #if defined(SIMDE_SSE2_NATIVE)
  3599. return SIMDE__M128D_C(_mm_unpackhi_pd(a.n, b.n));
  3600. #else
  3601. simde__m128d r;
  3602. SIMDE__VECTORIZE
  3603. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.f64[0])) / 2); i++) {
  3604. r.f64[(i * 2)] =
  3605. a.f64[i + ((sizeof(r) / sizeof(r.f64[0])) / 2)];
  3606. r.f64[(i * 2) + 1] =
  3607. b.f64[i + ((sizeof(r) / sizeof(r.f64[0])) / 2)];
  3608. }
  3609. return r;
  3610. #endif
  3611. }
  3612. SIMDE__FUNCTION_ATTRIBUTES
  3613. simde__m128i simde_mm_unpacklo_epi8(simde__m128i a, simde__m128i b)
  3614. {
  3615. #if defined(SIMDE_SSE2_NATIVE)
  3616. return SIMDE__M128I_C(_mm_unpacklo_epi8(a.n, b.n));
  3617. #elif defined(SIMDE_SSE2_NEON)
  3618. int8x8_t a1 = vreinterpret_s8_s16(vget_low_s16(a.neon_i16));
  3619. int8x8_t b1 = vreinterpret_s8_s16(vget_low_s16(b.neon_i16));
  3620. int8x8x2_t result = vzip_s8(a1, b1);
  3621. return SIMDE__M128I_NEON_C(i8,
  3622. vcombine_s8(result.val[0], result.val[1]));
  3623. #else
  3624. simde__m128i r;
  3625. SIMDE__VECTORIZE
  3626. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i8[0])) / 2); i++) {
  3627. r.i8[(i * 2)] = a.i8[i];
  3628. r.i8[(i * 2) + 1] = b.i8[i];
  3629. }
  3630. return r;
  3631. #endif
  3632. }
  3633. SIMDE__FUNCTION_ATTRIBUTES
  3634. simde__m128i simde_mm_unpacklo_epi16(simde__m128i a, simde__m128i b)
  3635. {
  3636. #if defined(SIMDE_SSE2_NATIVE)
  3637. return SIMDE__M128I_C(_mm_unpacklo_epi16(a.n, b.n));
  3638. #elif defined(SIMDE_SSE2_NEON)
  3639. int16x4_t a1 = vget_low_s16(a.neon_i16);
  3640. int16x4_t b1 = vget_low_s16(b.neon_i16);
  3641. int16x4x2_t result = vzip_s16(a1, b1);
  3642. return SIMDE__M128I_NEON_C(i16,
  3643. vcombine_s16(result.val[0], result.val[1]));
  3644. #else
  3645. simde__m128i r;
  3646. SIMDE__VECTORIZE
  3647. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i16[0])) / 2); i++) {
  3648. r.i16[(i * 2)] = a.i16[i];
  3649. r.i16[(i * 2) + 1] = b.i16[i];
  3650. }
  3651. return r;
  3652. #endif
  3653. }
  3654. SIMDE__FUNCTION_ATTRIBUTES
  3655. simde__m128i simde_mm_unpacklo_epi32(simde__m128i a, simde__m128i b)
  3656. {
  3657. #if defined(SIMDE_SSE2_NATIVE)
  3658. return SIMDE__M128I_C(_mm_unpacklo_epi32(a.n, b.n));
  3659. #elif defined(SIMDE_SSE2_NEON)
  3660. int32x2_t a1 = vget_low_s32(a.neon_i32);
  3661. int32x2_t b1 = vget_low_s32(b.neon_i32);
  3662. int32x2x2_t result = vzip_s32(a1, b1);
  3663. return SIMDE__M128I_NEON_C(i32,
  3664. vcombine_s32(result.val[0], result.val[1]));
  3665. #else
  3666. simde__m128i r;
  3667. SIMDE__VECTORIZE
  3668. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i32[0])) / 2); i++) {
  3669. r.i32[(i * 2)] = a.i32[i];
  3670. r.i32[(i * 2) + 1] = b.i32[i];
  3671. }
  3672. return r;
  3673. #endif
  3674. }
  3675. SIMDE__FUNCTION_ATTRIBUTES
  3676. simde__m128i simde_mm_unpacklo_epi64(simde__m128i a, simde__m128i b)
  3677. {
  3678. #if defined(SIMDE_SSE2_NATIVE)
  3679. return SIMDE__M128I_C(_mm_unpacklo_epi64(a.n, b.n));
  3680. #else
  3681. simde__m128i r;
  3682. SIMDE__VECTORIZE
  3683. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.i64[0])) / 2); i++) {
  3684. r.i64[(i * 2)] = a.i64[i];
  3685. r.i64[(i * 2) + 1] = b.i64[i];
  3686. }
  3687. return r;
  3688. #endif
  3689. }
  3690. SIMDE__FUNCTION_ATTRIBUTES
  3691. simde__m128d simde_mm_unpacklo_pd(simde__m128d a, simde__m128d b)
  3692. {
  3693. #if defined(SIMDE_SSE2_NATIVE)
  3694. return SIMDE__M128D_C(_mm_unpacklo_pd(a.n, b.n));
  3695. #else
  3696. simde__m128d r;
  3697. SIMDE__VECTORIZE
  3698. for (size_t i = 0; i < ((sizeof(r) / sizeof(r.f64[0])) / 2); i++) {
  3699. r.f64[(i * 2)] = a.f64[i];
  3700. r.f64[(i * 2) + 1] = b.f64[i];
  3701. }
  3702. return r;
  3703. #endif
  3704. }
  3705. SIMDE__FUNCTION_ATTRIBUTES
  3706. simde__m128d simde_mm_xor_pd(simde__m128d a, simde__m128d b)
  3707. {
  3708. #if defined(SIMDE_SSE2_NATIVE)
  3709. return SIMDE__M128D_C(_mm_xor_pd(a.n, b.n));
  3710. #else
  3711. simde__m128d r;
  3712. SIMDE__VECTORIZE
  3713. for (size_t i = 0; i < (sizeof(r.i64) / sizeof(r.i64[0])); i++) {
  3714. r.i64[i] = a.i64[i] ^ b.i64[i];
  3715. }
  3716. return r;
  3717. #endif
  3718. }
  3719. SIMDE__FUNCTION_ATTRIBUTES
  3720. simde__m128i simde_mm_xor_si128(simde__m128i a, simde__m128i b)
  3721. {
  3722. #if defined(SIMDE_SSE2_NATIVE)
  3723. return SIMDE__M128I_C(_mm_xor_si128(a.n, b.n));
  3724. #elif defined(SIMDE_SSE2_NEON)
  3725. return SIMDE__M128I_NEON_C(i32, veorq_s32(a.neon_i32, b.neon_i32));
  3726. #else
  3727. simde__m128i r;
  3728. SIMDE__VECTORIZE
  3729. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  3730. r.i32[i] = a.i32[i] ^ b.i32[i];
  3731. }
  3732. return r;
  3733. #endif
  3734. }
  3735. SIMDE__FUNCTION_ATTRIBUTES
  3736. simde__m128i simde_x_mm_not_si128(simde__m128i a)
  3737. {
  3738. #if defined(SIMDE_SSE2_NEON)
  3739. return SIMDE__M128I_NEON_C(i32, vmvnq_s32(a.neon_i32));
  3740. #else
  3741. simde__m128i r;
  3742. SIMDE__VECTORIZE
  3743. for (size_t i = 0; i < (sizeof(r.i32) / sizeof(r.i32[0])); i++) {
  3744. r.i32[i] = ~(a.i32[i]);
  3745. }
  3746. return r;
  3747. #endif
  3748. }
  3749. SIMDE__END_DECLS
  3750. #endif /* !defined(SIMDE__SSE2_H) */