ImproperDihedralTerm.h

Go to the documentation of this file.

#ifndef MMIMPROPERDIHEDRALTERM_H
#define MMIMPROPERDIHEDRALTERM_H

#include "../Topology/ImproperDihedralType.h"
#include "ImproperDihedral.h"
#include "InteractionTermBase.h"
#include <Utils/Math/AutomaticDifferentiation/Second3D.h>
#include <Utils/Typenames.h>

namespace Scine {

namespace Utils {
class AtomicSecondDerivativeCollection;
class Tensor33;
} // namespace Utils

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

  ImproperDihedralTerm(AtomIndex firstAtom, AtomIndex secondAtom, AtomIndex centralAtom, AtomIndex fourthAtom,
                       const ImproperDihedral& improperDihedral, const ImproperDihedralType& typeOfImproperDihedral);
  ~ImproperDihedralTerm();

  double evaluateImproperDihedralTerm(const Utils::PositionCollection& positions,
                                      Utils::AtomicSecondDerivativeCollection& derivatives) const;

  static double getTheta(const Eigen::Vector3d& A, const Eigen::Vector3d& B, const Eigen::Vector3d& G);

  ImproperDihedralType getTypeOfImproperDihedral() const;

  int getFirstAtom() const;

  int getSecondAtom() const;

  int getCentralAtom() const;

  int getFourthAtom() const;

 private:
  Utils::AutomaticDifferentiation::Second3D threeDimDer(const Utils::AutomaticDifferentiation::Second1D& energy,
                                                        const Utils::AutomaticDifferentiation::Second3D& alpha) const;
  void setSecondDerivative(Utils::AutomaticDifferentiation::Second3D& h, const Utils::Tensor33& tensor, double derFirst,
                           double derSecond) const;

  AtomIndex firstAtom_, secondAtom_, centralAtom_, fourthAtom_;
  ImproperDihedral improperDihedral_;
  ImproperDihedralType typeOfImproperDihedral_;
};

} // namespace MolecularMechanics
} // namespace Scine

#endif // MMIMPROPERDIHEDRALTERM_H