-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathPhaseHalfPi4F64Avx.h
129 lines (68 loc) · 3.35 KB
/
PhaseHalfPi4F64Avx.h
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
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
/*****************************************************************************
PhaseHalfPi4F64Avx.h
Author: Laurent de Soras, 2020
From the input signal, generates two signals with a pi/2 phase shift, using
AVX instruction set. Works on vectors of 4 double.
This object must be aligned on a 32-byte boundary!
Template parameters:
- NC: number of coefficients, > 0
--- 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***********************************************************************/
#pragma once
#if ! defined (hiir_PhaseHalfPi4F64Avx_HEADER_INCLUDED)
#define hiir_PhaseHalfPi4F64Avx_HEADER_INCLUDED
/*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
#include "hiir/def.h"
#include "hiir/StageDataF64Avx.h"
#include <immintrin.h>
#include <array>
namespace hiir
{
template <int NC>
class PhaseHalfPi4F64Avx
{
static_assert ((NC > 0), "Number of coefficient must be positive.");
/*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
public:
typedef double DataType;
static constexpr int _nbr_chn = 4;
static constexpr int NBR_COEFS = NC;
static constexpr double _delay = 0;
PhaseHalfPi4F64Avx () noexcept;
PhaseHalfPi4F64Avx (const PhaseHalfPi4F64Avx <NC> &other) = default;
PhaseHalfPi4F64Avx (PhaseHalfPi4F64Avx <NC> &&other) = default;
~PhaseHalfPi4F64Avx () = default;
PhaseHalfPi4F64Avx <NC> &
operator = (const PhaseHalfPi4F64Avx <NC> &other) = default;
PhaseHalfPi4F64Avx <NC> &
operator = (PhaseHalfPi4F64Avx <NC> &&other) = default;
void set_coefs (const double coef_arr []) noexcept;
hiir_FORCEINLINE void
process_sample (__m256d &out_0, __m256d &out_1, __m256d input) noexcept;
void process_block (double out_0_ptr [], double out_1_ptr [], const double in_ptr [], long nbr_spl) noexcept;
void clear_buffers () noexcept;
/*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
protected:
/*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
private:
static constexpr int _nbr_phases = 2;
typedef std::array <StageDataF64Avx, NBR_COEFS + 2> Filter; // Stages 0 and 1 contain only input memories
typedef std::array <Filter, _nbr_phases> FilterBiPhase;
FilterBiPhase _bifilter;
alignas (32) double
_prev [_nbr_chn];
int _phase; // 0 or 1
/*\\\ FORBIDDEN MEMBER FUNCTIONS \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
private:
bool operator == (const PhaseHalfPi4F64Avx &other) const = delete;
bool operator != (const PhaseHalfPi4F64Avx &other) const = delete;
}; // class PhaseHalfPi4F64Avx
} // namespace hiir
#include "hiir/PhaseHalfPi4F64Avx.hpp"
#endif // hiir_PhaseHalfPi4F64Avx_HEADER_INCLUDED
/*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/