combinedStatistics.h

00001 /* This file is part of the Palabos library.
00002  *
00003  * Copyright (C) 2011 FlowKit Sarl
00004  * Avenue de Chailly 23
00005  * 1012 Lausanne, Switzerland
00006  * E-mail contact: contact@flowkit.com
00007  *
00008  * The most recent release of Palabos can be downloaded at 
00009  * <http://www.palabos.org/>
00010  *
00011  * The library Palabos is free software: you can redistribute it and/or
00012  * modify it under the terms of the GNU Affero General Public License as
00013  * published by the Free Software Foundation, either version 3 of the
00014  * License, or (at your option) any later version.
00015  *
00016  * The library is distributed in the hope that it will be useful,
00017  * but WITHOUT ANY WARRANTY; without even the implied warranty of
00018  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00019  * GNU Affero General Public License for more details.
00020  *
00021  * You should have received a copy of the GNU Affero General Public License
00022  * along with this program.  If not, see <http://www.gnu.org/licenses/>.
00023 */
00024 
00028 #ifndef COMBINED_STATISTICS_H
00029 #define COMBINED_STATISTICS_H
00030 
00031 #include "core/globalDefs.h"
00032 #include "core/blockStatistics.h"
00033 #include <vector>
00034 
00035 namespace plb {
00036 
00037 class CombinedStatistics {
00038 public:
00039     virtual ~CombinedStatistics();
00040     virtual CombinedStatistics* clone() const =0;
00041     void combine (
00042             std::vector<BlockStatistics const*>& individualStatistics,
00043             BlockStatistics& result ) const;
00044 protected:
00045     virtual void reduceStatistics (
00046             std::vector<double>& averageObservables,
00047             std::vector<double>& sumWeights,
00048             std::vector<double>& sumObservables,
00049             std::vector<double>& maxObservables,
00050             std::vector<plint>& intSumObservables ) const =0;
00051 private:
00052     void computeLocalAverage (
00053             std::vector<BlockStatistics const*> const& individualStatistics,
00054             std::vector<double>& averageObservables,
00055             std::vector<double>& sumWeights ) const;
00056     void computeLocalSum (
00057             std::vector<BlockStatistics const*> const& individualStatistics,
00058             std::vector<double>& sumObservables ) const;
00059     void computeLocalMax (
00060             std::vector<BlockStatistics const*> const& individualStatistics,
00061             std::vector<double>& maxObservables ) const;
00062     void computeLocalIntSum (
00063             std::vector<BlockStatistics const*> const& individualStatistics,
00064             std::vector<plint>& intSumObservables ) const;
00065 };
00066 
00067 class SerialCombinedStatistics : public CombinedStatistics {
00068 public:
00069     virtual SerialCombinedStatistics* clone() const;
00070 protected:
00071     virtual void reduceStatistics (
00072             std::vector<double>& averageObservables,
00073             std::vector<double>& sumWeights,
00074             std::vector<double>& sumObservables,
00075             std::vector<double>& maxObservables,
00076             std::vector<plint>& intSumObservables ) const;
00077 };
00078 
00079 }  // namespace plb
00080 
00081 #endif