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.
- sequenzo/__init__.py +25 -1
- sequenzo/big_data/clara/clara.py +1 -1
- sequenzo/big_data/clara/utils/get_weighted_diss.c +156 -156
- sequenzo/big_data/clara/utils/get_weighted_diss.cp39-win_amd64.pyd +0 -0
- sequenzo/clustering/clustering_c_code.cp39-win_amd64.pyd +0 -0
- sequenzo/clustering/hierarchical_clustering.py +202 -8
- sequenzo/define_sequence_data.py +34 -2
- sequenzo/dissimilarity_measures/c_code.cp39-win_amd64.pyd +0 -0
- sequenzo/dissimilarity_measures/get_substitution_cost_matrix.py +1 -1
- sequenzo/dissimilarity_measures/src/DHDdistance.cpp +13 -37
- sequenzo/dissimilarity_measures/src/LCPdistance.cpp +13 -37
- sequenzo/dissimilarity_measures/src/OMdistance.cpp +12 -47
- sequenzo/dissimilarity_measures/src/OMspellDistance.cpp +103 -67
- sequenzo/dissimilarity_measures/src/dp_utils.h +160 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_arithmetic.hpp +41 -16
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_complex.hpp +4 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_details.hpp +7 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_logical.hpp +10 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_math.hpp +127 -43
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_memory.hpp +30 -2
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_swizzle.hpp +174 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/common/xsimd_common_trigo.hpp +14 -5
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx.hpp +111 -54
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx2.hpp +131 -9
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512bw.hpp +11 -113
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512dq.hpp +39 -7
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512f.hpp +336 -30
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512vbmi.hpp +9 -37
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_avx512vbmi2.hpp +58 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_common.hpp +1 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_common_fwd.hpp +35 -2
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_constants.hpp +3 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_emulated.hpp +17 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma3_avx.hpp +13 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma3_sse.hpp +18 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_fma4.hpp +13 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_isa.hpp +8 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_neon.hpp +363 -34
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_neon64.hpp +7 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_rvv.hpp +13 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_scalar.hpp +41 -4
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sse2.hpp +252 -16
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sse3.hpp +9 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_ssse3.hpp +12 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_sve.hpp +7 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_vsx.hpp +892 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/arch/xsimd_wasm.hpp +78 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_arch.hpp +3 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_config.hpp +13 -2
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_cpuid.hpp +5 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/config/xsimd_inline.hpp +5 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_all_registers.hpp +2 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_api.hpp +64 -1
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_batch.hpp +36 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_rvv_register.hpp +40 -31
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_traits.hpp +8 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/types/xsimd_vsx_register.hpp +77 -0
- sequenzo/dissimilarity_measures/src/xsimd/include/xsimd/xsimd.hpp +6 -0
- sequenzo/dissimilarity_measures/src/xsimd/test/test_basic_math.cpp +6 -0
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch.cpp +54 -2
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_bool.cpp +8 -0
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_cast.cpp +11 -4
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_complex.cpp +18 -0
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_int.cpp +8 -14
- sequenzo/dissimilarity_measures/src/xsimd/test/test_batch_manip.cpp +216 -173
- sequenzo/dissimilarity_measures/src/xsimd/test/test_load_store.cpp +6 -0
- sequenzo/dissimilarity_measures/src/xsimd/test/test_memory.cpp +1 -1
- sequenzo/dissimilarity_measures/src/xsimd/test/test_power.cpp +7 -4
- sequenzo/dissimilarity_measures/src/xsimd/test/test_select.cpp +6 -2
- sequenzo/dissimilarity_measures/src/xsimd/test/test_shuffle.cpp +32 -18
- sequenzo/dissimilarity_measures/src/xsimd/test/test_utils.hpp +21 -24
- sequenzo/dissimilarity_measures/src/xsimd/test/test_xsimd_api.cpp +69 -9
- sequenzo/dissimilarity_measures/utils/get_sm_trate_substitution_cost_matrix.c +156 -156
- sequenzo/dissimilarity_measures/utils/get_sm_trate_substitution_cost_matrix.cp39-win_amd64.pyd +0 -0
- sequenzo/dissimilarity_measures/utils/seqconc.c +156 -156
- sequenzo/dissimilarity_measures/utils/seqconc.cp39-win_amd64.pyd +0 -0
- sequenzo/dissimilarity_measures/utils/seqdss.c +156 -156
- sequenzo/dissimilarity_measures/utils/seqdss.cp39-win_amd64.pyd +0 -0
- sequenzo/dissimilarity_measures/utils/seqdur.c +156 -156
- sequenzo/dissimilarity_measures/utils/seqdur.cp39-win_amd64.pyd +0 -0
- sequenzo/dissimilarity_measures/utils/seqlength.c +156 -156
- sequenzo/dissimilarity_measures/utils/seqlength.cp39-win_amd64.pyd +0 -0
- sequenzo/sequence_characteristics/__init__.py +4 -0
- sequenzo/sequence_characteristics/complexity_index.py +17 -57
- sequenzo/sequence_characteristics/overall_cross_sectional_entropy.py +177 -111
- sequenzo/sequence_characteristics/plot_characteristics.py +30 -11
- sequenzo/sequence_characteristics/simple_characteristics.py +1 -0
- sequenzo/sequence_characteristics/state_frequencies_and_entropy_per_sequence.py +9 -3
- sequenzo/sequence_characteristics/turbulence.py +47 -67
- sequenzo/sequence_characteristics/variance_of_spell_durations.py +19 -9
- sequenzo/sequence_characteristics/within_sequence_entropy.py +5 -58
- sequenzo/visualization/plot_sequence_index.py +58 -35
- sequenzo/visualization/plot_state_distribution.py +57 -36
- sequenzo/with_event_history_analysis/__init__.py +35 -0
- sequenzo/with_event_history_analysis/sequence_analysis_multi_state_model.py +850 -0
- sequenzo/with_event_history_analysis/sequence_history_analysis.py +283 -0
- {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/METADATA +7 -6
- {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/RECORD +101 -94
- {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/WHEEL +0 -0
- {sequenzo-0.1.17.dist-info → sequenzo-0.1.18.dist-info}/licenses/LICENSE +0 -0
- {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
|
-
|
|
20
|
-
struct init_swizzle_base
|
|
21
|
+
namespace kernel
|
|
21
22
|
{
|
|
22
|
-
|
|
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
|
-
|
|
29
|
-
|
|
30
|
-
|
|
31
|
-
|
|
32
|
-
|
|
33
|
-
|
|
34
|
-
|
|
35
|
-
|
|
36
|
-
|
|
37
|
-
|
|
38
|
-
|
|
39
|
-
|
|
40
|
-
|
|
41
|
-
|
|
42
|
-
|
|
43
|
-
|
|
44
|
-
|
|
45
|
-
|
|
46
|
-
|
|
47
|
-
|
|
48
|
-
|
|
49
|
-
|
|
50
|
-
|
|
51
|
-
|
|
52
|
-
|
|
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
|
-
|
|
58
|
-
struct Reversor
|
|
66
|
+
namespace xsimd
|
|
59
67
|
{
|
|
60
|
-
|
|
68
|
+
template <template <class> class Pattern, class Vec>
|
|
69
|
+
void fill_pattern(Vec& dst, const Vec& src)
|
|
61
70
|
{
|
|
62
|
-
|
|
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
|
|
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
|
|
77
|
-
{
|
|
78
|
-
static constexpr T get(T i, T)
|
|
94
|
+
template <class T>
|
|
95
|
+
struct DupImag
|
|
79
96
|
{
|
|
80
|
-
|
|
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
|
|
86
|
-
{
|
|
87
|
-
|
|
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
|
|
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,
|
|
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,
|
|
121
|
-
|
|
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
|
|
225
|
+
static constexpr std::size_t N = B::size;
|
|
226
|
+
using vec_t = std::array<value_type, N>;
|
|
150
227
|
|
|
151
|
-
|
|
228
|
+
// Build the input [1,3,5,...,2N-1]
|
|
229
|
+
static vec_t make_lhs()
|
|
152
230
|
{
|
|
153
|
-
|
|
154
|
-
|
|
155
|
-
|
|
156
|
-
|
|
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
|
-
|
|
237
|
+
template <template <class> class Pattern>
|
|
238
|
+
void run()
|
|
166
239
|
{
|
|
167
|
-
|
|
168
|
-
|
|
169
|
-
|
|
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
|
-
|
|
173
|
-
|
|
244
|
+
auto b_lhs = batch_type::load_unaligned(lhs.data());
|
|
245
|
+
auto b_expect = batch_type::load_unaligned(expect.data());
|
|
174
246
|
|
|
175
|
-
|
|
176
|
-
|
|
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
|
|
256
|
+
void rotate_right()
|
|
180
257
|
{
|
|
181
|
-
|
|
182
|
-
|
|
183
|
-
|
|
184
|
-
|
|
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
|
-
|
|
202
|
-
|
|
203
|
-
|
|
204
|
-
|
|
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
|
-
|
|
222
|
-
|
|
223
|
-
|
|
224
|
-
|
|
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
|
-
|
|
243
|
-
|
|
244
|
-
|
|
245
|
-
|
|
246
|
-
|
|
247
|
-
|
|
248
|
-
|
|
249
|
-
|
|
250
|
-
|
|
251
|
-
|
|
252
|
-
|
|
253
|
-
|
|
254
|
-
|
|
255
|
-
|
|
256
|
-
|
|
257
|
-
|
|
258
|
-
|
|
259
|
-
SUBCASE("
|
|
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
|
-
#
|
|
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
|
|
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
|
|
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] =
|
|
40
|
-
rhs_input[i] =
|
|
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);
|