rayleigh_transient()
Description: Calculates pressure field generated by a transducer array in homogeneous media using the Rayleigh-Sommerfeld Integral.
Usage
pressure = rayleigh_transient(transducer_array, coordinate_grid, medium,time_struct, ndiv,\
input_function);
pressure = rayleigh_transient(transducer_array, coordinate_grid, medium,time_struct, ndiv, input_function,\
disp_flag);
pressure = rayleigh_transient(transducer_array, coordinate_grid, medium,time_struct, ndiv, input_function,\
disp_flag, nthreads);
Arguments
- transducer_array A transducer_array.
- coordinate_grid A coordinate grid struct like the ones created by set_coordinate_grid.
- medium A medium struct like the ones created by set_medium.
- time_struct A time struct generated by set_time_samples.
- ndiv The number of integral points to use.
- input_function The excitation function. Create using set_excitation_function.
- disp_flag Display flag, 1 = display, 0 = suppress.
- nthreads The number of threads to use in the calculation.
Output Parameters
- pressure A 3-d array representing the pressure at each point in space.
Notes
This function only calculates pressures in lossless media, so any attenuation parameter set by the user will be ignored.
Threading in FOCUS is implemented at the transducer level, meaning that single transducers will not benefit from this feature.
Specifying a number of threads larger than the number of CPU cores available will not result in a significant additional speed increase and may in fact result in slower speeds due to additional inter-thread communication overhead.