http://www.k-wave.org/forum/topic/attencompm-error Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:49: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/attencompm-error#post-8248 Wed, 21 Jul 2021 14:17:10 +0000 Bradley Treeby 8248@http://www.k-wave.org/forum/ <p>Hi Damien,</p> <p>Yes, this is a bug - thanks very much for reporting. I'll include the fix with next release. Fortunately, at least for the few examples I've tried, it doesn't have a huge impact as the values in the non-stationary convolution matrix are pretty small away from the diagonal anyway.</p> <p>Brad </p> http://www.k-wave.org/forum/topic/attencompm-error#post-8141 Tue, 27 Apr 2021 10:12:03 +0000 Damien GASTEAU 8141@http://www.k-wave.org/forum/ <p>Hello,</p> <p>I am currently working with the attenComp.m function and noticed on line 510 that the lower triangle coefficients are set to zeros twice.</p> <p>% zero out lower and upper triangles<br /> ones_mat = ones(N, N);<br /> tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;<br /> tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;</p> <p>Shouldn't it be something like:</p> <p>% zero out lower and upper triangles<br /> ones_mat = ones(N, N);<br /> tv_filter(tril(ones_mat, -ceil(N/2) + 1) == 1) = 0;<br /> tv_filter(triu(ones_mat, ceil(N/2) + 1) == 1) = 0; </p>