diff --git a/docs/examples/howto_propertylib/plot_propertylib.py b/docs/examples/howto_propertylib/plot_propertylib.py
index f8eca6a140cd6335eae785376ac6f4d866e9f9d4..b25ea84cc9177a1b9e06ccff19c48e86f0f3fecc 100644
--- a/docs/examples/howto_propertylib/plot_propertylib.py
+++ b/docs/examples/howto_propertylib/plot_propertylib.py
@@ -49,10 +49,6 @@ print(presets.displacement[1]([0.01, 0.02, 0.03]))
 
 print(presets.displacement.magnitude([0.03, 0.04]))
 
-# %%
-# Log of Magnitude of a 2D velocity vector from the Hydraulics collection:
-print(presets.velocity.log_magnitude(np.sqrt([50, 50])))
-
 # %%
 # Magnitude and trace of a 3D strain matrix:
 eps = np.array([1, 3, 9, 1, 2, 2]) * 1e-2
diff --git a/tests/test_propertylib.py b/tests/test_propertylib.py
index 4568e46d20bf605c06801557bc6c1975c779242c..e751aa2f9b69a0592c9c7be82c2140ca4f6748fd 100644
--- a/tests/test_propertylib.py
+++ b/tests/test_propertylib.py
@@ -47,12 +47,6 @@ class PhysicalPropertyTest(unittest.TestCase):
         self.equality(
             pp.velocity.magnitude, [[3, 4], [1, 0]], qty([5, 1], "m/s")
         )
-        self.equality(pp.velocity.log_magnitude, np.sqrt([50, 50]), qty(1, ""))
-        self.equality(
-            pp.velocity.log_magnitude,
-            [np.sqrt([50, 50]), [10, 0]],
-            qty([1, 1], ""),
-        )
 
     def test_displacement(self):
         """Test displacement property."""