Ignore:
File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/Actions/AnalysisAction/DipoleAngularCorrelationAction.cpp

    rbe21fa rf10b0c  
    5656
    5757  // 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());
    6060  World::getInstance().clearMoleculeSelection(); // TODO: This should be done in setTime or where molecules are re-done
    6161  World::getInstance().selectAllMolecules(MoleculeByFormula(DipoleFormula));
     
    7171      +toString(DipoleFormula)+" selects no atoms.");
    7272  range<size_t> timesteps = getMaximumTrajectoryBounds(atoms);
    73   ASSERT(params.timestepzero < timesteps.first,
     73  ASSERT(params.timestepzero.get() < timesteps.first,
    7474    "AnalysisDipoleAngularCorrelationAction::performCall() - time step zero "
    75     +toString(params.timestepzero)+" is beyond trajectory range ("
     75    +toString(params.timestepzero.get())+" is beyond trajectory range ("
    7676    +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) {
    7878    // calculate dipoles relative to zero orientation
    7979    DipoleAngularCorrelationMap *correlationmap = NULL;
     
    8787    // output correlation map
    8888    ofstream output;
    89     std::string filename = params.outputname.string()+"."+stepname+".dat";
     89    std::string filename = params.outputname.get().string()+"."+stepname+".dat";
    9090    output.open(filename.c_str());
    9191    OutputCorrelationMap<DipoleAngularCorrelationMap>(&output, correlationmap, OutputDipoleAngularCorrelation_Header, OutputDipoleAngularCorrelation_Value);
     
    9393
    9494    // 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() );
    9696
    9797    // free correlation map
     
    100100    // output binned map
    101101    ofstream binoutput;
    102     std::string binfilename = params.binoutputname.string()+"."+stepname+".dat";
     102    std::string binfilename = params.binoutputname.get().string()+"."+stepname+".dat";
    103103    binoutput.open(binfilename.c_str());
    104104    OutputCorrelationMap<BinPairMap> ( &binoutput, binmap, OutputCorrelation_Header, OutputCorrelation_Value );
Note: See TracChangeset for help on using the changeset viewer.