- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/AnalysisAction/DipoleAngularCorrelationAction.cpp
rbe21fa rf10b0c 56 56 57 57 // select atoms and obtain zero dipole orientation 58 Formula DipoleFormula(params.DipoleFormula );59 World::getInstance().setTime(params.timestepzero );58 Formula DipoleFormula(params.DipoleFormula.get()); 59 World::getInstance().setTime(params.timestepzero.get()); 60 60 World::getInstance().clearMoleculeSelection(); // TODO: This should be done in setTime or where molecules are re-done 61 61 World::getInstance().selectAllMolecules(MoleculeByFormula(DipoleFormula)); … … 71 71 +toString(DipoleFormula)+" selects no atoms."); 72 72 range<size_t> timesteps = getMaximumTrajectoryBounds(atoms); 73 ASSERT(params.timestepzero < timesteps.first,73 ASSERT(params.timestepzero.get() < timesteps.first, 74 74 "AnalysisDipoleAngularCorrelationAction::performCall() - time step zero " 75 +toString(params.timestepzero )+" is beyond trajectory range ("75 +toString(params.timestepzero.get())+" is beyond trajectory range (" 76 76 +toString(timesteps.first)+") of some atoms."); 77 for (size_t step = params.timestepzero ; step < timesteps.first; ++step) {77 for (size_t step = params.timestepzero.get(); step < timesteps.first; ++step) { 78 78 // calculate dipoles relative to zero orientation 79 79 DipoleAngularCorrelationMap *correlationmap = NULL; … … 87 87 // output correlation map 88 88 ofstream output; 89 std::string filename = params.outputname. string()+"."+stepname+".dat";89 std::string filename = params.outputname.get().string()+"."+stepname+".dat"; 90 90 output.open(filename.c_str()); 91 91 OutputCorrelationMap<DipoleAngularCorrelationMap>(&output, correlationmap, OutputDipoleAngularCorrelation_Header, OutputDipoleAngularCorrelation_Value); … … 93 93 94 94 // bin map 95 BinPairMap *binmap = BinData( correlationmap, params.BinWidth , params.BinStart, params.BinEnd);95 BinPairMap *binmap = BinData( correlationmap, params.BinWidth.get(), params.BinStart.get(), params.BinEnd.get() ); 96 96 97 97 // free correlation map … … 100 100 // output binned map 101 101 ofstream binoutput; 102 std::string binfilename = params.binoutputname. string()+"."+stepname+".dat";102 std::string binfilename = params.binoutputname.get().string()+"."+stepname+".dat"; 103 103 binoutput.open(binfilename.c_str()); 104 104 OutputCorrelationMap<BinPairMap> ( &binoutput, binmap, OutputCorrelation_Header, OutputCorrelation_Value );
Note:
See TracChangeset
for help on using the changeset viewer.