diff --git a/test.py b/test.py
new file mode 100644
index 00000000..02dcf7ff
--- /dev/null
+++ b/test.py
@@ -0,0 +1,55 @@
+#!/usr/bin/env python
+import torch
+from tqdm import tqdm
+
+from gfn.gym import HyperGrid
+from gfn.gflownet import TBGFlowNet
+from gfn.modules import DiscretePolicyEstimator
+from gfn.samplers import Sampler
+from gfn.utils import NeuralNet
+
+
+torch.manual_seed(0)
+exploration_rate = 0.5
+learning_rate = 0.0005
+
+env = HyperGrid(ndim=5, height=2)
+
+module_PF = NeuralNet(
+    input_dim=env.preprocessor.output_dim,
+    output_dim=env.n_actions,
+)
+module_PB = NeuralNet(
+    input_dim=env.preprocessor.output_dim,
+    output_dim=env.n_actions - 1,
+    torso=module_PF.torso,
+)
+
+pf_estimator = DiscretePolicyEstimator(module_PF, env.n_actions, is_backward=False, preprocessor=env.preprocessor)
+pb_estimator = DiscretePolicyEstimator(module_PB, env.n_actions, is_backward=True, preprocessor=env.preprocessor)
+gflownet = TBGFlowNet(init_logZ=0., pf=pf_estimator, pb=pb_estimator, off_policy=True)
+sampler = Sampler(estimator=pf_estimator)
+
+# Policy parameters have their own LR.
+non_logz_params = [v for k, v in dict(gflownet.named_parameters()).items() if k != "logZ"]
+optimizer = torch.optim.Adam(non_logz_params, lr=1e-3)
+
+# Log Z gets dedicated learning rate (typically higher).
+logz_params = [dict(gflownet.named_parameters())["logZ"]]
+optimizer.add_param_group({"params": logz_params, "lr": 1e-1})
+
+n_iterations = int(1e4)
+batch_size = int(1e5)
+
+for i in (pbar := tqdm(range(n_iterations))):
+    trajectories = sampler.sample_trajectories(
+        env,
+        n_trajectories=batch_size,
+        off_policy=True,
+        epsilon=exploration_rate,
+    )
+    optimizer.zero_grad()
+    loss = gflownet.loss(env, trajectories)
+    loss.backward()
+    optimizer.step()
+    pbar.set_postfix({"loss": loss.item()})