| 140 |
|
|
| 141 |
|
Vector3d newR12 = i->second * r12; |
| 142 |
|
// x, y and z's possible values range -halfLen_ to halfLen_ |
| 143 |
< |
int xbin = (newR12.x()+ halfLen_) / deltaR_; |
| 143 |
> |
int xbin = (newR12.x() + halfLen_) / deltaR_; |
| 144 |
|
int ybin = (newR12.y() + halfLen_) / deltaR_; |
| 145 |
|
int zbin = (newR12.z() + halfLen_) / deltaR_; |
| 146 |
|
|
| 159 |
|
//rdfStream << "#selection1: (" << selectionScript1_ << ")\t"; |
| 160 |
|
//rdfStream << "selection2: (" << selectionScript2_ << ")\n"; |
| 161 |
|
//rdfStream << "#nRBins = " << nRBins_ << "\t maxLen = " << len_ << "deltaR = " << deltaR_ <<"\n"; |
| 162 |
< |
for (int i = 0; i < histogram_.size(); ++i) { |
| 163 |
< |
|
| 164 |
< |
for(int j = 0; j < histogram_[i].size(); ++j) { |
| 165 |
< |
|
| 166 |
< |
for(int k = 0;k < histogram_[i][j].size(); ++k) { |
| 167 |
< |
rdfStream.write(reinterpret_cast<char *>(&histogram_[i][j][k] ), sizeof(histogram_[i][j][k] )); |
| 162 |
> |
for (unsigned int i = 0; i < histogram_.size(); ++i) { |
| 163 |
> |
for(unsigned int j = 0; j < histogram_[i].size(); ++j) { |
| 164 |
> |
for(unsigned int k = 0;k < histogram_[i][j].size(); ++k) { |
| 165 |
> |
rdfStream.write(reinterpret_cast<char *>(&histogram_[i][j][k] ), |
| 166 |
> |
sizeof(histogram_[i][j][k] )); |
| 167 |
|
} |
| 168 |
|
} |
| 169 |
|
} |
