146 |
|
} |
147 |
|
|
148 |
|
void GofXyz::writeRdf() { |
149 |
< |
std::ofstream rdfStream(outputFilename_.c_str()); |
149 |
> |
std::ofstream rdfStream(outputFilename_.c_str(), std::ios::binary); |
150 |
|
if (rdfStream.is_open()) { |
151 |
< |
rdfStream << "#g(x, y, z)\n"; |
152 |
< |
rdfStream << "#selection1: (" << selectionScript1_ << ")\t"; |
153 |
< |
rdfStream << "selection2: (" << selectionScript2_ << ")\n"; |
154 |
< |
rdfStream << "#nRBins = " << nRBins_ << "\t maxLen = " << len_ << "deltaR = " << deltaR_ <<"\n"; |
151 |
> |
//rdfStream << "#g(x, y, z)\n"; |
152 |
> |
//rdfStream << "#selection1: (" << selectionScript1_ << ")\t"; |
153 |
> |
//rdfStream << "selection2: (" << selectionScript2_ << ")\n"; |
154 |
> |
//rdfStream << "#nRBins = " << nRBins_ << "\t maxLen = " << len_ << "deltaR = " << deltaR_ <<"\n"; |
155 |
|
for (int i = 0; i < histogram_.size(); ++i) { |
156 |
|
|
157 |
|
for(int j = 0; j < histogram_[i].size(); ++j) { |
158 |
|
|
159 |
|
for(int k = 0;k < histogram_[i].size(); ++k) { |
160 |
< |
|
161 |
< |
rdfStream << histogram_[i][j][k]/nProcessed_ << "\t"; |
160 |
> |
rdfStream.write(reinterpret_cast<char *>(&histogram_[i][j][k] ), sizeof(histogram_[i][j][k] )); |
161 |
|
} |
163 |
– |
rdfStream << "\n"; |
162 |
|
} |
163 |
|
} |
164 |
|
|