vkazei / fastfwi Goto Github PK

Possible solution
add opts.scale_model -- flag
modify fwiFunc so that it rescales the model if needed

Hi,

I'm trying to run the "marmousi_II_FWI.m" example, but it shows that:

Undefined function or variable 'F'.

Error in marmousi_II_FWI (line 140)
        DClean = F(m.Base,model);

Could you tell me how to solve this problem? Thanks!

Dear Dr Vladimir Kazei

I use your published program fastFWI to run some simulations. And compare received data with 2D analytical results, modeling frequency is 90 kHz. The first model I used is a 1 meter by 1 meter homogeneous model, with transducers deployed along a circle and each transducers acted as both transmitters and receivers

The real part of received data at transducer locations obtained by finite difference modeling (FDM) and 2D analytical green functions are shown as follows

It can be seen that FDM data and analytical results matched very well.

However, when I increase the size of model to 2 meter (x direction) by 3 meter (y direction). And use parallel array setup, as shown in the following figure (grid spacing 2mm)

The real part of FDM wavefield and analytical wavefield (obtained by green functions) are shown as follows

It shows that with the increase of model size, the analytical results (homogeneous model) and FDM results begin to differ significantly. Such discrepancy influence the imaging results, leading to the algorithm failure at some circumstances. But I cannot figure out the possible reasons.

What is the possible reason that causes discrepancy between the FDM modeling wavefield and analytical results? Is such discrepancy caused by modeling grid size or some other reasons?

Hope to hear from your answers!

未定义函数或变量 'lbfgsb_wrapper'。

出错 lbfgsb (line 175)
lbfgsb_wrapper( m, x, l, u, nbd, f, g,
factr, pgtol, wa, iwa, task,iprint,...

出错 fwiFunc (line 61)
[mk, farsh, info] = lbfgsb(fh, l, u, opts );

出错 BP_tooth_FWI (line 159)
fwiResult = fwiFunc(iFWI, m.Init, D,
model, opts);

when I run the matlab code, it inform me cannot find the function F in line 134, is there need some other package?
Ps: the file is marmousi_ii_FWI.m