1 files changed, 93 insertions, 0 deletions
diff --git a/src/core/latticeStatistics.h b/src/core/latticeStatistics.h
new file mode 100644
index 0000000..bb32bcd
--- /dev/null
+++ b/src/core/latticeStatistics.h
@@ -0,0 +1,93 @@
+/*  This file is part of the OpenLB library
+ *
+ *  Copyright (C) 2006, 2007 Jonas Latt
+ *  E-mail contact: info@openlb.net
+ *  The most recent release of OpenLB can be downloaded at
+ *  <http://www.openlb.net/>
+ *
+ *  This program 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 2
+ *  of the License, or (at your option) any later version.
+ *
+ *  This program 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 this program; if not, write to the Free
+ *  Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
+ *  Boston, MA  02110-1301, USA.
+*/
+
+/** \file
+ * Interface for post-processing steps -- header file.
+ */
+#ifndef LATTICE_STATISTICS_H
+#define LATTICE_STATISTICS_H
+
+#include <vector>
+#include "io/ostreamManager.h"
+
+namespace olb {
+
+/////////////////// Statistics Postprocessing ////////////////////////////
+
+template<typename T>
+class LatticeStatistics {
+public:
+  enum { avRho=0, avEnergy=1 } AverageT;
+  enum { maxU=0 } MaxT;
+public:
+  LatticeStatistics();
+  ~LatticeStatistics();
+  void reset();
+  void reset(T average_rho_, T average_energy_, T maxU_, size_t numCells_);
+
+  int subscribeAverage();
+  int subscribeSum();
+  int subscribeMin();
+  int subscribeMax();
+
+  void incrementStats(T rho, T uSqr);
+  void gatherAverage(int whichAverage, T value);
+  void gatherSum(int whichSum, T value);
+  void gatherMin(int whichMin, T value);
+  void gatherMax(int whichMax, T value);
+  void incrementStats();
+  T getAverageRho() const;
+  T getAverageEnergy() const;
+  T getMaxU() const;
+  size_t const& getNumCells() const;
+
+  T getAverage(int whichAverage) const;
+  T getSum(int whichSum) const;
+  T getMin(int whichMin) const;
+  T getMax(int whichMax) const;
+
+  std::vector<T>& getAverageVect();
+  std::vector<T>& getSumVect();
+  std::vector<T>& getMinVect();
+  std::vector<T>& getMaxVect();
+
+  void incrementTime();
+  void resetTime(size_t value=0);
+  size_t getTime() const;
+  void print(int iterationStep, T physicalTime=-1) const;
+  void initialize();
+private:
+  mutable OstreamManager clout;
+  // variables for internal computations
+  std::vector<T> tmpAv, tmpSum, tmpMin, tmpMax;
+  size_t tmpNumCells;
+  // variables containing the public result
+  std::vector<T> averageVect, sumVect, minVect, maxVect;
+  size_t numCells;
+  size_t latticeTime;
+  bool firstCall;
+};
+
+}  // namespace olb
+
+#endif