67 |
|
// RealType avaliablePhysMem = physMem - rssMem; |
68 |
|
|
69 |
|
int bytesPerStuntDouble = DataStorage::getBytesPerStuntDouble(storageLayout); |
70 |
< |
int bytesPerFrame = (nRigidBodies_ + nAtoms_) * bytesPerStuntDouble; |
70 |
> |
int bytesPerCutoffGroup = DataStorage::getBytesPerStuntDouble(DataStorage::dslPosition); |
71 |
> |
int bytesPerFrame = (nRigidBodies_ + nAtoms_) * bytesPerStuntDouble |
72 |
> |
+ nCutoffGroups_ * bytesPerCutoffGroup; |
73 |
|
|
74 |
|
// total number of frames that can fit in memory |
75 |
|
//RealType frameCapacity = avaliablePhysMem / bytesPerFrame; |
78 |
|
// number of frames in each block given the need to hold multiple blocks |
79 |
|
// in memory at the same time: |
80 |
|
nSnapshotPerBlock_ = int(frameCapacity) / blockCapacity_; |
81 |
+ |
if (nSnapshotPerBlock_ <= 0) { |
82 |
+ |
std::cerr << "not enough memory to hold two configs!" << std::endl; |
83 |
+ |
} |
84 |
|
reader_ = new DumpReader(info, filename); |
85 |
|
nframes_ = reader_->getNFrames(); |
86 |
|
int nblocks = nframes_ / nSnapshotPerBlock_; |
104 |
|
//std::cout << "Memory available for OpenMD:\t" << (unsigned long)avaliablePhysMem << " bytes" <<std::endl; |
105 |
|
std::cout << "Memory requested for OpenMD:\t" << (unsigned long)memSize_ << " bytes" <<std::endl; |
106 |
|
std::cout << " Bytes per StuntDouble:\t" << (unsigned long)bytesPerStuntDouble <<std::endl; |
107 |
+ |
std::cout << " Bytes per Cutoff Group:\t" << (unsigned long)bytesPerCutoffGroup <<std::endl; |
108 |
|
std::cout << " Bytes per Frame:\t" << (unsigned long)bytesPerFrame <<std::endl; |
109 |
|
std::cout << " Frame Capacity:\t" << (unsigned long)frameCapacity <<std::endl; |
110 |
|
std::cout << " Frames in trajectory:\t" << (unsigned long)nframes_ <<std::endl; |
181 |
|
//in case, unloadBlock called multiple times |
182 |
|
activeRefCount_[i - activeBlocks_.begin()] = 0; |
183 |
|
} |
184 |
+ |
|
185 |
+ |
if (activeRefCount_[i-activeBlocks_.begin()] == 0) { |
186 |
+ |
internalUnload(block); |
187 |
+ |
} |
188 |
|
|
189 |
|
unloadSuccess = true; |
190 |
|
} else { |
208 |
|
} |
209 |
|
|
210 |
|
void BlockSnapshotManager::internalUnload(int block) { |
211 |
+ |
std::cerr << "called internal unload for block "<< block << "\n"; |
212 |
|
for (int i = blocks_[block].first; i < blocks_[block].second; ++i) { |
213 |
|
delete snapshots_[i]; |
214 |
|
snapshots_[i] = NULL; |