computational_cell_2d.hpp

#ifndef COMPUTATIONAL_CELL_HPP
#define COMPUTATIONAL_CELL_HPP 1

/*#include <boost/container/small_vector.hpp>
typedef boost::container::small_vector<double,0> tvector;
typedef boost::container::small_vector<bool,0> svector;*/
#include <vector>
typedef std::vector<double> tvector;
typedef std::vector<bool> svector;
#include <string>
#include "../../tessellation/geometry.hpp"
#ifdef RICH_MPI
#include "../../misc/serializable.hpp"
#endif // RICH_MPI

class ComputationalCell
#ifdef RICH_MPI
        : public Serializable
#endif // RICH_MPI

{
public:

        double density;

        double pressure;

        Vector2D velocity;

        tvector tracers;

        svector stickers;

        ComputationalCell(ComputationalCell const& other);
        ComputationalCell(void);

        ComputationalCell& operator+=(ComputationalCell const& other);
        ComputationalCell& operator-=(ComputationalCell const& other);

        ComputationalCell& operator*=(double s);

        ComputationalCell& operator=(ComputationalCell const& other);

#ifdef RICH_MPI
  size_t getChunkSize(void) const;

  vector<double> serialize(void) const;

  void unserialize
  (const vector<double>& data);
#endif // RICH_MPI
};

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

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

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

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

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

void ComputationalCellAddMult(ComputationalCell &res, ComputationalCell const& other, double scalar);

void ReplaceComputationalCell(ComputationalCell &res, ComputationalCell const& other);

void ReplaceComputationalCellDiff(ComputationalCell &res, ComputationalCell const& first, ComputationalCell const& second);

class Slope
#ifdef RICH_MPI
        : public Serializable
#endif // RICH_MPI
{
public:
        ComputationalCell xderivative;

        ComputationalCell yderivative;

        Slope(ComputationalCell const& x, ComputationalCell const& y);
        Slope(void);
#ifdef RICH_MPI
        size_t getChunkSize(void) const;

        vector<double> serialize(void) const;

        void unserialize(const vector<double>& data);
#endif//RICH_MPI
};
class TracerStickerNames
{
public:
        std::vector<std::string> tracer_names;
        std::vector<std::string> sticker_names;

        TracerStickerNames(void);
        TracerStickerNames(TracerStickerNames const& other);
        TracerStickerNames
        (const std::vector<std::string>& tracers,
         const std::vector<std::string>& stickers);

  TracerStickerNames& operator=(const TracerStickerNames& other);

        ~TracerStickerNames(void);
};

#endif // COMPUTATIONAL_CELL_HPP