| 1 | /*
|
|---|
| 2 | * Project: MoleCuilder
|
|---|
| 3 | * Description: creates and alters molecular systems
|
|---|
| 4 | * Copyright (C) 2010 University of Bonn. All rights reserved.
|
|---|
| 5 | * Please see the LICENSE file or "Copyright notice" in builder.cpp for details.
|
|---|
| 6 | */
|
|---|
| 7 |
|
|---|
| 8 | /** \file builder.cpp
|
|---|
| 9 | *
|
|---|
| 10 | * date: Jan 1, 2007
|
|---|
| 11 | * author: heber
|
|---|
| 12 | *
|
|---|
| 13 | */
|
|---|
| 14 |
|
|---|
| 15 | /*! \page Copyright notice
|
|---|
| 16 | *
|
|---|
| 17 | * MoleCuilder - to create and alter molecular systems
|
|---|
| 18 | * Copyright (C) 2010 University Bonn. All rights reserved.
|
|---|
| 19 | *
|
|---|
| 20 | */
|
|---|
| 21 |
|
|---|
| 22 | /*! \mainpage MoleCuilder - a molecular set builder
|
|---|
| 23 | *
|
|---|
| 24 | * This introductory shall briefly make acquainted with the program, helping in installing and a first run.
|
|---|
| 25 | *
|
|---|
| 26 | * \section about About the Program
|
|---|
| 27 | *
|
|---|
| 28 | * MoleCuilder is a program, written entirely in C++, that enables the construction of a coordinate set for the
|
|---|
| 29 | * atoms making up an molecule. It allows for both building of simple molecules by adding atom-wise giving bond
|
|---|
| 30 | * angles and distances or absolute coordinates, but also using them as templates. Regions can be specified and
|
|---|
| 31 | * ordered to be filled with a molecule in a certain manner. Greater conglomerations of molecules can be tesselated
|
|---|
| 32 | * and recognized as a region themselves to be subsequently surrounded by other (surface solvated) molecules.
|
|---|
| 33 | * In the end, MoleCuilder allows the construction of arbitrary nano structures, whether they be crystalline or
|
|---|
| 34 | * amorphic in nature.
|
|---|
| 35 | *
|
|---|
| 36 | *
|
|---|
| 37 | * \section install Installation
|
|---|
| 38 | *
|
|---|
| 39 | * Installation should without problems succeed as follows:
|
|---|
| 40 | * -# ./configure (or: mkdir build;mkdir run;cd build; ../configure --bindir=../run)
|
|---|
| 41 | * -# make
|
|---|
| 42 | * -# make install
|
|---|
| 43 | *
|
|---|
| 44 | * Further useful commands are
|
|---|
| 45 | * -# make clean uninstall: deletes .o-files and removes executable from the given binary directory\n
|
|---|
| 46 | * -# make doxygen-doc: Creates these html pages out of the documented source
|
|---|
| 47 | * -# make check: Runs an extensive set of unit tests and a testsuite which also gives a good overview on the set of
|
|---|
| 48 | * functions.
|
|---|
| 49 | *
|
|---|
| 50 | * \section run Running
|
|---|
| 51 | *
|
|---|
| 52 | * The program can be executed by running: ./molecuilder
|
|---|
| 53 | *
|
|---|
| 54 | * MoleCuilder has three interfaces at your disposal:
|
|---|
| 55 | * -# Textmenu: A simple interactive console-based menu, where awaits your choices and inputs in order to set atoms
|
|---|
| 56 | * as you like
|
|---|
| 57 | * -# CommandLineUI: Every command can also be chained up as a sequence of actions on the command line to be executed
|
|---|
| 58 | * with any user interaction.
|
|---|
| 59 | * -# GraphicalUI: A graphical user interface that also display the molecular structure being built and lots of other
|
|---|
| 60 | * informations to ease the construction of bigger geometries.
|
|---|
| 61 | *
|
|---|
| 62 | * The supported output formats right now are:
|
|---|
| 63 | * -# mpqc: Configuration files of the Massively Parallel Quantum Chemistry package (Sandia labs)
|
|---|
| 64 | * -# pcp: Configuration files of the Parallel Car-Parrinello program (Institute for Numerical Simulation)
|
|---|
| 65 | * -# tremolo: Configuration files of TREMOLO (Institute for Numerical Simulation)
|
|---|
| 66 | * -# xyz: the most basic format for the 3d arrangement of atoms consisting of a list of element and 3 coordinates.
|
|---|
| 67 | *
|
|---|
| 68 | */
|
|---|
| 69 |
|
|---|
| 70 | // include config.h
|
|---|
| 71 | #ifdef HAVE_CONFIG_H
|
|---|
| 72 | #include <config.h>
|
|---|
| 73 | #endif
|
|---|
| 74 |
|
|---|
| 75 | #include "Helpers/MemDebug.hpp"
|
|---|
| 76 |
|
|---|
| 77 | #include "bondgraph.hpp"
|
|---|
| 78 | #include "config.hpp"
|
|---|
| 79 | #include "Helpers/Log.hpp"
|
|---|
| 80 | #include "molecule.hpp"
|
|---|
| 81 | #include "periodentafel.hpp"
|
|---|
| 82 | #include "tesselationhelpers.hpp"
|
|---|
| 83 | #include "UIElements/UIFactory.hpp"
|
|---|
| 84 | #include "UIElements/Menu/MenuDescription.hpp"
|
|---|
| 85 | #include "UIElements/TextUI/TextUIFactory.hpp"
|
|---|
| 86 | #include "UIElements/CommandLineUI/CommandLineUIFactory.hpp"
|
|---|
| 87 | #include "UIElements/CommandLineUI/CommandLineParser.hpp"
|
|---|
| 88 | #ifdef USE_GUI_QT
|
|---|
| 89 | #include "UIElements/Qt4/QtUIFactory.hpp"
|
|---|
| 90 | #endif
|
|---|
| 91 | #include "UIElements/MainWindow.hpp"
|
|---|
| 92 | #include "UIElements/Dialog.hpp"
|
|---|
| 93 | //#include "Menu/ActionMenuItem.hpp"
|
|---|
| 94 | #include "Helpers/Verbose.hpp"
|
|---|
| 95 | #include "World.hpp"
|
|---|
| 96 |
|
|---|
| 97 | #include "Actions/ActionRegistry.hpp"
|
|---|
| 98 | #include "Actions/ActionHistory.hpp"
|
|---|
| 99 |
|
|---|
| 100 | #include "Parser/ChangeTracker.hpp"
|
|---|
| 101 | #include "Parser/FormatParserStorage.hpp"
|
|---|
| 102 |
|
|---|
| 103 | #include "UIElements/UIFactory.hpp"
|
|---|
| 104 | #include "UIElements/TextUI/TextUIFactory.hpp"
|
|---|
| 105 | #include "UIElements/CommandLineUI/CommandLineUIFactory.hpp"
|
|---|
| 106 | #include "UIElements/MainWindow.hpp"
|
|---|
| 107 | #include "UIElements/Dialog.hpp"
|
|---|
| 108 |
|
|---|
| 109 | #include "version.h"
|
|---|
| 110 |
|
|---|
| 111 |
|
|---|
| 112 | /********************************************** Main routine **************************************/
|
|---|
| 113 |
|
|---|
| 114 | /** In this function all dynamicly allocated member variables to static/global
|
|---|
| 115 | * variables are added to the ignore list of Memory/MemDebug.
|
|---|
| 116 | *
|
|---|
| 117 | * Use this to prevent their listing in the Memory::getState() at the end of the
|
|---|
| 118 | * program. Check with valgrind that truely no memory leak occurs!
|
|---|
| 119 | */
|
|---|
| 120 | void AddStaticEntitiestoIgnoreList()
|
|---|
| 121 | {
|
|---|
| 122 | // zeroVec and unitVec are global variables (on the stack) but vectorContent
|
|---|
| 123 | // within is situated on the heap and has to be ignored
|
|---|
| 124 | Memory::ignore(zeroVec.get());
|
|---|
| 125 | Memory::ignore(unitVec[0].get());
|
|---|
| 126 | Memory::ignore(unitVec[1].get());
|
|---|
| 127 | Memory::ignore(unitVec[2].get());
|
|---|
| 128 | }
|
|---|
| 129 |
|
|---|
| 130 | /** Cleans all singleton instances in an orderly fashion.
|
|---|
| 131 | * C++ does not guarantee any specific sequence of removal of single instances
|
|---|
| 132 | * which have static/global variables. Some singletons depend on others hence we
|
|---|
| 133 | * acertain a specific ordering here, which is is used via the atexit() hook.
|
|---|
| 134 | */
|
|---|
| 135 | void cleanUp(){
|
|---|
| 136 | FormatParserStorage::purgeInstance();
|
|---|
| 137 | ChangeTracker::purgeInstance();
|
|---|
| 138 | World::purgeInstance();
|
|---|
| 139 | MenuDescription::purgeInstance();
|
|---|
| 140 | UIFactory::purgeInstance();
|
|---|
| 141 | ValueStorage::purgeInstance();
|
|---|
| 142 | CommandLineParser::purgeInstance();
|
|---|
| 143 | ActionRegistry::purgeInstance();
|
|---|
| 144 | OptionRegistry::purgeInstance();
|
|---|
| 145 | ActionHistory::purgeInstance();
|
|---|
| 146 | // we have to remove these two static as otherwise their boost::shared_ptrs are still present
|
|---|
| 147 | Action::removeStaticStateEntities();
|
|---|
| 148 | // put some static variables' dynamic contents on the Memory::ignore map to avoid their
|
|---|
| 149 | // admonishing lateron
|
|---|
| 150 | AddStaticEntitiestoIgnoreList();
|
|---|
| 151 | logger::purgeInstance();
|
|---|
| 152 | errorLogger::purgeInstance();
|
|---|
| 153 | #ifdef LOG_OBSERVER
|
|---|
| 154 | cout << observerLog().getLog();
|
|---|
| 155 | #endif
|
|---|
| 156 | Memory::getState();
|
|---|
| 157 | }
|
|---|
| 158 |
|
|---|
| 159 | void dumpMemory(){
|
|---|
| 160 | ofstream ost("molecuilder.memdump");
|
|---|
| 161 | Memory::dumpMemory(ost);
|
|---|
| 162 | }
|
|---|
| 163 |
|
|---|
| 164 | int main(int argc, char **argv)
|
|---|
| 165 | {
|
|---|
| 166 | // while we are non interactive, we want to abort from asserts
|
|---|
| 167 | ASSERT_DO(Assert::Abort);
|
|---|
| 168 | ASSERT_HOOK(dumpMemory);
|
|---|
| 169 | string line;
|
|---|
| 170 | char **Arguments = NULL;
|
|---|
| 171 | int ArgcSize = 0;
|
|---|
| 172 | int ExitFlag = 0;
|
|---|
| 173 | bool ArgumentsCopied = false;
|
|---|
| 174 | std::string BondGraphFileName("\n");
|
|---|
| 175 |
|
|---|
| 176 | // print version check and copyright notice
|
|---|
| 177 | cout << MOLECUILDERVERSION << endl;
|
|---|
| 178 | cout << "MoleCuilder comes with ABSOLUTELY NO WARRANTY; for details type" << endl;
|
|---|
| 179 | cout << "`molecuilder --warranty'." << endl;
|
|---|
| 180 | cout << "`MoleCuilder - to create and alter molecular systems." << endl;
|
|---|
| 181 | cout << "Copyright (C) 2010 University Bonn. All rights reserved." << endl;
|
|---|
| 182 |
|
|---|
| 183 | setVerbosity(0);
|
|---|
| 184 | // need to init the history before any action is created
|
|---|
| 185 | ActionHistory::init();
|
|---|
| 186 |
|
|---|
| 187 | // from this moment on, we need to be sure to deeinitialize in the correct order
|
|---|
| 188 | // this is handled by the cleanup function
|
|---|
| 189 | atexit(cleanUp);
|
|---|
| 190 |
|
|---|
| 191 | // Parse command line options and if present create respective UI
|
|---|
| 192 | {
|
|---|
| 193 | // construct bond graph
|
|---|
| 194 | if (World::getInstance().getConfig()->BG == NULL) {
|
|---|
| 195 | World::getInstance().getConfig()->BG = new BondGraph(World::getInstance().getConfig()->GetIsAngstroem());
|
|---|
| 196 | if (World::getInstance().getConfig()->BG->LoadBondLengthTable(BondGraphFileName)) {
|
|---|
| 197 | DoLog(0) && (Log() << Verbose(0) << "Bond length table loaded successfully." << endl);
|
|---|
| 198 | } else {
|
|---|
| 199 | DoeLog(1) && (eLog()<< Verbose(1) << "Bond length table loading failed." << endl);
|
|---|
| 200 | }
|
|---|
| 201 | }
|
|---|
| 202 | // handle remaining arguments by CommandLineParser
|
|---|
| 203 | if (argc>1) {
|
|---|
| 204 | CommandLineParser::getInstance().InitializeCommandArguments();
|
|---|
| 205 | CommandLineParser::getInstance().Run(argc,argv);
|
|---|
| 206 | DoLog(0) && (Log() << Verbose(0) << "Setting UI to CommandLine." << endl);
|
|---|
| 207 | UIFactory::registerFactory(new CommandLineUIFactory::description());
|
|---|
| 208 | UIFactory::makeUserInterface("CommandLine");
|
|---|
| 209 | } else {
|
|---|
| 210 | // In the interactive mode, we can leave the user the choice in case of error
|
|---|
| 211 | ASSERT_DO(Assert::Ask);
|
|---|
| 212 | #ifdef USE_GUI_QT
|
|---|
| 213 | DoLog(0) && (Log() << Verbose(0) << "Setting UI to Qt4." << endl);
|
|---|
| 214 | UIFactory::registerFactory(new QtUIFactory::description());
|
|---|
| 215 | UIFactory::makeUserInterface("Qt4");
|
|---|
| 216 | #else
|
|---|
| 217 | DoLog(0) && (Log() << Verbose(0) << "Setting UI to Text." << endl);
|
|---|
| 218 | cout << MOLECUILDERVERSION << endl;
|
|---|
| 219 | UIFactory::registerFactory(new TextUIFactory::description());
|
|---|
| 220 | UIFactory::makeUserInterface("Text");
|
|---|
| 221 | #endif
|
|---|
| 222 | }
|
|---|
| 223 | }
|
|---|
| 224 |
|
|---|
| 225 | {
|
|---|
| 226 | MainWindow *mainWindow = UIFactory::getInstance().makeMainWindow();
|
|---|
| 227 | mainWindow->display();
|
|---|
| 228 | delete mainWindow;
|
|---|
| 229 | }
|
|---|
| 230 |
|
|---|
| 231 | FormatParserStorage::getInstance().SaveAll();
|
|---|
| 232 | ChangeTracker::getInstance().saveStatus();
|
|---|
| 233 |
|
|---|
| 234 | // free the new argv
|
|---|
| 235 | if (ArgumentsCopied) {
|
|---|
| 236 | for (int i=0; i<ArgcSize;i++)
|
|---|
| 237 | delete[](Arguments[i]);
|
|---|
| 238 | delete[](Arguments);
|
|---|
| 239 | }
|
|---|
| 240 | //delete[](ConfigFileName);
|
|---|
| 241 |
|
|---|
| 242 | ExitFlag = World::getInstance().getExitFlag();
|
|---|
| 243 | return (ExitFlag == 1 ? 0 : ExitFlag);
|
|---|
| 244 | }
|
|---|
| 245 |
|
|---|
| 246 | /********************************************** E N D **************************************************/
|
|---|