Changes in src/Parser/PcpParser.cpp [84c494:5773a9]
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Parser/PcpParser.cpp
r84c494 r5773a9 20 20 #include "verbose.hpp" 21 21 #include "World.hpp" 22 #include "Matrix.hpp"23 #include "Box.hpp"24 22 25 23 /** Constructor of PcpParser. … … 90 88 IsAngstroem = 1; 91 89 RelativeCoord = 0; 90 StructOpt = 0; 92 91 MaxTypes = 0; 93 92 } … … 211 210 // Unit cell and magnetic field 212 211 ParseForParameter(verbose,FileBuffer, "BoxLength", 0, 3, 3, lower_trigrid, BoxLength, 1, critical); /* Lattice->RealBasis */ 213 double * cell_size = new double[6];212 double * const cell_size = World::getInstance().getDomain(); 214 213 cell_size[0] = BoxLength[0]; 215 214 cell_size[1] = BoxLength[3]; … … 218 217 cell_size[4] = BoxLength[7]; 219 218 cell_size[5] = BoxLength[8]; 220 World::getInstance().setDomain(cell_size);221 delete[] cell_size;222 219 //if (1) fprintf(stderr,"\n"); 223 220 … … 330 327 void PcpParser::save(std::ostream* file) 331 328 { 332 const Matrix &domain = World::getInstance().getDomain().getM();329 const double * const cell_size = World::getInstance().getDomain(); 333 330 class ThermoStatContainer *Thermostats = World::getInstance().getThermostats(); 334 331 if (!file->fail()) { … … 415 412 *file << endl; 416 413 *file << "BoxLength\t\t\t# (Length of a unit cell)" << endl; 417 *file << domain.at(0,0)<< "\t" << endl;418 *file << domain.at(1,0) << "\t" << domain.at(1,1)<< "\t" << endl;419 *file << domain.at(2,0) << "\t" << domain.at(2,1) << "\t" << domain.at(2,2)<< "\t" << endl;414 *file << cell_size[0] << "\t" << endl; 415 *file << cell_size[1] << "\t" << cell_size[2] << "\t" << endl; 416 *file << cell_size[3] << "\t" << cell_size[4] << "\t" << cell_size[5] << "\t" << endl; 420 417 // FIXME 421 418 *file << endl;
Note:
See TracChangeset
for help on using the changeset viewer.