---

colvarproxy_system.h

Go to the documentation of this file.

// -*- c++ -*-

// This file is part of the Collective Variables module (Colvars).
// The original version of Colvars and its updates are located at:
// https://github.com/Colvars/colvars
// Please update all Colvars source files before making any changes.
// If you wish to distribute your changes, please submit them to the
// Colvars repository at GitHub.

#ifndef COLVARPROXY_SYSTEM_H
#define COLVARPROXY_SYSTEM_H


class colvarproxy_system {

public:

  colvarproxy_system();

  virtual ~colvarproxy_system();

  std::string units;

  virtual int set_unit_system(std::string const &units, bool check_only);

  inline cvm::real angstrom_to_internal(cvm::real l) const
  {
    return l * angstrom_value_;
  }

  inline cvm::real internal_to_angstrom(cvm::real l) const
  {
    return l / angstrom_value_;
  }

  inline cvm::real boltzmann() const
  {
    return boltzmann_;
  }

  inline cvm::real target_temperature() const
  {
    return target_temperature_;
  }

  virtual int set_target_temperature(cvm::real T);

  virtual cvm::real dt();

  virtual cvm::real rand_gaussian(void);

  virtual void add_energy(cvm::real energy);

  virtual cvm::rvector position_distance(cvm::atom_pos const &pos1,
                                         cvm::atom_pos const &pos2) const;

  void update_pbc_lattice();

  void reset_pbc_lattice();

  virtual void request_total_force(bool yesno);

  virtual bool total_forces_enabled() const;

  virtual bool total_forces_same_step() const;

  virtual int get_molid(int &molid);

  virtual int get_alch_lambda(cvm::real* lambda);

  void set_alch_lambda(cvm::real lambda);

  virtual int send_alch_lambda();

  virtual int get_dE_dlambda(cvm::real* dE_dlambda);

  virtual int apply_force_dE_dlambda(cvm::real* force);

  virtual int get_d2E_dlambda2(cvm::real* d2E_dlambda2);

  cvm::real indirect_lambda_biasing_force;

  virtual cvm::real get_accelMD_factor() const {
    cvm::error("Error: accessing the reweighting factor of accelerated MD  "
               "is not yet implemented in the MD engine.\n",
               COLVARS_NOT_IMPLEMENTED);
    return 1.0;
  }
  virtual bool accelMD_enabled() const {
    return false;
  }

protected:

  cvm::real cached_alch_lambda;

  bool cached_alch_lambda_changed;

  cvm::real boltzmann_;

  cvm::real target_temperature_;

  cvm::real angstrom_value_;

  cvm::real kcal_mol_value_;

  bool total_force_requested;

  enum Boundaries_type {
    boundaries_non_periodic,
    boundaries_pbc_ortho,
    boundaries_pbc_triclinic,
    boundaries_unsupported
  };

  Boundaries_type boundaries_type;

  cvm::rvector unit_cell_x, unit_cell_y, unit_cell_z;

  cvm::rvector reciprocal_cell_x, reciprocal_cell_y, reciprocal_cell_z;
};

#endif

---