source: src/Fragmentation/Exporters/HydrogenPool.hpp@ 101d2d

Action_Thermostats Add_AtomRandomPerturbation Add_SelectAtomByNameAction Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_StructOpt_integration_tests AutomationFragmentation_failures Candidate_v1.6.1 ChangeBugEmailaddress ChemicalSpaceEvaluator Docu_Python_wait EmpiricalPotential_contain_HomologyGraph_documentation Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph Fix_Verbose_Codepatterns ForceAnnealing_oldresults ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion Gui_displays_atomic_force_velocity IndependentFragmentGrids_IntegrationTest JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks StoppableMakroAction TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps
Last change on this file since 101d2d was 7294dc, checked in by Frederik Heber <frederik.heber@…>, 8 years ago

FIX: Fixing the use of the trajectories that were changed to maps.

  • as the atom's trajectories are now stored inside a map, i.e. not every step necessarily needs to be present, there will still some places in the code that were meant for the old vector storing.
  • Atom::UpdateStep() updates the current step, not the future step.
  • AtomInfo::append/removeTrajectoryStep() now no longer ASSERT that the step is not the current one. Due to the trajectory in map concept (with current it simply sought by lower_bound()) this is no longer necessary.
  • HydrogenPool no longer adds all trajectory steps till current time.
  • World has commented out code for also copying the current atom positions in case a new time step is added, but this is not used as so far we do not recognize whether it is a new time step ... and also it is not really necessary.
  • TESTS: Removed XFAILs from tests of previous commit. They are all working again.
  • Property mode set to 100644
File size: 2.2 KB
Line 
1/*
2 * HydrogenPool.hpp
3 *
4 * Created on: Mar 3, 2013
5 * Author: heber
6 */
7
8#ifndef HYDROGENPOOL_HPP_
9#define HYDROGENPOOL_HPP_
10
11// include config.h
12#ifdef HAVE_CONFIG_H
13#include <config.h>
14#endif
15
16#include <deque>
17#include <map>
18
19#include "types.hpp"
20
21class atom;
22
23/** This class contains a pool of hydrogen meant for saturation.
24 *
25 * These hydrogen atoms do not trigger any updates, they are mostly
26 * invisible to the system.
27 *
28 * If the pool is lacking hydrogens, it requests more from the World.
29 *
30 * The idea is that we always create a single fragment, saturate it with
31 * hydrogen atoms from the pool, then store it to file or convert to a
32 * job. Directly afterwards, the fragment is free'd and the hydrogen may
33 * go back into the pool.
34 */
35class HydrogenPool
36{
37 //!> grant unit test access to private parts
38 friend class HydrogenPoolTest;
39public:
40 HydrogenPool();
41 ~HydrogenPool();
42
43 /** Delivers a hydrogen atom for use.
44 *
45 * \return point to the hydrogen atom
46 */
47 atom * leaseHydrogen();
48
49 /** Releases a used hydrogen again.
50 *
51 * \param _atom hydrogen atom to release
52 */
53 void releaseHydrogen(atom * _atom);
54
55 /** Releases a used hydrogen again.
56 *
57 * \param _atom hydrogen atom to release
58 */
59 void releaseHydrogen(atomId_t _atom);
60
61 //!> enum to define the atomic number of atoms contained in the pool
62 enum { HYDROGEN = 1 };
63
64private:
65 /** Requests a hydrogen atom from the World and places it into the queue.
66 *
67 */
68 void requestHydrogenIntoPool();
69
70 /** Helper function that removes all present hydrogen atoms again.
71 *
72 */
73 void cleanup();
74
75private:
76 //!> typedef for the deque of available hydrogens.
77 typedef std::deque<atom *> HydrogenQueue_t;
78 //!> typedef for the map of all hydrogens currently in use in a fragment.
79 typedef std::map<atomId_t, atom *> HydrogenInUse_t;
80
81 //!> enum to define a warning constant for the limit of hydrogens
82 enum { WARNINGTHRESHOLD = 100 };
83
84 //!> counter for the total number of hydrogens requested so far.
85 int HydrogenCount;
86 //!> Queue of available hydrogen atoms
87 HydrogenQueue_t HydrogenQueue;
88 //!> Map for looking up ids of hydrogen in use for reinsertion into queue
89 HydrogenInUse_t HydrogenInUse;
90};
91
92#endif /* HYDROGENPOOL_HPP_ */
Note: See TracBrowser for help on using the repository browser.