MassAccumulator.h

Go to the documentation of this file.

#pragma once
 
 
#include <vector>
 
#include "tarch/logging/Log.h"
#include "tarch/multicore/BooleanSemaphore.h"
 
 
namespace applications {
  namespace exahype2 {
    namespace euler {
      namespace sphericalaccretion {
        class MassAccumulator;
      } // namespace sphericalaccretion
    }   // namespace euler
  }     // namespace exahype2
} // namespace applications
 
 
class applications::exahype2::euler::sphericalaccretion::MassAccumulator {
private:
  static tarch::logging::Log _log;
 
  tarch::multicore::BooleanSemaphore _semaphore;
 
  const double _shellWidth;
  const double _threshold;
 
  std::vector<double> _currentMass;
 
  std::vector<double> _previousMassAccumulated;
 
  double _maxDistanceThatMakesContribution;
 
  int getBucket(double radius) const;
 
public:
  MassAccumulator(double shellWidth = 0.01, double threshold = 0.01);
  void addMass(double mass, double radius);
 
  double getMass_piecewiseConstantInterpolation(double radius) const;
  double getMass_linearInterpolation(double radius) const;
 
  void finishAccumulation();
 
  double getMaxRelevantRadius() const;
 
  double getTotalMass() const;
};