- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
src/Actions/AnalysisAction/PairCorrelationAction.cpp
rbe21fa r7b1a19 44 44 /** =========== define the function ====================== */ 45 45 Action::state_ptr AnalysisPairCorrelationAction::performCall() { 46 int ranges[3] = {1, 1, 1};47 46 ofstream output; 48 47 ofstream binoutput; … … 51 50 52 51 // execute action 53 output.open(params.outputname. string().c_str());54 binoutput.open(params.binoutputname. string().c_str());52 output.open(params.outputname.get().string().c_str()); 53 binoutput.open(params.binoutputname.get().string().c_str()); 55 54 PairCorrelationMap *correlationmap = NULL; 56 ASSERT(params.elements. size() == 2,55 ASSERT(params.elements.get().size() == 2, 57 56 "AnalysisPairCorrelationAction::performCall() - Exactly two elements are required for pair correlation."); 58 std::vector<const element *>::const_iterator elemiter = params.elements. begin();57 std::vector<const element *>::const_iterator elemiter = params.elements.get().begin(); 59 58 const World::AtomComposite atoms_first = World::getInstance().getAllAtoms(AtomByType(*(elemiter++))); 60 59 const World::AtomComposite atoms_second = World::getInstance().getAllAtoms(AtomByType(*(elemiter++))); 61 ASSERT(elemiter == params.elements. end(),60 ASSERT(elemiter == params.elements.get().end(), 62 61 "AnalysisPairCorrelationAction::performCall() - Exactly two elements are required for pair correlation."); 63 double max_distance = params.BinEnd ;64 if (params.BinEnd <= 0.) {62 double max_distance = params.BinEnd.get(); 63 if (params.BinEnd.get() <= 0.) { 65 64 // find max distance within box from diagonal 66 65 const RealSpaceMatrix &M = World::getInstance().getDomain().getM(); … … 69 68 correlationmap = PairCorrelation(atoms_first, atoms_second, max_distance); 70 69 OutputCorrelationMap<PairCorrelationMap>(&output, correlationmap, OutputPairCorrelation_Header, OutputPairCorrelation_Value); 71 binmap = BinData( correlationmap, params.BinWidth , params.BinStart, params.BinEnd);70 binmap = BinData( correlationmap, params.BinWidth.get(), params.BinStart.get(), params.BinEnd.get() ); 72 71 OutputCorrelationMap<BinPairMap> ( &binoutput, binmap, OutputCorrelation_Header, OutputCorrelation_Value ); 73 72 delete(binmap);
Note:
See TracChangeset
for help on using the changeset viewer.