Changeset 1fd675 for src/Actions/TesselationAction
- Timestamp:
- Aug 28, 2010, 12:57:56 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, 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:
- 966e12
- Parents:
- 53d01c
- git-author:
- Frederik Heber <heber@…> (08/26/10 21:05:34)
- git-committer:
- Frederik Heber <heber@…> (08/28/10 00:57:56)
- Location:
- src/Actions/TesselationAction
- Files:
-
- 2 added
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/TesselationAction/ConvexEnvelopeAction.cpp
r53d01c r1fd675 20 20 #include "Helpers/MemDebug.hpp" 21 21 22 #include "Actions/TesselationAction/ConvexEnvelopeAction.hpp"23 #include "Actions/ActionRegistry.hpp"24 22 #include "boundary.hpp" 25 23 #include "config.hpp" … … 36 34 using namespace std; 37 35 38 #include "UIElements/UIFactory.hpp" 39 #include "UIElements/Dialog.hpp" 40 #include "Actions/ValueStorage.hpp" 36 #include "Actions/TesselationAction/ConvexEnvelopeAction.hpp" 41 37 42 /****** TesselationConvexEnvelopeAction *****/ 43 44 // memento to remember the state when undoing 45 46 //class TesselationConvexEnvelopeState : public ActionState { 47 //public: 48 // TesselationConvexEnvelopeState(molecule* _mol,std::string _lastName) : 49 // mol(_mol), 50 // lastName(_lastName) 51 // {} 52 // molecule* mol; 53 // std::string lastName; 54 //}; 55 56 const char TesselationConvexEnvelopeAction::NAME[] = "convex-envelope"; 57 58 TesselationConvexEnvelopeAction::TesselationConvexEnvelopeAction() : 59 Action(NAME) 60 {} 61 62 TesselationConvexEnvelopeAction::~TesselationConvexEnvelopeAction() 63 {} 64 65 void TesselationConvexEnvelope(std::string &filenameConvex, std::string &filenameNonConvex) { 66 ValueStorage::getInstance().setCurrentValue("convex-file", filenameConvex); 67 ValueStorage::getInstance().setCurrentValue("nonconvex-file", filenameConvex); 68 ActionRegistry::getInstance().getActionByName(TesselationConvexEnvelopeAction::NAME)->call(Action::NonInteractive); 69 }; 70 71 void TesselationConvexEnvelopeAction::getParametersfromValueStorage() 72 {}; 73 74 Dialog* TesselationConvexEnvelopeAction::fillDialog(Dialog *dialog) { 75 ASSERT(dialog,"No Dialog given when filling action dialog"); 76 77 dialog->queryEmpty(NAME, ValueStorage::getInstance().getDescription(NAME)); 78 dialog->queryString("convex-file", ValueStorage::getInstance().getDescription("convex-file")); 79 dialog->queryString("nonconvex-file", ValueStorage::getInstance().getDescription("nonconvex-file")); 80 81 return dialog; 82 } 83 38 // and construct the stuff 39 #include "ConvexEnvelopeAction.def" 40 #include "Action_impl_pre.hpp" 41 /** =========== define the function ====================== */ 84 42 Action::state_ptr TesselationConvexEnvelopeAction::performCall() { 85 string filenameConvex;86 string filenameNonConvex;87 43 molecule * mol = NULL; 88 44 bool Success = true; 89 45 config *configuration = World::getInstance().getConfig(); 90 46 91 ValueStorage::getInstance().queryCurrentValue("convex-file", filenameConvex);92 ValueStorage::getInstance().queryCurrentValue("nonconvex-file", filenameNonConvex);47 // obtain information 48 getParametersfromValueStorage(); 93 49 94 50 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { … … 97 53 const LinkedCell *LCList = NULL; 98 54 DoLog(0) && (Log() << Verbose(0) << "Evaluating volume of the convex envelope."); 99 DoLog(1) && (Log() << Verbose(1) << "Storing tecplot convex data in " << filenameConvex << "." << endl);100 DoLog(1) && (Log() << Verbose(1) << "Storing tecplot non-convex data in " << filenameNonConvex << "." << endl);55 DoLog(1) && (Log() << Verbose(1) << "Storing tecplot convex data in " << params.filenameConvex << "." << endl); 56 DoLog(1) && (Log() << Verbose(1) << "Storing tecplot non-convex data in " << params.filenameNonConvex << "." << endl); 101 57 LCList = new LinkedCell(mol, 100.); 102 58 //Boundaries *BoundaryPoints = NULL; 103 59 //FindConvexBorder(mol, BoundaryPoints, TesselStruct, LCList, argv[argptr]); 104 60 // TODO: Beide Funktionen sollten streams anstelle des Filenamen benutzen, besser fuer unit tests 105 FindNonConvexBorder(mol, TesselStruct, LCList, 50., filenameNonConvex.c_str());61 FindNonConvexBorder(mol, TesselStruct, LCList, 50., params.filenameNonConvex.c_str()); 106 62 //RemoveAllBoundaryPoints(TesselStruct, mol, argv[argptr]); 107 const double volumedifference = ConvexizeNonconvexEnvelope(TesselStruct, mol, filenameConvex.c_str());63 const double volumedifference = ConvexizeNonconvexEnvelope(TesselStruct, mol, params.filenameConvex.c_str()); 108 64 const double clustervolume = VolumeOfConvexEnvelope(TesselStruct, configuration); 109 65 DoLog(0) && (Log() << Verbose(0) << "The tesselated volume area is " << clustervolume << " " << (configuration->GetIsAngstroem() ? "angstrom" : "atomiclength") << "^3." << endl); … … 143 99 return NAME; 144 100 } 101 /** =========== end of function ====================== */ -
src/Actions/TesselationAction/ConvexEnvelopeAction.hpp
r53d01c r1fd675 11 11 #include "Actions/Action.hpp" 12 12 13 14 class TesselationListClass; 15 16 void TesselationConvexEnvelope(std::string &filenameConvex, std::string &filenameNonConvex); 17 18 class TesselationConvexEnvelopeAction : public Action { 19 friend void TesselationConvexEnvelope(std::string &filenameConvex, std::string &filenameNonConvex); 20 21 public: 22 TesselationConvexEnvelopeAction(); 23 virtual ~TesselationConvexEnvelopeAction(); 24 25 bool canUndo(); 26 bool shouldUndo(); 27 28 virtual const std::string getName(); 29 protected: 30 virtual Dialog * fillDialog(Dialog *dialog); 31 private: 32 virtual void getParametersfromValueStorage(); 33 virtual Action::state_ptr performCall(); 34 virtual Action::state_ptr performUndo(Action::state_ptr); 35 virtual Action::state_ptr performRedo(Action::state_ptr); 36 37 static const char NAME[]; 38 }; 13 #include "ConvexEnvelopeAction.def" 14 #include "Action_impl_header.hpp" 39 15 40 16 #endif // CONVEXENVELOPEACTION_HPP -
src/Actions/TesselationAction/NonConvexEnvelopeAction.cpp
r53d01c r1fd675 20 20 #include "Helpers/MemDebug.hpp" 21 21 22 #include "Actions/TesselationAction/NonConvexEnvelopeAction.hpp"23 #include "Actions/ActionRegistry.hpp"24 22 #include "boundary.hpp" 25 23 #include "linkedcell.hpp" … … 35 33 using namespace std; 36 34 37 #include "UIElements/UIFactory.hpp" 38 #include "UIElements/Dialog.hpp" 39 #include "Actions/ValueStorage.hpp" 35 #include "Actions/TesselationAction/NonConvexEnvelopeAction.hpp" 40 36 41 /****** TesselationNonConvexEnvelopeAction *****/ 42 43 // memento to remember the state when undoing 44 45 //class TesselationNonConvexEnvelopeState : public ActionState { 46 //public: 47 // TesselationNonConvexEnvelopeState(molecule* _mol,std::string _lastName) : 48 // mol(_mol), 49 // lastName(_lastName) 50 // {} 51 // molecule* mol; 52 // std::string lastName; 53 //}; 54 55 const char TesselationNonConvexEnvelopeAction::NAME[] = "nonconvex-envelope"; 56 57 TesselationNonConvexEnvelopeAction::TesselationNonConvexEnvelopeAction() : 58 Action(NAME) 59 {} 60 61 TesselationNonConvexEnvelopeAction::~TesselationNonConvexEnvelopeAction() 62 {} 63 64 void TesselationNonConvexEnvelope(double radius, std::string &filename) { 65 ValueStorage::getInstance().setCurrentValue(TesselationNonConvexEnvelopeAction::NAME, radius); 66 ValueStorage::getInstance().setCurrentValue("nonconvex-file", filename); 67 ActionRegistry::getInstance().getActionByName(TesselationNonConvexEnvelopeAction::NAME)->call(Action::NonInteractive); 68 }; 69 70 void TesselationNonConvexEnvelopeAction::getParametersfromValueStorage() 71 {}; 72 73 Dialog* TesselationNonConvexEnvelopeAction::fillDialog(Dialog *dialog) { 74 ASSERT(dialog,"No Dialog given when filling action dialog"); 75 76 dialog->queryDouble(NAME, ValueStorage::getInstance().getDescription(NAME)); 77 dialog->queryString("nonconvex-file", ValueStorage::getInstance().getDescription("nonconvex-file")); 78 79 return dialog; 80 } 81 37 // and construct the stuff 38 #include "NonConvexEnvelopeAction.def" 39 #include "Action_impl_pre.hpp" 40 /** =========== define the function ====================== */ 82 41 Action::state_ptr TesselationNonConvexEnvelopeAction::performCall() { 83 string filename;84 42 molecule * Boundary = NULL; 85 double SphereRadius = 0;86 43 bool Success = false; 87 44 clock_t start,end; 88 45 89 ValueStorage::getInstance().queryCurrentValue(NAME, SphereRadius);90 ValueStorage::getInstance().queryCurrentValue("nonconvex-file", filename);46 // obtain information 47 getParametersfromValueStorage(); 91 48 92 49 for (World::MoleculeSelectionIterator iter = World::getInstance().beginMoleculeSelection(); iter != World::getInstance().endMoleculeSelection(); ++iter) { … … 95 52 const LinkedCell *LCList = NULL; 96 53 DoLog(0) && (Log() << Verbose(0) << "Evaluating non-convex envelope of molecule." << Boundary->getId() << endl); 97 DoLog(1) && (Log() << Verbose(1) << "Using rolling ball of radius " << SphereRadius << " and storing tecplot data in " <<filename << "." << endl);54 DoLog(1) && (Log() << Verbose(1) << "Using rolling ball of radius " << params.SphereRadius << " and storing tecplot data in " << params.filename << "." << endl); 98 55 DoLog(1) && (Log() << Verbose(1) << "Specified molecule has " << Boundary->getAtomCount() << " atoms." << endl); 99 56 start = clock(); 100 LCList = new LinkedCell(Boundary, SphereRadius*2.);101 Success = FindNonConvexBorder(Boundary, T, LCList, SphereRadius,filename.c_str());102 //FindDistributionOfEllipsoids(T, &LCList, N, number, filename.c_str());57 LCList = new LinkedCell(Boundary, params.SphereRadius*2.); 58 Success = FindNonConvexBorder(Boundary, T, LCList, params.SphereRadius, params.filename.c_str()); 59 //FindDistributionOfEllipsoids(T, &LCList, N, number, params.filename.c_str()); 103 60 end = clock(); 104 61 DoLog(0) && (Log() << Verbose(0) << "Clocks for this operation: " << (end-start) << ", time: " << ((double)(end-start)/CLOCKS_PER_SEC) << "s." << endl); … … 137 94 return NAME; 138 95 } 96 /** =========== end of function ====================== */ -
src/Actions/TesselationAction/NonConvexEnvelopeAction.hpp
r53d01c r1fd675 11 11 #include "Actions/Action.hpp" 12 12 13 14 class TesselationListClass; 15 16 void TesselationNonConvexEnvelope(double radius, std::string &filename); 17 18 class TesselationNonConvexEnvelopeAction : public Action { 19 friend void TesselationNonConvexEnvelope(double radius, std::string &filename); 20 21 public: 22 TesselationNonConvexEnvelopeAction(); 23 virtual ~TesselationNonConvexEnvelopeAction(); 24 25 bool canUndo(); 26 bool shouldUndo(); 27 28 virtual const std::string getName(); 29 protected: 30 virtual Dialog * fillDialog(Dialog *dialog); 31 private: 32 virtual void getParametersfromValueStorage(); 33 virtual Action::state_ptr performCall(); 34 virtual Action::state_ptr performUndo(Action::state_ptr); 35 virtual Action::state_ptr performRedo(Action::state_ptr); 36 37 static const char NAME[]; 38 }; 39 13 #include "NonConvexEnvelopeAction.def" 14 #include "Action_impl_header.hpp" 40 15 41 16 #endif // NONCONVEXENVELOPEACTION_HPP
Note:
See TracChangeset
for help on using the changeset viewer.