k-Wave User Forum » User Favorites: fantasyhaha

k-Wave User Forum » User Favorites: fantasyhaha fantasyhaha Support for the k-Wave MATLAB toolbox en-US Tue, 12 May 2026 22:02:04 +0000 http://bbpress.org/?v=1.0.2 <![CDATA[Search]]> q http://www.k-wave.org/forum/search.php Bradley Treeby on "Problem in PA simulation" http://www.k-wave.org/forum/topic/problem-in-pa-simulation#post-8654 Tue, 08 Nov 2022 10:32:34 +0000 Bradley Treeby 8654@http://www.k-wave.org/forum/ Thanks fantasyhaha. Can you send more details about the differences in the output between the linux and windows simulations? fantasyhaha on "Problem in PA simulation" http://www.k-wave.org/forum/topic/problem-in-pa-simulation#post-8631 Fri, 14 Oct 2022 00:47:10 +0000 fantasyhaha 8631@http://www.k-wave.org/forum/ Thank you sir. I have checked the code and I can show you the test code as follow: clearvars; close all % ========================================================================= % SIMULATION % ========================================================================= % create the computational grid PML_size = 20; % size of the PML in grid points Nx = 446; % number of grid points in the x direction Ny = Nx; % number of grid points in the y direction Nz = Nx; % number of grid points in the z direction dx = 0.2e-3; % grid point spacing in the x direction [m] dy = dx; % grid point spacing in the y direction [m] dz = dx; % grid point spacing in the z direction [m] kgrid = kWaveGrid(Nx, dx, Ny, dy, Nz, dz); % define the properties of the propagation medium medium.sound_speed = 1500; % [m/s] % create initial pressure distribution using makeBall % ball_magnitude = 10; % [Pa] ball_x_pos = Nx/2; % [grid points] ball_y_pos = Ny/2; % [grid points] ball_z_pos = Nz/2; % [grid points] % ball_radius = 30; % [grid points] % p0_binary = ball_magnitude * makeBall(Nx, Ny, Nz, ball_x_pos, ball_y_pos, ball_z_pos, ball_radius); sz=80; p0_init=zeros(Nx,Ny,Nz); p0_init((Nx/2-sz/2+1):(Nx/2+sz/2+1),(Nx/2-sz/2+1):(Nx/2+sz/2+1),(Nx/2-sz/2+1):(Nx/2+sz/2+1))=10; load('brain_up_2.mat') brain=imresize3(V,[Nx Ny Nz]); brain_BW=imbinarize(brain); p0_init=p0_init.*brain_BW; % smooth the initial pressure distribution and restore the magnitude p0 = smooth(p0_init, true); % assign to the source structure source.p0 = p0; % define a Cartesian spherical sensor sensor_radius = 40e-3; % [m] center_pos = [0.5, 0.5, 0.5]*Nx*dx; % [m] num_sensor_points = 1024; sensor_mask = makeCartSphere(sensor_radius, num_sensor_points, center_pos, true); % assign to the sensor structure sensor.mask = sensor_mask; % create the time array cfl=0.3; % kgrid.makeTime(medium.sound_speed); kgrid.makeTime(medium.sound_speed,cfl); time_limit=(Nx+sz)*dx/medium.sound_speed; kgrid.t_array(kgrid.t_array>time_limit)=[]; % set the input arguements input_args = {'PMLSize', PML_size, 'PMLInside', false, 'PlotPML', true, ... 'Smooth', false, 'DataCast', 'single', 'CartInterp', 'nearest'}; % % plot the initial pressure and sensor surface in voxel form % voxelPlot(double(p0_init | cart2grid(kgrid, sensor_mask))); % view([60, 20]); dt_stability_limit = checkStability(kgrid, medium); % run the simulation sensor_data = kspaceFirstOrder3D(kgrid, medium, source, sensor, input_args{:}); %% bencox on "Problem in PA simulation" http://www.k-wave.org/forum/topic/problem-in-pa-simulation#post-8622 Fri, 07 Oct 2022 13:35:13 +0000 bencox 8622@http://www.k-wave.org/forum/ Hi fantasyhaha, That sounds strange. You run the same simulation script on the same version of matlab and it's different if you're running Linux or Windows? Please can you send more details about the differences in the output. Best wishes Ben fantasyhaha on "Problem in PA simulation" http://www.k-wave.org/forum/topic/problem-in-pa-simulation#post-8615 Sun, 02 Oct 2022 07:34:37 +0000 fantasyhaha 8615@http://www.k-wave.org/forum/ This problem seems to happen on the server with a Linux system(Ubuntu 22.04). When I test it on a high-performance PC with Windows 10 with the same Matlab version, it's stable. fantasyhaha on "Problem in PA simulation" http://www.k-wave.org/forum/topic/problem-in-pa-simulation#post-8613 Thu, 29 Sep 2022 02:05:05 +0000 fantasyhaha 8613@http://www.k-wave.org/forum/ Hello developers, We meet a problem in simulation in a relatively large space. Parameters are as follows: PML_size = 10; % size of the PML in grid points Nx = 420; % number of grid points in the x direction Ny = Nx; % number of grid points in the y direction Nz = Nx; % number of grid points in the z direction dx = 0.2e-3; % grid point spacing in the x direction [m] dy = dx; % grid point spacing in the y direction [m] dz = dx; % grid point spacing in the z direction [m] kgrid = kWaveGrid(Nx, dx, Ny, dy, Nz, dz); % define the properties of the propagation medium medium.sound_speed = 1500; % [m/s] sz=128; p0_init=zeros(Nx,Ny,Nz); p0_init((Nx/2-sz/2+1):(Nx/2+sz/2+1),(Nx/2-sz/2+1):(Nx/2+sz/2+1),(Nx/2-sz/2+1):(Nx/2+sz/2+1))=10; % smooth the initial pressure distribution and restore the magnitude p0 = smooth(p0_init, true); % assign to the source structure source.p0 = p0; % define a Cartesian spherical sensor sensor_radius = 40e-3; % [m] center_pos = [0, 0, 0]; % [m] num_sensor_points = 1024; sensor_mask = makeCartSphere(sensor_radius, num_sensor_points, center_pos, true); % assign to the sensor structure sensor.mask = sensor_mask; % create the time array cfl=0.15; % kgrid.makeTime(medium.sound_speed); kgrid.makeTime(medium.sound_speed); % kgrid.t_array(3470:end)=[]; % set the input arguements input_args = {'PMLSize', PML_size, 'PMLInside', false, 'PlotPML', true, ... 'Smooth', false, 'DataCast', 'single', 'CartInterp', 'nearest'}; % run the simulation sensor_data = kspaceFirstOrder3D(kgrid, medium, source, sensor, input_args{:}); We meet some problems during the simulation and the signal then the 'sensor_data' becomes continuous oscillation with very high amplitude.(I hope I can share you a figure to show the results) This happened before the wave reach our transducer or PML. We have already tested different 'cfl' values but can not solve this problem.