http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:02:11 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time#post-6876 Fri, 07 Jun 2019 12:44:13 +0000 Jiri Jaros 6876@http://www.k-wave.org/forum/ <p>Hi laughingrice,</p> <p>first of all you can use the parameter -s to start sampling at a timestep you need (not the whole timespan of the simulation).</p> <p>second, there's no decimation in the C++/CUDA code, however, it should not be so difficult to implement it by our own if you know a little bit about C++. If so, you can open KSpaceFirstOrder3DSolver.cpp and find a routine storeSensorData(). Now, if you want to sample every other timestep, just add a modulo operation here. The other change will be in IndexOutputStream::create() or CuboidOtputStream and you change the number of timesteps at this line</p> <p> // Derive dataset dimension sizes<br /> DimensionSizes datasetSize(nSampledElementsPerStep,<br /> (mReduceOp == ReduceOperator::kNone) ?<br /> params.getNt() - params.getSamplingStartTimeIndex() : 1,<br /> 1);</p> <p>You you want to do averaging over multiple time steps, you would have to reimplement the sample() method. </p> http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time#post-6845 Thu, 09 May 2019 19:23:52 +0000 laughingrice 6845@http://www.k-wave.org/forum/ <p>Also, you mention having to use the MATLAB interface, I'm running simulations in 2D right now, I don't believe that there is a non-matlab interface to that </p> http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time#post-6844 Thu, 09 May 2019 15:11:27 +0000 laughingrice 6844@http://www.k-wave.org/forum/ <p>Thanks Brad</p> <p>I am using CUDA so it will help a bit, but things are borderline in CPU memory as well (I'm still looking at around 40GB). Where are things processed inside the k-wave code? Wondering it it's easy enough to just decimate by a factor of two with no filtering (i.e save every second step). It usually works well enough for RTM imaging. </p> http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time#post-6843 Wed, 08 May 2019 21:31:57 +0000 Bradley Treeby 6843@http://www.k-wave.org/forum/ <p>Are you using the native CUDA code? This should stream the data from the GPU memory by default. There’s also the ‘StreamToDisk’ optional input if for some reason you need to use the MATLAB version.</p> <p>Unfortunately there’s currently no way to decimate the data, although it is on the development list.</p> <p>Brad </p> http://www.k-wave.org/forum/topic/subsampling-the-recorded-pressure-field-in-time#post-6842 Wed, 08 May 2019 21:25:00 +0000 laughingrice 6842@http://www.k-wave.org/forum/ <p>When setting sensor.record = {'p'} to record the pressure field at given grid positions, is it possible to subsample the recorded data in time? (preferably with a FIR decimation filter applied)</p> <p>I need to record the time dependent pressure field solution over most of the domain, but doing that at the full time resolution is both not needed (the nyquist frequency is significantly lower than the simulation frequency) and impossible due to memory requirements (about 48GB on a 10GB GPU).</p> <p>I couldn't find any way to do that in the documentation. I also tried to figure out if there is a way to checkpoint the simulation to offload data, but that doesn't seem possible either.</p> <p>Thanks </p>