Skip to content

Commit 20e6a8b

Browse files
committed
Fix construct_f for fixed point case
Signed-off-by: ErikQQY <2283984853@qq.com>
1 parent 61355db commit 20e6a8b

File tree

1 file changed

+3
-1
lines changed

1 file changed

+3
-1
lines changed

ext/NonlinearSolveSIAMFANLEquationsExt.jl

Lines changed: 3 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -67,7 +67,7 @@ function SciMLBase.__solve(prob::NonlinearProblem, alg::SIAMFANLEquationsJL, arg
6767
end
6868

6969
f!, u = NonlinearSolve.__construct_f(prob; alias_u0,
70-
can_handle_arbitrary_dims = Val(true), make_fixed_point = Val(true))
70+
can_handle_arbitrary_dims = Val(true))
7171

7272
# Allocate ahead for function
7373
N = length(u)
@@ -108,6 +108,8 @@ function SciMLBase.__solve(prob::NonlinearProblem, alg::SIAMFANLEquationsJL, arg
108108
sol = ptcsol(f!, u, FS, FPS; atol, rtol, maxit = maxiters,
109109
delta0 = delta, printerr = ShT)
110110
elseif method == :anderson
111+
f!, u = NonlinearSolve.__construct_f(prob; alias_u0,
112+
can_handle_arbitrary_dims = Val(true), make_fixed_point = Val(true))
111113
sol = aasol(f!, u, m, zeros(T, N, 2*m+4), atol = atol, rtol = rtol,
112114
maxit = maxiters, beta = beta)
113115
end

0 commit comments

Comments
 (0)