fix mutation and test

Author: ChrisRackauckas

src/concrete_solve.jl

@@ -1288,7 +1288,7 @@ function DiffEqBase._concrete_solve_adjoint(
     end
 
     forward, reverse = Enzyme.autodiff_thunk(splitmode, Enzyme.Const{typeof(f)}, Enzyme.Duplicated, Enzyme.Duplicated{typeof(u0)}, Enzyme.Duplicated{typeof(p)})
-    tape, result, shadow_result = forward(Enzyme.Const(f), Enzyme.Duplicated(u0, du0), Enzyme.Duplicated(p, dp))
+    tape, result, shadow_result = forward(Enzyme.Const(f), Enzyme.Duplicated(copy(u0), du0), Enzyme.Duplicated(copy(p), dp))
 
     function enzyme_sensitivity_backpass(Δ)
         reverse(Const(f), Duplicated(u0, du0), Duplicated(p, dp), Δ, tape)

test/concrete_solve_derivatives.jl

@@ -354,6 +354,7 @@ du09, dp9 = Zygote.gradient(
 @test ū0≈du06 rtol=1e-12
 @test_broken ū0≈du07 rtol=1e-12
 @test_broken ū0≈du08 rtol=1e-12
+@test ū0≈du09 rtol=1e-12
 @test adj≈dp1' rtol=1e-12
 @test adj≈dp2' rtol=1e-12
 @test adj≈dp3' rtol=1e-12