`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)