Changeset 8d5db8 for src/Fragmentation/Summation/SetValues/Fragment.hpp

Timestamp:

Oct 2, 2016, 1:49:57 PM (9 years ago)

Author:

Frederik Heber <heber@…>

Branches:

Fix_FitPotential_needs_atomicnumbers

Children:

a13e21

Parents:

f6f761

git-author:

Frederik Heber <heber@…> (10/02/16 13:24:41)

git-committer:

Frederik Heber <heber@…> (10/02/16 13:49:57)

Message:

MPQCData additionally stores the atomic number per nuclei.

• in case of DoSampleValenceOnly the charge does not represent the atomic number and moreover it is no longer unique. However, we need this unique association for fitting potentials to the fragment results.
• TESTFIX: needed to adapt HomologyContainerUnitTest which uses FragmentStub.
• TESTFIX: needed to adapt FragmentUnitTest.

File:

• src/Fragmentation/Summation/SetValues/Fragment.hpp (modified) (4 diffs)

src/Fragmentation/Summation/SetValues/Fragment.hpp

@@ -34 +34 @@
   typedef double charge_t;
   typedef std::vector< charge_t > charges_t;
-  typedef std::pair< position_t, charge_t> nucleus_t;
+  typedef unsigned int atomicNumber_t; // typedef is copied
+  typedef std::vector< atomicNumber_t > atomicnumbers_t;
+  typedef std::pair< position_t, std::pair<atomicNumber_t, charge_t> > nucleus_t;
   typedef std::vector< nucleus_t > nuclei_t;
 
@@ -52 +54 @@
    *
    * @param _positions given positions
+   * @param _atomicnumbers given atomic numbers
    * @param _charges given charges
    */
-  Fragment(const positions_t &_positions, const charges_t &_charges);
+  Fragment(
+      const positions_t &_positions,
+      const atomicnumbers_t &_atomicnumbers,
+      const charges_t &_charges);
 
   /** Adding another fragment onto this one.
@@ -86 +92 @@
   positions_t getPositions() const;
 
+  /** Getter for all stored atomic numbers.
+   *
+   * @return vector of atomic numbers
+   */
+  atomicnumbers_t getAtomicNumbers() const;
+
   /** Getter for all stored charges.
    *
@@ -107 +119 @@
    *
    * @param position position of nucleus to create
+   * @param atomicnumber atomic number of nucleus to create
    * @param charge charge of nucleus to create
    * @return nucleus with given \a position and \a charge
    */
-  static nucleus_t createNucleus(const position_t &position, const double charge);
+  static nucleus_t createNucleus(
+      const position_t &position,
+      const atomicNumber_t &atomicnumber,
+      const double charge);
 
   /** Helper function to check whether two positions are equal.