Hi,
K-wave is a good works, but some details trouble me, or it is a bug.
codes:
for p_index = 1:length(source.p(:,1))
source.p(p_index, :) = source.p(p_index, :) .* (dt./(c(p_index).*kgrid.dx));
end
in line 616 of the kspaceFirstOrder2D, the aim of the author is :
% compute the scale parameter seperately for each source position
% based on the sound speed at that position
But I think that the p_index to the c matrix can't reach the correctly positon of the source.
For example, when the source located in a circle line, and the struction of the source.p only a two dimension matrix, but the c is a 256*256 matrix, the code c(p_index) maybe not right?
In addtionally, would your can explain the equation of line 616 clearly to me? Thanks very much!
Sincerely,
ECE DUKE,Guoping Chen
22th, Dec., 2011