Set jacobian of Fun to ZeroOperator (#909)

LucasAschenbach · web-flow · commit 615ad0385f17 · 2023-08-28T16:14:34.000+04:00
* jacobian of Fun is ZeroOperator

* Decrease maxiterations for independent variable newton test
diff --git a/src/Extras/autodifferentiation.jl b/src/Extras/autodifferentiation.jl
@@ -61,7 +61,7 @@ last(d::DualFun) = DualFun(last(d.f),Evaluation(rangespace(d.J),rightendpoint)*d
 
 jacobian(d::DualFun) = d.J
 jacobian(a::Number) = zero(a)
-jacobian(f::Fun) = Operator(I,space(f))
+jacobian(f::Fun) = zero(Operator(I,space(f)))
 
 promote_rule(::Type{DF},::Type{T}) where {DF<:DualFun,T<:Number}=DualFun
 convert(::Type{DualFun},b::Number) = DualFun(b,0)
diff --git a/test/ExtrasTest.jl b/test/ExtrasTest.jl
@@ -123,8 +123,7 @@ include(joinpath(@__DIR__, "testutils.jl"))
                 u2' - u2*u2;
             ]
 
-        # note takes a few more iterations to converge to accuracy
-        u1,u2 = newton(N_ind, [u1,u2], maxiterations=25)
+        u1,u2 = newton(N_ind, [u1,u2])
 
         u1_exact = -1 / (x - 2)
         u2_exact = -1 / (x + 2)

Original file line number	Diff line number	Diff line change
`@@ -123,8 +123,7 @@ include(joinpath(@__DIR__, "testutils.jl"))`
`123`	`123`	`u2' - u2*u2;`
`124`	`124`	`]`
`125`	`125`
`126`		`- # note takes a few more iterations to converge to accuracy`
`127`		`- u1,u2 = newton(N_ind, [u1,u2], maxiterations=25)`
	`126`	`+ u1,u2 = newton(N_ind, [u1,u2])`
`128`	`127`
`129`	`128`	`u1_exact = -1 / (x - 2)`
`130`	`129`	`u2_exact = -1 / (x + 2)`