diff --git a/test/enzyme-linalg/inv.jl b/test/enzyme-linalg/inv.jl
deleted file mode 100644
index 35bd074fb..000000000
--- a/test/enzyme-linalg/inv.jl
+++ /dev/null
@@ -1,33 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using Enzyme, EnzymeTestUtils
-using Random
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-is_ci = get(ENV, "CI", "false") == "true"
-TDs = is_ci ? (Duplicated,) : (Const, Duplicated)
-
-@timedtestset "Enzyme - LinearAlgebra (inv):" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-        atol = default_tol(T)
-        rtol = default_tol(T)
-        @testset "inv: TD $TD" for TD in TDs
-            D1 = randn(T, V[1] ← V[1])
-            EnzymeTestUtils.test_reverse(inv, TD, (D1, TD); atol, rtol)
-            EnzymeTestUtils.test_forward(inv, TD, (D1, TD); atol, rtol)
-            if !is_ci
-                D2 = randn(T, V[1] ⊗ V[2] ← V[1] ⊗ V[2])
-                EnzymeTestUtils.test_reverse(inv, TD, (D2, TD); atol, rtol)
-                EnzymeTestUtils.test_forward(inv, TD, (D2, TD); atol, rtol)
-                # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-                @static if VERSION ≥ v"1.11.0-rc"
-                    D3 = randn(T, V[1] ⊗ V[2] ⊗ V[3] ← V[1] ⊗ V[2] ⊗ V[3])
-                    EnzymeTestUtils.test_reverse(inv, TD, (D3, TD); atol, rtol)
-                    EnzymeTestUtils.test_forward(inv, TD, (D3, TD); atol, rtol)
-                end
-            end
-        end
-    end
-end
diff --git a/test/enzyme-linalg/mul.jl b/test/enzyme-linalg/mul.jl
deleted file mode 100644
index 1250df342..000000000
--- a/test/enzyme-linalg/mul.jl
+++ /dev/null
@@ -1,70 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using VectorInterface
-using Enzyme, EnzymeTestUtils
-using Random
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-is_ci = get(ENV, "CI", "false") == "true"
-
-@timedtestset verbose = true "Enzyme - LinearAlgebra (mul):" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-        atol = default_tol(T)
-        rtol = default_tol(T)
-
-        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-        @static if VERSION < v"1.11.0-rc"
-            C = randn(T, V[1] ⊗ V[2] ← (V[4] ⊗ V[5])')
-        else
-            C = randn(T, V[1] ⊗ V[2] ← (V[3] ⊗ V[4] ⊗ V[5])')
-        end
-        A = randn(T, codomain(C) ← V[5]' ⊗ V[4]')
-        B = randn(T, domain(A) ← domain(C))
-        zero_αβs = ((Zero(), Zero()), (randn(T), Zero()), (Zero(), randn(T)))
-        αβs = !is_ci ? vcat(zero_αβs..., (randn(T), randn(T))) : ((randn(T), randn(T)),)
-        for TC in (Duplicated,), TA in (Duplicated,), TB in (Duplicated,)
-            for (α, β) in αβs
-                rTαs = if α === Zero()
-                    (Const,)
-                elseif !is_ci
-                    (Active, Const)
-                else
-                    (Active,)
-                end
-                rTβs = if β === Zero()
-                    (Const,)
-                elseif !is_ci
-                    (Active, Const)
-                else
-                    (Active,)
-                end
-                for Tα in rTαs, Tβ in rTβs
-                    EnzymeTestUtils.test_reverse(mul!, TC, (C, TC), (A, TA), (B, TB), (α, Tα), (β, Tβ); atol, rtol, testset_name = "mul! reverse Tα $Tα, Tβ $Tβ")
-                end
-                fTαs = if α === Zero()
-                    (Const,)
-                elseif !is_ci
-                    (Duplicated, Const)
-                else
-                    (Duplicated,)
-                end
-                fTβs = if β === Zero()
-                    (Const,)
-                elseif !is_ci
-                    (Duplicated, Const)
-                else
-                    (Duplicated,)
-                end
-                for Tα in fTαs, Tβ in fTβs
-                    EnzymeTestUtils.test_forward(mul!, TC, (C, TC), (A, TA), (B, TB), (α, Tα), (β, Tβ); atol, rtol, testset_name = "mul! forward Tα $Tα, Tβ $Tβ")
-                end
-            end
-            if !is_ci
-                EnzymeTestUtils.test_reverse(mul!, TC, (C, TC), (A, TA), (B, TB); atol, rtol, testset_name = "mul! reverse no α no β")
-                EnzymeTestUtils.test_forward(mul!, TC, (C, TC), (A, TA), (B, TB); atol, rtol, testset_name = "mul! forward no α no β")
-            end
-        end
-    end
-end
diff --git a/test/enzyme-linalg/norm.jl b/test/enzyme-linalg/norm.jl
deleted file mode 100644
index 504e51884..000000000
--- a/test/enzyme-linalg/norm.jl
+++ /dev/null
@@ -1,32 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using Enzyme, EnzymeTestUtils
-using Random
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-is_ci = get(ENV, "CI", "false") == "true"
-rRTs = is_ci ? (Active,) : (Const, Active)
-fRTs = is_ci ? (Duplicated,) : (Const, Duplicated)
-
-@timedtestset verbose = true "Enzyme - LinearAlgebra (norm):" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T), TC $TC" for V in spacelist, T in eltypes, TC in (Const, Duplicated)
-        atol = default_tol(T)
-        rtol = default_tol(T)
-        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-        @static if VERSION < v"1.11.0-rc"
-            C = randn(T, V[1] ⊗ V[2] ← (V[4] ⊗ V[5])')
-        else
-            C = randn(T, V[1] ⊗ V[2] ← (V[3] ⊗ V[4] ⊗ V[5])')
-        end
-        for RT in rRTs
-            EnzymeTestUtils.test_reverse(norm, RT, (C, TC), (2, Const); atol, rtol)
-            EnzymeTestUtils.test_reverse(norm, RT, (C', TC), (2, Const); atol, rtol)
-        end
-        for RT in fRTs
-            EnzymeTestUtils.test_forward(norm, RT, (C, TC), (2, Const); atol, rtol)
-            EnzymeTestUtils.test_forward(norm, RT, (C', TC), (2, Const); atol, rtol)
-        end
-    end
-end
diff --git a/test/enzyme-linalg/tr.jl b/test/enzyme-linalg/tr.jl
deleted file mode 100644
index f42e02ee1..000000000
--- a/test/enzyme-linalg/tr.jl
+++ /dev/null
@@ -1,39 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using Enzyme, EnzymeTestUtils
-using Random
-
-is_ci = get(ENV, "CI", "false") == "true"
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-rRTs = is_ci ? (Active,) : (Const, Active)
-fRTs = is_ci ? (Duplicated,) : (Const, Duplicated)
-TDs = is_ci ? (Duplicated,) : (Const, Duplicated)
-
-@timedtestset verbose = true "Enzyme - LinearAlgebra (tr):" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-        atol = default_tol(T)
-        rtol = default_tol(T)
-        D1 = randn(T, V[1] ← V[1])
-        D2 = randn(T, V[1] ⊗ V[2] ← V[1] ⊗ V[2])
-        D3 = randn(T, V[1] ⊗ V[2] ⊗ V[3] ← V[1] ⊗ V[2] ⊗ V[3])
-        @testset "tr reverse: RT $RT, TD $TD" for RT in rRTs, TD in TDs
-            EnzymeTestUtils.test_reverse(tr, RT, (D1, TD); atol, rtol)
-            EnzymeTestUtils.test_reverse(tr, RT, (D2, TD); atol, rtol)
-            # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-            @static if VERSION ≥ v"1.11.0-rc"
-                EnzymeTestUtils.test_reverse(tr, RT, (D3, TD); atol, rtol)
-            end
-        end
-        @testset "tr forward: RT $RT, TD $TD" for RT in fRTs, TD in TDs
-            EnzymeTestUtils.test_forward(tr, RT, (D1, TD); atol, rtol)
-            EnzymeTestUtils.test_forward(tr, RT, (D2, TD); atol, rtol)
-            # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-            @static if VERSION ≥ v"1.11.0-rc"
-                EnzymeTestUtils.test_forward(tr, RT, (D3, TD); atol, rtol)
-            end
-        end
-    end
-end
diff --git a/test/enzyme-vectorinterface/add.jl b/test/enzyme-vectorinterface/add.jl
deleted file mode 100644
index 5d44b2830..000000000
--- a/test/enzyme-vectorinterface/add.jl
+++ /dev/null
@@ -1,51 +0,0 @@
-using Test, TestExtras
-using TensorKit, Enzyme, EnzymeTestUtils
-using TensorOperations
-using Random
-
-#spacelist = ad_spacelist(fast_tests)
-spacelist = [ad_spacelist(fast_tests)[1]]
-eltypes = (Float64, ComplexF64)
-
-@testset "Enzyme - VectorInterface (add!) $(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-    atol = default_tol(T)
-    rtol = default_tol(T)
-
-    α = randn(T)
-    β = randn(T)
-
-    # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-    if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
-        CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
-    else
-        CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
-    end
-    C = randn(T, CV)
-    A = randn(T, CV)
-    for TC in (Duplicated,), TA in (Duplicated,)
-        C = randn(T, CV)
-        A = randn(T, CV)
-        EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA); atol, rtol, testset_name = "add! reverse TC $TC TA $TA no α no β")
-        EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA); atol, rtol, testset_name = "add! forward TC $TC TA $TA no α no β")
-        for Tα in (Active, Const)
-            C = randn(T, CV)
-            A = randn(T, CV)
-            EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA), (α, Tα); atol, rtol, testset_name = "add! reverse TC $TC TA $TA Tα $Tα no β")
-            for Tβ in (Active, Const)
-                C = randn(T, CV)
-                A = randn(T, CV)
-                EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA), (α, Tα), (β, Tβ); atol, rtol, testset_name = "add! reverse TC $TC TA $TA Tα $Tα Tβ $Tβ")
-            end
-        end
-        for Tα in (Duplicated, Const)
-            C = randn(T, CV)
-            A = randn(T, CV)
-            EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA), (α, Tα); atol, rtol, testset_name = "add! forward TC $TC TA $TA Tα $Tα no β")
-            for Tβ in (Duplicated, Const)
-                C = randn(T, CV)
-                A = randn(T, CV)
-                EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA), (α, Tα), (β, Tβ); atol, rtol, testset_name = "add! forward TC $TC TA $TA Tα $Tα Tβ $Tβ")
-            end
-        end
-    end
-end
diff --git a/test/enzyme-vectorinterface/inner.jl b/test/enzyme-vectorinterface/inner.jl
deleted file mode 100644
index 5df4a8017..000000000
--- a/test/enzyme-vectorinterface/inner.jl
+++ /dev/null
@@ -1,31 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using TensorOperations
-using Enzyme, EnzymeTestUtils
-using Random, FiniteDifferences
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-@testset "Enzyme - VectorInterface" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-        @testset for TC in (Duplicated,), TA in (Duplicated,), f in (identity, adjoint)
-            atol = default_tol(T)
-            rtol = default_tol(T)
-            # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-            if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
-                CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
-            else
-                CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
-            end
-            C = randn(T, CV)
-            A = randn(T, CV)
-            for RT in (Active, Const)
-                EnzymeTestUtils.test_reverse(inner, RT, (f(C), TC), (f(A), TA); atol, rtol)
-            end
-            for RT in (Duplicated, Const)
-                EnzymeTestUtils.test_forward(inner, RT, (f(C), TC), (f(A), TA); atol, rtol)
-            end
-        end
-    end
-end
diff --git a/test/enzyme-vectorinterface/scale.jl b/test/enzyme-vectorinterface/scale.jl
deleted file mode 100644
index 8d3f109d5..000000000
--- a/test/enzyme-vectorinterface/scale.jl
+++ /dev/null
@@ -1,46 +0,0 @@
-using Test, TestExtras
-using TensorKit
-using TensorOperations
-using Enzyme, EnzymeTestUtils
-using Random
-
-spacelist = ad_spacelist(fast_tests)
-eltypes = (Float64, ComplexF64)
-
-@testset "Enzyme - VectorInterface (scale!)" begin
-    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
-        atol = default_tol(T)
-        rtol = default_tol(T)
-        α = randn(T)
-        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
-        if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
-            CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
-        else
-            CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
-        end
-        @testset for TC in (Duplicated,)
-            for Tα in (Active, Const)
-                C = randn(T, CV)
-                EnzymeTestUtils.test_reverse(scale!, TC, (C, TC), (α, Tα); atol, rtol)
-                C = randn(T, CV)
-                EnzymeTestUtils.test_reverse(scale!, TC, (C', TC), (α, Tα); atol, rtol)
-                @testset for TA in (Duplicated,), (fc, fa) in ((identity, identity), (adjoint, adjoint))
-                    C = randn(T, CV)
-                    A = randn(T, CV)
-                    EnzymeTestUtils.test_reverse(scale!, TC, (fc(C), TC), (fa(A), TA), (α, Tα); atol, rtol)
-                end
-            end
-            for Tα in (Duplicated, Const)
-                C = randn(T, CV)
-                EnzymeTestUtils.test_forward(scale!, TC, (C, TC), (α, Tα); atol, rtol)
-                C = randn(T, CV)
-                EnzymeTestUtils.test_forward(scale!, TC, (C', TC), (α, Tα); atol, rtol)
-                @testset for TA in (Duplicated,), (fc, fa) in ((identity, identity), (adjoint, adjoint))
-                    C = randn(T, CV)
-                    A = randn(T, CV)
-                    EnzymeTestUtils.test_forward(scale!, TC, (fc(C), TC), (fa(A), TA), (α, Tα); atol, rtol)
-                end
-            end
-        end
-    end
-end
diff --git a/test/enzyme/linalg.jl b/test/enzyme/linalg.jl
new file mode 100644
index 000000000..ef7dec4d3
--- /dev/null
+++ b/test/enzyme/linalg.jl
@@ -0,0 +1,141 @@
+using Test, TestExtras
+using TensorKit
+using VectorInterface
+using Enzyme, EnzymeTestUtils
+using Random
+
+spacelist = ad_spacelist(fast_tests)
+eltypes = (Float64, ComplexF64)
+
+is_ci = get(ENV, "CI", "false") == "true"
+
+@timedtestset verbose = true "Enzyme - LinearAlgebra (mul):" begin
+    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+        atol = default_tol(T)
+        rtol = default_tol(T)
+
+        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+        @static if VERSION < v"1.11.0-rc"
+            C = randn(T, V[1] ⊗ V[2] ← (V[4] ⊗ V[5])')
+        else
+            C = randn(T, V[1] ⊗ V[2] ← (V[3] ⊗ V[4] ⊗ V[5])')
+        end
+        A = randn(T, codomain(C) ← V[5]' ⊗ V[4]')
+        B = randn(T, domain(A) ← domain(C))
+        zero_αβs = ((Zero(), Zero()), (randn(T), Zero()), (Zero(), randn(T)))
+        αβs = !is_ci ? vcat(zero_αβs..., (randn(T), randn(T))) : ((randn(T), randn(T)),)
+        for TC in (Duplicated,), TA in (Duplicated,), TB in (Duplicated,)
+            for (α, β) in αβs
+                rTαs = []
+                α === Zero() || push!(rTαs, Active)
+                (α === Zero() || !is_ci) && push!(rTαs, Const)
+                rTβs = if β === Zero()
+                    (Const,)
+                elseif !is_ci
+                    (Active, Const)
+                else
+                    (Active,)
+                end
+                for Tα in rTαs, Tβ in rTβs
+                    EnzymeTestUtils.test_reverse(mul!, TC, (C, TC), (A, TA), (B, TB), (α, Tα), (β, Tβ); atol, rtol, testset_name = "mul! reverse Tα $Tα, Tβ $Tβ")
+                end
+                fTαs = if α === Zero()
+                    (Const,)
+                elseif !is_ci
+                    (Duplicated, Const)
+                else
+                    (Duplicated,)
+                end
+                fTβs = if β === Zero()
+                    (Const,)
+                elseif !is_ci
+                    (Duplicated, Const)
+                else
+                    (Duplicated,)
+                end
+                for Tα in fTαs, Tβ in fTβs
+                    EnzymeTestUtils.test_forward(mul!, TC, (C, TC), (A, TA), (B, TB), (α, Tα), (β, Tβ); atol, rtol, testset_name = "mul! forward Tα $Tα, Tβ $Tβ")
+                end
+            end
+            if !is_ci
+                EnzymeTestUtils.test_reverse(mul!, TC, (C, TC), (A, TA), (B, TB); atol, rtol, testset_name = "mul! reverse no α no β")
+                EnzymeTestUtils.test_forward(mul!, TC, (C, TC), (A, TA), (B, TB); atol, rtol, testset_name = "mul! forward no α no β")
+            end
+        end
+    end
+end
+
+rRTs = is_ci ? (Active,) : (Const, Active)
+fRTs = is_ci ? (Duplicated,) : (Const, Duplicated)
+
+@timedtestset verbose = true "Enzyme - LinearAlgebra (norm):" begin
+    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T), TC $TC" for V in spacelist, T in eltypes, TC in (Const, Duplicated)
+        atol = default_tol(T)
+        rtol = default_tol(T)
+        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+        @static if VERSION < v"1.11.0-rc"
+            C = randn(T, V[1] ⊗ V[2] ← (V[4] ⊗ V[5])')
+        else
+            C = randn(T, V[1] ⊗ V[2] ← (V[3] ⊗ V[4] ⊗ V[5])')
+        end
+        for RT in rRTs
+            EnzymeTestUtils.test_reverse(norm, RT, (C, TC), (2, Const); atol, rtol)
+            EnzymeTestUtils.test_reverse(norm, RT, (C', TC), (2, Const); atol, rtol)
+        end
+        for RT in fRTs
+            EnzymeTestUtils.test_forward(norm, RT, (C, TC), (2, Const); atol, rtol)
+            EnzymeTestUtils.test_forward(norm, RT, (C', TC), (2, Const); atol, rtol)
+        end
+    end
+end
+
+TDs = is_ci ? (Duplicated,) : (Const, Duplicated)
+
+@timedtestset "Enzyme - LinearAlgebra (inv):" begin
+    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+        atol = default_tol(T)
+        rtol = default_tol(T)
+        @testset "inv: TD $TD" for TD in TDs
+            D1 = randn(T, V[1] ← V[1])
+            EnzymeTestUtils.test_reverse(inv, TD, (D1, TD); atol, rtol)
+            EnzymeTestUtils.test_forward(inv, TD, (D1, TD); atol, rtol)
+            if !is_ci
+                D2 = randn(T, V[1] ⊗ V[2] ← V[1] ⊗ V[2])
+                EnzymeTestUtils.test_reverse(inv, TD, (D2, TD); atol, rtol)
+                EnzymeTestUtils.test_forward(inv, TD, (D2, TD); atol, rtol)
+                # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+                @static if VERSION ≥ v"1.11.0-rc"
+                    D3 = randn(T, V[1] ⊗ V[2] ⊗ V[3] ← V[1] ⊗ V[2] ⊗ V[3])
+                    EnzymeTestUtils.test_reverse(inv, TD, (D3, TD); atol, rtol)
+                    EnzymeTestUtils.test_forward(inv, TD, (D3, TD); atol, rtol)
+                end
+            end
+        end
+    end
+end
+
+@timedtestset verbose = true "Enzyme - LinearAlgebra (tr):" begin
+    @timedtestset "$(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+        atol = default_tol(T)
+        rtol = default_tol(T)
+        D1 = randn(T, V[1] ← V[1])
+        D2 = randn(T, V[1] ⊗ V[2] ← V[1] ⊗ V[2])
+        D3 = randn(T, V[1] ⊗ V[2] ⊗ V[3] ← V[1] ⊗ V[2] ⊗ V[3])
+        @testset "tr reverse: RT $RT, TD $TD" for RT in rRTs, TD in TDs
+            EnzymeTestUtils.test_reverse(tr, RT, (D1, TD); atol, rtol)
+            EnzymeTestUtils.test_reverse(tr, RT, (D2, TD); atol, rtol)
+            # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+            @static if VERSION ≥ v"1.11.0-rc"
+                EnzymeTestUtils.test_reverse(tr, RT, (D3, TD); atol, rtol)
+            end
+        end
+        @testset "tr forward: RT $RT, TD $TD" for RT in fRTs, TD in TDs
+            EnzymeTestUtils.test_forward(tr, RT, (D1, TD); atol, rtol)
+            EnzymeTestUtils.test_forward(tr, RT, (D2, TD); atol, rtol)
+            # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+            @static if VERSION ≥ v"1.11.0-rc"
+                EnzymeTestUtils.test_forward(tr, RT, (D3, TD); atol, rtol)
+            end
+        end
+    end
+end
diff --git a/test/enzyme/vectorinterface.jl b/test/enzyme/vectorinterface.jl
new file mode 100644
index 000000000..1d2e07266
--- /dev/null
+++ b/test/enzyme/vectorinterface.jl
@@ -0,0 +1,107 @@
+using Test, TestExtras
+using TensorKit, Enzyme, EnzymeTestUtils
+using TensorOperations
+using Random
+
+spacelist = ad_spacelist(fast_tests)
+eltypes = (Float64, ComplexF64)
+
+@testset "Enzyme - VectorInterface (add!) $(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+    atol = default_tol(T)
+    rtol = default_tol(T)
+
+    α = randn(T)
+    β = randn(T)
+
+    # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+    if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
+        CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
+    else
+        CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
+    end
+    C = randn(T, CV)
+    A = randn(T, CV)
+    for TC in (Duplicated,), TA in (Duplicated,)
+        C = randn(T, CV)
+        A = randn(T, CV)
+        EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA); atol, rtol, testset_name = "add! reverse TC $TC TA $TA no α no β")
+        EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA); atol, rtol, testset_name = "add! forward TC $TC TA $TA no α no β")
+        for Tα in (Active, Const)
+            C = randn(T, CV)
+            A = randn(T, CV)
+            EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA), (α, Tα); atol, rtol, testset_name = "add! reverse TC $TC TA $TA Tα $Tα no β")
+            for Tβ in (Active, Const)
+                C = randn(T, CV)
+                A = randn(T, CV)
+                EnzymeTestUtils.test_reverse(add!, TC, (C, TC), (A, TA), (α, Tα), (β, Tβ); atol, rtol, testset_name = "add! reverse TC $TC TA $TA Tα $Tα Tβ $Tβ")
+            end
+        end
+        for Tα in (Duplicated, Const)
+            C = randn(T, CV)
+            A = randn(T, CV)
+            EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA), (α, Tα); atol, rtol, testset_name = "add! forward TC $TC TA $TA Tα $Tα no β")
+            for Tβ in (Duplicated, Const)
+                C = randn(T, CV)
+                A = randn(T, CV)
+                EnzymeTestUtils.test_forward(add!, TC, (C, TC), (A, TA), (α, Tα), (β, Tβ); atol, rtol, testset_name = "add! forward TC $TC TA $TA Tα $Tα Tβ $Tβ")
+            end
+        end
+    end
+end
+
+@testset "Enzyme - VectorInterface (scale!) $(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+    atol = default_tol(T)
+    rtol = default_tol(T)
+    α = randn(T)
+    # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+    if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
+        CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
+    else
+        CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
+    end
+    @testset for TC in (Duplicated,)
+        for Tα in (Active, Const)
+            C = randn(T, CV)
+            EnzymeTestUtils.test_reverse(scale!, TC, (C, TC), (α, Tα); atol, rtol)
+            C = randn(T, CV)
+            EnzymeTestUtils.test_reverse(scale!, TC, (C', TC), (α, Tα); atol, rtol)
+            @testset for TA in (Duplicated,), (fc, fa) in ((identity, identity), (adjoint, adjoint))
+                C = randn(T, CV)
+                A = randn(T, CV)
+                EnzymeTestUtils.test_reverse(scale!, TC, (fc(C), TC), (fa(A), TA), (α, Tα); atol, rtol)
+            end
+        end
+        for Tα in (Duplicated, Const)
+            C = randn(T, CV)
+            EnzymeTestUtils.test_forward(scale!, TC, (C, TC), (α, Tα); atol, rtol)
+            C = randn(T, CV)
+            EnzymeTestUtils.test_forward(scale!, TC, (C', TC), (α, Tα); atol, rtol)
+            @testset for TA in (Duplicated,), (fc, fa) in ((identity, identity), (adjoint, adjoint))
+                C = randn(T, CV)
+                A = randn(T, CV)
+                EnzymeTestUtils.test_forward(scale!, TC, (fc(C), TC), (fa(A), TA), (α, Tα); atol, rtol)
+            end
+        end
+    end
+end
+
+@testset "Enzyme - VectorInterface (inner) $(TensorKit.type_repr(sectortype(eltype(V)))) ($T)" for V in spacelist, T in eltypes
+    @testset for TC in (Duplicated,), TA in (Duplicated,), f in (identity, adjoint)
+        atol = default_tol(T)
+        rtol = default_tol(T)
+        # see https://github.com/QuantumKitHub/TensorKit.jl/issues/457
+        if VERSION < v"1.11.0-rc" && sectortype(eltype(V)) == Trivial
+            CV = V[1] ⊗ V[2] ← V[4] ⊗ V[5]
+        else
+            CV = V[1] ⊗ V[2] ← V[3] ⊗ V[4] ⊗ V[5]
+        end
+        C = randn(T, CV)
+        A = randn(T, CV)
+        for RT in (Active, Const)
+            EnzymeTestUtils.test_reverse(inner, RT, (f(C), TC), (f(A), TA); atol, rtol)
+        end
+        for RT in (Duplicated, Const)
+            EnzymeTestUtils.test_forward(inner, RT, (f(C), TC), (f(A), TA); atol, rtol)
+        end
+    end
+end