I've been trying to match the k-wave simulation to a standard finite difference acoustic simulation with a Ricker wavelet forcing function.

i.e solve u_tt - c^2 u_xx = f

As Ricker is the second order time derivative of a Gaussian, the source pressure is a Gaussian. I tried using that for source.p which did not produce the same results. Looking at the governing equations defined in eq 2.8 in the user manual, what is defined as the pressure source term matches the mass source in that equation (d\rho/dt), and using that (the normalized time derivative of the Gaussian pressure source) does match the result of the finite difference simulations.

I wanted to verify what does source.p actually define, a pressure source term as in the function help file and section 3.4 of the manual or a mass source term (flow of mass) as defined in section 2.2 of the manual.

Thanks