Chris@16
|
1 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
2 // mean.hpp
|
Chris@16
|
3 //
|
Chris@16
|
4 // Copyright 2005 Eric Niebler. Distributed under the Boost
|
Chris@16
|
5 // Software License, Version 1.0. (See accompanying file
|
Chris@16
|
6 // LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)
|
Chris@16
|
7
|
Chris@16
|
8 #ifndef BOOST_ACCUMULATORS_STATISTICS_MEAN_HPP_EAN_28_10_2005
|
Chris@16
|
9 #define BOOST_ACCUMULATORS_STATISTICS_MEAN_HPP_EAN_28_10_2005
|
Chris@16
|
10
|
Chris@16
|
11 #include <boost/mpl/placeholders.hpp>
|
Chris@16
|
12 #include <boost/accumulators/framework/accumulator_base.hpp>
|
Chris@16
|
13 #include <boost/accumulators/framework/extractor.hpp>
|
Chris@16
|
14 #include <boost/accumulators/numeric/functional.hpp>
|
Chris@16
|
15 #include <boost/accumulators/framework/parameters/sample.hpp>
|
Chris@16
|
16 #include <boost/accumulators/framework/depends_on.hpp>
|
Chris@16
|
17 #include <boost/accumulators/statistics_fwd.hpp>
|
Chris@16
|
18 #include <boost/accumulators/statistics/count.hpp>
|
Chris@16
|
19 #include <boost/accumulators/statistics/sum.hpp>
|
Chris@16
|
20
|
Chris@16
|
21 namespace boost { namespace accumulators
|
Chris@16
|
22 {
|
Chris@16
|
23
|
Chris@16
|
24 namespace impl
|
Chris@16
|
25 {
|
Chris@16
|
26 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
27 // mean_impl
|
Chris@16
|
28 // lazy, by default
|
Chris@16
|
29 template<typename Sample, typename SumFeature>
|
Chris@16
|
30 struct mean_impl
|
Chris@16
|
31 : accumulator_base
|
Chris@16
|
32 {
|
Chris@16
|
33 // for boost::result_of
|
Chris@16
|
34 typedef typename numeric::functional::fdiv<Sample, std::size_t>::result_type result_type;
|
Chris@16
|
35
|
Chris@16
|
36 mean_impl(dont_care) {}
|
Chris@16
|
37
|
Chris@16
|
38 template<typename Args>
|
Chris@16
|
39 result_type result(Args const &args) const
|
Chris@16
|
40 {
|
Chris@16
|
41 extractor<SumFeature> sum;
|
Chris@16
|
42 return numeric::fdiv(sum(args), count(args));
|
Chris@16
|
43 }
|
Chris@16
|
44 };
|
Chris@16
|
45
|
Chris@16
|
46 template<typename Sample, typename Tag>
|
Chris@16
|
47 struct immediate_mean_impl
|
Chris@16
|
48 : accumulator_base
|
Chris@16
|
49 {
|
Chris@16
|
50 // for boost::result_of
|
Chris@16
|
51 typedef typename numeric::functional::fdiv<Sample, std::size_t>::result_type result_type;
|
Chris@16
|
52
|
Chris@16
|
53 template<typename Args>
|
Chris@16
|
54 immediate_mean_impl(Args const &args)
|
Chris@16
|
55 : mean(numeric::fdiv(args[sample | Sample()], numeric::one<std::size_t>::value))
|
Chris@16
|
56 {
|
Chris@16
|
57 }
|
Chris@16
|
58
|
Chris@16
|
59 template<typename Args>
|
Chris@16
|
60 void operator ()(Args const &args)
|
Chris@16
|
61 {
|
Chris@16
|
62 std::size_t cnt = count(args);
|
Chris@16
|
63 this->mean = numeric::fdiv(
|
Chris@16
|
64 (this->mean * (cnt - 1)) + args[parameter::keyword<Tag>::get()]
|
Chris@16
|
65 , cnt
|
Chris@16
|
66 );
|
Chris@16
|
67 }
|
Chris@16
|
68
|
Chris@16
|
69 result_type result(dont_care) const
|
Chris@16
|
70 {
|
Chris@16
|
71 return this->mean;
|
Chris@16
|
72 }
|
Chris@16
|
73
|
Chris@16
|
74 private:
|
Chris@16
|
75 result_type mean;
|
Chris@16
|
76 };
|
Chris@16
|
77
|
Chris@16
|
78 } // namespace impl
|
Chris@16
|
79
|
Chris@16
|
80 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
81 // tag::mean
|
Chris@16
|
82 // tag::immediate_mean
|
Chris@16
|
83 // tag::mean_of_weights
|
Chris@16
|
84 // tag::immediate_mean_of_weights
|
Chris@16
|
85 // tag::mean_of_variates
|
Chris@16
|
86 // tag::immediate_mean_of_variates
|
Chris@16
|
87 //
|
Chris@16
|
88 namespace tag
|
Chris@16
|
89 {
|
Chris@16
|
90 struct mean
|
Chris@16
|
91 : depends_on<count, sum>
|
Chris@16
|
92 {
|
Chris@16
|
93 /// INTERNAL ONLY
|
Chris@16
|
94 ///
|
Chris@16
|
95 typedef accumulators::impl::mean_impl<mpl::_1, sum> impl;
|
Chris@16
|
96 };
|
Chris@16
|
97 struct immediate_mean
|
Chris@16
|
98 : depends_on<count>
|
Chris@16
|
99 {
|
Chris@16
|
100 /// INTERNAL ONLY
|
Chris@16
|
101 ///
|
Chris@16
|
102 typedef accumulators::impl::immediate_mean_impl<mpl::_1, tag::sample> impl;
|
Chris@16
|
103 };
|
Chris@16
|
104 struct mean_of_weights
|
Chris@16
|
105 : depends_on<count, sum_of_weights>
|
Chris@16
|
106 {
|
Chris@16
|
107 typedef mpl::true_ is_weight_accumulator;
|
Chris@16
|
108 /// INTERNAL ONLY
|
Chris@16
|
109 ///
|
Chris@16
|
110 typedef accumulators::impl::mean_impl<mpl::_2, sum_of_weights> impl;
|
Chris@16
|
111 };
|
Chris@16
|
112 struct immediate_mean_of_weights
|
Chris@16
|
113 : depends_on<count>
|
Chris@16
|
114 {
|
Chris@16
|
115 typedef mpl::true_ is_weight_accumulator;
|
Chris@16
|
116 /// INTERNAL ONLY
|
Chris@16
|
117 ///
|
Chris@16
|
118 typedef accumulators::impl::immediate_mean_impl<mpl::_2, tag::weight> impl;
|
Chris@16
|
119 };
|
Chris@16
|
120 template<typename VariateType, typename VariateTag>
|
Chris@16
|
121 struct mean_of_variates
|
Chris@16
|
122 : depends_on<count, sum_of_variates<VariateType, VariateTag> >
|
Chris@16
|
123 {
|
Chris@16
|
124 /// INTERNAL ONLY
|
Chris@16
|
125 ///
|
Chris@16
|
126 typedef mpl::always<accumulators::impl::mean_impl<VariateType, sum_of_variates<VariateType, VariateTag> > > impl;
|
Chris@16
|
127 };
|
Chris@16
|
128 template<typename VariateType, typename VariateTag>
|
Chris@16
|
129 struct immediate_mean_of_variates
|
Chris@16
|
130 : depends_on<count>
|
Chris@16
|
131 {
|
Chris@16
|
132 /// INTERNAL ONLY
|
Chris@16
|
133 ///
|
Chris@16
|
134 typedef mpl::always<accumulators::impl::immediate_mean_impl<VariateType, VariateTag> > impl;
|
Chris@16
|
135 };
|
Chris@16
|
136 }
|
Chris@16
|
137
|
Chris@16
|
138 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
139 // extract::mean
|
Chris@16
|
140 // extract::mean_of_weights
|
Chris@16
|
141 // extract::mean_of_variates
|
Chris@16
|
142 //
|
Chris@16
|
143 namespace extract
|
Chris@16
|
144 {
|
Chris@16
|
145 extractor<tag::mean> const mean = {};
|
Chris@16
|
146 extractor<tag::mean_of_weights> const mean_of_weights = {};
|
Chris@16
|
147 BOOST_ACCUMULATORS_DEFINE_EXTRACTOR(tag, mean_of_variates, (typename)(typename))
|
Chris@16
|
148
|
Chris@16
|
149 BOOST_ACCUMULATORS_IGNORE_GLOBAL(mean)
|
Chris@16
|
150 BOOST_ACCUMULATORS_IGNORE_GLOBAL(mean_of_weights)
|
Chris@16
|
151 }
|
Chris@16
|
152
|
Chris@16
|
153 using extract::mean;
|
Chris@16
|
154 using extract::mean_of_weights;
|
Chris@16
|
155 using extract::mean_of_variates;
|
Chris@16
|
156
|
Chris@16
|
157 // mean(lazy) -> mean
|
Chris@16
|
158 template<>
|
Chris@16
|
159 struct as_feature<tag::mean(lazy)>
|
Chris@16
|
160 {
|
Chris@16
|
161 typedef tag::mean type;
|
Chris@16
|
162 };
|
Chris@16
|
163
|
Chris@16
|
164 // mean(immediate) -> immediate_mean
|
Chris@16
|
165 template<>
|
Chris@16
|
166 struct as_feature<tag::mean(immediate)>
|
Chris@16
|
167 {
|
Chris@16
|
168 typedef tag::immediate_mean type;
|
Chris@16
|
169 };
|
Chris@16
|
170
|
Chris@16
|
171 // mean_of_weights(lazy) -> mean_of_weights
|
Chris@16
|
172 template<>
|
Chris@16
|
173 struct as_feature<tag::mean_of_weights(lazy)>
|
Chris@16
|
174 {
|
Chris@16
|
175 typedef tag::mean_of_weights type;
|
Chris@16
|
176 };
|
Chris@16
|
177
|
Chris@16
|
178 // mean_of_weights(immediate) -> immediate_mean_of_weights
|
Chris@16
|
179 template<>
|
Chris@16
|
180 struct as_feature<tag::mean_of_weights(immediate)>
|
Chris@16
|
181 {
|
Chris@16
|
182 typedef tag::immediate_mean_of_weights type;
|
Chris@16
|
183 };
|
Chris@16
|
184
|
Chris@16
|
185 // mean_of_variates<VariateType, VariateTag>(lazy) -> mean_of_variates<VariateType, VariateTag>
|
Chris@16
|
186 template<typename VariateType, typename VariateTag>
|
Chris@16
|
187 struct as_feature<tag::mean_of_variates<VariateType, VariateTag>(lazy)>
|
Chris@16
|
188 {
|
Chris@16
|
189 typedef tag::mean_of_variates<VariateType, VariateTag> type;
|
Chris@16
|
190 };
|
Chris@16
|
191
|
Chris@16
|
192 // mean_of_variates<VariateType, VariateTag>(immediate) -> immediate_mean_of_variates<VariateType, VariateTag>
|
Chris@16
|
193 template<typename VariateType, typename VariateTag>
|
Chris@16
|
194 struct as_feature<tag::mean_of_variates<VariateType, VariateTag>(immediate)>
|
Chris@16
|
195 {
|
Chris@16
|
196 typedef tag::immediate_mean_of_variates<VariateType, VariateTag> type;
|
Chris@16
|
197 };
|
Chris@16
|
198
|
Chris@16
|
199 // for the purposes of feature-based dependency resolution,
|
Chris@16
|
200 // immediate_mean provides the same feature as mean
|
Chris@16
|
201 template<>
|
Chris@16
|
202 struct feature_of<tag::immediate_mean>
|
Chris@16
|
203 : feature_of<tag::mean>
|
Chris@16
|
204 {
|
Chris@16
|
205 };
|
Chris@16
|
206
|
Chris@16
|
207 // for the purposes of feature-based dependency resolution,
|
Chris@16
|
208 // immediate_mean provides the same feature as mean
|
Chris@16
|
209 template<>
|
Chris@16
|
210 struct feature_of<tag::immediate_mean_of_weights>
|
Chris@16
|
211 : feature_of<tag::mean_of_weights>
|
Chris@16
|
212 {
|
Chris@16
|
213 };
|
Chris@16
|
214
|
Chris@16
|
215 // for the purposes of feature-based dependency resolution,
|
Chris@16
|
216 // immediate_mean provides the same feature as mean
|
Chris@16
|
217 template<typename VariateType, typename VariateTag>
|
Chris@16
|
218 struct feature_of<tag::immediate_mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
219 : feature_of<tag::mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
220 {
|
Chris@16
|
221 };
|
Chris@16
|
222
|
Chris@16
|
223 // So that mean can be automatically substituted with
|
Chris@16
|
224 // weighted_mean when the weight parameter is non-void.
|
Chris@16
|
225 template<>
|
Chris@16
|
226 struct as_weighted_feature<tag::mean>
|
Chris@16
|
227 {
|
Chris@16
|
228 typedef tag::weighted_mean type;
|
Chris@16
|
229 };
|
Chris@16
|
230
|
Chris@16
|
231 template<>
|
Chris@16
|
232 struct feature_of<tag::weighted_mean>
|
Chris@16
|
233 : feature_of<tag::mean>
|
Chris@16
|
234 {};
|
Chris@16
|
235
|
Chris@16
|
236 // So that immediate_mean can be automatically substituted with
|
Chris@16
|
237 // immediate_weighted_mean when the weight parameter is non-void.
|
Chris@16
|
238 template<>
|
Chris@16
|
239 struct as_weighted_feature<tag::immediate_mean>
|
Chris@16
|
240 {
|
Chris@16
|
241 typedef tag::immediate_weighted_mean type;
|
Chris@16
|
242 };
|
Chris@16
|
243
|
Chris@16
|
244 template<>
|
Chris@16
|
245 struct feature_of<tag::immediate_weighted_mean>
|
Chris@16
|
246 : feature_of<tag::immediate_mean>
|
Chris@16
|
247 {};
|
Chris@16
|
248
|
Chris@16
|
249 // So that mean_of_weights<> can be automatically substituted with
|
Chris@16
|
250 // weighted_mean_of_variates<> when the weight parameter is non-void.
|
Chris@16
|
251 template<typename VariateType, typename VariateTag>
|
Chris@16
|
252 struct as_weighted_feature<tag::mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
253 {
|
Chris@16
|
254 typedef tag::weighted_mean_of_variates<VariateType, VariateTag> type;
|
Chris@16
|
255 };
|
Chris@16
|
256
|
Chris@16
|
257 template<typename VariateType, typename VariateTag>
|
Chris@16
|
258 struct feature_of<tag::weighted_mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
259 : feature_of<tag::mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
260 {
|
Chris@16
|
261 };
|
Chris@16
|
262
|
Chris@16
|
263 // So that immediate_mean_of_weights<> can be automatically substituted with
|
Chris@16
|
264 // immediate_weighted_mean_of_variates<> when the weight parameter is non-void.
|
Chris@16
|
265 template<typename VariateType, typename VariateTag>
|
Chris@16
|
266 struct as_weighted_feature<tag::immediate_mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
267 {
|
Chris@16
|
268 typedef tag::immediate_weighted_mean_of_variates<VariateType, VariateTag> type;
|
Chris@16
|
269 };
|
Chris@16
|
270
|
Chris@16
|
271 template<typename VariateType, typename VariateTag>
|
Chris@16
|
272 struct feature_of<tag::immediate_weighted_mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
273 : feature_of<tag::immediate_mean_of_variates<VariateType, VariateTag> >
|
Chris@16
|
274 {
|
Chris@16
|
275 };
|
Chris@16
|
276
|
Chris@16
|
277 ////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
278 //// droppable_accumulator<mean_impl>
|
Chris@16
|
279 //// need to specialize droppable lazy mean to cache the result at the
|
Chris@16
|
280 //// point the accumulator is dropped.
|
Chris@16
|
281 ///// INTERNAL ONLY
|
Chris@16
|
282 /////
|
Chris@16
|
283 //template<typename Sample, typename SumFeature>
|
Chris@16
|
284 //struct droppable_accumulator<impl::mean_impl<Sample, SumFeature> >
|
Chris@16
|
285 // : droppable_accumulator_base<
|
Chris@16
|
286 // with_cached_result<impl::mean_impl<Sample, SumFeature> >
|
Chris@16
|
287 // >
|
Chris@16
|
288 //{
|
Chris@16
|
289 // template<typename Args>
|
Chris@16
|
290 // droppable_accumulator(Args const &args)
|
Chris@16
|
291 // : droppable_accumulator::base(args)
|
Chris@16
|
292 // {
|
Chris@16
|
293 // }
|
Chris@16
|
294 //};
|
Chris@16
|
295
|
Chris@16
|
296 }} // namespace boost::accumulators
|
Chris@16
|
297
|
Chris@16
|
298 #endif
|