Changeset bfe2c2 for src/Actions
- Timestamp:
- Jul 22, 2010, 9:11:25 PM (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, 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:
- 22e780
- Parents:
- 9a9b2a
- Location:
- src/Actions/WorldAction
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/WorldAction/RepeatBoxAction.cpp
r9a9b2a rbfe2c2 25 25 #include "UIElements/UIFactory.hpp" 26 26 #include "UIElements/Dialog.hpp" 27 #include " Actions/MapOfActions.hpp"27 #include "UIElements/ValueStorage.hpp" 28 28 #include "Descriptors/MoleculeDescriptor.hpp" 29 29 #include "Descriptors/MoleculePtrDescriptor.hpp" … … 38 38 {} 39 39 40 Dialog * WorldRepeatBoxAction::createDialog() { 41 Dialog *dialog = UIFactory::getInstance().makeDialog(); 42 43 dialog->queryVector(NAME, false, ValueStorage::getInstance().getDescription(NAME)); 44 45 return dialog; 46 } 47 40 48 Action::state_ptr WorldRepeatBoxAction::performCall() { 41 Dialog *dialog = UIFactory::getInstance().makeDialog();42 49 Vector Repeater; 43 50 int count; … … 48 55 MoleculeListClass *molecules = World::getInstance().getMolecules(); 49 56 50 dialog->queryVector(NAME, &Repeater, false, MapOfActions::getInstance().getDescription(NAME));51 //dialog->queryMolecule("molecule-by-id", &mol,MapOfActions::getInstance().getDescription("molecule-by-id")); 57 ValueStorage::getInstance().queryCurrentValue(NAME, Repeater); 58 52 59 vector<molecule *> AllMolecules; 53 60 if (mol != NULL) { … … 59 66 } 60 67 61 if(dialog->display()) { 62 (cout << "Repeating box " << Repeater << " times for (x,y,z) axis." << endl); 63 Matrix M = World::getInstance().getDomain().getM(); 64 Matrix newM = M; 65 Vector x,y; 66 int n[NDIM]; 67 Matrix repMat; 68 for (int axis = 0; axis < NDIM; axis++) { 69 Repeater[axis] = floor(Repeater[axis]); 70 if (Repeater[axis] < 1) { 71 DoeLog(1) && (eLog()<< Verbose(1) << "Repetition factor must be greater than 1!" << endl); 72 Repeater[axis] = 1; 73 } 74 repMat.at(axis,axis) = Repeater[axis]; 68 (cout << "Repeating box " << Repeater << " times for (x,y,z) axis." << endl); 69 Matrix M = World::getInstance().getDomain().getM(); 70 Matrix newM = M; 71 Vector x,y; 72 int n[NDIM]; 73 Matrix repMat; 74 for (int axis = 0; axis < NDIM; axis++) { 75 Repeater[axis] = floor(Repeater[axis]); 76 if (Repeater[axis] < 1) { 77 DoeLog(1) && (eLog()<< Verbose(1) << "Repetition factor must be greater than 1!" << endl); 78 Repeater[axis] = 1; 75 79 } 76 newM *= repMat; 77 World::getInstance().setDomain(newM); 80 repMat.at(axis,axis) = Repeater[axis]; 81 } 82 newM *= repMat; 83 World::getInstance().setDomain(newM); 78 84 79 molecule *newmol = NULL; 80 Vector ** vectors = NULL; 81 for (n[0] = 0; n[0] < Repeater[0]; n[0]++) { 82 y[0] = n[0]; 83 for (n[1] = 0; n[1] < Repeater[1]; n[1]++) { 84 y[1] = n[1]; 85 for (n[2] = 0; n[2] < Repeater[2]; n[2]++) { 86 y[2] = n[2]; 87 if ((n[0] == 0) && (n[1] == 0) && (n[2] == 0)) 88 continue; 89 for (vector<molecule *>::iterator MolRunner = AllMolecules.begin(); MolRunner != AllMolecules.end(); ++MolRunner) { 90 mol = *MolRunner; 91 DoLog(1) && (Log() << Verbose(1) << "Current mol is " << mol->name << "." << endl); 92 count = mol->getAtomCount(); // is changed becausing of adding, thus has to be stored away beforehand 93 if (count != 0) { // if there is more than none 94 Elements = new const element *[count]; 95 vectors = new Vector *[count]; 96 j = 0; 97 for(molecule::iterator AtomRunner = mol->begin(); AtomRunner != mol->end(); ++AtomRunner) { 98 Elements[j] = (*AtomRunner)->type; 99 vectors[j] = &(*AtomRunner)->x; 100 j++; 101 } 102 if (count != j) 103 DoeLog(1) && (eLog()<< Verbose(1) << "AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl); 104 x = y; 105 x *= M; 106 newmol = World::getInstance().createMolecule(); 107 molecules->insert(newmol); 108 for (int k=count;k--;) { // go through every atom of the original cell 109 Walker = World::getInstance().createAtom(); // create a new body 110 Walker->x = (*vectors[k]) + x; 111 Walker->type = Elements[k]; // insert original element 112 cout << "new atom is " << *Walker << endl; 113 newmol->AddAtom(Walker); // and add to the molecule (which increments ElementsInMolecule, AtomCount, ...) 114 } 115 // free memory 116 delete[](Elements); 117 delete[](vectors); 118 } else { 119 DoLog(1) && (Log() << Verbose(1) << "\t ... is empty." << endl); 85 molecule *newmol = NULL; 86 Vector ** vectors = NULL; 87 for (n[0] = 0; n[0] < Repeater[0]; n[0]++) { 88 y[0] = n[0]; 89 for (n[1] = 0; n[1] < Repeater[1]; n[1]++) { 90 y[1] = n[1]; 91 for (n[2] = 0; n[2] < Repeater[2]; n[2]++) { 92 y[2] = n[2]; 93 if ((n[0] == 0) && (n[1] == 0) && (n[2] == 0)) 94 continue; 95 for (vector<molecule *>::iterator MolRunner = AllMolecules.begin(); MolRunner != AllMolecules.end(); ++MolRunner) { 96 mol = *MolRunner; 97 DoLog(1) && (Log() << Verbose(1) << "Current mol is " << mol->name << "." << endl); 98 count = mol->getAtomCount(); // is changed becausing of adding, thus has to be stored away beforehand 99 if (count != 0) { // if there is more than none 100 Elements = new const element *[count]; 101 vectors = new Vector *[count]; 102 j = 0; 103 for(molecule::iterator AtomRunner = mol->begin(); AtomRunner != mol->end(); ++AtomRunner) { 104 Elements[j] = (*AtomRunner)->type; 105 vectors[j] = &(*AtomRunner)->x; 106 j++; 120 107 } 108 if (count != j) 109 DoeLog(1) && (eLog()<< Verbose(1) << "AtomCount " << count << " is not equal to number of atoms in molecule " << j << "!" << endl); 110 x = y; 111 x *= M; 112 newmol = World::getInstance().createMolecule(); 113 molecules->insert(newmol); 114 for (int k=count;k--;) { // go through every atom of the original cell 115 Walker = World::getInstance().createAtom(); // create a new body 116 Walker->x = (*vectors[k]) + x; 117 Walker->type = Elements[k]; // insert original element 118 cout << "new atom is " << *Walker << endl; 119 newmol->AddAtom(Walker); // and add to the molecule (which increments ElementsInMolecule, AtomCount, ...) 120 } 121 // free memory 122 delete[](Elements); 123 delete[](vectors); 124 } else { 125 DoLog(1) && (Log() << Verbose(1) << "\t ... is empty." << endl); 121 126 } 122 127 } 123 128 } 124 129 } 125 delete dialog;126 return Action::success;127 } else {128 delete dialog;129 return Action::failure;130 130 } 131 return Action::success; 131 132 } 132 133 -
src/Actions/WorldAction/RepeatBoxAction.hpp
r9a9b2a rbfe2c2 21 21 virtual const std::string getName(); 22 22 private: 23 virtual Dialog * createDialog(); 23 24 virtual Action::state_ptr performCall(); 24 25 virtual Action::state_ptr performUndo(Action::state_ptr);
Note:
See TracChangeset
for help on using the changeset viewer.