http://www.k-wave.org/forum/topic/error-in-hd5-read-when-using-a-plane-wave-source-with-c-version Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:01:47 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/error-in-hd5-read-when-using-a-plane-wave-source-with-c-version#post-5654 Wed, 31 Aug 2016 00:13:55 +0000 scottalm 5654@http://www.k-wave.org/forum/ <p>The problem ended up being how I was defining the source (along the time instead of the spatial dimension):</p> <p>source.uz = reshape(source.uz, [], 1);</p> <p>Should have been:</p> <p>source.uz = reshape(source.uz, 1, []);</p> <p>The Matlab code will execute this, but the C++ code won't (since it is improperly constructed). It would be nice if the Matlab code would do a check for dimension sizes, but at least there is a forum post if someone comes across this problem. </p> http://www.k-wave.org/forum/topic/error-in-hd5-read-when-using-a-plane-wave-source-with-c-version#post-5604 Wed, 27 Jul 2016 09:36:29 +0000 Bradley Treeby 5604@http://www.k-wave.org/forum/ <p>Hi scottalm,</p> <p>If you add a second time point to your source, e.g., something like</p> <p><code>source.uz(:, 2) = 0</code>;</p> <p>Does it work? I would need to check, but it's possible the C++ code explicitly checks for a 2D input signal.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/error-in-hd5-read-when-using-a-plane-wave-source-with-c-version#post-5595 Tue, 26 Jul 2016 10:12:56 +0000 Jiri Jaros 5595@http://www.k-wave.org/forum/ <p>Hi scottalm,</p> <p>I think the matlab script is not correct - you mixed the source spatial definition (position and size) with the driving signal definition.</p> <p>source.u_mask - this should be the spatial definition of the transducer<br /> source.uz - this is the driving signal</p> <p>I think you tried to define the shape of the source in the uz component of the driving signal.<br /> Jiri </p> http://www.k-wave.org/forum/topic/error-in-hd5-read-when-using-a-plane-wave-source-with-c-version#post-5592 Fri, 22 Jul 2016 00:40:08 +0000 scottalm 5592@http://www.k-wave.org/forum/ <p>I am using the SSE3 Linux binaries for this work, but I have seen the same error in the Windows AVX2 binaries. </p> <p>I'm not sure if this is a real bug or just me missing something somewhere, but it seems that when I try to use a plane wave source with the C++ version of k-wave I get an error reading in the HD5 file. </p> <p>As an example, running the following works:</p> <pre><code>% create the computational grid Nx = 108; Ny = 108; Nz = 108; dx = .05; dy = .05; dz = .05; kgrid = makeGrid(Nx, dx, Ny, dy, Nz, dz); model = ones(Nx, Ny, Nz); % density rho_v=[1e3]; % speed of sound c_v=[1500]; medium.sound_speed = c_v(model); % [m/s] medium.density = rho_v(model); % [kg/m^3] medium.alpha_coeff = 0; % [dB/(MHz^2 cm)] medium.alpha_power = 0.999; sensor.mask = ones(Nx, Ny, Nz); % record the maximum pressure sensor.record = {&#39;p_max&#39;}; source.u_mask = ones(Nx, Ny, Nz); source.uz = zeros(Nx, Ny, Nz); source.uz(:, :, 3) = 1; source.uz = smooth(kgrid, source.uz, true); source.uz = reshape(source.uz, [], 1); input_args = {&#39;PlotScale&#39;, &#39;auto&#39;, ... &#39;PMLInside&#39;, false, &#39;PlotPML&#39;, false}; sensor_data = kspaceFirstOrder3D(kgrid, medium, source, sensor, input_args{:})</code></pre> <p>But when I replace <code>kspaceFirstOrder3D</code> with <code>kspaceFirstOrder3DC</code>, it errors out giving me</p> <pre><code>K-Wave panic: Data loading was not successful! Error: Dataset &quot;uz_source_input&quot; has wrong dimension sizes!</code></pre> <p>The size is appropriate for Matlab, but doesn't seem to work for C++. Am I doing something incorrectly? </p>