http://www.k-wave.org/forum/topic/question-regarding-multiple-sensorrecord-within-one-simulation Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 23:06:40 +0000 http://bbpress.org/?v=1.0.2 q http://www.k-wave.org/forum/search.php http://www.k-wave.org/forum/topic/question-regarding-multiple-sensorrecord-within-one-simulation#post-8581 Thu, 14 Jul 2022 08:53:22 +0000 federico.rossano 8581@http://www.k-wave.org/forum/ <p>Hi Ben,</p> <p>Thank you a lot for your answer. Luckly, in the meantime, I was able to run large simulations without running out of RAM just dividing the sensor.record matrix in different parts and then using 'cat' Matlab function to concatenate all the "result matrices parts".<br /> I will attach my code here:</p> <p>parts=n; %number of split in the domain along z<br /> for i=1:parts %domain splitted in half along z<br /> x_1=1;<br /> x_2=Nx;<br /> y_1=1;<br /> y_2=Ny;<br /> z_1=floor((Nz*i)/parts)-(floor(Nz/parts))+1;<br /> z_2=ceil(Nz/parts)+ceil((Nz*(i-1))/parts);<br /> if z_2&gt;ceil(Nz*i/parts)<br /> z_2=ceil(Nz*i/parts);<br /> end<br /> sensor_cuboid=[x_1,y_1,z_1,x_2,y_2,z_2].';<br /> sensor.mask=sensor_cuboid;<br /> sensor.record={'p_max'};<br /> %Starting simulation<br /> reset(gpuDevice(1)); %reset GPU memory<br /> input_args={ 'DataCast','single','PMLInside', false};<br /> sensor_data = kspaceFirstOrder3D(kgrid, medium, source, sensor, input_args{:},'PlotSim', false);<br /> data=struct('P_max',sensor_data.p_max,'Nx',Nx,'Ny',Ny,'Nz',Nz,'dx',dx,'dy',dy,'dz',dz,'signal',signal);<br /> save(['your_folder\data_', num2str(i), '.mat'],'data','-v7.3','-nocompression');<br /> clear data<br /> clear sensor.record<br /> clear sensor.mask<br /> end</p> <p>Surprisingly the division of the domain led to less computation time.<br /> The script cannot divide the domain in exact parts and may be an overlap in the results. That will be easy to correct, just deleting the exact same slices running a for cycle.</p> <p>Best regards,<br /> Federico </p> http://www.k-wave.org/forum/topic/question-regarding-multiple-sensorrecord-within-one-simulation#post-8574 Wed, 13 Jul 2022 17:56:58 +0000 bencox 8574@http://www.k-wave.org/forum/ <p>Hi Federico, </p> <p>You can only define one <code>sensor.mask</code> region using the 'cuboid corners' approach, but you can define <code>sensor.mask</code> as a binary source mask with any set of points. If you want to record very many time series, then the output will get quite large. It may be worth thinking about whether you need the full time series, or whether you could define another quantity to be recorded, eg. by using <code>sensor.record = {&#39;p_max&#39;}</code>.</p> <p>Best wishes<br /> Ben </p> http://www.k-wave.org/forum/topic/question-regarding-multiple-sensorrecord-within-one-simulation#post-8557 Wed, 29 Jun 2022 09:46:31 +0000 federico.rossano 8557@http://www.k-wave.org/forum/ <p>Hi,<br /> I need to run simulations on very large 3D domains using a frequency of 1 MHz, this translates in high RAM occupation (&gt;128 GB) when the algorithm kSpaceFirstOrder3D tries to create the sensor.record matrix. To overcome this problem, my strategy is to define the sensor.record matrix as small as possible exactly where I need to observe the simulation results. To compute the whole domain solution, I was thinking about create multiple cuboids iteratively and save the solution within them once the computation is completed. My question is all about it, is there a way to run a single kSpaceFirstOrder3D and then store the solution in every matrix one at time?<br /> Or should I run kSpaceFirstOrder3D iteratively (for cycle) with different sensor.record matrix every time, saving the results and deleting the variables from the workspace? The second option will of course take a longer time because I will need to run the simulations N times to cover all my 3D domain with different sensor.record matrices.<br /> Thank you in advance, you are doing such a great work with this toolbox!<br /> Federico </p>