Changes in / [b8d1aeb:7ba324]

Location:

src

Files:

• Actions/small_actions.cpp (modified) (1 diff)
• Makefile.am (modified) (1 diff)
• analysis_correlation.cpp (modified) (4 diffs)
• analysis_correlation.hpp (modified) (2 diffs)
• atom_bondedparticle.cpp (modified) (2 diffs)
• atom_bondedparticle.hpp (modified) (1 diff)
• boundary.cpp (modified) (2 diffs)
• builder.cpp (modified) (4 diffs)
• molecule.hpp (modified) (1 diff)
• molecule_graph.cpp (modified) (1 diff)
• unittests/Makefile.am (modified) (1 diff)

src/Actions/small_actions.cpp

@@ -53 +53 @@
   return NAME;
 }
-

src/Makefile.am

@@ -34 +34 @@
 QTUI_DEFS = 
 
-SOURCE = ${ANALYSISSOURCE} ${ATOMSOURCE} ${PATTERNSOURCE} bond.cpp bondgraph.cpp boundary.cpp config.cpp element.cpp ellipsoid.cpp errorlogger.cpp graph.cpp helpers.cpp info.cpp leastsquaremin.cpp linkedcell.cpp log.cpp logger.cpp memoryusageobserver.cpp moleculelist.cpp molecule.cpp molecule_dynamics.cpp molecule_fragmentation.cpp molecule_geometry.cpp molecule_graph.cpp molecule_pointcloud.cpp parser.cpp periodentafel.cpp tesselation.cpp tesselationhelpers.cpp vector.cpp verbose.cpp 
-HEADER = ${ANALYSISHEADER} ${ATOMHEADER} ${PATTERNHEADER} bond.hpp bondgraph.hpp boundary.hpp config.hpp defs.hpp element.hpp ellipsoid.hpp errorlogger.hpp graph.hpp helpers.hpp info.hpp leastsquaremin.hpp linkedcell.hpp lists.hpp log.hpp logger.hpp memoryallocator.hpp memoryusageobserver.hpp molecule.hpp molecule_template.hpp parser.hpp periodentafel.hpp stackclass.hpp tesselation.hpp tesselationhelpers.hpp vector.hpp verbose.hpp
+SOURCE = ${ANALYSISSOURCE} ${ATOMSOURCE} ${PATTERNSOURCE} bond.cpp bondgraph.cpp boundary.cpp config.cpp element.cpp ellipsoid.cpp errorlogger.cpp graph.cpp helpers.cpp info.cpp leastsquaremin.cpp linkedcell.cpp log.cpp logger.cpp memoryusageobserver.cpp moleculelist.cpp molecule.cpp molecule_dynamics.cpp molecule_fragmentation.cpp molecule_geometry.cpp molecule_graph.cpp molecule_pointcloud.cpp parser.cpp periodentafel.cpp tesselation.cpp tesselationhelpers.cpp vector.cpp verbose.cpp menu.cpp
+HEADER = ${ANALYSISHEADER} ${ATOMHEADER} ${PATTERNHEADER} bond.hpp bondgraph.hpp boundary.hpp config.hpp defs.hpp element.hpp ellipsoid.hpp errorlogger.hpp graph.hpp helpers.hpp info.hpp leastsquaremin.hpp linkedcell.hpp lists.hpp log.hpp logger.hpp memoryallocator.hpp memoryusageobserver.hpp molecule.hpp molecule_template.hpp parser.hpp periodentafel.hpp stackclass.hpp tesselation.hpp tesselationhelpers.hpp vector.hpp verbose.hpp menu.hpp
 
 BOOST_LIB = $(BOOST_LDFLAGS) $(BOOST_MPL_LIB)

src/analysis_correlation.cpp

@@ -370 +370 @@
 {
   Info FunctionInfo(__func__);
-  *file << "# BinStart\tCount" << endl;
+  *file << "BinStart\tCount" << endl;
   for (BinPairMap::const_iterator runner = map->begin(); runner != map->end(); ++runner) {
     *file << runner->first << "\t" << runner->second << endl;
@@ -383 +383 @@
 {
   Info FunctionInfo(__func__);
-  *file << "# BinStart\tAtom1\tAtom2" << endl;
+  *file << "BinStart\tAtom1\tAtom2" << endl;
   for (PairCorrelationMap::const_iterator runner = map->begin(); runner != map->end(); ++runner) {
     *file << runner->first << "\t" << *(runner->second.first) << "\t" << *(runner->second.second) << endl;
@@ -396 +396 @@
 {
   Info FunctionInfo(__func__);
-  *file << "# BinStart\tAtom::x[i]-point.x[i]" << endl;
+  *file << "BinStart\tAtom::x[i]-point.x[i]" << endl;
   for (CorrelationToPointMap::const_iterator runner = map->begin(); runner != map->end(); ++runner) {
     *file << runner->first;
@@ -412 +412 @@
 {
   Info FunctionInfo(__func__);
-  *file << "# BinStart\tTriangle" << endl;
+  *file << "BinStart\tTriangle" << endl;
   for (CorrelationToSurfaceMap::const_iterator runner = map->begin(); runner != map->end(); ++runner) {
     *file << runner->first << "\t" << *(runner->second.first) << "\t" << *(runner->second.second) << endl;

src/analysis_correlation.hpp

@@ -90 +90 @@
 /** Puts given correlation data into bins of given size (histogramming).
  * Note that BinStart and BinEnd are desired to fill the complete range, even where Bins are zero. If this is
- * not desired, give equal BinStart and BinEnd and the map will contain only Bins where the count is one or greater.
+ * not desired, give equal BinStart and BinEnd and the map will contain only Bins where the count is one or greater. If a
+ * certain start value is desired, give BinStart and a BinEnd that is smaller than BinStart, then only BinEnd will be
+ * calculated automatically, i.e. BinStart = 0. and BinEnd = -1. .
  * Also note that the range is given inclusive, i.e. [ BinStart, BinEnd ].
  * \param *map map of doubles to count
@@ -114 +116 @@
   if (BinStart == BinEnd) { // if same, find range ourselves
     GetMinMax( map, start, end);
-  } else { // if not, initialise range to zero
+  } else if (BinEnd < BinStart) { // if BinEnd smaller, just look for new max
+    GetMinMax( map, start, end);
+    start = BinStart;
+  } else { // else take both values
     start = BinStart;
     end = BinEnd;

src/atom_bondedparticle.cpp

@@ -56 +56 @@
  * \param *AdjacencyFile output stream
  */
-void BondedParticle::OutputAdjacency(ofstream *AdjacencyFile) const
+void BondedParticle::OutputAdjacency(ofstream * const AdjacencyFile) const
 {
   *AdjacencyFile << nr << "\t";
@@ -62 +62 @@
     *AdjacencyFile << (*Runner)->GetOtherAtom(this)->nr << "\t";
   *AdjacencyFile << endl;
+};
+
+/** Output of atom::nr along each bond partner per line.
+ * Only bonds are printed where atom::nr is smaller than the one of the bond partner.
+ * \param *AdjacencyFile output stream
+ */
+void BondedParticle::OutputBonds(ofstream * const BondFile) const
+{
+  for (BondList::const_iterator Runner = ListOfBonds.begin(); Runner != ListOfBonds.end(); (++Runner))
+    if (nr < (*Runner)->GetOtherAtom(this)->nr)
+      *BondFile << nr << "\t" << (*Runner)->GetOtherAtom(this)->nr << "\n";
 };
 

src/atom_bondedparticle.hpp

@@ -44 +44 @@
   int CorrectBondDegree();
   void OutputBondOfAtom() const;
-  void OutputAdjacency(ofstream *AdjacencyFile) const;
+  void OutputAdjacency(ofstream * const AdjacencyFile) const;
+  void OutputBonds(ofstream * const BondFile) const;
   void OutputOrder(ofstream *file) const;
 

src/boundary.cpp

@@ -19 +19 @@
 
 #include<gsl/gsl_poly.h>
+#include<time.h>
 
 // ========================================== F U N C T I O N S =================================
@@ -840 +841 @@
   Log() << Verbose(1) << "INFO: Grid steps are " << N[0] << ", " << N[1] << ", " << N[2] << "." << endl;
 
+  // initialize seed of random number generator to current time
+  srand ( time(NULL) );
+
   // go over [0,1]^3 filler grid
   for (n[0] = 0; n[0] < N[0]; n[0]++)

src/builder.cpp

@@ -267 +267 @@
         if (i >= 2) {
           first->x.LSQdistance((const Vector **)atoms, i);
-
           first->x.Output();
           first->type = periode->AskElement();  // give type
@@ -1310 +1309 @@
             Log() << Verbose(0) << "\t-h/-H/-?\tGive this help screen." << endl;
             Log() << Verbose(0) << "\t-I\t Dissect current system of molecules into a set of disconnected (subgraphs of) molecules." << endl;
+            Log() << Verbose(0) << "\t-j\t<path> Store all bonds to file." << endl;
+            Log() << Verbose(0) << "\t-J\t<path> Store adjacency per atom to file." << endl;
             Log() << Verbose(0) << "\t-L <step0> <step1> <prefix>\tStore a linear interpolation between two configurations <step0> and <step1> into single config files with prefix <prefix> and as Trajectories into the current config file." << endl;
             Log() << Verbose(0) << "\t-m <0/1>\tCalculate (0)/ Align in(1) PAS with greatest EV along z axis." << endl;
@@ -1596 +1597 @@
                 int ranges[NDIM] = {1,1,1};
                 CorrelationToSurfaceMap *surfacemap = PeriodicCorrelationToSurface( molecules, elemental, TesselStruct, LCList, ranges );
-                //OutputCorrelationToSurface(&output, surfacemap);
-                BinPairMap *binmap = BinData( surfacemap, 0.5, 0., 0. );
+                OutputCorrelationToSurface(&output, surfacemap);
+                BinPairMap *binmap = BinData( surfacemap, 0.5, 0., 20. );
                 OutputCorrelation ( &binoutput, binmap );
                 output.close();
@@ -1681 +1682 @@
               }
               break;
+
+            case 'J':
+              if (ExitFlag == 0) ExitFlag = 1;
+              if ((argptr >= argc) || (argv[argptr][0] == '-')) {
+                ExitFlag =255;
+                eLog() << Verbose(0) << "Missing path of adjacency file: -j <path>" << endl;
+                performCriticalExit();
+              } else {
+                Log() << Verbose(0) << "Storing adjacency to path " << argv[argptr] << "." << endl;
+                configuration.BG->ConstructBondGraph(mol);
+                mol->StoreAdjacencyToFile(argv[argptr]);
+                argptr+=1;
+              }
+              break;
+
+            case 'j':
+              if (ExitFlag == 0) ExitFlag = 1;
+              if ((argptr >= argc) || (argv[argptr][0] == '-')) {
+                ExitFlag =255;
+                eLog() << Verbose(0) << "Missing path of bonds file: -j <path>" << endl;
+                performCriticalExit();
+              } else {
+                Log() << Verbose(0) << "Storing bonds to path " << argv[argptr] << "." << endl;
+                configuration.BG->ConstructBondGraph(mol);
+                mol->StoreBondsToFile(argv[argptr]);
+                argptr+=1;
+              }
+              break;
+
             case 'N':
               if (ExitFlag == 0) ExitFlag = 1;

src/molecule.hpp

@@ -277 +277 @@
   int FragmentMolecule(int Order, config *configuration);
   bool CheckOrderAtSite(bool *AtomMask, Graph *GlobalKeySetList, int Order, int *MinimumRingSize, char *path = NULL);
+  bool StoreBondsToFile(char *path);
   bool StoreAdjacencyToFile(char *path);
   bool CheckAdjacencyFileAgainstMolecule(char *path, atom **ListOfAtoms);

src/molecule_graph.cpp

@@ -992 +992 @@
   Log() << Verbose(1) << "Saving adjacency list ... ";
   if (AdjacencyFile != NULL) {
+    AdjacencyFile << "m\tn" << endl;
     ActOnAllAtoms(&atom::OutputAdjacency, &AdjacencyFile);
     AdjacencyFile.close();
+    Log() << Verbose(1) << "done." << endl;
+  } else {
+    Log() << Verbose(1) << "failed to open file " << line.str() << "." << endl;
+    status = false;
+  }
+
+  return status;
+}
+;
+
+/** Storing the bond structure of a molecule to file.
+ * Simply stores Atom::nr and then the Atom::nr of all bond partners, one per line.
+ * \param *out output stream for debugging
+ * \param *path path to file
+ * \return true - file written successfully, false - writing failed
+ */
+bool molecule::StoreBondsToFile(char *path)
+{
+  ofstream BondFile;
+  stringstream line;
+  bool status = true;
+
+  line << path << "/" << FRAGMENTPREFIX << ADJACENCYFILE;
+  BondFile.open(line.str().c_str(), ios::out);
+  Log() << Verbose(1) << "Saving adjacency list ... ";
+  if (BondFile != NULL) {
+    BondFile << "m\tn" << endl;
+    ActOnAllAtoms(&atom::OutputBonds, &BondFile);
+    BondFile.close();
     Log() << Verbose(1) << "done." << endl;
   } else {

src/unittests/Makefile.am

@@ -30 +30 @@
   ObserverTest \
   ${MENUTESTS}
-    
+
   
 check_PROGRAMS = $(TESTS)