LSR1Operator(1) fails the secant equation

Aperta
  • Qualcosa non va?
#225 6 commenti 0 reazioni 0 assegnatari Vedi su GitHub

@paraynaud ci sta già lavorando.

Dal 24/6/2022.

  • #235 di @paraynaud — aperta

Valutazione

Difficoltà
3/5
Tempo stimato
1-2 giorni
Idoneità per principianti
45/100
Tipo di issue
Bug
Chiarezza
Abbastanza chiara
Stato di attività
Ferma
Stack tecnologico
julia
Ambito
tooling

Direzione di ricerca

Inizia eseguendo la riproduzione in Julia dell’issue e individua l’implementazione di LSR1Operator, in particolare il comportamento attivato da push!(B, s, y) quando n=1. Confronta il suo risultato con i casi n=2 e LBFGSOperator; il lavoro è completato quando aggiornamenti unidimensionali ripetuti soddisfano l’equazione secante con un residuo paragonabile al risultato n=2 riportato.

Scritto dal modello di indicizzazione a partire dal testo della issue.

Descrizione

Hello all,

I did test the LSR1 operators for myself, and it seems that there is a problem with the LSR1Operator of dimension 1.
In my work, I generate automatically several LSR1Operators, and some of them are of dimension 1.

Here is an example to replicate the issue:

n = 1
B = LSR1Operator(n)
r = 0.0
for i = 1:100
  global r
  s = rand(n)
  y = rand(n)
  push!(B, s, y)  
  r += norm(B * s - y) #secant equation
end
r #166.39...

If you replace n=1 by n=2 then r = 1e-14.

The problem does not occur with the LBFGSOperators.

Lingua principale
Julia
Stelle
214
Fork
46
Nessuna PR unita negli ultimi 30g

Preparare l'ambiente

Come iniziare

  1. Leggi tutta la issue e poi la guida ai contributi del progetto.
  2. Commenta sulla issue per dire che te ne occupi tu — evita che due persone facciano lo stesso lavoro.
  3. Fai un fork del repository e lavora su un branch.
  4. Apri una pull request che faccia riferimento al numero della issue.

Altre issue di JuliaSmoothOptimizers/LinearOperators.jl

Tutte le issue di JuliaSmoothOptimizers/LinearOperators.jl

Issue simili

Altre issue su Julia

Ricevi le nuove issue nella tua casella

Un breve riepilogo di issue GitHub adatte ai principianti.