ThermostatConfiguration.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_CONFIGURATIONS_THERMOSTATCONFIGURATION_H_
#define _MOLECULARDYNAMICS_COUPLING_CONFIGURATIONS_THERMOSTATCONFIGURATION_H_
 
#include "coupling/paralleltopology/ParallelTopologyFactory.h"
#include "tarch/configuration/Configuration.h"
#include "tarch/configuration/ParseConfiguration.h"
#include "tarch/la/Vector.h"
#include <iostream>
 
namespace coupling {
namespace configurations {
class ThermostatConfiguration;
}
} // namespace coupling

class coupling::configurations::ThermostatConfiguration : public tarch::configuration::Configuration {
public:
  enum ThermostatRegion { onlyOutestLayer, outerLayers, all, nowhere, none };
  ThermostatConfiguration() : _type{none}, _isValid(true) {}
 
  virtual ~ThermostatConfiguration() {}
 
  void parseSubtag(tinyxml2::XMLElement* node) {
    std::string value;
    tarch::configuration::ParseConfiguration::readStringMandatory(value, node, "type");
    if (value == "onlyOutestLayer") {
      _type = onlyOutestLayer;
    } else if (value == "outerLayers") {
      _type = outerLayers;
      tarch::configuration::ParseConfiguration::readIntMandatory(_cells, node, "number-layers");
      if (_cells < 0) {
        std::cout << "ERROR coupling::ThermostatConfiguration: Wrong number of cells to use!" << std::endl;
        _isValid = false;
        exit(EXIT_FAILURE);
      }
    } else if (value == "all") {
      _type = all;
    } else if (value == "nowhere") {
      _type = nowhere;
    } else {
      std::cout << "ERROR coupling::ThermostatConfiguration: Wrong type!" << std::endl;
      _isValid = false;
      exit(EXIT_FAILURE);
    }
  }

  std::string getTag() const { return "thermostat"; }

  bool isValid() const { return _isValid; }
 
  ThermostatRegion getThermostatRegionType() const { return _type; }
 
  unsigned int getCells2Use() const { return _cells; }
 
protected:
  ThermostatConfiguration(ThermostatRegion type) : _type{type}, _isValid(true) {}
 
private:
  ThermostatRegion _type;
  int _cells{0};
  bool _isValid;
};
 
#endif // _MOLECULARDYNAMICS_COUPLING_CONFIGURATIONS_THERMOSTATCONFIGURATION_H_