source: src/Actions/AnalysisAction/PointCorrelationAction.cpp@ d02e07

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

Split PairCorrelation into Pair, Point and SurfaceCorrelation.

  • also converted to new createDialog().
  • Property mode set to 100644
File size: 4.2 KB
Line 
1/*
2 * PointCorrelationAction.cpp
3 *
4 * Created on: May 9, 2010
5 * Author: heber
6 */
7
8#include "Helpers/MemDebug.hpp"
9
10#include "Actions/AnalysisAction/PointCorrelationAction.hpp"
11#include "analysis_correlation.hpp"
12#include "boundary.hpp"
13#include "linkedcell.hpp"
14#include "verbose.hpp"
15#include "log.hpp"
16#include "element.hpp"
17#include "molecule.hpp"
18#include "periodentafel.hpp"
19#include "vector.hpp"
20#include "World.hpp"
21
22#include <iostream>
23#include <string>
24
25using namespace std;
26
27#include "UIElements/UIFactory.hpp"
28#include "UIElements/Dialog.hpp"
29#include "Actions/MapOfActions.hpp"
30
31const char AnalysisPointCorrelationAction::NAME[] = "point-correlation";
32
33AnalysisPointCorrelationAction::AnalysisPointCorrelationAction() :
34 Action(NAME)
35{}
36
37AnalysisPointCorrelationAction::~AnalysisPointCorrelationAction()
38{}
39
40Dialog* AnalysisPointCorrelationAction::createDialog() {
41 Dialog *dialog = UIFactory::getInstance().makeDialog();
42 int ranges[3] = {1, 1, 1};
43 double BinEnd = 0.;
44 double BinStart = 0.;
45 double BinWidth = 0.;
46 molecule *Boundary = NULL;
47 string outputname;
48 string binoutputname;
49 bool periodic;
50 ofstream output;
51 ofstream binoutput;
52 std::vector< element *> elements;
53 string type;
54 Vector Point;
55 BinPairMap *binmap = NULL;
56
57 dialog->queryVector("position", &Point, false, MapOfActions::getInstance().getDescription("position"));
58 dialog->queryElement("elements", &elements, MapOfActions::getInstance().getDescription("elements"));
59 dialog->queryDouble("bin-start", &BinStart, MapOfActions::getInstance().getDescription("bin-start"));
60 dialog->queryDouble("bin-width", &BinWidth, MapOfActions::getInstance().getDescription("bin-width"));
61 dialog->queryDouble("bin-end", &BinEnd, MapOfActions::getInstance().getDescription("bin-end"));
62 dialog->queryString("output-file", &outputname, MapOfActions::getInstance().getDescription("output-file"));
63 dialog->queryString("bin-output-file", &binoutputname, MapOfActions::getInstance().getDescription("bin-output-file"));
64 dialog->queryBoolean("periodic", &periodic, MapOfActions::getInstance().getDescription("periodic"));
65
66 return dialog;
67}
68
69Action::state_ptr AnalysisPointCorrelationAction::performCall() {
70 int ranges[3] = {1, 1, 1};
71 double BinEnd = 0.;
72 double BinStart = 0.;
73 double BinWidth = 0.;
74 molecule *Boundary = NULL;
75 string outputname;
76 string binoutputname;
77 bool periodic;
78 ofstream output;
79 ofstream binoutput;
80 std::vector< element *> elements;
81 string type;
82 Vector Point;
83 BinPairMap *binmap = NULL;
84 MoleculeListClass *molecules = World::getInstance().getMolecules();
85
86 // obtain information
87 MapOfActions::getInstance().queryCurrentValue("position", Point);
88 MapOfActions::getInstance().queryCurrentValue("elements", elements);
89 MapOfActions::getInstance().queryCurrentValue("bin-start", BinStart);
90 MapOfActions::getInstance().queryCurrentValue("bin-width", BinWidth);
91 MapOfActions::getInstance().queryCurrentValue("bin-end", BinEnd);
92 MapOfActions::getInstance().queryCurrentValue("output-file", outputname);
93 MapOfActions::getInstance().queryCurrentValue("bin-output-file", binoutputname);
94 MapOfActions::getInstance().queryCurrentValue("periodic", periodic);
95
96 // execute action
97 output.open(outputname.c_str());
98 binoutput.open(binoutputname.c_str());
99 cout << "Point to correlate to is " << Point << endl;
100 CorrelationToPointMap *correlationmap = NULL;
101 if (periodic)
102 correlationmap = PeriodicCorrelationToPoint(molecules, elements, &Point, ranges);
103 else
104 correlationmap = CorrelationToPoint(molecules, elements, &Point);
105 OutputCorrelationToPoint(&output, correlationmap);
106 binmap = BinData( correlationmap, BinWidth, BinStart, BinEnd );
107 OutputCorrelation ( &binoutput, binmap );
108 delete(binmap);
109 delete(correlationmap);
110 output.close();
111 binoutput.close();
112 return Action::success;
113}
114
115Action::state_ptr AnalysisPointCorrelationAction::performUndo(Action::state_ptr _state) {
116 return Action::success;
117}
118
119Action::state_ptr AnalysisPointCorrelationAction::performRedo(Action::state_ptr _state){
120 return Action::success;
121}
122
123bool AnalysisPointCorrelationAction::canUndo() {
124 return true;
125}
126
127bool AnalysisPointCorrelationAction::shouldUndo() {
128 return true;
129}
130
131const string AnalysisPointCorrelationAction::getName() {
132 return NAME;
133}
Note: See TracBrowser for help on using the repository browser.