source: src/atom_trajectoryparticle.hpp@ 986ed3

Action_Thermostats Add_AtomRandomPerturbation Add_FitFragmentPartialChargesAction Add_RotateAroundBondAction Add_SelectAtomByNameAction Added_ParseSaveFragmentResults AddingActions_SaveParseParticleParameters Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_ParticleName_to_Atom Adding_StructOpt_integration_tests AtomFragments Automaking_mpqc_open AutomationFragmentation_failures Candidate_v1.5.4 Candidate_v1.6.0 Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator CombiningParticlePotentialParsing Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError Fix_BoundInBox_CenterInBox_MoleculeActions Fix_ChargeSampling_PBC Fix_ChronosMutex Fix_FitPartialCharges Fix_FitPotential_needs_atomicnumbers Fix_ForceAnnealing Fix_IndependentFragmentGrids Fix_ParseParticles Fix_ParseParticles_split_forward_backward_Actions Fix_PopActions Fix_QtFragmentList_sorted_selection Fix_Restrictedkeyset_FragmentMolecule Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns Fix_fitting_potentials Fixes ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion FragmentAction_writes_AtomFragments FragmentMolecule_checks_bonddegrees GeometryObjects Gui_Fixes Gui_displays_atomic_force_velocity ImplicitCharges IndependentFragmentGrids IndependentFragmentGrids_IndividualZeroInstances IndependentFragmentGrids_IntegrationTest IndependentFragmentGrids_Sole_NN_Calculation JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool JobMarket_unresolvable_hostname_fix MoreRobust_FragmentAutomation ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PdbParser_setsAtomName PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks Rewrite_FitPartialCharges RotateToPrincipalAxisSystem_UndoRedo SaturateAtoms_findBestMatching SaturateAtoms_singleDegree StoppableMakroAction Subpackage_CodePatterns Subpackage_JobMarket Subpackage_LinearAlgebra Subpackage_levmar Subpackage_mpqc_open Subpackage_vmg Switchable_LogView ThirdParty_MPQC_rebuilt_buildsystem TrajectoryDependenant_MaxOrder TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps TremoloParser_setsAtomName Ubuntu_1604_changes stable
Last change on this file since 986ed3 was ef7d30, checked in by Frederik Heber <heber@…>, 15 years ago

MEMFIX: ForceMatrix in VelocityVerletIntegration had a stupid offset.
  • Property mode set to 100644
File size: 2.2 KB
Line 
1/*
2 * atom_trajectoryparticle.hpp
3 *
4 * Created on: Oct 19, 2009
5 * Author: heber
6 */
7
8#ifndef ATOM_TRAJECTORYPARTICLE_HPP_
9#define ATOM_TRAJECTORYPARTICLE_HPP_
10
11using namespace std;
12
13/*********************************************** includes ***********************************/
14
15// include config.h
16#ifdef HAVE_CONFIG_H
17#include <config.h>
18#endif
19
20#include <fstream>
21
22#include <gsl/gsl_randist.h>
23
24#include "atom_atominfo.hpp"
25#include "atom_particleinfo.hpp"
26#include "atom_trajectoryparticleinfo.hpp"
27
28/****************************************** forward declarations *****************************/
29
30class atom;
31class config;
32class ForceMatrix;
33class Vector;
34
35/********************************************** declarations *******************************/
36
37
38class TrajectoryParticle : public TrajectoryParticleInfo, public virtual AtomInfo, public virtual ParticleInfo {
39public:
40
41 TrajectoryParticle();
42 virtual ~TrajectoryParticle();
43
44 // constraint potential and dynamics stuff
45 void AddKineticToTemperature(double *temperature, int step) const;
46 void EvaluateConstrainedForce(int startstep, int endstep, atom **PermutationMap, ForceMatrix *Force) const;
47 void CorrectVelocity(double *ActualTemp, int Step, Vector *CoGVelocity);
48
49 // trajectory stuff
50 void ResizeTrajectory(int MaxSteps);
51 void CopyStepOnStep(int dest, int src);
52 void VelocityVerletUpdate(int MDSteps, config *configuration, ForceMatrix *Force, const size_t offset);
53 void SumUpKineticEnergy( int Step, double *TotalMass, Vector *TotalVelocity ) const;
54
55 // thermostats
56 void Thermostat_Woodcock(double ScaleTempFactor, int Step, double *ekin);
57 void Thermostat_Gaussian_init(int Step, double *G, double *E);
58 void Thermostat_Gaussian_least_constraint(int Step, double G_over_E, double *ekin, config *configuration);
59 void Thermostat_Langevin(int Step, gsl_rng * r, double *ekin, config *configuration);
60 void Thermostat_Berendsen(int Step, double ScaleTempFactor, double *ekin, config *configuration);
61 void Thermostat_NoseHoover_init(int Step, double *delta_alpha);
62 void Thermostat_NoseHoover_scale(int Step, double *ekin, config *configuration);
63
64private:
65
66};
67
68
69#endif /* ATOM_TRAJECTORYPARTICLE_HPP_ */
Note: See TracBrowser for help on using the repository browser.