http://www.k-wave.org/forum/topic/save-only-uz_non_staggered Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:28:59 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4989 Thu, 19 Feb 2015 19:41:46 +0000 Bradley Treeby 4989@http://www.k-wave.org/forum/ <p>Hi Anthony,</p> <p>Yes we did consider this for the calculation of I_avg or Q for CW fields. It would involve the user having to define the length of the acoustic period over which to average, and preferably defining the time discretisation to use an integer number of points per period. However, it would save a lot of memory as you say. Any other suggestions most welcome!</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4980 Thu, 05 Feb 2015 09:07:09 +0000 Anthony 4980@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>Thanks for this important remark (and the pointer to the kspaceFirstOrder_saveIntensity.m function) :-) By the way, I hadn't seen it anywhere in the doc, maybe it would be worth adding it or mentioning it in the release notes, where you speak about "the additional input options".</p> <p>If you want to perform only heat deposition calculation, maybe you could only unstagger the last period as you basically only want to have a mean deposited power over one period. It would save a lot of memory...</p> <p>Anthony </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4979 Thu, 05 Feb 2015 00:07:33 +0000 Bradley Treeby 4979@http://www.k-wave.org/forum/ <p>Hi Anthony,</p> <p>You are correct - the output <code>&#39;u_non_staggered&#39;</code> is only non-staggered in space, it is still staggered in time. To calculate intensity within the MATLAB code, this is then passed to the <code>timeShift</code> function before being multiplied by the pressure. You can take a look inside the private function <code>kspaceFirstOrder_saveIntensity.m</code> to see exactly what is done. Once we figure out an efficient and general way to do it, we'd like to put the calculate of intensity (and heat deposition) back inside the C++ code. I agree it would be more convenient this way.</p> <p>Brad. </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4971 Wed, 04 Feb 2015 16:49:13 +0000 Anthony 4971@http://www.k-wave.org/forum/ <p>Hi Jiri,</p> <p>Thanks a lot for your reply: the solution you describe seems to work according to my first tests.</p> <p>I understand the problem of memory and it is definitely a good reason not to allow for computing the intensity directly in the simulation code. </p> <p>You answer just raised a question in my mind: I thought that "u_unstaggered" was unstaggered both in space and time but your post let me think that it is only unstaggered in space, isnt'it? </p> <p>Actually, the reason why I only want to save uz is that I compute the integral of Iz over each XY plane to check the evolution of the power along the propagation (so I can check the increase in absorption due non linear effects around the focus for example...). </p> <p>Cheers,<br /> Anthony </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4970 Tue, 03 Feb 2015 20:28:18 +0000 Jiri Jaros 4970@http://www.k-wave.org/forum/ <p>To answer the second question "Calculate intensity in the C code", I have to reply that it is impractical due to memory limitations.</p> <p>As you may know, p and u are staggered both in space and time to improve accuracy (Brad could explain that in detail). In k-Wave 1.0, we used a linear interpolation to get u at the same grid as p, however it had to low accuracy. </p> <p>In k-wave 1.1 we introduced fft-based shifts that improved accuracy. The down side is that to do so in space, you need all spatial data for ux, uy and uz even if you sample 1 single grid point! Thus to sample even a tiny region you need 6 more matrices! (speaking about a domain 512^3, we need additional 3GBs of RAM). Let me say tat This is the easy part.</p> <p>Now it starts to be tricky, to do the same shift in TIME we need the whole time history - over the sampling period. I talk here about tens to hundreds of GBs of RAM...:-)</p> <p>The only solution for now is to progressively store the non-staggered velocity to file, then terminate the simulation, load the data in a different order (1 gridpoint at time, but at all time steps), run the FFT shift and store the data back to file. - this is not implemented in the C code now, and you have to do that in Matlab (I agree that's not very convenient).</p> <p>Cheers<br /> Jiri </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4969 Tue, 03 Feb 2015 20:08:39 +0000 Jiri Jaros 4969@http://www.k-wave.org/forum/ <p>Hi Hypermoi,</p> <p>To only store uz_non_staggered you need to comment out a couple of lines at different places in the C code.</p> <p>This is not a clear solution</p> <p>1. MatrixContainer.cpp<br /> You have to inspect the `TMatrixContainer::AddMatricesIntoContainer() routine and remove these lines<br /> <pre><code>MatrixContainer[ux_shifted].SetAllValues(....); MatrixContainer[uy_shifted].SetAllValues(....);</code></pre> <p>This removes the matrices from the Matrix container (prevents allocation of memory)</p> <p>2. <code>TOutputStreamContainer::AddStreamsIntoContainer()</code><br /> Here, remove these two lines<br /> <pre><code>OutputStreamContainer[ux_shifted_sensor_raw] = CreateNewOutputStream(...) OutputStreamContainer[uy_shifted_sensor_raw] = CreateNewOutputStream(...)</code></pre> <p>which will stop sampling the quantities. </p> <p>3. You need to disable calculating the x and y components</p> <p>Inspect <code>TKSpaceFirstOrder3DSolver::Calculate_shifted_velocity()</code> and remove the code for X and Y</p> <p>If it does not work, let me know and I'll send you a patch.</p> <p>I might also add a command line parameter to let you choose which components to sample, however we thought it would not be of any use :) </p> http://www.k-wave.org/forum/topic/save-only-uz_non_staggered#post-4967 Tue, 03 Feb 2015 10:32:20 +0000 Anthony 4967@http://www.k-wave.org/forum/ <p>Hi Brad,</p> <p>I would like to save only uz_non_staggered in my result file because don't need the other components and u_non_staggered is simply too big. </p> <p>I guess there is juste one or two lines to comment in the source code but I can't figure out precisely what to do. My hypothesis is that it should work if I comment the lines 705 to 714 in MatrixContainer.cpp but I am not familiar enough with the methods you use to write the HDF5 files...</p> <p>Best regards,<br /> Anthony</p> <p>PS: by the way, it was convenient to be able to save directly intensity in k-wave 1.0 (instead of pressure and velocity), I would love to see it back in 1.2 :-) </p>