Changes in src/molecule.cpp [35b698:906822]

File:

• src/molecule.cpp (modified) (5 diffs)

src/molecule.cpp

@@ -79 +79 @@
 void molecule::setName(const std::string _name){
   OBSERVE;
-  cout << "Set name of molecule " << getId() << " to " << _name << endl;
   strncpy(name,_name.c_str(),MAXSTRINGSIZE);
 }
@@ -155 +154 @@
 molecule::const_iterator molecule::erase( atom * key )
 {
-  cout << "trying to erase atom" << endl;
   molecule::const_iterator iter = find(key);
   if (iter != end()){
@@ -741 +739 @@
   else
     length = strlen(molname) - strlen(endname);
-  cout << "Set name of molecule " << getId() << " to " << molname << endl;
   strncpy(name, molname, length);
   name[length]='\0';
@@ -883 +880 @@
         ElementNo[i] = current++;
     }
-    ActOnAllAtoms( &atom::OutputArrayIndexed, (ostream * const) output, (const int *)ElementNo, (int *)AtomNo, (const char *) NULL );
+    ActOnAllAtoms( &atom::OutputArrayIndexed, output, (const int *)ElementNo, (int *)AtomNo, (const char *) NULL );
     return true;
   }
@@ -1006 +1003 @@
   for(int i=MAX_ELEMENTS;i--;)
     ElementCount += (ElementsInMolecule[i] != 0 ? 1 : 0);
+};
+
+
+/** Counts necessary number of valence electrons and returns number and SpinType.
+ * \param configuration containing everything
+ */
+void molecule::CalculateOrbitals(class config &configuration)
+{
+  configuration.MaxPsiDouble = configuration.PsiMaxNoDown = configuration.PsiMaxNoUp = configuration.PsiType = 0;
+  for(int i=MAX_ELEMENTS;i--;) {
+    if (ElementsInMolecule[i] != 0) {
+      //Log() << Verbose(0) << "CalculateOrbitals: " << elemente->FindElement(i)->name << " has a valence of " << (int)elemente->FindElement(i)->Valence << " and there are " << ElementsInMolecule[i] << " of it." << endl;
+      configuration.MaxPsiDouble += ElementsInMolecule[i]*((int)elemente->FindElement(i)->Valence);
+    }
+  }
+  configuration.PsiMaxNoDown = configuration.MaxPsiDouble/2 + (configuration.MaxPsiDouble % 2);
+  configuration.PsiMaxNoUp = configuration.MaxPsiDouble/2;
+  configuration.MaxPsiDouble /= 2;
+  configuration.PsiType = (configuration.PsiMaxNoDown == configuration.PsiMaxNoUp) ? 0 : 1;
+  if ((configuration.PsiType == 1) && (configuration.ProcPEPsi < 2) && ((configuration.PsiMaxNoDown != 1) || (configuration.PsiMaxNoUp != 0))) {
+    configuration.ProcPEGamma /= 2;
+    configuration.ProcPEPsi *= 2;
+  } else {
+    configuration.ProcPEGamma *= configuration.ProcPEPsi;
+    configuration.ProcPEPsi = 1;
+  }
+  cout << configuration.PsiMaxNoDown << ">" << configuration.PsiMaxNoUp << endl;
+  if (configuration.PsiMaxNoDown > configuration.PsiMaxNoUp) {
+    configuration.InitMaxMinStopStep = configuration.MaxMinStopStep = configuration.PsiMaxNoDown;
+    cout << configuration.PsiMaxNoDown << " " << configuration.InitMaxMinStopStep << endl;
+  } else {
+    configuration.InitMaxMinStopStep = configuration.MaxMinStopStep = configuration.PsiMaxNoUp;
+    cout << configuration.PsiMaxNoUp << " " << configuration.InitMaxMinStopStep << endl;
+  }
 };