k-Wave

1. bob2024
   

   Member
   Joined: Aug '24
   Posts: 3

   What does this castZeros() function mean here? There is no such function in MATLAB.
   If you call it you will get error: Unrecognized function or variable 'castZeros'.

   `
   % create storage and scaling variables - all variables are saved as
   % fields of a structure called sensor_data

   % if only p is being stored (i.e., if no user input is given for
   % sensor.record), then sensor_data.p is copied to sensor_data at the
   % end of the simulation

   % time history of the acoustic pressure
   if flags.record_p || flags.record_I || flags.record_I_avg
   sensor_data.p = castZeros([num_sensor_points, num_recorded_time_points]);
   end

   % maximum pressure
   if flags.record_p_max
   sensor_data.p_max = castZeros([num_sensor_points, 1]);
   end

   ......
   `

   Posted 1 week ago #

2. bob2024
   

   Member
   Joined: Aug '24
   Posts: 3

   Sorry, my bad. The function is actually defined in kspaceFirstOrder_inputChecking.m

   
   switch data_cast
   case 'off'
   castZeros = @(sz) zeros(sz);
   precision = 'double';
   case 'single'
   castZeros = @(sz) zeros(sz, 'single');
   precision = 'single';
   case 'gsingle'
   castZeros = @(sz) gzeros(sz, 'single');
   precision = 'single';
   case 'gdouble'
   castZeros = @(sz) gzeros(sz, 'double');
   precision = 'double';
   case 'gpuArray'
   if strcmp(data_cast_prepend, 'single')
   if verLessThan('distcomp', '6.2')

   % original syntax to initialise array of zeros
   castZeros = @(sz) parallel.gpu.GPUArray.zeros(sz, 'single');

   else

   % current syntax to initialise array of zeros
   castZeros = @(sz) gpuArray.zeros(sz, 'single');

   end
   precision = 'single';
   else
   if verLessThan('distcomp', '6.2')

   % original syntax to initialise array of zeros
   castZeros = @(sz) parallel.gpu.GPUArray.zeros(sz, 'double');

   else

   % current syntax to initialise array of zeros
   castZeros = @(sz) gpuArray.zeros(sz, 'double');

   end
   precision = 'double';
   end
   case 'kWaveGPUsingle'
   castZeros = @(sz) zeros(sz, GPUsingle);
   precision = 'single';
   case 'kWaveGPUdouble'
   castZeros = @(sz) zeros(sz, GPUdouble);
   precision = 'double';
   otherwise
   error('Unknown ''DataCast'' option');
   end
   

   Posted 1 week ago #

RSS feed for this topic

Reply

You must log in to post.