OHTSTJob.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 "eon/Eigen.h"
#include "eon/Job.h"
#include "eon/Matter.h"
#include "eon/Parameters.h"
 
namespace eonc {

class OHTSTJob : public Job {
 
public:
  OHTSTJob(std::unique_ptr<Parameters> parameters)
      : Job(std::move(parameters)) {}
  ~OHTSTJob(void) = default;
  std::vector<std::string> run(void);
 
private:
  struct PlaneAverages {
    double fn{0.0};      
    VectorXd rotNorm;    
    VectorXd rotRaw;     
    VectorXd pos;        
  };
 
  PlaneAverages samplePlane(Matter &matter, const VectorXd &gamma,
                            const VectorXd &normal);

  double reactantQRatio(Matter &matter, const VectorXd &gammaR,
                        const VectorXd &normal);

  void drawThermalVelocities(VectorXd &vel, const VectorXd *normal);

  bool symmetryReflect(const VectorXd &xR, VectorXd &x, VectorXd &v,
                       const VectorXd &xOld, const VectorXd *normal);
 
  std::vector<VectorXd> m_symDirs; 
  VectorXd m_symXR;                
 
  VectorXd m_masses3N;        
  double m_dt{0.0};           
  double m_kbt{0.0};          
  double m_andersenProb{0.0}; 
  long m_seedState{12345};    
  double uniformDraw();       
  double gaussDraw();         
  bool m_gaussHave{false};
  double m_gaussSpare{0.0};
};
 
} // namespace eonc
 
using eonc::OHTSTJob;