BasinHoppingSaddleSearch.h

Go to the documentation of this file.

/*
** This file is part of eOn.
**
** SPDX-License-Identifier: BSD-3-Clause
**
** Copyright (c) 2010--present, eOn Development Team
** All rights reserved.
**
** Repo:
** https://github.com/TheochemUI/eOn
*/
#pragma once
#include "Eigen.h"
#include "EonLogger.h"
#include "Matter.h"
#include "MinModeSaddleSearch.h"
#include "SaddleSearchMethod.h"
#include <vector>
 
namespace eonc {
 
class BasinHoppingSaddleSearch : public SaddleSearchMethod {
public:
  BasinHoppingSaddleSearch(std::shared_ptr<Matter> reactant,
                           std::shared_ptr<Matter> displacement,
                           std::shared_ptr<Potential> potPassed,
                           const Parameters &parametersPassed)
      : SaddleSearchMethod(potPassed, parametersPassed),
        reactant{std::make_shared<Matter>(*reactant)},
        saddle{displacement} {
    eigenvector.resize(reactant->numberOfAtoms(), 3);
    eigenvector.setZero();
  }
  ~BasinHoppingSaddleSearch() = default;
 
  int run(void);
  double getEigenvalue();
  AtomMatrix getEigenvector();
  std::string_view describeStatus(int status) const override {
    return MinModeSaddleSearch::statusMessage(status);
  }
  int getStatus() const override { return status; }
 
  double eigenvalue{0.0};
  AtomMatrix eigenvector;
 
  std::shared_ptr<Matter> reactant;
  std::shared_ptr<Matter> saddle;
  std::shared_ptr<Matter> product;
 
  int status{0};
 
private:
  eonc::log::Scoped log;
};
 
} // namespace eonc
 
using eonc::BasinHoppingSaddleSearch;