http://www.k-wave.org/forum/topic/b-mode-imaging-matlab-vs-omp Support for the k-Wave MATLAB toolbox en-US Wed, 13 May 2026 01:13:48 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/b-mode-imaging-matlab-vs-omp#post-7480 Tue, 05 May 2020 17:09:44 +0000 Cyrill 7480@http://www.k-wave.org/forum/ <p>Thanks Brad for the help! With your advise I tracked down that all files were indeed identical (including sensor_data) and the problem was that I had not updated the transducer class on the readout:</p> <p>for angle_index = 1:number_scan_lines<br /> % update the current steering angle<br /> transducer.steering_angle = steering_angles(angle_index);<br /> filename = sprintf('output_orig_%d.h5',angle_index);<br /> sensor_data = h5read(filename,'/p');<br /> sensor_data = transducer.combine_sensor_data(sensor_data);<br /> scan_lines(angle_index, :) = transducer.scan_line(sensor_data);<br /> end</p> <p>Everything works now - thanks!! </p> http://www.k-wave.org/forum/topic/b-mode-imaging-matlab-vs-omp#post-7478 Mon, 04 May 2020 16:36:59 +0000 Bradley Treeby 7478@http://www.k-wave.org/forum/ <p>Hi Cyrill,</p> <p>Can you try the following:</p> <p>1. When you call <code>kspaceFirstOrder3DG</code>, set the optional input <code>&#39;DataPath&#39;</code> to save the input file somewhere you can access easily, and set <code>&#39;DeleteData&#39;, false</code>.</p> <p>2. Use <code>h5compare</code> to compare the input file created by <code>kspaceFirstOrder3DG</code> and the file you create manually.</p> <p>They should be identical (both use <code>kspaceFirstOrder3D</code> to create the input file).</p> <p>Also, have you set values for <code>sensor.record_start_index</code>? If so, this should be passed to the code on the cluster using the <code>-s</code> command line input.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/b-mode-imaging-matlab-vs-omp#post-7466 Fri, 01 May 2020 16:22:10 +0000 Cyrill 7466@http://www.k-wave.org/forum/ <p>I have been building on the example for B-mode Images Using A Phased Array. However, when using the C++ code (either CPU or CUDA) in MATLAB I get inconsistent results with the OMP results. Unfortunately, I don’t get either code to misbehave and even after hours in the forum I’m struggling to make any headway – but at least I learned a lot!!</p> <p>Running the C++ code directly from MATLAB:</p> <p>for angle_index = 1:number_scan_lines<br /> transducer.steering_angle = steering_angles(angle_index);<br /> sensor_data = kspaceFirstOrder3DG(kgrid, medium, transducer, transducer, input_args{:});<br /> scan_lines(angle_index, :) = transducer.scan_line(sensor_data);<br /> end</p> <p>Running the OMP code:</p> <p>for angle_index = 1:number_scan_lines<br /> transducer.steering_angle = steering_angles(angle_index);<br /> filename=sprintf('input_orig_%d.h5',angle_index);<br /> kspaceFirstOrder3D(kgrid, medium, transducer, transducer, input_args{:},'SaveToDisk',filename);<br /> end</p> <p>Run on cluster: kspaceFirstOrder-OMP -i input_orig_##.h5 -o output_orig_##.h5</p> <p>for angle_index = 1:number_scan_lines<br /> filename = sprintf('output_orig_%d.h5',angle_index);<br /> sensor_data = h5read(filename,'/p');<br /> sensor_data = transducer.combine_sensor_data(sensor_data);<br /> scan_lines(angle_index, :) = transducer.scan_line(sensor_data);<br /> end </p>