source: src/Parser/XyzParser.cpp@ e97a44

Action_Thermostats Add_AtomRandomPerturbation Add_FitFragmentPartialChargesAction Add_RotateAroundBondAction Add_SelectAtomByNameAction Added_ParseSaveFragmentResults AddingActions_SaveParseParticleParameters Adding_Graph_to_ChangeBondActions Adding_MD_integration_tests Adding_ParticleName_to_Atom Adding_StructOpt_integration_tests AtomFragments Automaking_mpqc_open AutomationFragmentation_failures Candidate_v1.5.4 Candidate_v1.6.0 Candidate_v1.6.1 ChangeBugEmailaddress ChangingTestPorts ChemicalSpaceEvaluator CombiningParticlePotentialParsing Combining_Subpackages Debian_Package_split Debian_package_split_molecuildergui_only Disabling_MemDebug Docu_Python_wait EmpiricalPotential_contain_HomologyGraph EmpiricalPotential_contain_HomologyGraph_documentation Enable_parallel_make_install Enhance_userguide Enhanced_StructuralOptimization Enhanced_StructuralOptimization_continued Example_ManyWaysToTranslateAtom Exclude_Hydrogens_annealWithBondGraph FitPartialCharges_GlobalError Fix_BoundInBox_CenterInBox_MoleculeActions Fix_ChargeSampling_PBC Fix_ChronosMutex Fix_FitPartialCharges Fix_FitPotential_needs_atomicnumbers Fix_ForceAnnealing Fix_IndependentFragmentGrids Fix_ParseParticles Fix_ParseParticles_split_forward_backward_Actions Fix_PopActions Fix_QtFragmentList_sorted_selection Fix_Restrictedkeyset_FragmentMolecule Fix_StatusMsg Fix_StepWorldTime_single_argument Fix_Verbose_Codepatterns Fix_fitting_potentials Fixes ForceAnnealing_goodresults ForceAnnealing_oldresults ForceAnnealing_tocheck ForceAnnealing_with_BondGraph ForceAnnealing_with_BondGraph_continued ForceAnnealing_with_BondGraph_continued_betteresults ForceAnnealing_with_BondGraph_contraction-expansion FragmentAction_writes_AtomFragments FragmentMolecule_checks_bonddegrees GeometryObjects Gui_Fixes Gui_displays_atomic_force_velocity ImplicitCharges IndependentFragmentGrids IndependentFragmentGrids_IndividualZeroInstances IndependentFragmentGrids_IntegrationTest IndependentFragmentGrids_Sole_NN_Calculation JobMarket_RobustOnKillsSegFaults JobMarket_StableWorkerPool JobMarket_unresolvable_hostname_fix MoreRobust_FragmentAutomation ODR_violation_mpqc_open PartialCharges_OrthogonalSummation PdbParser_setsAtomName PythonUI_with_named_parameters QtGui_reactivate_TimeChanged_changes Recreated_GuiChecks Rewrite_FitPartialCharges RotateToPrincipalAxisSystem_UndoRedo SaturateAtoms_findBestMatching SaturateAtoms_singleDegree StoppableMakroAction Subpackage_CodePatterns Subpackage_JobMarket Subpackage_LinearAlgebra Subpackage_levmar Subpackage_mpqc_open Subpackage_vmg Switchable_LogView ThirdParty_MPQC_rebuilt_buildsystem TrajectoryDependenant_MaxOrder TremoloParser_IncreasedPrecision TremoloParser_MultipleTimesteps TremoloParser_setsAtomName Ubuntu_1604_changes stable
Last change on this file since e97a44 was e97a44, checked in by Frederik Heber <heber@…>, 14 years ago

Added verbose message on save to all derived FormatParser and ChangeTracker.

  • Property mode set to 100644
File size: 2.0 KB
Line 
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
18using namespace std;
19
20/**
21 * Constructor.
22 */
23XyzParser::XyzParser() {
24 comment = "";
25}
26
27/**
28 * Destructor.
29 */
30XyzParser::~XyzParser() {
31}
32
33/**
34 * Loads an XYZ file into the World.
35 *
36 * \param XYZ file
37 */
38void 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 */
60void 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}
Note: See TracBrowser for help on using the repository browser.