Method updateStoragesAndSizes

src/Numerical/ODEVector.cpp:35–45

}

void ODEVector::updateStoragesAndSizes(std::vector<real*> newStorages,
                                       std::vector<std::size_t> newSizes) {
  storages = std::move(newStorages);
  sizes = std::move(newSizes);
  offsets.clear();
  std::size_t curOffset = 0;
  for (const unsigned long size : this->sizes) {
    offsets.push_back(curOffset);
    curOffset += size;
  }
}

std::pair<real*, size_t> ODEVector::getSubvector(size_t storageIdx) {
  return {storages[storageIdx], sizes[storageIdx]};