Layered ASA
Layered Medium
The layered medium structure is a typical medium structure with the extra fields of zint, zstart, and zend. To generate a set of layered media, a function, make_tissue_struct, has been provided.
layer_struct=make_tissue_struct(z,num_layers,layer_breaks,layer_media)
To use this function, you need to pass it the z vector you'll be using for the calculation, the number of layers, the location of the layers, and the layers themselves.
Calling the layered ASA function
Calling the function is done in the same way as a normal ASA call.
Pressure=layereasa(source,z,layer_struct,nfft,delta,type,f0,1,1)
If you are using nonlinear media in the layered ASA calculation, the last two arguments control the number of harmonics to use in the calculation, and the number of harmonics to return in the final answer, respectively.
NOTE: This code is still in beta, and results have not been verified.