QUELLTYPB=4 fails to do FWI in AC problem

[ADharaUTEXAS123007]

When I do SEISMO=2, QUELLTYPB=4, The LBFGS method fails to find the step length.

[daniel-koehn]

Does this happen in case of the Marmousi-2 benchmark problem or did you use another model? I tested the parameter combination (SEISMO=2,QUELLTYPB=4, LBFGS) for the Marmousi-2 model extensively, as discussed in this thread:

#23 (comment)

Additionally, you have to set GRAD_FORM=2 in the DENISE parameter file to use a specific type of acoustic FWI gradients. In this case, the updates of the density model also seem to be quite crappy. So I would recommend to update only the P-wave velocity model as discussed in the thread.

Best regards,

Daniel

[ADharaUTEXAS123007]

Thanks Daniel.
This worked as expected.

But I want to bring to your attention in file psv.c, line 180
/* adjoint explosive source */
if ((QUELLTYPB >= 4) && (mode == 1))
{
psource(nt, (*wavePSV).psxx, (*wavePSV).psyy, (*acq).srcpos_loc_back, (*seisPSVfwi).sectionpdiff, nsrc_loc, 1);
}

it is still nsrc_loc instead of ntr. Is this as expected?

[ADharaUTEXAS123007]

closing this issue.