improved_center_gravity.hpp

Go to the documentation of this file.

#ifndef CENTERGRAVITY_HPP
#define CENTERGRAVITY_HPP 1

#include "ConservativeForce.hpp"

class ImprovedCenterGravity: public Acceleration
{
public:
  ImprovedCenterGravity(double M,double Rmin,const Vector2D& center,double Rsoft=0);

  Vector2D operator()
  (const Tessellation& tess,
   const vector<ComputationalCell>& cells,
   const vector<Extensive>& fluxes,
   const double time,
   const int point,
   TracerStickerNames const& tracerstickernames) const;

  Vector2D const& get_center(void) const;

private:
  const double M_;
  const double Rmin_;
  const Vector2D center_;
  const double Rsoft_;
};

#endif // CENTERGRAVITY_HPP