Changeset 9ee38b for src/Actions/AnalysisAction/PointCorrelationAction.cpp
- Timestamp:
- Aug 28, 2010, 12:52:58 AM (15 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, Candidate_v1.7.0, 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:
- b4fa106
- Parents:
- 0b2ce9
- git-author:
- Frederik Heber <heber@…> (08/25/10 17:20:37)
- git-committer:
- Frederik Heber <heber@…> (08/28/10 00:52:58)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/AnalysisAction/PointCorrelationAction.cpp
r0b2ce9 r9ee38b 20 20 #include "Helpers/MemDebug.hpp" 21 21 22 #include "Actions/AnalysisAction/PointCorrelationAction.hpp"23 #include "Actions/ActionRegistry.hpp"24 22 #include "analysis_correlation.hpp" 25 23 #include "boundary.hpp" … … 39 37 using namespace std; 40 38 41 #include "UIElements/UIFactory.hpp" 42 #include "UIElements/Dialog.hpp" 43 #include "Actions/ValueStorage.hpp" 39 #include "Actions/AnalysisAction/PointCorrelationAction.hpp" 44 40 45 const char AnalysisPointCorrelationAction::NAME[] = "point-correlation";46 41 47 AnalysisPointCorrelationAction::AnalysisPointCorrelationAction() : 48 Action(NAME) 49 {} 42 // and construct the stuff 43 #include "PointCorrelationAction.def" 44 #include "Action_impl_pre.hpp" 50 45 51 AnalysisPointCorrelationAction::~AnalysisPointCorrelationAction() 52 {} 53 54 void AnalysisPointCorrelation(std::vector< element *> &elements, Vector &position, double BinStart, double BinWidth, double BinEnd, string &outputname, string &binoutputname, bool periodic) { 55 ValueStorage::getInstance().setCurrentValue("elements", elements); 56 ValueStorage::getInstance().setCurrentValue("position", position); 57 ValueStorage::getInstance().setCurrentValue("bin-start", BinStart); 58 ValueStorage::getInstance().setCurrentValue("bin-width", BinWidth); 59 ValueStorage::getInstance().setCurrentValue("bin-end", BinEnd); 60 ValueStorage::getInstance().setCurrentValue("output-file", outputname); 61 ValueStorage::getInstance().setCurrentValue("bin-output-file", binoutputname); 62 ValueStorage::getInstance().setCurrentValue("periodic", periodic); 63 ActionRegistry::getInstance().getActionByName(AnalysisPointCorrelationAction::NAME)->call(Action::NonInteractive); 64 }; 65 66 void AnalysisPointCorrelationAction::getParametersfromValueStorage() 67 {}; 68 69 Dialog* AnalysisPointCorrelationAction::fillDialog(Dialog *dialog) { 70 ASSERT(dialog,"No Dialog given when filling action dialog"); 71 72 dialog->queryVector("position", false, ValueStorage::getInstance().getDescription("position")); 73 dialog->queryElements("elements", ValueStorage::getInstance().getDescription("elements")); 74 dialog->queryDouble("bin-start", ValueStorage::getInstance().getDescription("bin-start")); 75 dialog->queryDouble("bin-width", ValueStorage::getInstance().getDescription("bin-width")); 76 dialog->queryDouble("bin-end", ValueStorage::getInstance().getDescription("bin-end")); 77 dialog->queryString("output-file", ValueStorage::getInstance().getDescription("output-file")); 78 dialog->queryString("bin-output-file", ValueStorage::getInstance().getDescription("bin-output-file")); 79 dialog->queryBoolean("periodic", ValueStorage::getInstance().getDescription("periodic")); 80 81 return dialog; 82 } 83 46 /** =========== define the function ====================== */ 84 47 Action::state_ptr AnalysisPointCorrelationAction::performCall() { 85 48 int ranges[3] = {1, 1, 1}; 86 double BinEnd = 0.;87 double BinStart = 0.;88 double BinWidth = 0.;89 string outputname;90 string binoutputname;91 bool periodic;92 49 ofstream output; 93 50 ofstream binoutput; 94 std::vector<const element *> elements;95 51 string type; 96 Vector Point;97 52 BinPairMap *binmap = NULL; 98 53 99 54 // obtain information 100 ValueStorage::getInstance().queryCurrentValue("position", Point); 101 ValueStorage::getInstance().queryCurrentValue("elements", elements); 102 ValueStorage::getInstance().queryCurrentValue("bin-start", BinStart); 103 ValueStorage::getInstance().queryCurrentValue("bin-width", BinWidth); 104 ValueStorage::getInstance().queryCurrentValue("bin-end", BinEnd); 105 ValueStorage::getInstance().queryCurrentValue("output-file", outputname); 106 ValueStorage::getInstance().queryCurrentValue("bin-output-file", binoutputname); 107 ValueStorage::getInstance().queryCurrentValue("periodic", periodic); 55 getParametersfromValueStorage(); 108 56 109 57 // execute action 110 output.open( outputname.c_str());111 binoutput.open( binoutputname.c_str());112 cout << "Point to correlate to is " << Point << endl;58 output.open(params.outputname.c_str()); 59 binoutput.open(params.binoutputname.c_str()); 60 cout << "Point to correlate to is " << params.Point << endl; 113 61 CorrelationToPointMap *correlationmap = NULL; 114 for(std::vector<const element *>::iterator iter = elements.begin(); iter !=elements.end(); ++iter)62 for(std::vector<const element *>::iterator iter = params.elements.begin(); iter != params.elements.end(); ++iter) 115 63 cout << "element is " << (*iter)->getSymbol() << endl; 116 64 std::vector<molecule*> molecules = World::getInstance().getSelectedMolecules(); 117 if (p eriodic)118 correlationmap = PeriodicCorrelationToPoint(molecules, elements, &Point, ranges);65 if (params.periodic) 66 correlationmap = PeriodicCorrelationToPoint(molecules, params.elements, ¶ms.Point, ranges); 119 67 else 120 correlationmap = CorrelationToPoint(molecules, elements, &Point);68 correlationmap = CorrelationToPoint(molecules, params.elements, ¶ms.Point); 121 69 OutputCorrelationToPoint(&output, correlationmap); 122 binmap = BinData( correlationmap, BinWidth, BinStart,BinEnd );70 binmap = BinData( correlationmap, params.BinWidth, params.BinStart, params.BinEnd ); 123 71 OutputCorrelation ( &binoutput, binmap ); 124 72 delete(binmap); … … 148 96 return NAME; 149 97 } 98 /** =========== end of function ====================== */
Note:
See TracChangeset
for help on using the changeset viewer.
