AN EFFICIENT ALGORITHM FOR CHOOSING SCATTERING DIRECTIONS IN MONTE-CARLO WORK WITH ARBITRARY PHASE FUNCTIONS

被引:13
作者
BARKSTROM, BR
机构
[1] Atmospheric Sciences Division, NASA Langley Research Center, Hampton
关键词
D O I
10.1016/0022-4073(94)00104-F
中图分类号
O43 [光学];
学科分类号
070207 ; 0803 ;
摘要
This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(-5).
引用
收藏
页码:23 / 38
页数:16
相关论文
共 21 条
[1]  
Abramowitz M., 1965, HDB MATH FUNCTIONS
[2]   SCATTERING OF SOLAR-RADIATION AS A FUNCTION OF CLOUD DIMENSIONS AND ORIENTATION [J].
AIDA, M .
JOURNAL OF QUANTITATIVE SPECTROSCOPY & RADIATIVE TRANSFER, 1977, 17 (03) :303-310
[3]  
[Anonymous], 1986, NUMERICAL RECIPES
[4]  
DANIELSO.RE, 1969, J ATMOS SCI, V26, P1078, DOI 10.1175/1520-0469(1969)026<1078:TTOVRT>2.0.CO
[5]  
2
[6]  
Deirmendjian D., 1970, Electromagnetic scattering on spherical polydispersions
[7]  
Hobson E., 1931, THEORY SPHERICAL ELL
[8]  
Kattawar G. W., 1973, Journal of Physical Oceanography, V3, P353, DOI 10.1175/1520-0485(1973)003<0353:MCCOTP>2.0.CO
[9]  
2
[10]  
Knuth D. E., 1969, ART COMPUTER PROGRAM, V2