-
Notifications
You must be signed in to change notification settings - Fork 549
Expand file tree
/
Copy pathmean.cpp
More file actions
98 lines (86 loc) · 3.15 KB
/
mean.cpp
File metadata and controls
98 lines (86 loc) · 3.15 KB
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
/*******************************************************
* Copyright (c) 2014, ArrayFire
* All rights reserved.
*
* This file is distributed under 3-clause BSD license.
* The complete license agreement can be obtained at:
* http://arrayfire.com/licenses/BSD-3-Clause
********************************************************/
#include <af/algorithm.h>
#include <af/array.h>
#include <af/dim4.hpp>
#include <af/statistics.h>
#include "common.hpp"
#include "error.hpp"
#include "half.hpp"
#ifdef AF_CUDA
#include <cuda_fp16.h>
#include <traits.hpp>
#endif
namespace af {
array mean(const array& in, const dim_t dim) {
af_array temp = 0;
AF_THROW(af_mean(&temp, in.get(), getFNSD(dim, in.dims())));
return array(temp);
}
array mean(const array& in, const array& weights, const dim_t dim) {
af_array temp = 0;
AF_THROW(af_mean_weighted(&temp, in.get(), weights.get(),
getFNSD(dim, in.dims())));
return array(temp);
}
#define INSTANTIATE_MEAN(T) \
template<> \
AFAPI T mean(const array& in) { \
double ret_val; \
AF_THROW(af_mean_all(&ret_val, NULL, in.get())); \
return cast<T>(ret_val); \
} \
template<> \
AFAPI T mean(const array& in, const array& wts) { \
double ret_val; \
AF_THROW(af_mean_all_weighted(&ret_val, NULL, in.get(), wts.get())); \
return cast<T>(ret_val); \
}
template<>
AFAPI af_cfloat mean(const array& in) {
double real, imag;
AF_THROW(af_mean_all(&real, &imag, in.get()));
return {static_cast<float>(real), static_cast<float>(imag)};
}
template<>
AFAPI af_cdouble mean(const array& in) {
double real, imag;
AF_THROW(af_mean_all(&real, &imag, in.get()));
return {real, imag};
}
template<>
AFAPI af_cfloat mean(const array& in, const array& weights) {
double real, imag;
AF_THROW(af_mean_all_weighted(&real, &imag, in.get(), weights.get()));
return {static_cast<float>(real), static_cast<float>(imag)};
}
template<>
AFAPI af_cdouble mean(const array& in, const array& weights) {
double real, imag;
AF_THROW(af_mean_all_weighted(&real, &imag, in.get(), weights.get()));
return {real, imag};
}
INSTANTIATE_MEAN(float);
INSTANTIATE_MEAN(double);
INSTANTIATE_MEAN(int);
INSTANTIATE_MEAN(unsigned int);
INSTANTIATE_MEAN(char);
INSTANTIATE_MEAN(signed char);
INSTANTIATE_MEAN(unsigned char);
INSTANTIATE_MEAN(long long);
INSTANTIATE_MEAN(unsigned long long);
INSTANTIATE_MEAN(short);
INSTANTIATE_MEAN(unsigned short);
INSTANTIATE_MEAN(af_half);
INSTANTIATE_MEAN(half_float::half); // Add support for public API
#ifdef AF_CUDA
INSTANTIATE_MEAN(__half);
#endif
#undef INSTANTIATE_MEAN
} // namespace af