Changeset de29ad6 for src/Box.cpp

Timestamp:

Dec 28, 2011, 3:24:23 PM (13 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:

68c923

Parents:

bde4a6

git-author:

Frederik Heber <heber@…> (12/28/11 14:58:52)

git-committer:

Frederik Heber <heber@…> (12/28/11 15:24:23)

Message:

Some optimization to speed up Subgraph dissections.

• Replaced std::list by std::vector in VectorSet<>::minDistSquared() and Box::internal_explode() as it caused lots of dynamic allocation. This has been the main cause of the slowdown of Box::internal_explode().
• Box has internalized vector<int>s coords and index of internal_explode() to avoid dynamic allocation for them as well. This is worth it as it has a heavily used function.
• Made Box internal_list non-static. There is only one box anyway.
• changes caused in BoxUnitTest and Box cstor's and dstor.
• speedup of Subgraph

• enhanced verbosity of BondGraph::CreateAdjacency().

File:

• src/Box.cpp (modified) (11 diffs)

src/Box.cpp

@@ -37 +37 @@
 
 
-using namespace std;
-
-VECTORSET(std::list) Box::internal_list;
-
 Box::Box() :
     M(new RealSpaceMatrix()),
     Minv(new RealSpaceMatrix())
 {
+  internal_list.reserve(pow(3,3));
+  coords.reserve(NDIM);
+  index.reserve(NDIM);
   M->setIdentity();
   Minv->setIdentity();
@@ -55 +54 @@
   M(new RealSpaceMatrix(*src.M)),
   Minv(new RealSpaceMatrix(*src.Minv))
-{}
+{
+  internal_list.reserve(pow(3,3));
+  coords.reserve(NDIM);
+  index.reserve(NDIM);
+}
 
 Box::Box(RealSpaceMatrix _M) :
@@ -61 +64 @@
     Minv(new RealSpaceMatrix())
 {
+  internal_list.reserve(pow(3,3));
+  coords.reserve(NDIM);
+  index.reserve(NDIM);
   ASSERT(M->determinant()!=0,"Matrix in Box construction was not invertible");
   *Minv = M->invert();
@@ -135 +141 @@
 
 
-VECTORSET(std::list) Box::explode(const Vector &point,int n) const{
+VECTORSET(std::vector) Box::explode(const Vector &point,int n) const{
   ASSERT(isInside(point),"Exploded point not inside Box");
   internal_explode(point, n);
-  VECTORSET(std::list) res(internal_list);
+  VECTORSET(std::vector) res(internal_list);
   return res;
 }
@@ -144 +150 @@
 void Box::internal_explode(const Vector &point,int n) const{
   internal_list.clear();
+  size_t list_index = 0;
 
   Vector translater = translateOut(point);
@@ -150 +157 @@
   // count the number of coordinates we need to do
   int dims = 0; // number of dimensions that are not ignored
-  vector<int> coords;
-  vector<int> index;
+  coords.clear();
+  index.clear();
   for(int i=0;i<NDIM;++i){
     if(conditions[i]==Ignore){
@@ -161 +168 @@
     dims++;
   } // there are max vectors in total we need to create
+  internal_list.resize(pow(2*n+1, dims));
 
   if(!dims){
     // all boundaries are ignored
-    internal_list.push_back(point);
+    internal_list[list_index++] = point;
     return;
   }
@@ -206 +214 @@
     // add back all ignored coordinates (not handled in above loop)
     helper+=mask;
-    internal_list.push_back(translateIn(helper));
+    ASSERT(list_index < internal_list.size(),
+        "Box::internal_explode() - we have estimated the number of vectors wrong: "
+        +toString(list_index) +" >= "+toString(internal_list.size())+".");
+    internal_list[list_index++] = translateIn(helper);
     // set the new indexes
     int pos=0;
@@ -221 +232 @@
 }
 
-VECTORSET(std::list) Box::explode(const Vector &point) const{
+VECTORSET(std::vector) Box::explode(const Vector &point) const{
   ASSERT(isInside(point),"Exploded point not inside Box");
   return explode(point,1);
@@ -268 +279 @@
 }
 
-const vector<pair<Plane,Plane> >  Box::getBoundingPlanes() const
-{
-  vector<pair<Plane,Plane> > res;
+const std::vector<std::pair<Plane,Plane> >  Box::getBoundingPlanes() const
+{
+  std::vector<std::pair<Plane,Plane> > res;
   for(int i=0;i<NDIM;++i){
     Vector base1,base2,base3;
@@ -314 +325 @@
 }
 
-ostream & operator << (ostream& ost, const Box &m)
+std::ostream & operator << (std::ostream& ost, const Box &m)
 {
   ost << m.getM();