-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathAverage.h
95 lines (85 loc) · 3.17 KB
/
Average.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
// (c) 2016 Neeraj Sharma <[email protected]>.
// This file is part of SignalBlocks.
//
// SignalBlocks is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// SignalBlocks is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with SignalBlocks. If not, see <http://www.gnu.org/licenses/>.
//
#ifndef SIGNALBLOCKS_AVERAGE_H
#define SIGNALBLOCKS_AVERAGE_H
#include "../../common/Port.h"
#include <cassert>
#include <stdint.h>
namespace signalblocks {
///
///
/// @caution No under/overflow checking is performed.
///
/** Average computes an average of the inputs.
* Computes an average of the inputs (over scalar, vector or matrix).
*
* For scalar the average is the value itself
* and for vector & matrix the output is the
* average of all the elements.
*
* The output of this block is always scalar.
*/
template<class T>
class Average
: public Port<1, 1, T> {
public:
Average(std::string name)
: Port<1, 1, T>(std::move(name), "An averaging block.") {
}
protected: // Port interface
virtual void Process(int sourceIndex, const T& data, const TimeTick& startTime) {
// average of single value is the value itself
this->LeakData(0, data, startTime);
}
virtual void Process(
int sourceIndex, std::unique_ptr<T[]> data, int len, const TimeTick& startTime) {
assert(sourceIndex == 0);
assert(len > 0);
T* raw_data = data.get();
T data_sum = raw_data[0];
for (int i = 1; i < len; ++i) {
data_sum += raw_data[i];
}
T result = data_sum / len;
this->LeakData(0, result, startTime);
}
virtual void ProcessMatrix(int sourceIndex,
std::unique_ptr<T[]> data,
const std::vector<int>& dims,
const TimeTick& startTime) {
assert(sourceIndex == 0);
assert(! dims.empty());
// Matrix are stored in row-major format. See IPort.h for details.
int len = 1;
//printf("dims=[");
for (size_t i = 0; i < (dims.size() - 1); ++i) {
//printf("%d, ", dims[i]);
len *= dims[i];
}
//printf("%d]\n", dims[dims.size()-1]);
assert(len > 0);
T* raw_data = data.get();
T data_sum = raw_data[0];
for (int i = 1; i < len; ++i) {
data_sum += raw_data[i];
}
T result = data_sum / len;
this->LeakData(0, result, startTime);
}
};
}
#endif // SIGNALBLOCKS_AVERAGE_H