diff --git a/src/replay/__init__.py b/src/replay/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..0a917e157a16a99ff22784e6d3572be4957efe65
--- /dev/null
+++ b/src/replay/__init__.py
@@ -0,0 +1,7 @@
+# This init file exposes isotropic and anisotropic replay methods
+# for easy importing in other parts of the code.
+
+from .isotropic_replay import isotropic_replay
+from .anisotropic_replay import anisotropic_replay
+
+__all__ = ["isotropic_replay", "anisotropic_replay"]