Changeset c5e75f3 for src/Potentials

Timestamp:

Nov 4, 2016, 9:37:50 AM (8 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Action_Thermostats, Add_AtomRandomPerturbation, Add_RotateAroundBondAction, Add_SelectAtomByNameAction, Adding_Graph_to_ChangeBondActions, Adding_MD_integration_tests, Adding_StructOpt_integration_tests, Automaking_mpqc_open, AutomationFragmentation_failures, Candidate_v1.6.0, Candidate_v1.6.1, ChangeBugEmailaddress, ChangingTestPorts, ChemicalSpaceEvaluator, 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_ChronosMutex, Fix_StatusMsg, Fix_StepWorldTime_single_argument, Fix_Verbose_Codepatterns, ForceAnnealing_goodresults, ForceAnnealing_oldresults, ForceAnnealing_tocheck, ForceAnnealing_with_BondGraph, ForceAnnealing_with_BondGraph_continued, ForceAnnealing_with_BondGraph_continued_betteresults, ForceAnnealing_with_BondGraph_contraction-expansion, GeometryObjects, Gui_displays_atomic_force_velocity, IndependentFragmentGrids_IntegrationTest, JobMarket_RobustOnKillsSegFaults, JobMarket_StableWorkerPool, JobMarket_unresolvable_hostname_fix, ODR_violation_mpqc_open, PartialCharges_OrthogonalSummation, PythonUI_with_named_parameters, QtGui_reactivate_TimeChanged_changes, Recreated_GuiChecks, RotateToPrincipalAxisSystem_UndoRedo, StoppableMakroAction, Subpackage_CodePatterns, Subpackage_JobMarket, Subpackage_LinearAlgebra, Subpackage_levmar, Subpackage_mpqc_open, Subpackage_vmg, ThirdParty_MPQC_rebuilt_buildsystem, TremoloParser_IncreasedPrecision, TremoloParser_MultipleTimesteps, Ubuntu_1604_changes, stable

Children:

d5ca1a

Parents:

c7aac9

git-author:

Frederik Heber <heber@…> (10/03/16 10:35:00)

git-committer:

Frederik Heber <heber@…> (11/04/16 09:37:50)

Message:

ParticleType_t is now an unsigned int.

• this allows direct use of atomicNumber_t and also signed integers dont' make much sense for a particle type, i.e. an id.

Location:

src/Potentials

Files:

• CompoundPotential.cpp (modified) (2 diffs)
• PotentialTrainer.cpp (modified) (2 diffs)
• SerializablePotential.hpp (modified) (1 diff)

src/Potentials/CompoundPotential.cpp

@@ -70 +70 @@
         potential->getParticleTypes();
 
-    // create charges
-    Fragment::charges_t charges;
-    charges.resize(types.size());
-    std::transform(types.begin(), types.end(),
-        charges.begin(), boost::lambda::_1);
     // convert into count map
-    Extractors::elementcounts_t counts_per_charge =
-        Extractors::_detail::getElementCounts(charges);
-//    ASSERT( !counts_per_charge.empty(),
-//        "getFirstGraphwithSpecifiedElements() - charge counts are empty?");
-    LOG(2, "DEBUG: counts_per_charge is " << counts_per_charge << ".");
+    Extractors::elementcounts_t counts_per_element =
+        Extractors::_detail::getElementCounts(types);
+//    ASSERT( !counts_per_element.empty(),
+//        "getFirstGraphwithSpecifiedElements() - element counts are empty?");
+    LOG(2, "DEBUG: counts_per_element is " << counts_per_element << ".");
 
     // check whether graph contains suitable types
-    Extractors::elementcounts_t::const_iterator countiter = counts_per_charge.begin();
-    for (; countiter != counts_per_charge.end(); ++countiter)
+    Extractors::elementcounts_t::const_iterator countiter = counts_per_element.begin();
+    for (; countiter != counts_per_element.end(); ++countiter)
       if (!graph.hasGreaterEqualTimesAtomicNumber(
           static_cast<size_t>(countiter->first),
@@ -91 +86 @@
         break;
     // if we have a match for every count, store model
-    if( countiter == counts_per_charge.end()) {
+    if( countiter == counts_per_element.end()) {
       LOG(1, "INFO: Potential " << potentialiter->first << " matches with fragment.");
       models.push_back(static_cast<FunctionModel*>(potential));

src/Potentials/PotentialTrainer.cpp

@@ -193 +193 @@
   ASSERT( !types.empty(),
       "getFirstGraphwithSpecifiedElements() - charges is empty?");
-  // create charges
-  Fragment::charges_t charges;
-  charges.resize(types.size());
-  std::transform(types.begin(), types.end(),
-      charges.begin(), boost::lambda::_1);
+
   // convert into count map
-  Extractors::elementcounts_t counts_per_charge =
-      Extractors::_detail::getElementCounts(charges);
-  ASSERT( !counts_per_charge.empty(),
-      "getFirstGraphwithSpecifiedElements() - charge counts are empty?");
-  LOG(1, "DEBUG: counts_per_charge is " << counts_per_charge << ".");
+  Extractors::elementcounts_t counts_per_element =
+      Extractors::_detail::getElementCounts(types);
+  ASSERT( !counts_per_element.empty(),
+      "getFirstGraphwithSpecifiedElements() - element counts are empty?");
+  LOG(1, "DEBUG: counts_per_element is " << counts_per_element << ".");
   // we want to check each (unique) key only once
   HomologyContainer::const_key_iterator olditer = homologies.key_end();
@@ -216 +212 @@
     // check whether we have the same set of atomic numbers
     const HomologyGraph::nodes_t &nodes = (*iter).getNodes();
-    Extractors::elementcounts_t nodes_counts_per_charge;
+    Extractors::elementcounts_t nodes_counts_per_element;
     for (HomologyGraph::nodes_t::const_iterator nodeiter = nodes.begin();
         nodeiter != nodes.end(); ++nodeiter) {
       const Extractors::element_t elem = nodeiter->first.getAtomicNumber();
       const std::pair<Extractors::elementcounts_t::iterator, bool> inserter =
-          nodes_counts_per_charge.insert( std::make_pair(elem, (Extractors::count_t)nodeiter->second ) );
+          nodes_counts_per_element.insert( std::make_pair(elem, (Extractors::count_t)nodeiter->second ) );
       if (!inserter.second)
         inserter.first->second += (Extractors::count_t)nodeiter->second;
     }
-    LOG(1, "DEBUG: Node (" << *iter << ")'s counts_per_charge is " << nodes_counts_per_charge << ".");
-    if (counts_per_charge == nodes_counts_per_charge)
+    LOG(1, "DEBUG: Node (" << *iter << ")'s counts_per_element is " << nodes_counts_per_element << ".");
+    if (counts_per_element == nodes_counts_per_element)
       return *iter;
   }

src/Potentials/SerializablePotential.hpp

@@ -39 +39 @@
 public:
   //!> typedef for particle designation
-  typedef int ParticleType_t;
+  typedef unsigned int ParticleType_t;
   //!> typedef for a vector of particle designations
   typedef std::vector<ParticleType_t> ParticleTypes_t;