commit d09048f5d4035f89afb9adb1eea955436e608645
parent 5e8d57ee99513ac85ec72f8792ff23a449fc28f2
Author: Vincent Forest <vincent.forest@meso-star.com>
Date: Wed, 4 Jan 2023 22:15:14 +0100
Edition mineure des sources noweb
Supprime les espaces en fin de ligne, les lignes vides et quelques
typos.
Diffstat:
1 file changed, 10 insertions(+), 8 deletions(-)
diff --git a/src/sgs_compute_sensitivity_translation.nw b/src/sgs_compute_sensitivity_translation.nw
@@ -775,14 +775,13 @@ capteur pour différente valeur du paramètre géométrique $\PI$, soit pour
différentes hauteurs de la paroi spéculaire. Le calcul du poids associé au
calcul du flux est décrit en annexe \ref{flux}.
-
-
\appendix
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Annexe CL de sensibilité
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Détails de la condition à la limite de sensibilité}
+
Nous récupérons ici la condition à la limite pour une paroi spéculaire donnée
dans \citep{papier_sensib} $s = s(\vec{x},\vec{\omega},\PI)$:
\begin{equation}
@@ -905,7 +904,7 @@ Ce qui permet d'obtenir :
\partial_{1,\vec{u}} I
\end{equation}
et de résoudre la dérivée en estimant $\partial_{1,\vec{u}}I$ le long de la
-surface et en utilisant l'ETR pour résoudre $partial{1,\vec{\omega}}I$:
+surface et en utilisant l'ETR pour résoudre $\partial{1,\vec{\omega}}I$:
\begin{equation}
\partial_{1,\vec{\omega}}I = \mathcal{C} \lbrack I \rbrack
\end{equation}
@@ -951,8 +950,8 @@ paroi de droite :
\item soit la source est transportée directement depuis la paroi de droite
jusqu'au récepteur,
\item soit la source est transportée en direction de la paroi spéculaire et est
-ensuite réfléchie jusqu'à atteindre le récepteur.
-\end{itemize}
+ensuite réfléchie jusqu'à atteindre le récepteur.
+\end{itemize}
Dans notre configuration le paramètre géométrique $\PI$ n'a d'influence que sur
la hauteur de la paroi spéculaire. Ainsi, la contribution radiative de la
source émise par la paroi de droite, transportée directement vers le
@@ -965,7 +964,7 @@ En pratique cela signifie que, en vue du calcul des différences finies, le
calcul par Monte-Carlo du flux au récepteur n'est pas entièrement nécessaire.
Nous pouvons nous contenter d'estimer l'unique partie du flux qui sera
perturbée par une variation de $\PI$, soit les contributions portées par les
-chemins réfléchis par la paroi spéculaire.
+chemins réfléchis par la paroi spéculaire.
En terme d'algorithme nous pouvons donc réutiliser les chemins déjà
échantillonnés pour le problème couplé puisque leur statistique correspond
exactement à celle d'un chemin émis par la paroi de droite et réfléchi par la
@@ -973,7 +972,7 @@ paroi du haut.
Le poids ([[weight_flux_part_spec]]), qui correspond à la partie du flux venant
de la réflection sur la paroi spéculaire, est donc calculé au même moment que
-celui de la sensibilité.
+celui de la sensibilité.
<<calcul du poids>>=
weight_flux_part_spec = Sb * rho * PI * get_Sr_area(scn);
@@ -990,6 +989,9 @@ seront comptés comme une contribution nulle.
\paragraph{TODO} Initialisation du poids du flux au même endroit que
l'initialisation de celui de la sensib.
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+% Annexe fonction de calcul
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
\section{Fonction de calcul}
\paragraph{}
@@ -1101,7 +1103,7 @@ struct smc_estimator_status status = SMC_ESTIMATOR_STATUS_NULL;
\paragraph{} Le lecteur attentif aura remarqué que l'intégrateur utilise la
fonction [[run_realisation]] et non directement la fonction [[réalisation]]
développée dans ce document (voir [[<<exécuter l'intégration Monte Carlo>>]]).
-[[run_réalisation]] est une fonction intermédiaire, qui ne fait qu'appeler la
+[[run_réalisation]] est une fonction intermédiaire qui ne fait qu'appeler la
fonction [[realisation]]. Une fonction parfaitement dispensable sauf à la
bibliothèque \texttt{Star-MonteCarlo} qui nous impose la signature de la
fonction à utiliser, c'est à dire le type des paramètres d'entrées et de
