JuliaNLSolvers / NLsolve.jl

Julia solvers for systems of nonlinear equations and mixed complementarity problems

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

`nlsolve` allocates a Jacobian even when using `anderson`

antoine-levitt opened this issue · comments

The basic nlsolve(f,x0,method=:anderson) allocates a dense jacobian. This breaks large systems where a dense jacobian doesn't fit into memory.

While we're at it: NLsolve doesn't support NonDifferentiable.

Workaround: trick OnceDifferentiable into thinking you have a custom jacobian: OnceDifferentiable(F!,identity,x0,x0,[])

Thanks, ill fix it unless you do :)

I'd rather you do it if it's easy for you :-) I'm not very comfortable around the NDifferentiable code.

Ill do it :)

#203 . It’s just a quick “do something that works”, but I’ll come back to it. But it works I think ;)

Many thanks, and merry christmas :-)

You too! (have a merry christmas that is)