http://www.k-wave.org/forum/topic/error-with-gpumat Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 22:00:14 +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-with-gpumat#post-27 Wed, 21 Jul 2010 07:47:28 +0000 bISHOP 27@http://www.k-wave.org/forum/ <p>Thank you. Interpolating to nearest has worked.</p> <p>I think I´m going to be a regular user of this forum with exclusive attention.</p> <p>Great job!</p> <p>Daniel. </p> http://www.k-wave.org/forum/topic/error-with-gpumat#post-26 Tue, 20 Jul 2010 13:07:34 +0000 Bradley Treeby 26@http://www.k-wave.org/forum/ <p>Hi Daniel,</p> <p>Thanks for your feedback and post. Unfortunately, the function which calculates the pressure at sensor points which are defined in Cartesian space using linear interpolation (<code>gridDataFast</code>) does not currently support execution on the GPU. Hopefully this will be fixed by the next release. To get around this in the meantime, you can either set the interpolation mode to nearest neighbour, e.g., </p> <p><code>sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, &#39;DataCast&#39;, &#39;GPUsingle&#39;, &#39;CartInterp&#39;, &#39;nearest&#39;);</code></p> <p>or you can use a binary sensor mask as in the second photoacoustic example. </p> http://www.k-wave.org/forum/topic/error-with-gpumat#post-25 Tue, 20 Jul 2010 12:12:34 +0000 bISHOP 25@http://www.k-wave.org/forum/ <p>Hi.</p> <p>I´m sorry about being first to post... an error, but I´m glad to help with this fantastic tool you have developed. I´ve followed the steps to install CUDA, GPUmat and k-Wave and everything seemed to be allright until I tried the homogeneus propagation medium (the first one, at <a href="http://www.k-wave.org/documentation/example_ps_homogeneous_medium.php" rel="nofollow">http://www.k-wave.org/documentation/example_ps_homogeneous_medium.php</a> ) with 'Datacast', 'GPUsingle' in the kspaceFirstOrder2D call. MATLAB is not happy. MATLAB disagrees with my code this way:</p> <pre><code>Running k-space simulation... dt: 20ns, t_end: 12.06us, time steps: 604 input grid size: 128 by 128 pixels (12.8 by 12.8mm) maximum supported frequency: 7.5MHz smoothing p0 distribution... calculating Delaunay triangulation... casting variables to kWaveGPUsingle type... precomputation completed in 2.2222s starting time loop... ??? Error using ==&gt; times Operations between Matlab arrays and GPU arrays are not supported. (ERROR code 5) Error in ==&gt; gridDataFast at 96 zi = sum(z(tri) .* w,2); Error in ==&gt; kspaceFirstOrder2D at 1194 sensor_data(:, t_index) = gridDataFast(kgrid.z, kgrid.x, p, sensor_z, sensor_x, del_tri); Error in ==&gt; ej211 at 40 sensor_data = kspaceFirstOrder2D(kgrid, medium, source, sensor, &#39;DataCast&#39;, &#39;GPUsingle&#39;);</code></pre> <p>The only modifications are two addpath lines at the beggining and the DataCast parameter. </p> <p>I´m using a Intel Core2 Quad (Q6600, 2.4Ghz) with 2GB RAM, Windows XP SP3, MATLAB 7.7.0 (R2008b), CUDA 3.1, GPUmat 0.251 and k-Wave Version B.0.2</p> <p>Thank you and I hope it is an error of mine.</p> <p>Daniel. </p>