search

JuliaNLSolvers / NLsolve.jl

Julia solvers for systems of nonlinear equations and mixed complementarity problems
Other
330 stars 66 forks source link

Add linsolve choice. #201

Closed pkofod closed 5 years ago

pkofod commented 5 years ago

Seems to work fine with gmres! and idrs!. Other examples from other (?) packages are welcome as well.

ChrisRackauckas commented 5 years ago

Is it possible to swap out the Jacobian type? I would like to see a test that does gmres with a matrix-free representation of the Jacobian given by a LinearMaps.jl FunctionMap.

pkofod commented 5 years ago

Is it possible to swap out the Jacobian type? I would like to see a test that does gmres with a matrix-free representation of the Jacobian given by a LinearMaps.jl FunctionMap.

That should be possible, yes. Let me see if I can find the appropriate types you're talking about and come up with an example (if you have a link handy that shows how it works I'd be happy to see it)

pkofod commented 5 years ago

That should be possible, yes. Let me see if I can find the appropriate types you're talking about and come up with an example (if you have a link handy that shows how it works I'd be happy to see it)

Yeah it probably isn't ready right now. The problem is that some internals assume .= will work, and I suspect you're asking for some sort of immutable storage type, righT?

ChrisRackauckas commented 5 years ago

oh okay. But things like sparse work?

pkofod commented 5 years ago

If you’re willing to reuse the object, sure. What sparse problem could benefit from this? Are you thinking calling some nonjulia library?

pkofod commented 5 years ago

This works so I'm putting it in. I'm going to collect some examples and maybe make some corrections to better handle immutable jacobians before tagging a version.