http://www.k-wave.org/forum/topic/3d-case-with-non-equal-dimensions-fails-to-run Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 22:02:51 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/3d-case-with-non-equal-dimensions-fails-to-run#post-63 Wed, 22 Sep 2010 18:42:04 +0000 Bradley Treeby 63@http://www.k-wave.org/forum/ <p>Hi jkiv,</p> <p>Thanks for your feedback, I'm glad you got your example to run. Please let us know if you have any more questions.</p> <p>Thanks,</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/3d-case-with-non-equal-dimensions-fails-to-run#post-62 Wed, 22 Sep 2010 17:21:20 +0000 jkiv 62@http://www.k-wave.org/forum/ <p>Good day,</p> <p>I'm trying to run a 3D example whose dimensions are non-equal. I'm trying 64x32x16 but during sound-speed/density smoothing the following term fails (in smooth.m):</p> <p><code>f_sm = abs(ifftn(fftshift(fftshift(fftn(f)).*Hann_filt)));</code></p> <p>The issue is <code>Hann_filt</code> is 16x64x32 whereas <code>f</code> is 64x32x16 and thus the <code>.*</code> fails.</p> <p>I might be setting up my data wrong and I'd like assurance from the developers that smooth.m is correct for the 3D case. This is how I'm setting up my data:</p> <pre><code>% physical dimensions X = 1e-2; % m Y = X; Z = X; % resolution Nx = 64; % points Ny = 32; Nz = 16; % dimensions per pixel dx = X/Nx; dy = Y/Ny; dz = Z/Nz; kgrid = makeGrid(Nx, dx, Ny, dy, Nz, dz);</code></pre> <p>and my c, rho, etc. are made to be 64x32x16 (e.g. <code>1500*ones(64,32,16)</code>).</p> <p>Any ideas?<br /> Thanks in advance!</p> <p>EDIT: I re-visited the provided examples for 3D (which I originally dismissed because they were equal dimensions) and the c, rho, etc. are made with <code>Nz x Nx x Ny</code> and not <code>Nx x Ny x Nz</code>. I've done the same and it remedied the situation. Thanks. </p>