diff --git a/src/gauging.jl b/src/gauging.jl
index 2ad9d0f4..49384ca6 100644
--- a/src/gauging.jl
+++ b/src/gauging.jl
@@ -128,7 +128,7 @@ function VidalITensorNetwork(
   ψ::ITensorNetwork;
   (cache!)=nothing,
   update_cache=isnothing(cache!),
-  cache_update_kwargs=default_cache_update_kwargs(cache!),
+  cache_update_kwargs=default_cache_update_kwargs(Algorithm("bp")),
   kwargs...,
 )
   if isnothing(cache!)
diff --git a/test/test_forms.jl b/test/test_forms.jl
index a58822e5..7e6ada8b 100644
--- a/test/test_forms.jl
+++ b/test/test_forms.jl
@@ -62,7 +62,7 @@ using Test: @test, @testset
   @test underlying_graph(ket_network(qf)) == underlying_graph(ψket)
   @test underlying_graph(operator_network(qf)) == underlying_graph(A)
 
-  ∂qf_∂v = only(environment(qf, state_vertices(qf, [v])))
+  ∂qf_∂v = only(environment(qf, state_vertices(qf, [v]); alg="exact"))
   @test (∂qf_∂v) * (qf[ket_vertex(qf, v)] * qf[bra_vertex(qf, v)]) ≈ contract(qf)
 
   ∂qf_∂v_bp = environment(qf, state_vertices(qf, [v]); alg="bp", update_cache=false)