JuliaNLSolvers / NLsolve.jl

Julia solvers for systems of nonlinear equations and mixed complementarity problems

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Field Initial Point is showing content of zero

BeastyBlacksmith opened this issue · comments

When I take a look at the ouput object of nlsolve the fields initial_x and zero show the same data, namely the content of zero.

Could you be more specific? I mean, what solver etc.

That showed up with the default setting. So it is method = :trust_region, don't think the other values matter but for completeness: iterations = 1_000, ftol = 1e-10, xtol = 0.0
and

versioninfo()
Julia Version 1.0.3
Commit 099e826241 (2018-12-18 01:34 UTC)
Platform Info:
  OS: Linux (x86_64-pc-linux-gnu)
  CPU: Intel(R) Core(TM) i7-3770 CPU @ 3.40GHz
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-6.0.0 (ORCJIT, ivybridge)
Environment:
  JULIA_EDITOR = atom -a
  JULIA_NUM_THREADS = 4

Okay... Are they the same array or do they just hold the same values? Are you sure it's not just that it converges in 0 iterations? Try changing an element of res.zero and see if the initial point field changes as well.

They just hold the same values. I can see that it iterates because I know the true initial values and also see the number of function evaluations in the output be around 50.

o.initial_x === o.zero
false

but

all( o.initial_x .== o.zero )
true

I'll look into it, but I don't have time right now. Thanks.