http://www.k-wave.org/forum/topic/correction-to-pml-implementation Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:32:31 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/correction-to-pml-implementation#post-402 Fri, 23 Mar 2012 22:53:58 +0000 bencox 402@http://www.k-wave.org/forum/ <p>Well spotted! And thanks for writing in. Glad you're enjoying it. Keep the comments and bug reports coming please!</p> <p>Ben </p> http://www.k-wave.org/forum/topic/correction-to-pml-implementation#post-398 Fri, 23 Mar 2012 11:14:11 +0000 Jonathan Hargreaves 398@http://www.k-wave.org/forum/ <p>Hello! This release is very well timed as I was about to write to you with a bug report but you've fixed it. Here's the info anyway for the benefit of anyone else experiencing the same issue:</p> <p>I’d been comparing the code to Tabei et al’s 2002 JASA paper and realised there was a mismatch, specifically that the following lines of code don’t match their equation 26:</p> <p>kspaceFirstOrder3D.m version B.0.4 lines 1032 – 1035:<br /> <code>rhox = bsxfun(@times, pml_x, rhox - dt.*rho0 .* duxdx);</code><br /> <code>rhoy = bsxfun(@times, pml_y, rhoy - dt.*rho0 .* duydy);</code><br /> <code>rhoz = bsxfun(@times, pml_z, rhoz - dt.*rho0 .* duzdz);</code></p> <p>I wasn’t getting the PML attenuation they reported so I suspected there may be a problem. And low and behold in the new version you've fixed it! :-)</p> <p>kspaceFirstOrder3D.m version B.0.5 lines 812 – 814:<br /> <code>rhox = bsxfun(@times, pml_x, bsxfun(@times, pml_x, rhox) - dt.*rho0 .* duxdx);</code><br /> <code>rhoy = bsxfun(@times, pml_y, bsxfun(@times, pml_y, rhoy) - dt.*rho0 .* duydy);</code><br /> <code>rhoz = bsxfun(@times, pml_z, bsxfun(@times, pml_z, rhoz) - dt.*rho0 .* duzdz);</code></p> <p>Many thanks for your hard work in creating this excellent resource and your generosity in sharing and disseminating it.</p> <p>Jon </p>