Mcm

La classe mcm se concentre sur les aspects communs à l'ensemble des algorithmes de Monte Carlo pour l'étude des phénomènes de transport linéaires. Cette classe peut donc être employée bien au delà du transfert radiatif. On définit le vecteur des observables (correspondant aux grandeurs étudiées) et le vecteur des paramètres (dont l'influence est importante soit pour un besoin de compréhension de mécanismes, soit lors de la recherche d'un optimum où l'étude du couplage avec d'autres phénomènes). L'utilisateur définit lui même entièrement l'algorithme effectuant (à partir d'un générateur aléatoire uniforme sur $[0,1]$) la succession des générations aléatoires nécessaires et calculant les poids correspondant à chaque observable, ainsi qu'à la sensibilité de chaque observable à chaque paramètre. Mais ensuite, sur cette seule base, mcm assure la mise en oeuvre répétée de ces générations aléatoires en parallèle, le contrôle de l'indépendance statistique, le calcul des incertitudes statistiques correspondant à un nombre donné de réalisations, le lancement de séquences supplémentaires si la précision souhaitée n'est pas atteinte, et la gestion des entrées-sorties.

Les sources sont téléchargeables ici