source: src/Parser/XyzParser.cpp@ b5c53d

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 b5c53d was b5c53d, checked in by Frederik Heber <heber@…>, 14 years ago

Merge branch 'StructureRefactoring' into stable

Conflicts:

src/Actions/AtomAction/AddAction.cpp
src/Actions/AtomAction/ChangeElementAction.cpp
src/Parser/XyzParser.cpp
src/analysis_correlation.cpp
src/atom.cpp
src/config.cpp
src/molecule.cpp

  • AtomInfo::element were privatized in stable and element::symbol, ::name in StructureRefactoring (overlapped in various lines).
  • Property mode set to 100644
File size: 2.2 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 "molecule.hpp"
16#include "element.hpp"
17#include "periodentafel.hpp"
18
19using namespace std;
20
21/**
22 * Constructor.
23 */
24XyzParser::XyzParser() :
25 comment("")
26{}
27
28/**
29 * Destructor.
30 */
31XyzParser::~XyzParser() {
32}
33
34/**
35 * Loads an XYZ file into the World.
36 *
37 * \param XYZ file
38 */
39void XyzParser::load(istream* file) {
40 atom* newAtom = NULL;
41 molecule* newmol = NULL;
42 int numberOfAtoms;
43 char commentBuffer[512], type[3];
44 double tmp;
45
46 // the first line tells number of atoms, the second line is always a comment
47 *file >> numberOfAtoms >> ws;
48 file->getline(commentBuffer, 512);
49 comment = commentBuffer;
50
51 newmol = World::getInstance().createMolecule();
52 newmol->ActiveFlag = true;
53 World::getInstance().getMolecules()->insert(newmol);
54 for (int i = 0; i < numberOfAtoms; i++) {
55 newAtom = World::getInstance().createAtom();
56 *file >> type;
57 for (int j=0;j<NDIM;j++) {
58 *file >> tmp;
59 newAtom->set(j, tmp);
60 }
61 newAtom->setType(World::getInstance().getPeriode()->FindElement(type));
62 newmol->AddAtom(newAtom);
63 }
64}
65
66/**
67 * Saves the current state of the World into the given XYZ file.
68 *
69 * \param XYZ file
70 */
71void XyzParser::save(ostream* file) {
72 DoLog(0) && (Log() << Verbose(0) << "Saving changes to xyz." << std::endl);
73 if (comment == "") {
74 time_t now = time((time_t *)NULL); // Get the system time and put it into 'now' as 'calender time'
75 comment = "Created by molecuilder on ";
76 // ctime ends in \n\0, we have to cut away the newline
77 std::string time(ctime(&now));
78 size_t pos = time.find('\n');
79 if (pos != 0)
80 comment += time.substr(0,pos);
81 else
82 comment += time;
83 }
84 *file << World::getInstance().numAtoms() << endl << "\t" << comment << endl;
85
86 vector<atom*> atoms = World::getInstance().getAllAtoms();
87 for(vector<atom*>::iterator it = atoms.begin(); it != atoms.end(); it++) {
88 *file << noshowpoint << (*it)->getType()->getSymbol() << "\t" << (*it)->at(0) << "\t" << (*it)->at(1) << "\t" << (*it)->at(2) << endl;
89 }
90}
Note: See TracBrowser for help on using the repository browser.