ElectrostaticTerm.h

Go to the documentation of this file.

#ifndef MOLECULARMECHANICS_ELECTROSTATICTERM_H
#define MOLECULARMECHANICS_ELECTROSTATICTERM_H

#include "Electrostatic.h"
#include "InteractionTermBase.h"
#include <Utils/Typenames.h>
#include <memory>
#include <vector>

namespace Scine {

namespace Utils {
class FullSecondDerivativeCollection;
} // namespace Utils

namespace MolecularMechanics {
class ElectrostaticTerm : public InteractionTermBase {
 public:
  using AtomIndex = int;

  ElectrostaticTerm(AtomIndex firstAtom, AtomIndex secondAtom, const Electrostatic& electrostatic,
                    std::shared_ptr<double> cutoffRadius);

  double evaluateElectrostaticTerm(const Utils::PositionCollection& positions, Utils::FullSecondDerivativeCollection& derivatives,
                                   const std::vector<double>& atomicCharges, const double& scalingFactorForEachCharge) const;

  int getFirstAtom() const;
  int getSecondAtom() const;

 private:
  AtomIndex firstAtom_, secondAtom_;
  Electrostatic electrostatic_;
  // Pointer so that it is updated whenever the settings are updated.
  // Reference does not work because then the class loses its assignment operator.
  std::shared_ptr<double> cutoffRadius_;
};

} // namespace MolecularMechanics
} // namespace Scine

#endif // MOLECULARMECHANICS_ELECTROSTATICTERM_H