diff --git a/src/napari_stress/_measurements/toolbox.py b/src/napari_stress/_measurements/toolbox.py
index 7df89a25..7007f26c 100644
--- a/src/napari_stress/_measurements/toolbox.py
+++ b/src/napari_stress/_measurements/toolbox.py
@@ -499,9 +499,13 @@ def comprehensive_analysis(
         gamma=gamma,
     )
 
-    max_min_anisotropy = 2 * gamma * (
-        Expander_ellipsoid.properties["maximum_mean_curvature"]
-        - Expander_ellipsoid.properties["minimum_mean_curvature"]
+    max_min_anisotropy = (
+        2
+        * gamma
+        * (
+            Expander_ellipsoid.properties["maximum_mean_curvature"]
+            - Expander_ellipsoid.properties["minimum_mean_curvature"]
+        )
     )
 
     result = measurements.tissue_stress_tensor(
diff --git a/src/napari_stress/_tests/test_measurements.py b/src/napari_stress/_tests/test_measurements.py
index 63a57319..8ac93767 100644
--- a/src/napari_stress/_tests/test_measurements.py
+++ b/src/napari_stress/_tests/test_measurements.py
@@ -439,21 +439,21 @@ def test_comprehensive_stress_toolbox_4d(make_napari_viewer):
     # Check against legacy measurements: load legacy data
     legacy_results_file = Path(__file__).parent / "stress_tissue_anisotropy.csv"
     df_legacy = pd.read_csv(legacy_results_file, header=None).T
-    df_legacy.columns = ['tissue_stress_anisotropy']
+    df_legacy.columns = ["tissue_stress_anisotropy"]
 
-    results_layer = viewer.layers['Result of lebedev quadrature on ellipsoid']
+    results_layer = viewer.layers["Result of lebedev quadrature on ellipsoid"]
     result_tissue_stress = results_layer.metadata[
-        types._METADATAKEY_STRESS_TISSUE_ANISO]
-    
-    difference_abs = result_tissue_stress - df_legacy['tissue_stress_anisotropy'].values
-    difference_rel = difference_abs / df_legacy['tissue_stress_anisotropy'].values
+        types._METADATAKEY_STRESS_TISSUE_ANISO
+    ]
+
+    difference_abs = result_tissue_stress - df_legacy["tissue_stress_anisotropy"].values
+    difference_rel = difference_abs / df_legacy["tissue_stress_anisotropy"].values
 
     # make sure that the results are within limits of agreement
     assert np.all(np.abs(difference_rel).mean() < 0.1)
     assert np.all(np.abs(difference_abs).mean() < 0.1)
 
 
-
 def test_curvature(make_napari_viewer):
     from napari_stress._spherical_harmonics.spherical_harmonics_napari import (
         perform_lebedev_quadrature,
@@ -687,6 +687,7 @@ def test_stresses():
     measurements.maximal_tissue_anisotropy(ellipsoid, gamma=gamma)
 
 
-if __name__ == '__main__':
+if __name__ == "__main__":
     import napari
-    test_comprehensive_stress_toolbox_4d(napari.Viewer)
\ No newline at end of file
+
+    test_comprehensive_stress_toolbox_4d(napari.Viewer)