Precomputed quantities useful for the diffusion approximation in high-scattering opacity regions. We follow [77] Note that r_d in that manuscript should be (distance diffused)/(time elapsed) i.e. the paper is missing a normalization of r_d by delta_t. The upper bound of the integral in Eq (30) should just be sqrt(3*tau/4) All quantities in this struct are fixed; we could also just hard-code them to save us the on-the-fly calculation, or have a single instance of this struct in the global cache. More...
#include <Scattering.hpp>
Public Attributes | |
| const double | MinimumOpacityForDiffusion = 3.0 |
| const double | MaximumOpacityForCorrection = 11.0 |
| const std::array< double, 21 > | BvsRhoForScattering |
| Definition of the vector B_i, equation (35) More... | |
| std::array< double, 1001 > | ScatteringRofP |
| Storage for the function r(P) implicitly defined by Eq (29) ScatteringRofP[i] = r(0.001*i) Calculation performed in constructor. | |
| std::array< double, 101 > | OpacityDependentCorrectionToRofP |
| Storage for the opacity dependent correction on the right-hand side of Eq (30) We use 101 points between the min and max opacities defined above. Calculation performed in constructor. | |
Precomputed quantities useful for the diffusion approximation in high-scattering opacity regions. We follow [77] Note that r_d in that manuscript should be (distance diffused)/(time elapsed) i.e. the paper is missing a normalization of r_d by delta_t. The upper bound of the integral in Eq (30) should just be sqrt(3*tau/4) All quantities in this struct are fixed; we could also just hard-code them to save us the on-the-fly calculation, or have a single instance of this struct in the global cache.
| const std::array<double, 21> Particles::MonteCarlo::DiffusionMonteCarloParameters::BvsRhoForScattering |
Definition of the vector B_i, equation (35)