diff --git a/mcdc/kernel.py b/mcdc/kernel.py
index 47b35cef..3aff6b2c 100644
--- a/mcdc/kernel.py
+++ b/mcdc/kernel.py
@@ -3744,6 +3744,8 @@ def modified_gram_schmidt(V, u):
     v2 = w2 / np.linalg.norm(w2)
     V = np.append(V, v2, axis=1)
     return V
+
+
 # =============================================================================
 # Variance Deconvolution
 # =============================================================================
diff --git a/mcdc/main.py b/mcdc/main.py
index db368b6b..da3e7369 100644
--- a/mcdc/main.py
+++ b/mcdc/main.py
@@ -716,18 +716,31 @@ def generate_hdf5(mcdc):
                     "iqmc/tally/flux", data=np.squeeze(T["iqmc"]["score"]["flux"])
                 )
                 f.create_dataset(
-                    "iqmc/tally/fission_source", data=T["iqmc"]["score"]["fission-source"]
+                    "iqmc/tally/fission_source",
+                    data=T["iqmc"]["score"]["fission-source"],
                 )
                 f.create_dataset(
                     "iqmc/tally/fission_power", data=T["iqmc"]["score"]["fission-power"]
                 )
                 f.create_dataset("iqmc/tally/source_constant", data=T["iqmc"]["source"])
-                f.create_dataset("iqmc/tally/source_x", data=T["iqmc"]["score"]["tilt-x"])
-                f.create_dataset("iqmc/tally/source_y", data=T["iqmc"]["score"]["tilt-y"])
-                f.create_dataset("iqmc/tally/source_z", data=T["iqmc"]["score"]["tilt-z"])
-                f.create_dataset("iqmc/tally/source_xy", data=T["iqmc"]["score"]["tilt-xy"])
-                f.create_dataset("iqmc/tally/source_xz", data=T["iqmc"]["score"]["tilt-xz"])
-                f.create_dataset("iqmc/tally/source_yz", data=T["iqmc"]["score"]["tilt-yz"])
+                f.create_dataset(
+                    "iqmc/tally/source_x", data=T["iqmc"]["score"]["tilt-x"]
+                )
+                f.create_dataset(
+                    "iqmc/tally/source_y", data=T["iqmc"]["score"]["tilt-y"]
+                )
+                f.create_dataset(
+                    "iqmc/tally/source_z", data=T["iqmc"]["score"]["tilt-z"]
+                )
+                f.create_dataset(
+                    "iqmc/tally/source_xy", data=T["iqmc"]["score"]["tilt-xy"]
+                )
+                f.create_dataset(
+                    "iqmc/tally/source_xz", data=T["iqmc"]["score"]["tilt-xz"]
+                )
+                f.create_dataset(
+                    "iqmc/tally/source_yz", data=T["iqmc"]["score"]["tilt-yz"]
+                )
                 # iteration data
                 f.create_dataset("iqmc/itteration_count", data=T["iqmc"]["itt"])
                 f.create_dataset("iqmc/final_residual", data=T["iqmc"]["res"])
diff --git a/test/regression/iqmc_cooper2/input.py b/test/regression/iqmc_cooper2/input.py
index efb22de1..1de33bda 100644
--- a/test/regression/iqmc_cooper2/input.py
+++ b/test/regression/iqmc_cooper2/input.py
@@ -55,7 +55,7 @@
     maxitt=maxit,
     tol=tol,
     generator=generator,
-    score=['tilt-x', 'tilt-y', 'tilt-xy']
+    score=["tilt-x", "tilt-y", "tilt-xy"],
 )
 
 # =============================================================================
diff --git a/test/regression/iqmc_kornreich_pi/input.py b/test/regression/iqmc_kornreich_pi/input.py
index 6ea5aebf..50efbebf 100644
--- a/test/regression/iqmc_kornreich_pi/input.py
+++ b/test/regression/iqmc_kornreich_pi/input.py
@@ -57,7 +57,7 @@
     tol=tol,
     generator=generator,
     eigenmode_solver=solver,
-    score=['tilt-x']
+    score=["tilt-x"],
 )
 # Setting
 mcdc.setting(N_particle=N)
diff --git a/test/regression/iqmc_reed/input.py b/test/regression/iqmc_reed/input.py
index 05fd4c84..866fb7a9 100644
--- a/test/regression/iqmc_reed/input.py
+++ b/test/regression/iqmc_reed/input.py
@@ -60,7 +60,7 @@
     tol=tol,
     generator=generator,
     fixed_source_solver=solver,
-    score=['tilt-x']
+    score=["tilt-x"],
 )
 
 # =============================================================================