diff --git a/docs/src/tutorials/linear_methods/discrete_time.md b/docs/src/tutorials/linear_methods/discrete_time.md
index b11f49b0c..a618d1556 100644
--- a/docs/src/tutorials/linear_methods/discrete_time.md
+++ b/docs/src/tutorials/linear_methods/discrete_time.md
@@ -47,7 +47,7 @@ V = vertices_list(X0) .|> Singleton |> UnionSetArray
 plot!(c)
 plot!(V)
 plot!(M(pi/4) * c) # rotate 45 degrees
-plot!(M(pi/4) * V)
+plot!(linear_map(M(pi/4), V))
 
 xlims!(0.0, 1.8) # hide
 ylims!(-0.4, 1.4) # hide
@@ -60,9 +60,9 @@ fig = DisplayAs.Text(DisplayAs.PNG(fig))  # hide
 X = sample(X0, 25) .|> Singleton |> UnionSetArray
 
 plot!(X)
-plot!(M(pi/4) * X, c=:blue)     # rotate 45 degrees
-plot!(M(pi/2) * X, c=:green)    # rotate 90 degrees
-plot!(M(4pi/3) * X, c=:orange)  # rotate 180 degrees
+plot!(linear_map(M(pi/4), X), c=:blue)     # rotate 45 degrees
+plot!(linear_map(M(pi/2), X), c=:green)    # rotate 90 degrees
+plot!(linear_map(M(4pi/3), X), c=:orange)  # rotate 180 degrees
 
 plot!(X0, c=:white)
 plot!(M(pi/4) * X0, c=:white)
diff --git a/src/Flowpipes/setops.jl b/src/Flowpipes/setops.jl
index 1b24f3b8d..1a72588d6 100644
--- a/src/Flowpipes/setops.jl
+++ b/src/Flowpipes/setops.jl
@@ -217,10 +217,6 @@ function _split_symmetric_box(D::Int, partition::Vector{Int})
     return convert.(IA.IntervalBox, Sp)
 end
 
-function Base.:(*)(M::AbstractMatrix, X::UnionSetArray{N,<:AbstractSingleton{N}}) where {N}
-    return UnionSetArray([linear_map(M, p) for p in array(X)])
-end
-
 # ==================================
 # Zonotope order reduction methods
 # ==================================