Added n_cut parameter. If the number of fluctuations is above n_cut, only n_cut fluctuations are rendered and the amplitudes are suppressed by sqrt( n_flucs / n_cut ).
For the case of apertures, I just take the average number of fluctuations in the apertures and check if its above n_cut. I don't think this should matter much since we are are only cancelling at large n.
For reference, for a 0.25 arcsec aperture, 10^-20 eV renders 60,000 fluctuations and 10^-19.5 renders 600,000 fluctuations. So a reasonable value for n_cut should be something like order 10,000.
Added n_cut parameter. If the number of fluctuations is above n_cut, only n_cut fluctuations are rendered and the amplitudes are suppressed by sqrt( n_flucs / n_cut ).
For the case of apertures, I just take the average number of fluctuations in the apertures and check if its above n_cut. I don't think this should matter much since we are are only cancelling at large n.
For reference, for a 0.25 arcsec aperture, 10^-20 eV renders 60,000 fluctuations and 10^-19.5 renders 600,000 fluctuations. So a reasonable value for n_cut should be something like order 10,000.