hydrodynamic_variables.hpp

Go to the documentation of this file.

#ifndef HYDRODYNAMIC_VARIABLES_HPP
#define HYDRODYNAMIC_VARIABLES_HPP 1

#include "../../tessellation/geometry.hpp"

class Conserved
{
public:

  Conserved(void);

  Conserved(double mass,
            Vector2D const& momentum,
            double energy);

  double Mass;

  Vector2D Momentum;

  double Energy;

  Conserved& operator+=(Conserved const& f);

  Conserved& operator-=(Conserved const& c);

  Conserved& operator=(Conserved const&other);

  Conserved(const Conserved& other);
};

class Primitive
{
public:

  enum {NUMBER_PRIMITIVE_VARIABLES=6};

  Primitive(void);

  Primitive(double density_i,
            double pressure_i,
            Vector2D const& velocity_i,
            double energy_i,
            double sound_speed_i);

  double Density;

  double Pressure;

  Vector2D Velocity;

  double Energy;

  double SoundSpeed;

  int GetVarNo(void) const;

  Primitive& operator+=(Primitive const& p);

  double& operator[](int index);

  double operator[](int index) const;

  Primitive& operator=(Primitive const&other);

  Primitive(const Primitive& other);
};

bool primitive_has_nan(Primitive const& p);

Primitive operator+(Primitive const& p1, Primitive const& p2);

Primitive operator-(Primitive const& p1, Primitive const& p2);

Primitive operator/(Primitive const& p, double s);

Primitive operator*(Primitive const& p, double s);

Primitive operator*(double s, Primitive const& p);

Conserved operator*(double d, Conserved const& c);

Conserved operator+(Conserved const& c1, Conserved const& c2);

Conserved operator-(Conserved const& c1, Conserved const& c2);

Conserved operator/(Conserved const& c, double d);

double TotalEnergyDensity(Primitive const& p);

Conserved Primitive2Conserved(Primitive const& p);

Conserved Primitive2Conserved(Primitive const& p, double volume);

Conserved Primitive2Flux(Primitive const& p, Vector2D const& n);

#endif // HYDRODYNAMIC_VARIABLES_HPP