source: src/Actions/FragmentationAction/FragmentationAction.cpp@ d0fbec

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 d0fbec was 35b698, checked in by Frederik Heber <heber@…>, 15 years ago

BIG CHANGE: config::load and config::save in ParseCommandLineOptions() and main() replaced with FormatParser replacements.

Fragmentation:

  • FIX: MoleculeFillWithMoleculeAction: filler atoms have to be removed before the system can be stored to file.
  • FIX: PcpParser::load() - has to put the molecule also into World's MoleculeListClass (otherwise the name cannot be set right after loading)
  • new Libparser.a
  • all sources from PARSER subdir are compiled into libparser such that only ParserUnitTest is recompiled.

Testfixes:

  • testsuite-fragmentation - changes to due to different -f calling syntax.
  • most of the xyz files had to be replaced due to a single whitespace at the end of each entry: Domain/6, Simple_configuration/2, Simple_configuration/3, Simple_configuration/4, Simple_configuration/5, Simple_configuration/8
  • in many cases were the number orbitals (and thus MaxMinStopStep) wrong: Filling/1, Simple_configuration/4, Simple_configuration/5

Signed-off-by: Frederik Heber <heber@…>

  • Property mode set to 100644
File size: 3.2 KB
Line 
1/*
2 * FragmentationAction.cpp
3 *
4 * Created on: May 9, 2010
5 * Author: heber
6 */
7
8#include "Helpers/MemDebug.hpp"
9
10#include "Actions/FragmentationAction/FragmentationAction.hpp"
11#include "atom.hpp"
12#include "bondgraph.hpp"
13#include "config.hpp"
14#include "log.hpp"
15#include "molecule.hpp"
16#include "Descriptors/MoleculeDescriptor.hpp"
17#include "stackclass.hpp"
18#include "World.hpp"
19
20#include <iostream>
21#include <string>
22
23using namespace std;
24
25#include "UIElements/UIFactory.hpp"
26#include "UIElements/Dialog.hpp"
27#include "Actions/MapOfActions.hpp"
28
29const char FragmentationFragmentationAction::NAME[] = "fragment-mol";
30
31FragmentationFragmentationAction::FragmentationFragmentationAction() :
32 Action(NAME)
33{}
34
35FragmentationFragmentationAction::~FragmentationFragmentationAction()
36{}
37
38Action::state_ptr FragmentationFragmentationAction::performCall() {
39 Dialog *dialog = UIFactory::getInstance().makeDialog();
40 clock_t start,end;
41 molecule *mol = NULL;
42 double distance = -1.;
43 int order = 0;
44 std::string path;
45 config *configuration = World::getInstance().getConfig();
46 int ExitFlag = 0;
47
48 cout << "pre-dialog"<< endl;
49 dialog->queryString(NAME, &path, MapOfActions::getInstance().getDescription(NAME));
50 dialog->queryMolecule("molecule-by-id", &mol, MapOfActions::getInstance().getDescription("molecule-by-id"));
51 dialog->queryDouble("distance", &distance, MapOfActions::getInstance().getDescription("distance"));
52 dialog->queryInt("order", &order, MapOfActions::getInstance().getDescription("order"));
53
54 if(dialog->display()) {
55 cout << "POST-dialog"<< endl;
56 ASSERT(mol != NULL, "No molecule has been picked for fragmentation.");
57 DoLog(0) && (Log() << Verbose(0) << "Fragmenting molecule with bond distance " << distance << " angstroem, order of " << order << "." << endl);
58 DoLog(0) && (Log() << Verbose(0) << "Creating connection matrix..." << endl);
59 start = clock();
60 mol->CreateAdjacencyList(distance, configuration->GetIsAngstroem(), &BondGraph::CovalentMinMaxDistance, NULL);
61 DoLog(0) && (Log() << Verbose(0) << "Fragmenting molecule with current connection matrix ..." << endl);
62 if (mol->hasBondStructure()) {
63 ExitFlag = mol->FragmentMolecule(order, path);
64 }
65 World::getInstance().setExitFlag(ExitFlag);
66 end = clock();
67 DoLog(0) && (Log() << Verbose(0) << "Clocks for this operation: " << (end-start) << ", time: " << ((double)(end-start)/CLOCKS_PER_SEC) << "s." << endl);
68 delete dialog;
69 return Action::success;
70 } else {
71 delete dialog;
72 return Action::failure;
73 }
74}
75
76Action::state_ptr FragmentationFragmentationAction::performUndo(Action::state_ptr _state) {
77// ParserLoadXyzState *state = assert_cast<ParserLoadXyzState*>(_state.get());
78
79 return Action::failure;
80// string newName = state->mol->getName();
81// state->mol->setName(state->lastName);
82//
83// return Action::state_ptr(new ParserLoadXyzState(state->mol,newName));
84}
85
86Action::state_ptr FragmentationFragmentationAction::performRedo(Action::state_ptr _state){
87 return Action::failure;
88}
89
90bool FragmentationFragmentationAction::canUndo() {
91 return false;
92}
93
94bool FragmentationFragmentationAction::shouldUndo() {
95 return false;
96}
97
98const string FragmentationFragmentationAction::getName() {
99 return NAME;
100}
Note: See TracBrowser for help on using the repository browser.