sequenzo 0.1.17__cp39-cp39-win_amd64.whl → 0.1.18__cp39-cp39-win_amd64.whl

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.

Potentially problematic release.


This version of sequenzo might be problematic. Click here for more details.

Files changed (101) hide show
  1. sequenzo/__init__.py +25 -1
  2. sequenzo/big_data/clara/clara.py +1 -1
  3. sequenzo/big_data/clara/utils/get_weighted_diss.c +156 -156
  4. sequenzo/big_data/clara/utils/get_weighted_diss.cp39-win_amd64.pyd +0 -0
  5. sequenzo/clustering/clustering_c_code.cp39-win_amd64.pyd +0 -0
  6. sequenzo/clustering/hierarchical_clustering.py +202 -8
  7. sequenzo/define_sequence_data.py +34 -2
  8. sequenzo/dissimilarity_measures/c_code.cp39-win_amd64.pyd +0 -0
  9. sequenzo/dissimilarity_measures/get_substitution_cost_matrix.py +1 -1
  10. sequenzo/dissimilarity_measures/src/DHDdistance.cpp +13 -37
  11. sequenzo/dissimilarity_measures/src/LCPdistance.cpp +13 -37
  12. sequenzo/dissimilarity_measures/src/OMdistance.cpp +12 -47
  13. sequenzo/dissimilarity_measures/src/OMspellDistance.cpp +103 -67
  14. sequenzo/dissimilarity_measures/src/dp_utils.h +160 -0
  15. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_arithmetic.hpp +41 -16
  16. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_complex.hpp +4 -0
  17. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_details.hpp +7 -0
  18. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_logical.hpp +10 -0
  19. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_math.hpp +127 -43
  20. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_memory.hpp +30 -2
  21. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_swizzle.hpp +174 -0
  22. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_trigo.hpp +14 -5
  23. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx.hpp +111 -54
  24. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx2.hpp +131 -9
  25. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512bw.hpp +11 -113
  26. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512dq.hpp +39 -7
  27. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512f.hpp +336 -30
  28. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512vbmi.hpp +9 -37
  29. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512vbmi2.hpp +58 -0
  30. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_common.hpp +1 -0
  31. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_common_fwd.hpp +35 -2
  32. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_constants.hpp +3 -1
  33. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_emulated.hpp +17 -0
  34. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma3_avx.hpp +13 -0
  35. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma3_sse.hpp +18 -0
  36. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma4.hpp +13 -0
  37. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_isa.hpp +8 -0
  38. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_neon.hpp +363 -34
  39. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_neon64.hpp +7 -0
  40. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_rvv.hpp +13 -0
  41. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_scalar.hpp +41 -4
  42. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sse2.hpp +252 -16
  43. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sse3.hpp +9 -0
  44. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_ssse3.hpp +12 -1
  45. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sve.hpp +7 -0
  46. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_vsx.hpp +892 -0
  47. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_wasm.hpp +78 -1
  48. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_arch.hpp +3 -1
  49. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_config.hpp +13 -2
  50. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_cpuid.hpp +5 -0
  51. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_inline.hpp +5 -1
  52. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_all_registers.hpp +2 -0
  53. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_api.hpp +64 -1
  54. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_batch.hpp +36 -0
  55. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_rvv_register.hpp +40 -31
  56. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_traits.hpp +8 -0
  57. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_vsx_register.hpp +77 -0
  58. sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/xsimd.hpp +6 -0
  59. sequenzo/dissimilarity_measures/src/xsimd/test/test_basic_math.cpp +6 -0
  60. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch.cpp +54 -2
  61. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_bool.cpp +8 -0
  62. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_cast.cpp +11 -4
  63. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_complex.cpp +18 -0
  64. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_int.cpp +8 -14
  65. sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_manip.cpp +216 -173
  66. sequenzo/dissimilarity_measures/src/xsimd/test/test_load_store.cpp +6 -0
  67. sequenzo/dissimilarity_measures/src/xsimd/test/test_memory.cpp +1 -1
  68. sequenzo/dissimilarity_measures/src/xsimd/test/test_power.cpp +7 -4
  69. sequenzo/dissimilarity_measures/src/xsimd/test/test_select.cpp +6 -2
  70. sequenzo/dissimilarity_measures/src/xsimd/test/test_shuffle.cpp +32 -18
  71. sequenzo/dissimilarity_measures/src/xsimd/test/test_utils.hpp +21 -24
  72. sequenzo/dissimilarity_measures/src/xsimd/test/test_xsimd_api.cpp +69 -9
  73. sequenzo/dissimilarity_measures/utils/get_sm_trate_substitution_cost_matrix.c +156 -156
  74. sequenzo/dissimilarity_measures/utils/get_sm_trate_substitution_cost_matrix.cp39-win_amd64.pyd +0 -0
  75. sequenzo/dissimilarity_measures/utils/seqconc.c +156 -156
  76. sequenzo/dissimilarity_measures/utils/seqconc.cp39-win_amd64.pyd +0 -0
  77. sequenzo/dissimilarity_measures/utils/seqdss.c +156 -156
  78. sequenzo/dissimilarity_measures/utils/seqdss.cp39-win_amd64.pyd +0 -0
  79. sequenzo/dissimilarity_measures/utils/seqdur.c +156 -156
  80. sequenzo/dissimilarity_measures/utils/seqdur.cp39-win_amd64.pyd +0 -0
  81. sequenzo/dissimilarity_measures/utils/seqlength.c +156 -156
  82. sequenzo/dissimilarity_measures/utils/seqlength.cp39-win_amd64.pyd +0 -0
  83. sequenzo/sequence_characteristics/__init__.py +4 -0
  84. sequenzo/sequence_characteristics/complexity_index.py +17 -57
  85. sequenzo/sequence_characteristics/overall_cross_sectional_entropy.py +177 -111
  86. sequenzo/sequence_characteristics/plot_characteristics.py +30 -11
  87. sequenzo/sequence_characteristics/simple_characteristics.py +1 -0
  88. sequenzo/sequence_characteristics/state_frequencies_and_entropy_per_sequence.py +9 -3
  89. sequenzo/sequence_characteristics/turbulence.py +47 -67
  90. sequenzo/sequence_characteristics/variance_of_spell_durations.py +19 -9
  91. sequenzo/sequence_characteristics/within_sequence_entropy.py +5 -58
  92. sequenzo/visualization/plot_sequence_index.py +58 -35
  93. sequenzo/visualization/plot_state_distribution.py +57 -36
  94. sequenzo/with_event_history_analysis/__init__.py +35 -0
  95. sequenzo/with_event_history_analysis/sequence_analysis_multi_state_model.py +850 -0
  96. sequenzo/with_event_history_analysis/sequence_history_analysis.py +283 -0
  97. {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/METADATA +7 -6
  98. {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/RECORD +101 -94
  99. {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/WHEEL +0 -0
  100. {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/licenses/LICENSE +0 -0
  101. {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/top_level.txt +0 -0
@@ -3,6 +3,7 @@
3
3
  * Martin Renou *
4
4
  * Copyright (c) QuantStack *
5
5
  * Copyright (c) Serge Guelton *
6
+ * Copyright (c) Marco Barbone *
6
7
  * *
7
8
  * Distributed under the terms of the BSD 3-Clause License. *
8
9
  * *
@@ -14,101 +15,182 @@
14
15
 
15
16
  #include "test_utils.hpp"
16
17
 
18
+ // Compile time tests for include/xsimd/arch/common/xsimd_common_swizzle.hpp
17
19
  namespace xsimd
18
20
  {
19
- template <typename T, std::size_t N>
20
- struct init_swizzle_base
21
+ namespace kernel
21
22
  {
22
- using swizzle_vector_type = std::array<T, N>;
23
- swizzle_vector_type lhs_in, exped_reverse, exped_fill, exped_dup, exped_ror, exped_rol;
24
-
25
- template <int... Indices>
26
- std::vector<swizzle_vector_type> create_swizzle_vectors()
23
+ namespace detail
27
24
  {
28
- std::vector<swizzle_vector_type> vects;
29
-
30
- /* Generate input data */
31
- for (std::size_t i = 0; i < N; ++i)
32
- {
33
- lhs_in[i] = static_cast<T>(2 * i + 1);
34
- }
35
- vects.push_back(std::move(lhs_in));
36
-
37
- /* Expected reversed data */
38
- for (std::size_t i = 0; i < N; ++i)
39
- {
40
- exped_reverse[i] = lhs_in[N - 1 - i];
41
- exped_fill[i] = lhs_in[N - 1];
42
- exped_dup[i] = lhs_in[2 * (i / 2)];
43
- exped_ror[i] = lhs_in[(i - 1) % N];
44
- exped_rol[i] = lhs_in[(i + 1) % N];
45
- }
46
- vects.push_back(std::move(exped_reverse));
47
- vects.push_back(std::move(exped_fill));
48
- vects.push_back(std::move(exped_dup));
49
- vects.push_back(std::move(exped_ror));
50
- vects.push_back(std::move(exped_rol));
51
-
52
- return vects;
25
+ // ────────────────────────────────────────────────────────────────────────
26
+ // compile-time tests (identity, all-different, dup-lo, dup-hi)
27
+ // 8-lane identity
28
+ static_assert(is_identity<std::uint32_t, 0, 1, 2, 3, 4, 5, 6, 7>(), "identity failed");
29
+ // 8-lane reverse is all-different but not identity
30
+ static_assert(is_all_different<std::uint32_t, 7, 6, 5, 4, 3, 2, 1, 0>(), "all-diff failed");
31
+ static_assert(!is_identity<std::uint32_t, 7, 6, 5, 4, 3, 2, 1, 0>(), "identity on reverse");
32
+ // 8-lane dup-lo (repeat 0..3 twice)
33
+ static_assert(is_dup_lo<std::uint32_t, 0, 1, 2, 3, 0, 1, 2, 3>(), "dup_lo failed");
34
+ static_assert(!is_dup_hi<std::uint32_t, 0, 1, 2, 3, 0, 1, 2, 3>(), "dup_hi on dup_lo");
35
+ // 8-lane dup-hi (repeat 4..7 twice)
36
+ static_assert(is_dup_hi<std::uint32_t, 4, 5, 6, 7, 4, 5, 6, 7>(), "dup_hi failed");
37
+ static_assert(!is_dup_lo<std::uint32_t, 4, 5, 6, 7, 4, 5, 6, 7>(), "dup_lo on dup_hi");
38
+ // ────────────────────────────────────────────────────────────────────────
39
+ // 4-lane identity
40
+ static_assert(is_identity<std::uint32_t, 0, 1, 2, 3>(), "4-lane identity failed");
41
+ // 4-lane reverse all-different but not identity
42
+ static_assert(is_all_different<std::uint32_t, 3, 2, 1, 0>(), "4-lane all-diff failed");
43
+ static_assert(!is_identity<std::uint32_t, 3, 2, 1, 0>(), "4-lane identity on reverse");
44
+ // 4-lane dup-lo (repeat 0..1 twice)
45
+ static_assert(is_dup_lo<std::uint32_t, 0, 1, 0, 1>(), "4-lane dup_lo failed");
46
+ static_assert(!is_dup_hi<std::uint32_t, 0, 1, 0, 1>(), "4-lane dup_hi on dup_lo");
47
+ // 4-lane dup-hi (repeat 2..3 twice)
48
+ static_assert(is_dup_hi<std::uint32_t, 2, 3, 2, 3>(), "4-lane dup_hi failed");
49
+ static_assert(!is_dup_lo<std::uint32_t, 2, 3, 2, 3>(), "4-lane dup_lo on dup_hi");
50
+
51
+ static_assert(is_cross_lane<0, 1, 0, 1>(), "dup-lo only → crossing");
52
+ static_assert(is_cross_lane<2, 3, 2, 3>(), "dup-hi only → crossing");
53
+ static_assert(is_cross_lane<0, 3, 3, 3>(), "one low + rest high → crossing");
54
+ static_assert(!is_cross_lane<1, 0, 2, 3>(), "mixed low/high → no crossing");
55
+ static_assert(!is_cross_lane<0, 1, 2, 3>(), "mixed low/high → no crossing");
56
+
57
+ static_assert(no_duplicates_v<0, 1, 2, 3>(), "N=4: [0,1,2,3] → distinct");
58
+ static_assert(!no_duplicates_v<0, 1, 2, 2>(), "N=4: [0,1,2,2] → dup");
59
+
60
+ static_assert(no_duplicates_v<0, 1, 2, 3, 4, 5, 6, 7>(), "N=8: [0..7] → distinct");
61
+ static_assert(!no_duplicates_v<0, 1, 2, 3, 4, 5, 6, 0>(), "N=8: last repeats 0");
53
62
  }
54
- };
63
+ }
55
64
  }
56
65
 
57
- template <class T>
58
- struct Reversor
66
+ namespace xsimd
59
67
  {
60
- static constexpr T get(T i, T n)
68
+ template <template <class> class Pattern, class Vec>
69
+ void fill_pattern(Vec& dst, const Vec& src)
61
70
  {
62
- return n - 1 - i;
71
+ using size_type = typename Vec::size_type;
72
+ for (size_type i = 0; i < src.size(); ++i)
73
+ {
74
+ dst[i] = src[Pattern<size_type>::get(i, static_cast<size_type>(src.size()))];
75
+ }
63
76
  }
64
- };
65
77
 
66
- template <class T>
67
- struct Last
68
- {
69
- static constexpr T get(T, T n)
78
+ template <class T>
79
+ struct Reversor
70
80
  {
71
- return n - 1;
72
- }
73
- };
81
+ static constexpr T get(T i, T n) { return n - 1 - i; }
82
+ };
83
+ template <class T>
84
+ struct Last
85
+ {
86
+ static constexpr T get(T, T n) { return n - 1; }
87
+ };
88
+ template <class T>
89
+ struct DupReal
90
+ {
91
+ static constexpr T get(T i, T) { return (i & ~T { 1 }); }
92
+ };
74
93
 
75
- template <class T>
76
- struct Dup
77
- {
78
- static constexpr T get(T i, T)
94
+ template <class T>
95
+ struct DupImag
79
96
  {
80
- return 2 * (i / 2);
81
- }
82
- };
97
+ static constexpr T get(T i, T) { return (i & ~T { 1 }) + 1; }
98
+ };
99
+ template <class T>
100
+ struct SwapRI
101
+ {
102
+ static constexpr T get(T i, T)
103
+ {
104
+ return i ^ T { 1 };
105
+ }
106
+ };
107
+ template <class T>
108
+ struct Identity
109
+ {
110
+ static constexpr T get(T i, T) { return i; }
111
+ };
112
+ template <class T>
113
+ struct DupLowPair
114
+ {
115
+ static constexpr T get(T i, T) { return i / 2; }
116
+ };
117
+ template <class T>
118
+ struct DupHighPair
119
+ {
120
+ static constexpr T get(T i, T n) { return n / 2 + i / 2; }
121
+ };
83
122
 
84
- template <class T>
85
- struct as_index
86
- {
87
- using type = xsimd::as_unsigned_integer_t<T>;
88
- };
123
+ template <class T>
124
+ struct RotateRight1
125
+ {
126
+ static constexpr T get(T i, T n) { return (i + n - 1) % n; }
127
+ };
128
+ template <class T>
129
+ struct RotateLeft1
130
+ {
131
+ static constexpr T get(T i, T n) { return (i + 1) % n; }
132
+ };
89
133
 
90
- template <class T>
91
- struct as_index<std::complex<T>> : as_index<T>
92
- {
93
- };
134
+ template <class T>
135
+ struct ReversePairs
136
+ {
137
+ static constexpr T get(T i, T) { return (i & ~T { 1 }) | (1 - (i & T { 1 })); }
138
+ };
139
+ template <class T>
140
+ struct EvenThenOdd
141
+ {
142
+ static constexpr T get(T i, T n)
143
+ {
144
+ return (i < n / 2 ? 2 * i : 2 * (i - n / 2) + 1);
145
+ }
146
+ };
147
+ template <class T>
148
+ struct OddThenEven
149
+ {
150
+ static constexpr T get(T i, T n)
151
+ {
152
+ return (i < n / 2 ? 2 * i + 1 : 2 * (i - n / 2));
153
+ }
154
+ };
155
+ template <class T>
156
+ struct InterleavePairs
157
+ {
158
+ static constexpr T get(T i, T n)
159
+ {
160
+ return (i & 1) ? (i / 2 + n / 2) : (i / 2);
161
+ }
162
+ };
163
+ template <class T>
164
+ struct as_index
165
+ {
166
+ using type = xsimd::as_unsigned_integer_t<T>;
167
+ };
94
168
 
169
+ template <class T>
170
+ struct as_index<std::complex<T>> : as_index<T>
171
+ {
172
+ };
173
+ } // namespace xsimd
174
+
175
+ //------------------------------------------------------------------------------
176
+ // insert_test: unchanged from original
177
+ //------------------------------------------------------------------------------
95
178
  template <class B>
96
179
  struct insert_test
97
180
  {
98
181
  using batch_type = B;
99
182
  using value_type = typename B::value_type;
100
- static constexpr size_t size = B::size;
101
183
 
102
184
  void insert_first()
103
185
  {
104
186
  value_type fill_value = 0;
105
187
  value_type sentinel_value = 1;
106
188
  batch_type v(fill_value);
107
- batch_type w = insert(v, sentinel_value, ::xsimd::index<0>());
108
- std::array<value_type, batch_type::size> data;
189
+ batch_type w = xsimd::insert(v, sentinel_value, xsimd::index<0>());
190
+ std::array<value_type, batch_type::size> data {};
109
191
  w.store_unaligned(data.data());
110
192
  CHECK_SCALAR_EQ(data.front(), sentinel_value);
111
- for (size_t i = 1; i < batch_type::size; ++i)
193
+ for (std::size_t i = 1; i < batch_type::size; ++i)
112
194
  CHECK_SCALAR_EQ(data[i], fill_value);
113
195
  }
114
196
 
@@ -117,10 +199,11 @@ struct insert_test
117
199
  value_type fill_value = 0;
118
200
  value_type sentinel_value = 1;
119
201
  batch_type v(fill_value);
120
- batch_type w = insert(v, sentinel_value, ::xsimd::index<batch_type::size - 1>());
121
- std::array<value_type, batch_type::size> data;
202
+ batch_type w = xsimd::insert(v, sentinel_value,
203
+ xsimd::index<batch_type::size - 1>());
204
+ std::array<value_type, batch_type::size> data {};
122
205
  w.store_unaligned(data.data());
123
- for (size_t i = 0; i < batch_type::size - 1; ++i)
206
+ for (std::size_t i = 0; i < batch_type::size - 1; ++i)
124
207
  CHECK_SCALAR_EQ(data[i], fill_value);
125
208
  CHECK_SCALAR_EQ(data.back(), sentinel_value);
126
209
  }
@@ -129,15 +212,8 @@ struct insert_test
129
212
  TEST_CASE_TEMPLATE("[insert_test]", B, BATCH_TYPES)
130
213
  {
131
214
  insert_test<B> Test;
132
- SUBCASE("insert_first")
133
- {
134
- Test.insert_first();
135
- }
136
-
137
- SUBCASE("insert_last")
138
- {
139
- Test.insert_last();
140
- }
215
+ SUBCASE("insert_first") { Test.insert_first(); }
216
+ SUBCASE("insert_last") { Test.insert_last(); }
141
217
  }
142
218
 
143
219
  template <class B>
@@ -146,120 +222,87 @@ struct swizzle_test
146
222
  using batch_type = B;
147
223
  using value_type = typename B::value_type;
148
224
  using arch_type = typename B::arch_type;
149
- static constexpr size_t size = B::size;
225
+ static constexpr std::size_t N = B::size;
226
+ using vec_t = std::array<value_type, N>;
150
227
 
151
- void rotate_right()
228
+ // Build the input [1,3,5,...,2N-1]
229
+ static vec_t make_lhs()
152
230
  {
153
- xsimd::init_swizzle_base<value_type, size> swizzle_base;
154
- auto swizzle_vecs = swizzle_base.create_swizzle_vectors();
155
- auto v_lhs = swizzle_vecs[0];
156
- auto v_exped = swizzle_vecs[4];
157
-
158
- B b_lhs = B::load_unaligned(v_lhs.data());
159
- B b_exped = B::load_unaligned(v_exped.data());
160
-
161
- B b_res = xsimd::rotate_right<1>(b_lhs);
162
- CHECK_BATCH_EQ(b_res, b_exped);
231
+ vec_t v;
232
+ for (std::size_t i = 0; i < N; ++i)
233
+ v[i] = static_cast<value_type>(2 * i + 1);
234
+ return v;
163
235
  }
164
236
 
165
- void rotate_left()
237
+ template <template <class> class Pattern>
238
+ void run()
166
239
  {
167
- xsimd::init_swizzle_base<value_type, size> swizzle_base;
168
- auto swizzle_vecs = swizzle_base.create_swizzle_vectors();
169
- auto v_lhs = swizzle_vecs[0];
170
- auto v_exped = swizzle_vecs[5];
240
+ vec_t lhs = make_lhs();
241
+ vec_t expect = lhs;
242
+ xsimd::fill_pattern<Pattern>(expect, lhs);
171
243
 
172
- B b_lhs = B::load_unaligned(v_lhs.data());
173
- B b_exped = B::load_unaligned(v_exped.data());
244
+ auto b_lhs = batch_type::load_unaligned(lhs.data());
245
+ auto b_expect = batch_type::load_unaligned(expect.data());
174
246
 
175
- B b_res = xsimd::rotate_left<1>(b_lhs);
176
- CHECK_BATCH_EQ(b_res, b_exped);
247
+ using idx_t = typename xsimd::as_index<value_type>::type;
248
+ auto idx_batch = xsimd::make_batch_constant<idx_t, Pattern<idx_t>, arch_type>();
249
+
250
+ CHECK_BATCH_EQ(xsimd::swizzle(b_lhs, idx_batch), b_expect);
251
+ CHECK_BATCH_EQ(xsimd::swizzle(b_lhs,
252
+ static_cast<xsimd::batch<idx_t, arch_type>>(idx_batch)),
253
+ b_expect);
177
254
  }
178
255
 
179
- void swizzle_reverse()
256
+ void rotate_right()
180
257
  {
181
- xsimd::init_swizzle_base<value_type, size> swizzle_base;
182
- auto swizzle_vecs = swizzle_base.create_swizzle_vectors();
183
- auto v_lhs = swizzle_vecs[0];
184
- auto v_exped = swizzle_vecs[1];
185
-
186
- B b_lhs = B::load_unaligned(v_lhs.data());
187
- B b_exped = B::load_unaligned(v_exped.data());
188
-
189
- using index_type = typename as_index<value_type>::type;
190
- auto index_batch = xsimd::make_batch_constant<index_type, Reversor<index_type>, arch_type>();
191
-
192
- B b_res = xsimd::swizzle(b_lhs, index_batch);
193
- CHECK_BATCH_EQ(b_res, b_exped);
194
-
195
- B b_dyres = xsimd::swizzle(b_lhs, (xsimd::batch<index_type, arch_type>)index_batch);
196
- CHECK_BATCH_EQ(b_dyres, b_exped);
258
+ vec_t lhs = make_lhs(), expect;
259
+ std::rotate_copy(lhs.begin(), lhs.end() - 1, lhs.end(), expect.begin());
260
+ CHECK_BATCH_EQ(xsimd::rotate_right<1>(batch_type::load_unaligned(lhs.data())),
261
+ batch_type::load_unaligned(expect.data()));
197
262
  }
198
-
199
- void swizzle_fill()
263
+ void rotate_left()
200
264
  {
201
- xsimd::init_swizzle_base<value_type, size> swizzle_base;
202
- auto swizzle_vecs = swizzle_base.create_swizzle_vectors();
203
- auto v_lhs = swizzle_vecs[0];
204
- auto v_exped = swizzle_vecs[2];
205
-
206
- B b_lhs = B::load_unaligned(v_lhs.data());
207
- B b_exped = B::load_unaligned(v_exped.data());
208
-
209
- using index_type = typename as_index<value_type>::type;
210
- auto index_batch = xsimd::make_batch_constant<index_type, Last<index_type>, arch_type>();
211
-
212
- B b_res = xsimd::swizzle(b_lhs, index_batch);
213
- CHECK_BATCH_EQ(b_res, b_exped);
214
-
215
- B b_dyres = xsimd::swizzle(b_lhs, (xsimd::batch<index_type, arch_type>)index_batch);
216
- CHECK_BATCH_EQ(b_dyres, b_exped);
265
+ vec_t lhs = make_lhs(), expect;
266
+ std::rotate_copy(lhs.begin(), lhs.begin() + 1, lhs.end(), expect.begin());
267
+ CHECK_BATCH_EQ(xsimd::rotate_left<1>(batch_type::load_unaligned(lhs.data())),
268
+ batch_type::load_unaligned(expect.data()));
217
269
  }
218
-
219
- void swizzle_dup()
270
+ void rotate_left_inv()
220
271
  {
221
- xsimd::init_swizzle_base<value_type, size> swizzle_base;
222
- auto swizzle_vecs = swizzle_base.create_swizzle_vectors();
223
- auto v_lhs = swizzle_vecs[0];
224
- auto v_exped = swizzle_vecs[3];
225
-
226
- B b_lhs = B::load_unaligned(v_lhs.data());
227
- B b_exped = B::load_unaligned(v_exped.data());
228
-
229
- using index_type = typename as_index<value_type>::type;
230
- auto index_batch = xsimd::make_batch_constant<index_type, Dup<index_type>, arch_type>();
231
-
232
- B b_res = xsimd::swizzle(b_lhs, index_batch);
233
- CHECK_BATCH_EQ(b_res, b_exped);
234
-
235
- B b_dyres = xsimd::swizzle(b_lhs, (xsimd::batch<index_type, arch_type>)index_batch);
236
- CHECK_BATCH_EQ(b_dyres, b_exped);
272
+ vec_t lhs = make_lhs(), expect;
273
+ std::rotate_copy(lhs.begin(), lhs.end() - 1, lhs.end(), expect.begin());
274
+ CHECK_BATCH_EQ(xsimd::rotate_left<N - 1>(batch_type::load_unaligned(lhs.data())),
275
+ batch_type::load_unaligned(expect.data()));
237
276
  }
238
277
  };
239
278
 
279
+ // Macro to instantiate one SUBCASE per pattern
280
+ #define XSIMD_SWIZZLE_PATTERN_CASE(PAT) \
281
+ SUBCASE(#PAT) { swizzle_test<B>().template run<xsimd::PAT>(); }
282
+
240
283
  TEST_CASE_TEMPLATE("[swizzle]", B, BATCH_SWIZZLE_TYPES)
241
284
  {
242
- swizzle_test<B> Test;
243
- SUBCASE("reverse")
244
- {
245
- Test.swizzle_reverse();
246
- }
247
-
248
- SUBCASE("rotate")
249
- {
250
- Test.rotate_left();
251
- Test.rotate_right();
252
- }
253
-
254
- SUBCASE("fill")
255
- {
256
- Test.swizzle_fill();
257
- }
258
-
259
- SUBCASE("dup")
260
- {
261
- Test.swizzle_dup();
262
- }
285
+ // All existing patterns:
286
+ XSIMD_SWIZZLE_PATTERN_CASE(Reversor);
287
+ XSIMD_SWIZZLE_PATTERN_CASE(Last);
288
+ XSIMD_SWIZZLE_PATTERN_CASE(DupReal);
289
+ XSIMD_SWIZZLE_PATTERN_CASE(DupImag);
290
+ XSIMD_SWIZZLE_PATTERN_CASE(SwapRI);
291
+ XSIMD_SWIZZLE_PATTERN_CASE(Identity);
292
+ XSIMD_SWIZZLE_PATTERN_CASE(DupLowPair);
293
+ XSIMD_SWIZZLE_PATTERN_CASE(DupHighPair);
294
+ XSIMD_SWIZZLE_PATTERN_CASE(RotateRight1);
295
+ XSIMD_SWIZZLE_PATTERN_CASE(RotateLeft1);
296
+ XSIMD_SWIZZLE_PATTERN_CASE(ReversePairs);
297
+ XSIMD_SWIZZLE_PATTERN_CASE(EvenThenOdd);
298
+ XSIMD_SWIZZLE_PATTERN_CASE(OddThenEven);
299
+ XSIMD_SWIZZLE_PATTERN_CASE(InterleavePairs);
300
+ // Rotation checks:
301
+ SUBCASE("rotate_left") { swizzle_test<B>().rotate_left(); }
302
+ SUBCASE("rotate_left_inv") { swizzle_test<B>().rotate_left_inv(); }
303
+ SUBCASE("rotate_right") { swizzle_test<B>().rotate_right(); }
263
304
  }
264
305
 
265
- #endif
306
+ #undef XSIMD_SWIZZLE_PATTERN_CASE
307
+
308
+ #endif /* XSIMD_NO_SUPPORTED_ARCHITECTURE */
@@ -39,7 +39,9 @@ struct load_store_test
39
39
  using ulong_vector_type = std::vector<unsigned long, allocator<unsigned long>>;
40
40
  #endif
41
41
  using float_vector_type = std::vector<float, allocator<float>>;
42
+ #if !XSIMD_WITH_NEON || XSIMD_WITH_NEON64
42
43
  using double_vector_type = std::vector<double, allocator<double>>;
44
+ #endif
43
45
 
44
46
  int8_vector_type i8_vec;
45
47
  uint8_vector_type ui8_vec;
@@ -54,7 +56,9 @@ struct load_store_test
54
56
  ulong_vector_type ul_vec;
55
57
  #endif
56
58
  float_vector_type f_vec;
59
+ #if !XSIMD_WITH_NEON || XSIMD_WITH_NEON64
57
60
  double_vector_type d_vec;
61
+ #endif
58
62
 
59
63
  array_type expected;
60
64
 
@@ -73,7 +77,9 @@ struct load_store_test
73
77
  init_test_vector(ul_vec);
74
78
  #endif
75
79
  init_test_vector(f_vec);
80
+ #if !XSIMD_WITH_NEON || XSIMD_WITH_NEON64
76
81
  init_test_vector(d_vec);
82
+ #endif
77
83
  }
78
84
 
79
85
  void test_load()
@@ -46,7 +46,7 @@ TEST_CASE("[is_aligned]")
46
46
  float f[100];
47
47
  void* unaligned_f = static_cast<void*>(&f[0]);
48
48
  constexpr std::size_t alignment = xsimd::default_arch::alignment();
49
- std::size_t aligned_f_size;
49
+ std::size_t aligned_f_size = sizeof(f);
50
50
  void* aligned_f = std::align(alignment, sizeof(f), unaligned_f, aligned_f_size);
51
51
  CHECK_UNARY(xsimd::is_aligned(aligned_f));
52
52
 
@@ -13,6 +13,7 @@
13
13
  #ifndef XSIMD_NO_SUPPORTED_ARCHITECTURE
14
14
 
15
15
  #include "test_utils.hpp"
16
+ #include <iostream>
16
17
 
17
18
  template <class B>
18
19
  struct power_test
@@ -40,9 +41,9 @@ struct power_test
40
41
  for (size_t i = 0; i < nb_input; ++i)
41
42
  {
42
43
  zero_input[i] = 0;
43
- lhs_input[i] = value_type(i) / 4 + value_type(1.2) * std::sqrt(value_type(i + 0.25));
44
+ lhs_input[i] = value_type(i / 4 + 1.2 * std::sqrt(i + 0.25));
44
45
  zlhs_input[i] = lhs_input[i] * (i % 2);
45
- rhs_input[i] = value_type(10.2) / (i + 2) + value_type(0.25);
46
+ rhs_input[i] = value_type(10.2 / (i + 2) + 0.25);
46
47
  }
47
48
 
48
49
  expected.resize(nb_input);
@@ -103,6 +104,7 @@ struct power_test
103
104
  CHECK_EQ(diff, 0);
104
105
  #endif
105
106
  }
107
+ #ifndef __FAST_MATH__
106
108
  // pow 0^-x
107
109
  {
108
110
  std::transform(zero_input.cbegin(), zero_input.cend(), rhs_input.cbegin(), expected.begin(),
@@ -120,17 +122,18 @@ struct power_test
120
122
  INFO("pow(0, -x)");
121
123
  CHECK_EQ(diff, 0);
122
124
  }
125
+ #endif
123
126
  // ipow
124
127
  {
125
128
  long k = 0;
126
129
  std::transform(lhs_input.cbegin(), lhs_input.cend(), expected.begin(),
127
130
  [&k, this](const value_type& l)
128
- { auto arg = k / size - nb_input / size / 2; ++k; return std::pow(l, arg); });
131
+ { auto arg = k / size / 8000 - nb_input / size / 8000 / 2; ++k; return std::pow(l, arg); });
129
132
  batch_type lhs_in, out;
130
133
  for (size_t i = 0; i < nb_input; i += size)
131
134
  {
132
135
  detail::load_batch(lhs_in, lhs_input, i);
133
- out = pow(lhs_in, i / size - nb_input / size / 2);
136
+ out = pow(lhs_in, i / size / 8000 - nb_input / size / 8000 / 2);
134
137
  detail::store_batch(out, res, i);
135
138
  }
136
139
  size_t diff = detail::get_nb_diff(res, expected);
@@ -34,10 +34,14 @@ struct select_test
34
34
  nb_input = size * 10000;
35
35
  lhs_input.resize(nb_input);
36
36
  rhs_input.resize(nb_input);
37
+ auto clamp = [](double v)
38
+ {
39
+ return static_cast<value_type>(std::min(v, static_cast<double>(std::numeric_limits<value_type>::max())));
40
+ };
37
41
  for (size_t i = 0; i < nb_input; ++i)
38
42
  {
39
- lhs_input[i] = value_type(i) / 4 + value_type(1.2) * std::sqrt(value_type(i + 0.25));
40
- rhs_input[i] = value_type(10.2) / (i + 2) + value_type(0.25);
43
+ lhs_input[i] = clamp(i / 4 + 1.2 * std::sqrt(i + 0.25));
44
+ rhs_input[i] = clamp(10.2 / (i + 2) + 0.25);
41
45
  }
42
46
  expected.resize(nb_input);
43
47
  res.resize(nb_input);