Changeset bfe2c2 for src/Actions


Ignore:
Timestamp:
Jul 22, 2010, 9:11:25 PM (15 years ago)
Author:
Frederik Heber <heber@…>
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
Message:

converted RepeatBoxAction to new createDialog().

Location:
src/Actions/WorldAction
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • src/Actions/WorldAction/RepeatBoxAction.cpp

    r9a9b2a rbfe2c2  
    2525#include "UIElements/UIFactory.hpp"
    2626#include "UIElements/Dialog.hpp"
    27 #include "Actions/MapOfActions.hpp"
     27#include "UIElements/ValueStorage.hpp"
    2828#include "Descriptors/MoleculeDescriptor.hpp"
    2929#include "Descriptors/MoleculePtrDescriptor.hpp"
     
    3838{}
    3939
     40Dialog * WorldRepeatBoxAction::createDialog() {
     41  Dialog *dialog = UIFactory::getInstance().makeDialog();
     42
     43  dialog->queryVector(NAME, false, ValueStorage::getInstance().getDescription(NAME));
     44
     45  return dialog;
     46}
     47
    4048Action::state_ptr WorldRepeatBoxAction::performCall() {
    41   Dialog *dialog = UIFactory::getInstance().makeDialog();
    4249  Vector Repeater;
    4350  int count;
     
    4855  MoleculeListClass *molecules = World::getInstance().getMolecules();
    4956
    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
    5259  vector<molecule *> AllMolecules;
    5360  if (mol != NULL) {
     
    5966  }
    6067
    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;
    7579    }
    76     newM *= repMat;
    77     World::getInstance().setDomain(newM);
     80    repMat.at(axis,axis) = Repeater[axis];
     81  }
     82  newM *= repMat;
     83  World::getInstance().setDomain(newM);
    7884
    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++;
    120107            }
     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);
    121126          }
    122127        }
    123128      }
    124129    }
    125     delete dialog;
    126     return Action::success;
    127   } else {
    128     delete dialog;
    129     return Action::failure;
    130130  }
     131  return Action::success;
    131132}
    132133
  • src/Actions/WorldAction/RepeatBoxAction.hpp

    r9a9b2a rbfe2c2  
    2121  virtual const std::string getName();
    2222private:
     23  virtual Dialog * createDialog();
    2324  virtual Action::state_ptr performCall();
    2425  virtual Action::state_ptr performUndo(Action::state_ptr);
Note: See TracChangeset for help on using the changeset viewer.