Changeset f23d7d for molecuilder/src/builder.cpp

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()

File:

• molecuilder/src/builder.cpp (modified) (3 diffs)

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;