source: src/unittests/manipulateAtomsTest.cpp@ b54ac8

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 b54ac8 was b54ac8, checked in by Tillmann Crueger <crueger@…>, 15 years ago

Added templates that allow arbitrary calculations on atoms to be mapped to sets of Atoms

  • Property mode set to 100644
File size: 4.3 KB
Line 
1/*
2 * manipulateAtomsTest.cpp
3 *
4 * Created on: Feb 18, 2010
5 * Author: crueger
6 */
7
8#include "manipulateAtomsTest.hpp"
9
10#include <cppunit/CompilerOutputter.h>
11#include <cppunit/extensions/TestFactoryRegistry.h>
12#include <cppunit/ui/text/TestRunner.h>
13#include <iostream>
14#include <boost/bind.hpp>
15
16#include "Descriptors/AtomDescriptor.hpp"
17#include "Descriptors/AtomIdDescriptor.hpp"
18#include "Actions/ManipulateAtomsProcess.hpp"
19#include "Actions/ActionRegistry.hpp"
20
21#include "World.hpp"
22#include "atom.hpp"
23
24// Registers the fixture into the 'registry'
25CPPUNIT_TEST_SUITE_REGISTRATION( manipulateAtomsTest );
26
27// some stubs
28class AtomStub : public atom {
29public:
30 AtomStub(int _id) :
31 atom(),
32 id(_id),
33 manipulated(false)
34 {}
35
36 virtual int getId(){
37 return id;
38 }
39
40 virtual void doSomething(){
41 manipulated = true;
42 }
43
44 bool manipulated;
45private:
46 int id;
47};
48
49class countObserver : public Observer{
50public:
51 countObserver() :
52 count(0)
53 {}
54 virtual ~countObserver(){}
55
56 void update(Observable *){
57 count++;
58 }
59
60 void subjectKilled(Observable *)
61 {}
62
63 int count;
64};
65
66// set up and tear down
67void manipulateAtomsTest::setUp(){
68 World::get();
69 for(int i=0;i<ATOM_COUNT;++i){
70 atoms[i]= new AtomStub(i);
71 }
72}
73void manipulateAtomsTest::tearDown(){
74 World::destroy();
75 for(int i=0;i<ATOM_COUNT;++i){
76 delete atoms[i];
77 }
78 ActionRegistry::purgeRegistry();
79}
80
81// some helper functions
82bool hasAll(std::vector<atom*> atoms,int min, int max, std::set<int> excluded = std::set<int>()){
83 for(int i=min;i<max;++i){
84 if(!excluded.count(i)){
85 std::vector<atom*>::iterator iter;
86 bool res=false;
87 for(iter=atoms.begin();iter!=atoms.end();++iter){
88 res |= (*iter)->getId() == i;
89 }
90 if(!res) {
91 cout << "Atom " << i << " missing in returned list" << endl;
92 return false;
93 }
94 }
95 }
96 return true;
97}
98
99bool hasNoDuplicates(std::vector<atom*> atoms){
100 std::set<int> found;
101 std::vector<atom*>::iterator iter;
102 for(iter=atoms.begin();iter!=atoms.end();++iter){
103 int id = (*iter)->getId();
104 if(found.count(id))
105 return false;
106 found.insert(id);
107 }
108 return true;
109}
110
111void operation(atom* _atom){
112 AtomStub *atom = dynamic_cast<AtomStub*>(_atom);
113 assert(atom);
114 atom->doSomething();
115}
116
117
118void manipulateAtomsTest::testManipulateSimple(){
119 ManipulateAtomsProcess *proc = World::get()->manipulateAtoms(boost::bind(operation,_1),"FOO",AllAtoms());
120 proc->call();
121 std::vector<atom*> allAtoms = World::get()->getAllAtoms(AllAtoms());
122 std::vector<atom*>::iterator iter;
123 for(iter=allAtoms.begin();iter!=allAtoms.end();++iter){
124 AtomStub *atom;
125 atom = dynamic_cast<AtomStub*>(*iter);
126 assert(atom);
127 CPPUNIT_ASSERT(atom->manipulated);
128 }
129}
130
131void manipulateAtomsTest::testManipulateExcluded(){
132 ManipulateAtomsProcess *proc = World::get()->manipulateAtoms(boost::bind(operation,_1),"FOO",AllAtoms() && !AtomById(ATOM_COUNT/2));
133 proc->call();
134 std::vector<atom*> allAtoms = World::get()->getAllAtoms(AllAtoms());
135 std::vector<atom*>::iterator iter;
136 for(iter=allAtoms.begin();iter!=allAtoms.end();++iter){
137 AtomStub *atom;
138 atom = dynamic_cast<AtomStub*>(*iter);
139 assert(atom);
140 if(atom->getId()!=(int)ATOM_COUNT/2)
141 CPPUNIT_ASSERT(atom->manipulated);
142 else
143 CPPUNIT_ASSERT(!atom->manipulated);
144 }
145}
146
147void manipulateAtomsTest::testObserver(){
148 countObserver *obs = new countObserver();
149 World::get()->signOn(obs);
150 ManipulateAtomsProcess *proc = World::get()->manipulateAtoms(boost::bind(operation,_1),"FOO",AllAtoms() && !AtomById(ATOM_COUNT/2));
151 proc->call();
152
153 CPPUNIT_ASSERT_EQUAL(1,obs->count);
154 World::get()->signOff(obs);
155 delete obs;
156}
157
158/********************************************** Main routine **************************************/
159
160int main(int argc, char **argv)
161{
162 // Get the top level suite from the registry
163 CppUnit::Test *suite = CppUnit::TestFactoryRegistry::getRegistry().makeTest();
164
165 // Adds the test to the list of test to run
166 CppUnit::TextUi::TestRunner runner;
167 runner.addTest( suite );
168
169 // Change the default outputter to a compiler error format outputter
170 runner.setOutputter( new CppUnit::CompilerOutputter( &runner.result(),
171 std::cerr ) );
172 // Run the tests.
173 bool wasSucessful = runner.run();
174
175 // Return error code 1 if the one of test failed.
176 return wasSucessful ? 0 : 1;
177};
Note: See TracBrowser for help on using the repository browser.