2.5D g0 time domain - frequency domain and theory section - reference equation mismatch

[TorosSenan]

Dear SFS toolbox creators,

Thanks for such an excellent contribution. I am learning a lot by following the theory and the MATLAB scripts.

While digging in, I got confused with the g0 definition in MATLAB script driving_function_mono_wfs_ps.m, line 209 (sqrt(dref / (dref + ds)), which is also stated in the Eq. 9.11. However, [Sch16, Eq. 2.141] and [Sta97, 3.16], which are cited in the theory section of SFS (just above the Eq. 9.11), define the g0 for a line parallel to the linear secondary sources as sqrt(dref / (dref-ds)). In driving_function_imp_wfs_ps.m, the line 166 also defines g0 following the [Sta97, 3.16], sqrt(dref / (dref - ds)).

The Eq. 9.11 is followed by the statement "...which is the traditional formulation of a point source in WFS as given by eq. (2.27) in [Ver97]" which is exactly correct. In [Ver97, Eq. 2.29b], sqrt(dref / (dref-ds)) is used for focused sources, as it is also the case in sfs-toolbox MATLAB implementation both for frequency and time-domain calculations of focused sources.

So, which MATLAB script in sfs-toolbox (frequency or time domain) calculates the 2.5d operator of a point source correctly for a reference line? Or what am I missing?

Thank you again for your efforts.