-
Notifications
You must be signed in to change notification settings - Fork 1
/
StageProcF64Sse2.hpp
112 lines (63 loc) · 2.76 KB
/
StageProcF64Sse2.hpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
/*****************************************************************************
StageProcF64Sse2.hpp
Author: Laurent de Soras, 2020
--- Legal stuff ---
This program is free software. It comes without any warranty, to
the extent permitted by applicable law. You can redistribute it
and/or modify it under the terms of the Do What The Fuck You Want
To Public License, Version 2, as published by Sam Hocevar. See
http://www.wtfpl.net/ for more details.
*Tab=3***********************************************************************/
#if ! defined (hiir_StageProcF64Sse2_CODEHEADER_INCLUDED)
#define hiir_StageProcF64Sse2_CODEHEADER_INCLUDED
/*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
namespace hiir
{
/*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
template <int CUR>
void StageProcF64Sse2 <CUR>::process_sample_pos (__m128d &x, StageDataF64Sse2 *stage_arr)
{
StageProcF64Sse2 <CUR>::process_sample_pos_rec (x, stage_arr);
_mm_store_pd (stage_arr [CUR]._mem, x);
}
template <int CUR>
void StageProcF64Sse2 <CUR>::process_sample_neg (__m128d &x, StageDataF64Sse2 *stage_arr)
{
StageProcF64Sse2 <CUR>::process_sample_neg_rec (x, stage_arr);
_mm_store_pd (stage_arr [CUR]._mem, x);
}
template <int CUR>
void StageProcF64Sse2 <CUR>::process_sample_pos_rec (__m128d &x, StageDataF64Sse2 *stage_arr)
{
StageProcF64Sse2 <CUR - 1>::process_sample_pos_rec (x, stage_arr);
const auto tmp = _mm_load_pd (stage_arr [CUR - 1]._mem);
_mm_store_pd (stage_arr [CUR - 1]._mem, x);
x = _mm_sub_pd (x, _mm_load_pd (stage_arr [CUR]._mem ));
x = _mm_mul_pd (x, _mm_load_pd (stage_arr [CUR]._coef));
x = _mm_add_pd (x, tmp);
}
template <>
void StageProcF64Sse2 <0>::process_sample_pos_rec (__m128d & /*x*/, StageDataF64Sse2 * /*stage_arr*/)
{
// Nothing, stops the recursion
}
template <int CUR>
void StageProcF64Sse2 <CUR>::process_sample_neg_rec (__m128d &x, StageDataF64Sse2 *stage_arr)
{
StageProcF64Sse2 <CUR - 1>::process_sample_neg_rec (x, stage_arr);
const auto tmp = _mm_load_pd (stage_arr [CUR - 1]._mem);
_mm_store_pd (stage_arr [CUR - 1]._mem, x);
x = _mm_add_pd (x, _mm_load_pd (stage_arr [CUR]._mem ));
x = _mm_mul_pd (x, _mm_load_pd (stage_arr [CUR]._coef));
x = _mm_sub_pd (x, tmp);
}
template <>
void StageProcF64Sse2 <0>::process_sample_neg_rec (__m128d & /*x*/, StageDataF64Sse2 * /*stage_arr*/)
{
// Nothing, stops the recursion
}
/*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
/*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
} // namespace hiir
#endif // hiir_StageProcF64Sse2_CODEHEADER_INCLUDED
/*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/