@@ -65,7 +65,7 @@ cache4 = ForwardDiff.Dual{typeof(ForwardDiff.Tag(Nothing, eltype(x))), eltype(x)
6565@test auto_hesvecgrad! (dy, h, x, v, cache1, cache2)≈ ForwardDiff. hessian (g, x) * v rtol= 1e-2
6666@test auto_hesvecgrad (h, x, v)≈ ForwardDiff. hessian (g, x) * v rtol= 1e-2
6767
68- # JacVec
68+ # ## JacVec
6969
7070L = JacVec (f, x)
7171@test L * x ≈ auto_jacvec (f, x, x)
@@ -88,21 +88,7 @@ dy=rand(N);_dy=copy(dy);@test mul!(dy,L,v,a,b) ≈ a*num_jacvec(f,x,v) + b*_dy r
8888out = similar (v)
8989gmres! (out, L, v)
9090
91- #=
92- ff1 = ODEFunction(lorenz, jac_prototype = JacVec{Float64}(lorenz, u0))
93- ff2 = ODEFunction(lorenz, jac_prototype = JacVec{Float64}(lorenz, u0, autodiff=false))
94-
95- for ff in [ff1, ff2]
96- prob = ODEProblem(ff, u0, tspan)
97- @test solve(prob, TRBDF2()).retcode == :Success
98- @test solve(prob, TRBDF2(linsolve = KrylovJL_GMRES())).retcode == :Success
99- @test solve(prob, Exprb32()).retcode == :Success
100- @test solve(prob, Rosenbrock23()).retcode == :Success
101- @test solve(prob, Rosenbrock23(linsolve = KrylovJL_GMRES())).retcode == :Success
102- end
103- =#
104-
105- # HesVec
91+ # ## HesVec
10692
10793x = rand (N)
10894v = rand (N)
@@ -127,7 +113,7 @@ dy=rand(N);_dy=copy(dy);@test mul!(dy,L,v,a,b)≈a*numauto_hesvec(g,x,v)+b*_dy r
127113out = similar (v)
128114gmres! (out, L, v)
129115
130- # HesVecGrad
116+ # ## HesVecGrad
131117
132118x = rand (N)
133119v = rand (N)
0 commit comments