Changeset 35b698 for src/moleculelist.cpp

Timestamp:

Jun 23, 2010, 3:52:50 PM (15 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:

d0fbec

Parents:

dc0d21

git-author:

Frederik Heber <heber@…> (06/23/10 15:20:35)

git-committer:

Frederik Heber <heber@…> (06/23/10 15:52:50)

Message:

BIG CHANGE: config::load and config::save in ParseCommandLineOptions() and main() replaced with FormatParser replacements.

Fragmentation:

• FragmentMolecule() - configuration is not needed anymore, but a path (including prefix)
• class molecule:
  • lots of changes from char * to std:string (due to above): LinearInterpolationBetweenConfiguration(), CheckOrderAtSite(), StoreBondsToFile(), StoreAdjacencyToFile(), CheckAdjacencyFileAgainstMolecule(), CheckAdjacencyFileAgainstMolecule_Init(), ParseOrderAtSiteFromFile(), StoreOrderAtSiteFile(), StoreForcesFile(), ScanAdaptiveFileIntoMap(), StoreForcesFile(), AddHydrogenCorrection(), OutputConfigForListOfFragments()
  • in all above functions, also stringstreams have been replaced by simple strings and file != NULL by file.good()
  • no more CalculateOrbitals(), shifted to class PcpParser
• FragmentationAction: argument of -f is now the prefix, molecule is obtained by molecule-by-id
• no more strcpy needed in: MoleculeLinearInterpolationofTrajectoriesAction, MoleculeSaveAdjacencyAction, MoleculeSaveBondsAction
• oldmenu::FragmentAtoms() - needs to ask for path with prefix now
• config: RetrieveConfigPathAndName() and configpath removed (replaced by above path with prefix)
• changes due to this removal: SaveConfig(), config::Load(), config::LoadOld()

• FIX: MoleculeFillWithMoleculeAction: filler atoms have to be removed before the system can be stored to file.
• FIX: PcpParser::load() - has to put the molecule also into World's MoleculeListClass (otherwise the name cannot be set right after loading)

• ParseCommandLineOptions()
  • new parameter BondGraphFileName and creation of bond graph is done in main()
  • no more config::load, just parser->load(), ConfigFileName is set extra and we always set configPresent to present.
  • FIX: filename is only set when no molecule was allocated before
  • FIX: case 't' and 'f' had wrong argc checks
• main()
  • FormatParserStorage is set with mpqc, pcp and xyz
  • BondGraph is constructed from new string BondGrapghFileName

• new Libparser.a
• all sources from PARSER subdir are compiled into libparser such that only ParserUnitTest is recompiled.

Testfixes:

• testsuite-fragmentation - changes to due to different -f calling syntax.
• most of the xyz files had to be replaced due to a single whitespace at the end of each entry: Domain/6, Simple_configuration/2, Simple_configuration/3, Simple_configuration/4, Simple_configuration/5, Simple_configuration/8
• in many cases were the number orbitals (and thus MaxMinStopStep) wrong: Filling/1, Simple_configuration/4, Simple_configuration/5

Signed-off-by: Frederik Heber <heber@…>

File:

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

src/moleculelist.cpp

@@ -380 +380 @@
  * bonded to the same atom, then we add for this pair a correction term constructed from a Morse
  * potential function fit to QM calculations with respecting to the interatomic hydrogen distance.
- * \param *out output stream for debugging
- * \param *path path to file
- */
-bool MoleculeListClass::AddHydrogenCorrection(char *path)
+ * \param &path path to file
+ */
+bool MoleculeListClass::AddHydrogenCorrection(std::string &path)
 {
   bond *Binder = NULL;
@@ -401 +400 @@
   // 0a. find dimension of matrices with constants
   line = path;
-  line.append("/");
-  line += FRAGMENTPREFIX;
   line += "1";
   line += FITCONSTANTSUFFIX;
   input.open(line.c_str());
-  if (input == NULL) {
+  if (input.fail()) {
     DoLog(1) && (Log() << Verbose(1) << endl << "Unable to open " << line << ", is the directory correct?" << endl);
     return false;
@@ -570 +567 @@
 
 /** Store force indices, i.e. the connection between the nuclear index in the total molecule config and the respective atom in fragment config.
- * \param *out output stream for debugging
- * \param *path path to file
+ * \param &path path to file
  * \param *SortIndex Index to map from the BFS labeling to the sequence how of Ion_Type in the config
  * \return true - file written successfully, false - writing failed
  */
-bool MoleculeListClass::StoreForcesFile(char *path,
-    int *SortIndex)
+bool MoleculeListClass::StoreForcesFile(std::string &path, int *SortIndex)
 {
   bool status = true;
-  ofstream ForcesFile;
-  stringstream line;
+  string filename(path);
+  filename += FORCESFILE;
+  ofstream ForcesFile(filename.c_str());
   periodentafel *periode=World::getInstance().getPeriode();
 
   // open file for the force factors
   DoLog(1) && (Log() << Verbose(1) << "Saving  force factors ... ");
-  line << path << "/" << FRAGMENTPREFIX << FORCESFILE;
-  ForcesFile.open(line.str().c_str(), ios::out);
-  if (ForcesFile != NULL) {
+  if (!ForcesFile.fail()) {
     //Log() << Verbose(1) << "Final AtomicForcesList: ";
     //output << prefix << "Forces" << endl;
@@ -612 +606 @@
   } else {
     status = false;
-    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << line.str() << "." << endl);
+    DoLog(1) && (Log() << Verbose(1) << "failed to open file " << filename << "." << endl);
   }
   ForcesFile.close();
@@ -621 +615 @@
 /** Writes a config file for each molecule in the given \a **FragmentList.
  * \param *out output stream for debugging
- * \param *configuration standard configuration to attach atoms in fragment molecule to.
+ * \param &prefix path and prefix to the fragment config files
  * \param *SortIndex Index to map from the BFS labeling to the sequence how of Ion_Type in the config
  * \return true - success (each file was written), false - something went wrong.
  */
-bool MoleculeListClass::OutputConfigForListOfFragments(config *configuration, int *SortIndex)
+bool MoleculeListClass::OutputConfigForListOfFragments(std::string &prefix, int *SortIndex)
 {
   ofstream outputFragment;
-  char FragmentName[MAXSTRINGSIZE];
+  std::string FragmentName;
   char PathBackup[MAXSTRINGSIZE];
   bool result = true;
@@ -646 +640 @@
   for (MoleculeList::iterator ListRunner = ListOfMolecules.begin(); ListRunner != ListOfMolecules.end(); ListRunner++) {
     // save default path as it is changed for each fragment
-    path = configuration->GetDefaultPath();
+    path = World::getInstance().getConfig()->GetDefaultPath();
     if (path != NULL)
       strcpy(PathBackup, path);
@@ -659 +653 @@
     // output xyz file
     FragmentNumber = FixedDigitNumber(ListOfMolecules.size(), FragmentCounter++);
-    sprintf(FragmentName, "%s/%s%s.conf.xyz", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
-    outputFragment.open(FragmentName, ios::out);
+    FragmentName = prefix + FragmentNumber + ".conf.xyz";
+    outputFragment.open(FragmentName.c_str(), ios::out);
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as XYZ ...");
     if ((intermediateResult = (*ListRunner)->OutputXYZ(&outputFragment)))
@@ -683 +677 @@
     for (int k = 0; k < NDIM; k++) {
       j += k + 1;
-      BoxDimension[k] = 2.5 * (configuration->GetIsAngstroem() ? 1. : 1. / AtomicLengthToAngstroem);
+      BoxDimension[k] = 2.5 * (World::getInstance().getConfig()->GetIsAngstroem() ? 1. : 1. / AtomicLengthToAngstroem);
       cell_size[j] = BoxDimension[k] * 2.;
     }
@@ -690 +684 @@
     // also calculate necessary orbitals
     (*ListRunner)->CountElements(); // this is a bugfix, atoms should shoulds actually be added correctly to this fragment
-    (*ListRunner)->CalculateOrbitals(*configuration);
+    //(*ListRunner)->CalculateOrbitals(*World::getInstance().getConfig);
 
     // change path in config
-    //strcpy(PathBackup, configuration->configpath);
-    sprintf(FragmentName, "%s/%s%s/", PathBackup, FRAGMENTPREFIX, FragmentNumber);
-    configuration->SetDefaultPath(FragmentName);
+    FragmentName = PathBackup;
+    FragmentName += "/";
+    FragmentName += FRAGMENTPREFIX;
+    FragmentName += FragmentNumber;
+    FragmentName += "/";
+    World::getInstance().getConfig()->SetDefaultPath(FragmentName.c_str());
 
     // and save as config
-    sprintf(FragmentName, "%s/%s%s.conf", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
+    FragmentName = prefix + FragmentNumber + ".conf";
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as config ...");
-    if ((intermediateResult = configuration->Save(FragmentName, (*ListRunner)->elemente, (*ListRunner))))
+    if ((intermediateResult = World::getInstance().getConfig()->Save(FragmentName.c_str(), (*ListRunner)->elemente, (*ListRunner))))
       DoLog(0) && (Log() << Verbose(0) << " done." << endl);
     else
@@ -707 +704 @@
 
     // restore old config
-    configuration->SetDefaultPath(PathBackup);
+    World::getInstance().getConfig()->SetDefaultPath(PathBackup);
 
     // and save as mpqc input file
-    sprintf(FragmentName, "%s/%s%s.conf", configuration->configpath, FRAGMENTPREFIX, FragmentNumber);
+    FragmentName = prefix + FragmentNumber + ".conf";
     DoLog(2) && (Log() << Verbose(2) << "Saving bond fragment No. " << FragmentNumber << "/" << FragmentCounter - 1 << " as mpqc input ...");
-    if ((intermediateResult = configuration->SaveMPQC(FragmentName, (*ListRunner))))
+    if ((intermediateResult = World::getInstance().getConfig()->SaveMPQC(FragmentName.c_str(), (*ListRunner))))
       DoLog(2) && (Log() << Verbose(2) << " done." << endl);
     else
@@ -768 +765 @@
 
   // 1. dissect the molecule into connected subgraphs
-  if (!configuration->BG->ConstructBondGraph(mol)) {
-    World::getInstance().destroyMolecule(mol);
-    DoeLog(1) && (eLog()<< Verbose(1) << "There are no bonds." << endl);
+  if (configuration->BG != NULL) {
+    if (!configuration->BG->ConstructBondGraph(mol)) {
+      World::getInstance().destroyMolecule(mol);
+      DoeLog(1) && (eLog()<< Verbose(1) << "There are no bonds." << endl);
+      return;
+    }
+  } else {
+    DoeLog(1) && (eLog()<< Verbose(1) << "There is no BondGraph class present to create bonds." << endl);
     return;
   }