| 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 | }
|
|---|