http://www.k-wave.org/forum/topic/isppa-and-ispta-computation Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 22:29:20 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/isppa-and-ispta-computation#post-8668 Thu, 17 Nov 2022 16:05:04 +0000 federico.rossano 8668@http://www.k-wave.org/forum/ <p>Hi Ben,</p> <p>Thank you for your answer and for the amazing work.</p> <p>I would just only clarify the point about p_max. In the k-Wave User Manual page 38, I see:<br /> [...] The exceptions are the averaged quantities (‘p_max’, ‘p_rms’, ‘u_max, ‘p_rms, ‘I_avg’) and the quantities recorded over all the grid points within the domain (‘p_max_all’, ‘p_min_all’,‘p_final’, ‘u_max_all’, ‘u_min_all’, ‘u_final’). The averaged quantities return the maximum, average, or root-mean-squared values at each sensor point for the complete simulation and are indexed as (sensor_point_index). </p> <p>For "averaged quantities" you mean averaged on time? As for p_max, does k-Wave take the maximum pressure in each voxel for each time step and then the result is the average of these maximum pressures recorded over the time vector? Or it just takes the maximum pressure recorded at any time during the simulation?</p> <p>Thank you so much again,<br /> Besti wishes,<br /> Federico </p> http://www.k-wave.org/forum/topic/isppa-and-ispta-computation#post-8623 Fri, 07 Oct 2022 13:54:32 +0000 bencox 8623@http://www.k-wave.org/forum/ <p>Hi Federco,</p> <p>There are various things you could try. </p> <p>kspaceFirstOrder3D will directly output average intensity <code>I_avg</code> as well as the time-varying intensity <code>I</code> from which you could compute the average. You can ask it to do that by setting <code>sensor.record = {&#39;I_avg&#39;}</code> or <code>sensor.record = {&#39;I}</code>. </p> <p>Beware there may be a transient at the start of a signal which may have a higher peak than the steady state, so you may want to consider setting <code>sensor.record_start_index</code> to avoid a possible transient.</p> <p>If your wave is close enough to a single frequency plane wave in a steady state then you could also find the pressure amplitude and use the formula you gave above for the time-averaged intensity. You could use <code>p_max</code> if you just look at the steady state part of the propagation by avoiding the transient. (<code>p_max</code> is the maximum pressure at the defined sensor points, not everywhere.)</p> <p>If you are recording on very many points, eg. at all points in the domain, and that means the memory require is still very high, could you reduce the region over which you are recording? Or break it into subregions and run one simulation per subregion?</p> <p>Best wishes<br /> Ben </p> http://www.k-wave.org/forum/topic/isppa-and-ispta-computation#post-8616 Mon, 03 Oct 2022 10:47:24 +0000 federico.rossano 8616@http://www.k-wave.org/forum/ <p>Hi,</p> <p>I would like to ask if there is a strategy to calculate Isppa and Ispta values without integrating over the pulse duration the square of p(t). I have problems due to RAM memory with computations of both big domains and time vectors, so I am unable to compute sensor.p to achieve the I(t) through p^2/2*rho*c.</p> <p>I've already considered to start recording my signal in the steady state, but it still requires a lot of memory. Maybe I can calculate Isppa/Ispta with p_max or p_max_all?<br /> And what is the difference between the latter two?</p> <p>Thank you in advance for the help!<br /> Federico </p>