Skip to content

Instantly share code, notes, and snippets.

# from https://stackoverflow.com/questions/22712292/how-to-use-numpy-savez-in-a-loop-for-save-more-than-one-array/22716159#22716159
import numpy as np
import tempfile
class my_savez(object):
def __init__(self, file, compression=None, compresslevel=None):
# Import is postponed to here since zipfile depends on gzip, an optional
# component of the so-called standard library.
import zipfile
@golobor
golobor / bottle-brush-scaling-theory.ipynb
Created June 22, 2020 20:26
Bottle Brush P(s) theory for Gibcus et at 2018
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@golobor
golobor / hela-pairing.ipynb
Created October 14, 2022 15:18
region-wide-P(s)
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@golobor
golobor / 230306-symmetrize-stack.ipynb
Last active March 6, 2023 20:32
230306-symmetrize-stack
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
from contextlib import contextmanager
import functools
import numpy as np
import gsd
import gsd.hoomd
import hoomd
import hoomd.md
@golobor
golobor / h5cp.py
Created December 12, 2023 12:51
h5cp
import h5py
import argparse
def copy_dataset(source_file, destination_file, source_dataset, destination_dataset, overwrite):
with h5py.File(source_file, 'r') as source:
with h5py.File(destination_file, 'a') as destination:
if source_dataset not in source:
print(f"Dataset '{source_dataset}' does not exist in the source file.")
return