| 1 | /*
 | 
|---|
| 2 |  * XyzParser.cpp
 | 
|---|
| 3 |  *
 | 
|---|
| 4 |  *  Created on: Mar 2, 2010
 | 
|---|
| 5 |  *      Author: metzler
 | 
|---|
| 6 |  */
 | 
|---|
| 7 | 
 | 
|---|
| 8 | #include "Helpers/MemDebug.hpp"
 | 
|---|
| 9 | 
 | 
|---|
| 10 | #include "Helpers/Log.hpp"
 | 
|---|
| 11 | #include "Helpers/Verbose.hpp"
 | 
|---|
| 12 | #include "XyzParser.hpp"
 | 
|---|
| 13 | #include "World.hpp"
 | 
|---|
| 14 | #include "atom.hpp"
 | 
|---|
| 15 | #include "element.hpp"
 | 
|---|
| 16 | #include "periodentafel.hpp"
 | 
|---|
| 17 | 
 | 
|---|
| 18 | using namespace std;
 | 
|---|
| 19 | 
 | 
|---|
| 20 | /**
 | 
|---|
| 21 |  * Constructor.
 | 
|---|
| 22 |  */
 | 
|---|
| 23 | XyzParser::XyzParser() {
 | 
|---|
| 24 |   comment = "";
 | 
|---|
| 25 | }
 | 
|---|
| 26 | 
 | 
|---|
| 27 | /**
 | 
|---|
| 28 |  * Destructor.
 | 
|---|
| 29 |  */
 | 
|---|
| 30 | XyzParser::~XyzParser() {
 | 
|---|
| 31 | }
 | 
|---|
| 32 | 
 | 
|---|
| 33 | /**
 | 
|---|
| 34 |  * Loads an XYZ file into the World.
 | 
|---|
| 35 |  *
 | 
|---|
| 36 |  * \param XYZ file
 | 
|---|
| 37 |  */
 | 
|---|
| 38 | void XyzParser::load(istream* file) {
 | 
|---|
| 39 |   atom* newAtom;
 | 
|---|
| 40 |   int numberOfAtoms;
 | 
|---|
| 41 |   char commentBuffer[512], type[3];
 | 
|---|
| 42 | 
 | 
|---|
| 43 |   // the first line tells number of atoms, the second line is always a comment
 | 
|---|
| 44 |   *file >> numberOfAtoms >> ws;
 | 
|---|
| 45 |   file->getline(commentBuffer, 512);
 | 
|---|
| 46 |   comment = commentBuffer;
 | 
|---|
| 47 | 
 | 
|---|
| 48 |   for (int i = 0; i < numberOfAtoms; i++) {
 | 
|---|
| 49 |     newAtom = World::getInstance().createAtom();
 | 
|---|
| 50 |     *file >> type >> ws >> newAtom->x[0] >> ws >> newAtom->x[1] >> ws >> newAtom->x[2];
 | 
|---|
| 51 |     newAtom->setType(World::getInstance().getPeriode()->FindElement(type));
 | 
|---|
| 52 |   }
 | 
|---|
| 53 | }
 | 
|---|
| 54 | 
 | 
|---|
| 55 | /**
 | 
|---|
| 56 |  * Saves the current state of the World into the given XYZ file.
 | 
|---|
| 57 |  *
 | 
|---|
| 58 |  * \param XYZ file
 | 
|---|
| 59 |  */
 | 
|---|
| 60 | void XyzParser::save(ostream* file) {
 | 
|---|
| 61 |   DoLog(0) && (Log() << Verbose(0) << "Saving changes to xyz." << std::endl);
 | 
|---|
| 62 |   if (comment == "") {
 | 
|---|
| 63 |     time_t now = time((time_t *)NULL);   // Get the system time and put it into 'now' as 'calender time'
 | 
|---|
| 64 |     comment = "Created by molecuilder on ";
 | 
|---|
| 65 |     // ctime ends in \n\0, we have to cut away the newline
 | 
|---|
| 66 |     std::string time(ctime(&now));
 | 
|---|
| 67 |     size_t pos = time.find('\n');
 | 
|---|
| 68 |     if (pos != 0)
 | 
|---|
| 69 |       comment += time.substr(0,pos);
 | 
|---|
| 70 |     else
 | 
|---|
| 71 |       comment += time;
 | 
|---|
| 72 |   }
 | 
|---|
| 73 |   *file << World::getInstance().numAtoms() << endl << "\t" << comment << endl;
 | 
|---|
| 74 | 
 | 
|---|
| 75 |   vector<atom*> atoms = World::getInstance().getAllAtoms();
 | 
|---|
| 76 |   for(vector<atom*>::iterator it = atoms.begin(); it != atoms.end(); it++) {
 | 
|---|
| 77 |     *file << noshowpoint << (*it)->getType()->symbol << "\t" << (*it)->x[0] << "\t" << (*it)->x[1] << "\t" << (*it)->x[2] << endl;
 | 
|---|
| 78 |   }
 | 
|---|
| 79 | }
 | 
|---|