- Timestamp:
- May 20, 2014, 9:14:56 AM (11 years ago)
- Branches:
- 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
- Children:
- 74459a
- Parents:
- 975b83
- git-author:
- Frederik Heber <heber@…> (08/26/13 22:56:36)
- git-committer:
- Frederik Heber <heber@…> (05/20/14 09:14:56)
- Location:
- src
- Files:
-
- 1 added
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/ActionQueue.cpp
r975b83 r415ddd 42 42 #include "CodePatterns/Singleton_impl.hpp" 43 43 44 #include <boost/date_time/posix_time/posix_time.hpp> 45 #include <boost/version.hpp> 44 46 #include <string> 45 47 #include <sstream> … … 56 58 AR(new ActionRegistry()), 57 59 history(new ActionHistory), 58 CurrentAction(0) 60 CurrentAction(0), 61 run_thread(boost::bind(&ActionQueue::run, this)), 62 run_thread_isIdle(true) 59 63 {} 60 64 61 65 ActionQueue::~ActionQueue() 62 66 { 67 stop(); 68 63 69 // free all actions contained in actionqueue 64 70 for (ActionQueue_t::iterator iter = actionqueue.begin(); !actionqueue.empty(); iter = actionqueue.begin()) { … … 80 86 Action *newaction = _action->clone(state); 81 87 newaction->prepare(state); 88 mtx_queue.lock(); 82 89 actionqueue.push_back( newaction ); 83 try { 84 newaction->call(); 85 } catch(ActionFailureException &e) { 86 std::cerr << "Action " << *boost::get_error_info<ActionNameString>(e) << " has failed." << std::endl; 87 World::getInstance().setExitFlag(5); 88 } 90 { 91 boost::lock_guard<boost::mutex> lock(mtx_idle); 92 run_thread_isIdle = (CurrentAction == actionqueue.size()); 93 } 94 mtx_queue.unlock(); 89 95 } 90 96 91 97 void ActionQueue::insertAction(Action *_action, enum Action::QueryOptions state) 92 98 { 93 queueAction(_action, state); 99 Action *newaction = _action->clone(state); 100 newaction->prepare(state); 101 mtx_queue.lock(); 102 tempqueue.push_back( newaction ); 103 { 104 boost::lock_guard<boost::mutex> lock(mtx_idle); 105 run_thread_isIdle = !((CurrentAction != actionqueue.size()) || !tempqueue.empty()); 106 } 107 mtx_queue.unlock(); 108 } 109 110 void ActionQueue::run() 111 { 112 bool Interrupted = false; 113 do { 114 // sleep for some time and wait for queue to fill up again 115 try { 116 #if BOOST_VERSION < 105000 117 run_thread.sleep(boost::get_system_time() + boost::posix_time::milliseconds(100)); 118 #else 119 run_thread.sleep_for(100); 120 #endif 121 } catch(boost::thread_interrupted &e) { 122 LOG(2, "INFO: ActionQueue has received stop signal."); 123 Interrupted = true; 124 } 125 // LOG(1, "DEBUG: Start of ActionQueue's run() loop."); 126 // call all currently present Actions 127 mtx_queue.lock(); 128 insertTempQueue(); 129 bool status = (CurrentAction != actionqueue.size()); 130 mtx_queue.unlock(); 131 while (status) { 132 // boost::this_thread::disable_interruption di; 133 LOG(0, "Calling Action " << actionqueue[CurrentAction]->getName() << " ... "); 134 try { 135 actionqueue[CurrentAction]->call(); 136 } catch(ActionFailureException &e) { 137 std::cerr << "Action " << *boost::get_error_info<ActionNameString>(e) << " has failed." << std::endl; 138 World::getInstance().setExitFlag(5); 139 } 140 // access actionqueue, hence using mutex 141 mtx_queue.lock(); 142 // step on to next action and check for end 143 CurrentAction++; 144 // insert new actions (before [CurrentAction]) if they have been spawned 145 // we must have an extra vector for this, as we cannot change actionqueue 146 // while an action instance is "in-use" 147 insertTempQueue(); 148 status = (CurrentAction != actionqueue.size()); 149 mtx_queue.unlock(); 150 } 151 { 152 boost::lock_guard<boost::mutex> lock(mtx_idle); 153 run_thread_isIdle = !((CurrentAction != actionqueue.size()) || !tempqueue.empty()); 154 } 155 cond_idle.notify_one(); 156 // LOG(1, "DEBUG: End of ActionQueue's run() loop."); 157 } while (!Interrupted); 158 } 159 160 void ActionQueue::insertTempQueue() 161 { 162 if (!tempqueue.empty()) { 163 ActionQueue_t::iterator InsertionIter = actionqueue.begin(); 164 std::advance(InsertionIter, CurrentAction); 165 actionqueue.insert( InsertionIter, tempqueue.begin(), tempqueue.end() ); 166 tempqueue.clear(); 167 } 168 } 169 170 void ActionQueue::wait() 171 { 172 boost::unique_lock<boost::mutex> lock(mtx_idle); 173 while(!run_thread_isIdle) 174 { 175 cond_idle.wait(lock); 176 } 177 } 178 179 void ActionQueue::stop() 180 { 181 // notify actionqueue thread that we wish to terminate 182 run_thread.interrupt(); 183 // wait till it ends 184 run_thread.join(); 94 185 } 95 186 -
src/Actions/ActionQueue.hpp
r975b83 r415ddd 16 16 #include "CodePatterns/Singleton.hpp" 17 17 18 #include <boost/thread.hpp> 18 19 #include <vector> 19 20 20 21 #include "Actions/Action.hpp" 21 22 #include "Actions/ActionState.hpp" 23 24 void stopQueue(); 25 void waitQueue(); 22 26 23 27 namespace MoleCuilder { … … 113 117 void redoLast(); 114 118 115 /** Getter for the last executed action. 116 * 117 * \note this is necessary for Reactions to actually have a chance of getting 118 * the calculated value as Action's are always cloned. 119 * 120 * \return const ref to last action 121 */ 122 const Action &getLastAction() const { 123 return *(actionqueue.back()); 124 } 119 boost::thread &getRunThread() 120 { return run_thread; } 125 121 126 122 private: … … 140 136 void clear(); 141 137 138 /** Runs the ActionQueue. 139 * 140 */ 141 void run(); 142 143 friend void ::stopQueue(); 144 145 /** Stops the internal thread. 146 * 147 */ 148 void stop(); 149 150 friend void ::waitQueue(); 151 152 /** Wait till all currently queued actions are processed. 153 * 154 */ 155 void wait(); 156 142 157 /** Insert an action after CurrentAction. 143 158 * … … 148 163 void insertAction(Action *_action, enum Action::QueryOptions state); 149 164 165 /** Moves all action from tempqueue into real queue. 166 * 167 */ 168 void insertTempQueue(); 169 150 170 private: 151 171 /** Private cstor for ActionQueue. … … 172 192 //!> point to current action in actionqueue 173 193 size_t CurrentAction; 194 195 //!> internal temporary actionqueue of actions used by insertAction() 196 ActionQueue_t tempqueue; 197 198 //!> internal thread to call Actions 199 boost::thread run_thread; 200 201 //!> internal mutex to synchronize access to queue 202 boost::mutex mtx_queue; 203 204 //!> conditional variable notifying when run_thread is idling 205 boost::condition_variable cond_idle; 206 207 //!> flag indicating whether run_thread is idle or not 208 bool run_thread_isIdle; 209 210 //!> internal mutex to synchronize access to run_thread_isIdle 211 boost::mutex mtx_idle; 174 212 }; 175 213 -
src/Python/Makefile.am
r975b83 r415ddd 8 8 Python/getSelectedMolarMass.cpp \ 9 9 Python/PythonScripting.cpp \ 10 Python/reinit.cpp 10 Python/reinit.cpp \ 11 Python/wait.cpp 11 12 12 13 PYTHONHEADER = \ -
src/Python/PythonScripting_impl.hpp
r975b83 r415ddd 85 85 "Reinitializes the internal state of the python module as if it had been freshly imported,saves all input files beforehand." 86 86 ); 87 boost::python::def( 88 "wait", 89 MoleCuilder::detail::module_wait, 90 "Waits until all currently queued actions have been processed." 91 ); 87 92 boost::python::def< MoleCuilder::detail::doubleVec() >( 88 93 "getBoundingBox", -
src/UIElements/UIFactory.hpp
r975b83 r415ddd 77 77 static void makeUserInterface(std::string type); 78 78 static void registerFactory(factoryDescription *factoryDesc); 79 static bool isFactoryPresent() { return !factories.empty(); } 79 80 virtual std::string getUIName(){ return "none"; } 80 81 protected: -
src/builder.cpp
r975b83 r415ddd 49 49 doUI(); 50 50 51 waitQueue(); 52 stopQueue(); 51 53 if (UIFactory::getInstance().getUIName() == "CommandLine") 52 54 return saveAll(); -
src/builder_init.hpp
r975b83 r415ddd 20 20 void doUI(); 21 21 22 23 22 #endif /* BUILDER_INIT_HPP_ */ -
src/cleanUp.cpp
r975b83 r415ddd 92 92 void purgeStaticInstances() 93 93 { 94 // make sure that ActionQueue is already purged! 95 94 96 Chronos::purgeInstance(); 95 97 RandomNumberDistributionFactory::purgeInstance(); … … 105 107 UIFactory::purgeInstance(); 106 108 CommandLineParser::purgeInstance(); 107 MoleCuilder::ActionQueue::purgeInstance();108 MoleCuilder::OptionRegistry::purgeInstance();109 109 logger::purgeInstance(); 110 110 errorLogger::purgeInstance(); … … 112 112 } 113 113 114 void stopAndPurgeQueue() 115 { 116 // adding a wait here such that ActionQueue's run_thread is done any may stop 117 waitQueue(); 118 MoleCuilder::ActionQueue::purgeInstance(); 119 MoleCuilder::OptionRegistry::purgeInstance(); 120 } 121 122 void printTimings(const MoleCuilder::ActionQueue::ActionTokens_t &tokens); 123 114 124 /** Cleans all singleton instances in an orderly fashion. 115 125 * C++ does not guarantee any specific sequence of removal of single instances 116 126 * which have static/global variables. Some singletons depend on others hence we 117 * a certain a specific ordering here, which is is used via the atexit() hook.127 * ascertain a specific ordering here, which is is used via the atexit() hook. 118 128 */ 119 129 void cleanUp() 120 130 { 121 // give timings per Action 122 printTimings(); 131 // stop queue 132 const MoleCuilder::ActionQueue::ActionTokens_t tokens = MoleCuilder::ActionQueue::getInstance().getListOfActions(); 133 stopAndPurgeQueue(); 134 // give timings per Action, requires still present Chronos 135 printTimings(tokens); 136 123 137 // purge static instances from memory 124 138 purgeStaticInstances(); … … 138 152 /** We give a list of all times per action and a total time. 139 153 * 140 */ 141 void printTimings() 142 { 143 const MoleCuilder::ActionQueue::ActionTokens_t tokens = MoleCuilder::ActionQueue::getInstance().getListOfActions(); 154 * We require a list of Action tokens such that we do not access the ActionQueue. 155 * Timings are printed when Actions have all been performed, hence the queue should 156 * already be purged and we cannot access it anymore. 157 * 158 * \param tokens 159 */ 160 void printTimings(const MoleCuilder::ActionQueue::ActionTokens_t &tokens) 161 { 162 // const MoleCuilder::ActionQueue::ActionTokens_t tokens = MoleCuilder::ActionQueue::getInstance().getListOfActions(); 144 163 const Chronos &Chron = Chronos::getInstance(); 145 164 if (!DoLog(2)) { … … 198 217 return (ExitFlag == 1 ? 0 : ExitFlag); 199 218 } 219 220 /** Stops the queue such that all Actions are done. 221 * 222 */ 223 void stopQueue() 224 { 225 MoleCuilder::ActionQueue::getInstance().stop(); 226 } 227 228 /** Waits for the queue to idle. 229 * 230 */ 231 void waitQueue() 232 { 233 MoleCuilder::ActionQueue::getInstance().wait(); 234 } -
src/cleanUp.hpp
r975b83 r415ddd 17 17 void cleanUp(); 18 18 void dumpMemory(); 19 void printTimings();20 19 void purgeStaticInstances(); 21 20 int saveAll(); 21 void stopQueue(); 22 void waitQueue(); 22 23 23 24 #endif /* CLEANUP_HPP_ */ -
src/documentation/Makefile.am
r975b83 r415ddd 13 13 ${top_srcdir}/src/documentation/constructs/constructs.dox \ 14 14 ${top_srcdir}/src/documentation/constructs/descriptors.dox \ 15 ${top_srcdir}/src/documentation/constructs/filling.dox \ 15 16 ${top_srcdir}/src/documentation/constructs/fragmentation.dox \ 16 17 ${top_srcdir}/src/documentation/constructs/linearalgebra.dox \ … … 18 19 ${top_srcdir}/src/documentation/constructs/molecules.dox \ 19 20 ${top_srcdir}/src/documentation/constructs/observers_observables.dox \ 21 ${top_srcdir}/src/documentation/constructs/parameters.dox \ 20 22 ${top_srcdir}/src/documentation/constructs/parsers.dox \ 23 ${top_srcdir}/src/documentation/constructs/potentials.dox \ 24 ${top_srcdir}/src/documentation/constructs/qt-gui.dox \ 25 ${top_srcdir}/src/documentation/constructs/queries.dox \ 21 26 ${top_srcdir}/src/documentation/constructs/randomnumbers.dox \ 22 27 ${top_srcdir}/src/documentation/constructs/serialization.dox \ 28 ${top_srcdir}/src/documentation/constructs/shaperegistry.dox \ 23 29 ${top_srcdir}/src/documentation/constructs/shapes.dox \ 24 30 ${top_srcdir}/src/documentation/constructs/tesselation.dox \ 31 ${top_srcdir}/src/documentation/constructs/validators.dox \ 32 ${top_srcdir}/src/documentation/constructs/values.dox \ 25 33 ${top_srcdir}/src/documentation/constructs/world.dox \ 26 34 ${top_srcdir}/src/documentation/howtos/action.dox \ -
src/documentation/userinterfaces/python.dox
r975b83 r415ddd 20 20 * script. 21 21 * 22 * This is done in \b src/ Actions/pyMoleCuilder.cpp.22 * This is done in \b src/Python/PythonScripting.cpp. 23 23 * 24 24 * There again some preprocessor magic is happening. One the one hand we … … 43 43 * mol.SelectAtomById("0") 44 44 * mol.AtomRemove() 45 * mol.wait() 46 * mol.getSelectedMolarMass() 47 * mol.wait() 45 48 * \endcode 46 49 * which loads a file \b test.xyz into the (internal) World, selects the first 47 * atom and removes it. 50 * atom and removes it. Notice \b mol.wait() at the end. This might be necessary 51 * as actions are executed in a different thread than the python script itself. 52 * Hence, if you require values from molecuilder you have to make sure that 53 * all your actions have been processed by this second thread. That's what 54 * wait() is good for. It waits until action queue thread is idle. Then you 55 * can be sure that molecuilder has removed all atoms, performed all selections 56 * and any value you retrieve is up-to-date. 57 * 58 * Note that there are two \b wait()s present in the example. As the Actions 59 * are executed in another thread and the above commands just tell the MoleCuilder 60 * library (the ActionQueue to be precise) to enqueue the requested action, 61 * we have to wait (in the main thread) until the actions actually have been 62 * executed before we continue (i.e. when we need the new state where the 63 * atoms have been removed) and before we \b terminate! 48 64 * 49 65 * \section userinterfaces-python-running Running a test script … … 72 88 * \code 73 89 * pyMoleCuilder.WorldInput("test.xyz") 90 * pyMoleCuilder.wait() 74 91 * \endcode 75 92 * … … 78 95 * \code 79 96 * pyMoleCuilder.SelectionAllMolecules() 97 * pyMoleCuilder.wait() 80 98 * \endcode 81 99 * … … 84 102 * \code 85 103 * pyMoleCuilder.SelectAtomById("0") 104 * pyMoleCuilder.wait() 86 105 * \endcode 106 * 107 * \warning Again, take note of the added wait()s that ensure the all enqueued 108 * actions also have been executed. This is especially important in scripts as 109 * otherwise your script may deadlock. That's because ActionQueue's destructor 110 * waits for the thread that executes the actions to end, and in another thread 111 * we still want to access to ActionQueue whose instance is however locked as 112 * it is about the get destroyed. 87 113 * 88 114 * \subsection userinterfaces-python-notes-cleanup Cleaning up or reset state ... … … 115 141 * gives you the docu string on WorldInputAction. 116 142 * 143 * \subsection userinterfaces-python-notes-wait Waiting for the action queue 117 144 * 118 * \date 2013-03-18 145 * Note again that actions are executed in a different thread as the python 146 * script. Hence, we require synchronization at certain intervals where you 147 * require molecuilder to be up to speed. All commands you executed such 148 * as 149 * \code 150 * import pyMoleCuilder as mol 151 * mol.WorldInput("foo.xyz") 152 * mol.wait() 153 * \endcode 154 * just queue this specific input action but not execute it right away. That's 155 * left to the other thread. Hence, you need to wait() before: 156 * -# you access mol.get...() functions as these are not actions themselves. 157 * -# you need to have files written by molecuilder to be parsed in the python 158 * script. 159 * 160 * 161 * \date 2013-09-28 119 162 * 120 163 */
Note:
See TracChangeset
for help on using the changeset viewer.