Chris@16
|
1 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
2 // weighted_sum.hpp
|
Chris@16
|
3 //
|
Chris@16
|
4 // Copyright 2006 Eric Niebler, Olivier Gygi. 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_WEIGHTED_SUM_HPP_EAN_28_10_2005
|
Chris@16
|
9 #define BOOST_ACCUMULATORS_STATISTICS_WEIGHTED_SUM_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/parameters/weight.hpp>
|
Chris@16
|
17 #include <boost/accumulators/framework/accumulators/external_accumulator.hpp>
|
Chris@16
|
18 #include <boost/accumulators/framework/depends_on.hpp>
|
Chris@16
|
19 #include <boost/accumulators/statistics_fwd.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 // weighted_sum_impl
|
Chris@16
|
28 template<typename Sample, typename Weight, typename Tag>
|
Chris@16
|
29 struct weighted_sum_impl
|
Chris@16
|
30 : accumulator_base
|
Chris@16
|
31 {
|
Chris@16
|
32 typedef typename numeric::functional::multiplies<Sample, Weight>::result_type weighted_sample;
|
Chris@16
|
33
|
Chris@16
|
34 // for boost::result_of
|
Chris@16
|
35 typedef weighted_sample result_type;
|
Chris@16
|
36
|
Chris@16
|
37 template<typename Args>
|
Chris@16
|
38 weighted_sum_impl(Args const &args)
|
Chris@16
|
39 : weighted_sum_(
|
Chris@16
|
40 args[parameter::keyword<Tag>::get() | Sample()]
|
Chris@16
|
41 * numeric::one<Weight>::value
|
Chris@16
|
42 )
|
Chris@16
|
43 {
|
Chris@16
|
44 }
|
Chris@16
|
45
|
Chris@16
|
46 template<typename Args>
|
Chris@16
|
47 void operator ()(Args const &args)
|
Chris@16
|
48 {
|
Chris@16
|
49 // what about overflow?
|
Chris@16
|
50 this->weighted_sum_ += args[parameter::keyword<Tag>::get()] * args[weight];
|
Chris@16
|
51 }
|
Chris@16
|
52
|
Chris@16
|
53 result_type result(dont_care) const
|
Chris@16
|
54 {
|
Chris@16
|
55 return this->weighted_sum_;
|
Chris@16
|
56 }
|
Chris@16
|
57
|
Chris@16
|
58 private:
|
Chris@16
|
59
|
Chris@16
|
60 weighted_sample weighted_sum_;
|
Chris@16
|
61 };
|
Chris@16
|
62
|
Chris@16
|
63 } // namespace impl
|
Chris@16
|
64
|
Chris@16
|
65 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
66 // tag::weighted_sum
|
Chris@16
|
67 //
|
Chris@16
|
68 namespace tag
|
Chris@16
|
69 {
|
Chris@16
|
70 struct weighted_sum
|
Chris@16
|
71 : depends_on<>
|
Chris@16
|
72 {
|
Chris@16
|
73 /// INTERNAL ONLY
|
Chris@16
|
74 ///
|
Chris@16
|
75 typedef accumulators::impl::weighted_sum_impl<mpl::_1, mpl::_2, tag::sample> impl;
|
Chris@16
|
76 };
|
Chris@16
|
77
|
Chris@16
|
78 template<typename VariateType, typename VariateTag>
|
Chris@16
|
79 struct weighted_sum_of_variates
|
Chris@16
|
80 : depends_on<>
|
Chris@16
|
81 {
|
Chris@16
|
82 /// INTERNAL ONLY
|
Chris@16
|
83 ///
|
Chris@16
|
84 typedef accumulators::impl::weighted_sum_impl<VariateType, mpl::_2, VariateTag> impl;
|
Chris@16
|
85 };
|
Chris@16
|
86
|
Chris@16
|
87 struct abstract_weighted_sum_of_variates
|
Chris@16
|
88 : depends_on<>
|
Chris@16
|
89 {
|
Chris@16
|
90 };
|
Chris@16
|
91 }
|
Chris@16
|
92
|
Chris@16
|
93 ///////////////////////////////////////////////////////////////////////////////
|
Chris@16
|
94 // extract::weighted_sum
|
Chris@16
|
95 //
|
Chris@16
|
96 namespace extract
|
Chris@16
|
97 {
|
Chris@16
|
98 extractor<tag::weighted_sum> const weighted_sum = {};
|
Chris@16
|
99 extractor<tag::abstract_weighted_sum_of_variates> const weighted_sum_of_variates = {};
|
Chris@16
|
100
|
Chris@16
|
101 BOOST_ACCUMULATORS_IGNORE_GLOBAL(weighted_sum)
|
Chris@16
|
102 BOOST_ACCUMULATORS_IGNORE_GLOBAL(weighted_sum_of_variates)
|
Chris@16
|
103 }
|
Chris@16
|
104
|
Chris@16
|
105 using extract::weighted_sum;
|
Chris@16
|
106 using extract::weighted_sum_of_variates;
|
Chris@16
|
107
|
Chris@16
|
108 template<typename VariateType, typename VariateTag>
|
Chris@16
|
109 struct feature_of<tag::weighted_sum_of_variates<VariateType, VariateTag> >
|
Chris@16
|
110 : feature_of<tag::abstract_weighted_sum_of_variates>
|
Chris@16
|
111 {
|
Chris@16
|
112 };
|
Chris@16
|
113
|
Chris@16
|
114 }} // namespace boost::accumulators
|
Chris@16
|
115
|
Chris@16
|
116 #endif
|