source: src/Actions/GraphAction/DestroyAdjacencyAction.cpp@ b5b01e

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 b5b01e was b5b01e, checked in by Frederik Heber <heber@…>, 11 years ago

ActionState extracted into own header file, rename Action::state_ptr -> ActionState::ptr.

  • Property mode set to 100644
File size: 5.1 KB
Line 
1/*
2 * Project: MoleCuilder
3 * Description: creates and alters molecular systems
4 * Copyright (C) 2013 University of Bonn. All rights reserved.
5 * Copyright (C) 2013 Frederik Heber. All rights reserved.
6 *
7 *
8 * This file is part of MoleCuilder.
9 *
10 * MoleCuilder is free software: you can redistribute it and/or modify
11 * it under the terms of the GNU General Public License as published by
12 * the Free Software Foundation, either version 2 of the License, or
13 * (at your option) any later version.
14 *
15 * MoleCuilder is distributed in the hope that it will be useful,
16 * but WITHOUT ANY WARRANTY; without even the implied warranty of
17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18 * GNU General Public License for more details.
19 *
20 * You should have received a copy of the GNU General Public License
21 * along with MoleCuilder. If not, see <http://www.gnu.org/licenses/>.
22 */
23
24/*
25 * DestroyAdjacencyAction.cpp
26 *
27 * Created on: May 1, 2013
28 * Author: heber
29 */
30
31// include config.h
32#ifdef HAVE_CONFIG_H
33#include <config.h>
34#endif
35
36#include "CodePatterns/MemDebug.hpp"
37
38#include "CodePatterns/Assert.hpp"
39#include "CodePatterns/Log.hpp"
40
41#include "Atom/atom.hpp"
42#include "Descriptors/AtomIdDescriptor.hpp"
43#include "Descriptors/AtomSelectionDescriptor.hpp"
44#include "Graph/BondGraph.hpp"
45#include "World.hpp"
46#include "WorldTime.hpp"
47
48#include <iostream>
49#include <list>
50#include <string>
51
52#include "Actions/GraphAction/DestroyAdjacencyAction.hpp"
53
54using namespace MoleCuilder;
55
56// Storing undo state
57struct BondInfo_t {
58 atomId_t leftatom;
59 atomId_t rightatom;
60 size_t degree;
61};
62
63// and construct the stuff
64#include "DestroyAdjacencyAction.def"
65#include "Action_impl_pre.hpp"
66/** =========== define the function ====================== */
67ActionState::ptr GraphDestroyAdjacencyAction::performCall() {
68 BondGraph *BG = World::getInstance().getBondGraph();
69 ASSERT(BG != NULL, "GraphDestroyAdjacencyAction: BondGraph is NULL.");
70
71 // check precondition
72 World::AtomComposite Set = World::getInstance().getAllAtoms(AtomsBySelection());
73 if (Set.empty()) {
74 ELOG(2, "No atoms selected.");
75 return Action::failure;
76 }
77
78 // count all present bonds
79 std::vector<BondInfo_t> bonds;
80 {
81 size_t count_bonds = 0;
82 for (World::AtomComposite::const_iterator iter = Set.begin();
83 iter != Set.end();
84 ++iter) {
85 const atom * const Walker = *iter;
86 count_bonds += Walker->getListOfBonds().size();
87 }
88 bonds.reserve(count_bonds/2);
89 }
90 // Storing undo info
91 for (World::AtomComposite::const_iterator iter = Set.begin();
92 iter != Set.end();
93 ++iter) {
94 const atom * const Walker = *iter;
95 const BondList& ListOfBonds = Walker->getListOfBonds();
96 for (BondList::const_iterator bonditer = ListOfBonds.begin();
97 bonditer != ListOfBonds.end();
98 ++bonditer) {
99 const bond::ptr &CurrentBond = *bonditer;
100 // if both atoms are in selected set, we check ids otherwise not
101 if (((!World::getInstance().isSelected(CurrentBond->leftatom))
102 || (!World::getInstance().isSelected(CurrentBond->rightatom)))
103 || (CurrentBond->leftatom->getId() < CurrentBond->rightatom->getId())) {
104 BondInfo_t BondInfo;
105 BondInfo.leftatom = CurrentBond->leftatom->getId();
106 BondInfo.rightatom = CurrentBond->rightatom->getId();
107 BondInfo.degree = CurrentBond->getDegree();
108 bonds.push_back(BondInfo);
109 }
110 }
111 }
112 GraphDestroyAdjacencyState *UndoState = new GraphDestroyAdjacencyState(bonds, params);
113
114 LOG(0, "STATUS: Clearing the Adjacency.");
115 BG->cleanAdjacencyList(Set);
116
117 return ActionState::ptr(UndoState);
118}
119
120ActionState::ptr GraphDestroyAdjacencyAction::performUndo(ActionState::ptr _state) {
121 GraphDestroyAdjacencyState *state = assert_cast<GraphDestroyAdjacencyState*>(_state.get());
122
123 const size_t CurrentTime = WorldTime::getTime();
124 std::vector<BondInfo_t> &bonds = state->bonds;
125 for(std::vector<BondInfo_t>::const_iterator iter = bonds.begin();
126 iter != bonds.end(); ++iter) {
127 atom * const Walker = World::getInstance().getAtom(AtomById(iter->leftatom));
128 ASSERT( Walker != NULL,
129 "GraphDestroyAdjacencyAction::performUndo() - "+toString(iter->leftatom)+" missing.");
130 atom * const OtherWalker = World::getInstance().getAtom(AtomById(iter->rightatom));
131 ASSERT( OtherWalker != NULL,
132 "GraphDestroyAdjacencyAction::performUndo() - "+toString(iter->rightatom)+" missing.");
133 bond::ptr CurrentBond = Walker->addBond(CurrentTime, OtherWalker);
134 CurrentBond->setDegree(iter->degree);
135 }
136 return ActionState::ptr(_state);
137}
138
139ActionState::ptr GraphDestroyAdjacencyAction::performRedo(ActionState::ptr _state){
140 BondGraph *BG = World::getInstance().getBondGraph();
141 ASSERT(BG != NULL, "GraphDestroyAdjacencyAction: BondGraph is NULL.");
142
143 World::AtomComposite Set = World::getInstance().getAllAtoms(AtomsBySelection());
144 BG->cleanAdjacencyList(Set);
145
146 return ActionState::ptr(_state);
147}
148
149bool GraphDestroyAdjacencyAction::canUndo() {
150 return true;
151}
152
153bool GraphDestroyAdjacencyAction::shouldUndo() {
154 return true;
155}
156/** =========== end of function ====================== */
Note: See TracBrowser for help on using the repository browser.