conserved_3d.cpp

#include "conserved_3d.hpp"

using std::size_t;

Conserved3D::Conserved3D(void):
  mass(0), momentum(), energy(0), tracers() {}

Conserved3D::Conserved3D(double mass_i,
                         const Vector3D& momentum_i,
                         double energy_i):
  mass(mass_i), momentum(momentum_i), energy(energy_i), tracers() {}

Conserved3D::Conserved3D(double mass_i,
                         const Vector3D& momentum_i,
                         double energy_i,
                         const vector<double>& tracers_i):
  mass(mass_i), momentum(momentum_i), 
  energy(energy_i), tracers(tracers_i) {}

namespace {
  vector<double> operator*(double s, const vector<double>& v)
  {
    vector<double> res(v.size());
    for(size_t i=0;i<v.size();++i)
      res[i] = s*v[i];
    return res;
  }

  vector<double> operator/(const vector<double>& v,double s)
  {
    vector<double> res(v.size());
    for(size_t i=0;i<v.size();++i)
      res[i] = v[i]/s;
    return res;
  }
}

Conserved3D& Conserved3D::operator-=(const Conserved3D& diff)
{
  mass -= diff.mass;
  momentum -= diff.momentum;
  energy -= diff.energy;
  for(size_t i=0;i<tracers.size();++i)
    tracers[i] -= diff.tracers[i];
  return *this;
}

Conserved3D& Conserved3D::operator+=(const Conserved3D& diff)
{
  mass += diff.mass;
  momentum += diff.momentum;
  energy += diff.energy;
  for(size_t i=0;i<tracers.size();++i)
    tracers[i] += diff.tracers[i];
  return *this;
}

Conserved3D operator*(double s, const Conserved3D& c)
{
  return Conserved3D(s*c.mass,
                     s*c.momentum,
                     s*c.energy,
                     s*c.tracers);
}

Conserved3D operator/(const Conserved3D& c, double s)
{
  return Conserved3D(c.mass/s,
                     c.momentum/s,
                     c.energy/s,
                     c.tracers/s);
}