CouplingCell.h

// Copyright (C) 2015 Technische Universitaet Muenchen
// This file is part of the Mamico project. For conditions of distribution
// and use, please see the copyright notice in Mamico's main folder, or at
// www5.in.tum.de/mamico
#ifndef _MOLECULARDYNAMICS_COUPLING_DATA_STRUCTURES_COUPLINGCELL_H_
#define _MOLECULARDYNAMICS_COUPLING_DATA_STRUCTURES_COUPLINGCELL_H_
 
#include "tarch/la/Vector.h"
#include <cstdlib>
#include <fstream>
#include <iostream>
#include <string>
 
namespace coupling {
namespace datastructures {
 
template <unsigned int dim> class CouplingCell;
} // namespace datastructures
} // namespace coupling

template <unsigned int dim> class coupling::datastructures::CouplingCell {
public:
  CouplingCell()
      : _microscopicMass(0.0), _microscopicMomentum(0.0), _macroscopicMass(0.0), _macroscopicMomentum(0.0), _potentialEnergy(0.0), _temperature(0.0),
        _currentVelocity(0.0) {}

  virtual ~CouplingCell() {}

  void setMicroscopicMass(const double& mass) { _microscopicMass = mass; }
  const double& getMicroscopicMass() const { return _microscopicMass; }
  void setMicroscopicMomentum(const tarch::la::Vector<dim, double>& momentum) { _microscopicMomentum = momentum; }
  const tarch::la::Vector<dim, double>& getMicroscopicMomentum() const { return _microscopicMomentum; }

  void setMacroscopicMass(const double& mass) { _macroscopicMass = mass; }
  const double& getMacroscopicMass() const { return _macroscopicMass; }
  void setMacroscopicMomentum(const tarch::la::Vector<dim, double>& momentum) { _macroscopicMomentum = momentum; }
  const tarch::la::Vector<dim, double>& getMacroscopicMomentum() const { return _macroscopicMomentum; }

  const double& getPotentialEnergy() const { return _potentialEnergy; }
  void setPotentialEnergy(const double& potentialEnergy) { _potentialEnergy = potentialEnergy; }

  void addMicroscopicMass(const double& mass) { _microscopicMass += mass; }
  void addMicroscopicMomentum(const tarch::la::Vector<dim, double>& momentum) { _microscopicMomentum += momentum; }

  void addMacroscopicMass(const double& mass) { _macroscopicMass += mass; }
  void addMacroscopicMomentum(const tarch::la::Vector<dim, double>& momentum) { _macroscopicMomentum += momentum; }

  void setCurrentVelocity(const tarch::la::Vector<dim, double>& velocity) { _currentVelocity = velocity; }
  const tarch::la::Vector<dim, double>& getCurrentVelocity() const { return _currentVelocity; }

  void setTemperature(const double& temperature) { _temperature = temperature; }
  const double& getTemperature() const { return _temperature; }

  double _microscopicMass;
  tarch::la::Vector<dim, double> _microscopicMomentum;

  double _macroscopicMass;
  tarch::la::Vector<dim, double> _macroscopicMomentum;

  double _potentialEnergy;

  double _temperature;

  tarch::la::Vector<dim, double> _currentVelocity;
};
 
#endif // _MOLECULARDYNAMICS_COUPLING_DATA_STRUCTURES_COUPLINGCELL_H_