ConservativeForce.hpp

Go to the documentation of this file.

#ifndef CONSFORCE_HPP
#define CONSFORCE_HPP 1

#include "../SourceTerm.hpp"

class Acceleration
{
public:
        virtual Vector2D operator()
                (const Tessellation& tess,
                        const vector<ComputationalCell>& cells,
                        const vector<Extensive>& fluxes,
                        const double time,
                        const int point,
                        TracerStickerNames const& tracerstickernames) const = 0;

        virtual ~Acceleration(void);
};
class ConservativeForce : public SourceTerm
{
public:
        explicit ConservativeForce(const Acceleration& acc, bool mass_flux = false);

        ~ConservativeForce(void);

        vector<Extensive> operator()
                (const Tessellation& tess,
                        const PhysicalGeometry& pg,
                        const CacheData& cd,
                        const vector<ComputationalCell>& cells,
                        const vector<Extensive>& fluxes,
                        const vector<Vector2D>& point_velocities,
                        const double t,
                        TracerStickerNames const& tracerstickernames) const;

private:
        const Acceleration& acc_;
        const bool mass_flux_;

        ConservativeForce(const ConservativeForce& origin);
        ConservativeForce& operator=(const ConservativeForce& origin);
};

#endif // CONSFORCE_HPP