search

daniel-koehn / DENISE-Black-Edition

2D time-domain isotropic (visco)elastic FD modeling and full waveform inversion (FWI) code for P/SV-waves
GNU General Public License v2.0
121 stars 66 forks source link

Bug in acoustic l-bfgs optimization #9

Open daniel-koehn opened 7 years ago

daniel-koehn commented 7 years ago

Currently, acoustic FWI (MODE=1, PHYSICS=2) does only converge, when using PCG optimization (GRAD_METHOD=1). The l-bfgs optimization implementation (GRAD_METHOD=2) does not converge.

I suspect there is a bug when using an even number of parameter classes, because the l-bfgs implementation works in GERMAINE for the mono-parameter Vp inversion and in DENISE with 3-parameter Vp, Vs, Density inversion.