| 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 | 
  | 
 | 
| 72 | 
+ | 
      int bytesPerFrameData = Snapshot::getFrameDataSize(); | 
| 73 | 
+ | 
      int bytesPerFrame = (nRigidBodies_ + nAtoms_) * bytesPerStuntDouble  | 
| 74 | 
+ | 
        + nCutoffGroups_ * bytesPerCutoffGroup  | 
| 75 | 
+ | 
        + bytesPerFrameData; | 
| 76 | 
+ | 
 | 
| 77 | 
  | 
      // total number of frames that can fit in memory | 
| 78 | 
  | 
      //RealType frameCapacity = avaliablePhysMem / bytesPerFrame; | 
| 79 | 
  | 
      RealType frameCapacity = (RealType) memSize_ / (RealType) bytesPerFrame; | 
| 81 | 
  | 
      // number of frames in each block given the need to hold multiple blocks  | 
| 82 | 
  | 
      // in memory at the same time: | 
| 83 | 
  | 
      nSnapshotPerBlock_ = int(frameCapacity) / blockCapacity_; | 
| 84 | 
+ | 
      if (nSnapshotPerBlock_ <= 0) { | 
| 85 | 
+ | 
       std::cerr << "not enough memory to hold two configs!" << std::endl; | 
| 86 | 
+ | 
      } | 
| 87 | 
  | 
      reader_ = new DumpReader(info, filename); | 
| 88 | 
  | 
      nframes_ = reader_->getNFrames(); | 
| 89 | 
  | 
      int nblocks = nframes_ / nSnapshotPerBlock_; | 
| 99 | 
  | 
 | 
| 100 | 
  | 
      snapshots_.insert(snapshots_.begin(), nframes_, static_cast<Snapshot*>(NULL));    | 
| 101 | 
  | 
 | 
| 102 | 
< | 
      std::cout << "-----------------------------------------------------"<<std::endl; | 
| 102 | 
> | 
      std::cout << "-----------------------------------------------------" | 
| 103 | 
> | 
                << std::endl; | 
| 104 | 
  | 
      std::cout << "BlockSnapshotManager memory report:" << std::endl; | 
| 105 | 
  | 
      std::cout << "\n"; | 
| 106 | 
  | 
      // std::cout << "  Physical Memory available:\t" << (unsigned long)physMem <<  " bytes" <<std::endl; | 
| 107 | 
  | 
      //std::cout << "     Resident Memory in use:\t" << (unsigned long)rssMem << " bytes" <<std::endl; | 
| 108 | 
  | 
      //std::cout << "Memory available for OpenMD:\t" << (unsigned long)avaliablePhysMem << " bytes" <<std::endl; | 
| 109 | 
< | 
      std::cout << "Memory requested for OpenMD:\t" << (unsigned long)memSize_ << " bytes" <<std::endl; | 
| 110 | 
< | 
      std::cout << "      Bytes per StuntDouble:\t" << (unsigned long)bytesPerStuntDouble <<std::endl; | 
| 111 | 
< | 
      std::cout << "            Bytes per Frame:\t" << (unsigned long)bytesPerFrame <<std::endl; | 
| 112 | 
< | 
      std::cout << "             Frame Capacity:\t" << (unsigned long)frameCapacity <<std::endl; | 
| 113 | 
< | 
      std::cout << "       Frames in trajectory:\t" << (unsigned long)nframes_ <<std::endl; | 
| 114 | 
< | 
      std::cout << "        Snapshots per Block:\t" << (unsigned long)nSnapshotPerBlock_ <<std::endl; | 
| 115 | 
< | 
      std::cout << "     Total number of Blocks:\t" << (unsigned long)nblocks << std::endl; | 
| 116 | 
< | 
      std::cout << "-----------------------------------------------------"<<std::endl; | 
| 109 | 
> | 
      std::cout << "Memory requested for OpenMD:\t"  | 
| 110 | 
> | 
                << (unsigned long)memSize_ << " bytes" << std::endl; | 
| 111 | 
> | 
      std::cout << "        Bytes per FrameData:\t"  | 
| 112 | 
> | 
                << (unsigned long)bytesPerFrameData << std::endl; | 
| 113 | 
> | 
      std::cout << "      Bytes per StuntDouble:\t"  | 
| 114 | 
> | 
                << (unsigned long)bytesPerStuntDouble << std::endl; | 
| 115 | 
> | 
      std::cout << "     Bytes per Cutoff Group:\t"  | 
| 116 | 
> | 
                << (unsigned long)bytesPerCutoffGroup << std::endl; | 
| 117 | 
> | 
      std::cout << "            Bytes per Frame:\t"  | 
| 118 | 
> | 
                << (unsigned long)bytesPerFrame << std::endl; | 
| 119 | 
> | 
      std::cout << "             Frame Capacity:\t" | 
| 120 | 
> | 
                << (unsigned long)frameCapacity << std::endl; | 
| 121 | 
> | 
      std::cout << "       Frames in trajectory:\t"  | 
| 122 | 
> | 
                << (unsigned long)nframes_ << std::endl; | 
| 123 | 
> | 
      std::cout << "        Snapshots per Block:\t"  | 
| 124 | 
> | 
                << (unsigned long)nSnapshotPerBlock_ << std::endl; | 
| 125 | 
> | 
      std::cout << "     Total number of Blocks:\t"  | 
| 126 | 
> | 
                << (unsigned long)nblocks << std::endl; | 
| 127 | 
> | 
      std::cout << "-----------------------------------------------------" | 
| 128 | 
> | 
                << std::endl; | 
| 129 | 
  | 
     | 
| 130 | 
  | 
    } | 
| 131 | 
  | 
 | 
| 196 | 
  | 
        //in case, unloadBlock called multiple times | 
| 197 | 
  | 
        activeRefCount_[i - activeBlocks_.begin()]  = 0; | 
| 198 | 
  | 
      } | 
| 199 | 
+ | 
 | 
| 200 | 
+ | 
      if (activeRefCount_[i-activeBlocks_.begin()] == 0) { | 
| 201 | 
+ | 
        internalUnload(block); | 
| 202 | 
+ | 
      } | 
| 203 | 
  | 
         | 
| 204 | 
  | 
      unloadSuccess = true; | 
| 205 | 
  | 
    } else { | 
| 223 | 
  | 
  } | 
| 224 | 
  | 
 | 
| 225 | 
  | 
  void BlockSnapshotManager::internalUnload(int block) { | 
| 226 | 
+ | 
    std::cerr << "called internal unload for block "<< block << "\n"; | 
| 227 | 
  | 
    for (int i = blocks_[block].first; i < blocks_[block].second; ++i) { | 
| 228 | 
  | 
      delete snapshots_[i]; | 
| 229 | 
  | 
      snapshots_[i] = NULL; |