ARTnSaddleSearch.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 "EonLogger.h"
#include "Matter.h"
#include "MinModeSaddleSearch.h"
#include "Parameters.h"
#include "SaddleSearchMethod.h"
#include <limits>
 
namespace eonc {
 
class IARTnResource;

class ARTnSaddleSearch : public SaddleSearchMethod {
public:
  static constexpr int STATUS_GOOD = 0;
  static constexpr int STATUS_BAD_MAX_ITERATIONS =
      MinModeSaddleSearch::STATUS_BAD_MAX_ITERATIONS;
  static constexpr int STATUS_BAD_ARTN_ERROR = 22;
 
  ARTnSaddleSearch(std::shared_ptr<Matter> matterPassed,
                   std::shared_ptr<Potential> potPassed, AtomMatrix modeInitial,
                   const Parameters &paramsPassed);
  ~ARTnSaddleSearch() override;

  int run() override;
  int run(IARTnResource &res);
  double getEigenvalue() override;
  AtomMatrix getEigenvector() override;
  std::string_view describeStatus(int status) const override;
  int getStatus() const override { return status; }
  int getIterationCount() const override { return iteration; }
  int getForceCalls() const override { return forcecalls; }
 
private:
  std::shared_ptr<Matter> matter;
  double eigenvalue{std::numeric_limits<double>::quiet_NaN()};
  AtomMatrix eigenvector, mode;
  int status{0};
  int iteration{0};
  int forcecalls{0};
  eonc::log::Scoped log;
};
 
} // namespace eonc