Changeset f23d7d

Timestamp:

Feb 9, 2009, 11:17:38 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

2218dca

Parents:

97c751

Message:

Basis for MPQC can now be specified via command line switch

• -B is ParseCommandLineOptions added
• configuration::basis added, set to 3-21G in constructor
• used in configuration::SaveMPQC()

Location:

molecuilder/src

Files:

• builder.cpp (modified) (3 diffs)
• config.cpp (modified) (2 diffs)
• moleculelist.cpp (modified) (2 diffs)
• molecules.hpp (modified) (1 diff)

molecuilder/src/builder.cpp

@@ -697 +697 @@
   char filename[MAXSTRINGSIZE];
   ofstream output;
+  string basis("3-21G");
 
   cout << Verbose(0) << "Storing configuration ... " << endl;
@@ -806 +807 @@
             cout << "\t-A <source>\tCreate adjacency list from bonds parsed from 'dbond'-style file." <<endl;
             cout << "\t-b x1 x2 x3\tCenter atoms in domain with given edge lengths of (x1,x2,x3)." << endl;
+            cout << "\t-B <basis>\tgive gaussian basis for MPQC output." << endl;
             cout << "\t-c x1 x2 x3\tCenter atoms in domain with a minimum distance to boundary of (x1,x2,x3)." << endl;
             cout << "\t-D <bond distance>\tDepth-First-Search Analysis of the molecule, giving cycles and tree/back edges." << endl;
@@ -961 +963 @@
         if (config_present == present) {
           switch(argv[argptr-1][1]) {
+            case 'B':
+              if ((argptr >= argc) || (!IsValidNumber(argv[argptr])) || (argv[argptr][0] == '-')) {
+                ExitFlag = 255;
+                cerr << "Not enough or invalid arguments given for giving gaussian basis: -B <basis>" << endl;
+              } else {
+                cout << Verbose(1) << "Setting gaussian basis to " << argv[argptr] << "." << endl;
+                configuration.basis.copy(argv[argptr],strlen(argv[argptr]));
+                argptr+=1;
+              }
+              break;
             case 'D':
               ExitFlag = 1;

molecuilder/src/config.cpp

@@ -23 +23 @@
   configpath[0]='\0';
   configname[0]='\0';
+  basis="3-21G";
   
   FastParsing = false;
@@ -1131 +1132 @@
   *output << ")" << endl;
   *output << "basis<GaussianBasisSet>: (" << endl;
-  *output << "\tname = \"3-21G\"" << endl;
+  *output << "\tname = \"" << basis << "\"" << endl;
   *output << "\tmolecule = $:molecule" << endl;
   *output << ")" << endl;

molecuilder/src/moleculelist.cpp

@@ -393 +393 @@
   int FragmentCounter = 0;
   ofstream output;
+  string basis("3-21G");
   
   // store the fragments as config and as xyz
@@ -404 +405 @@
 
     // correct periodic
-    ListOfMolecules[i]->ScanForPeriodicCorrection(out);
+    //ListOfMolecules[i]->ScanForPeriodicCorrection(out);
 
     // output xyz file

molecuilder/src/molecules.hpp

@@ -387 +387 @@
     bool FastParsing;
     double Deltat;
+    string basis;
 
     private: