2 Copyright (C) 2008-2015 EDF R&D
4 This file is part of SALOME ADAO module.
6 This library is free software; you can redistribute it and/or
7 modify it under the terms of the GNU Lesser General Public
8 License as published by the Free Software Foundation; either
9 version 2.1 of the License, or (at your option) any later version.
11 This library is distributed in the hope that it will be useful,
12 but WITHOUT ANY WARRANTY; without even the implied warranty of
13 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
14 Lesser General Public License for more details.
16 You should have received a copy of the GNU Lesser General Public
17 License along with this library; if not, write to the Free Software
18 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
20 See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
22 Author: Jean-Philippe Argaud, jean-philippe.argaud@edf.fr, EDF R&D
24 .. index:: single: 3DVAR
25 .. _section_ref_algorithm_3DVAR:
27 Algorithme de calcul "*3DVAR*"
28 ------------------------------
33 Cet algorithme réalise une estimation d'état par minimisation variationnelle de
34 la fonctionnelle :math:`J` d'écart classique en assimilation de données
37 .. math:: J(\mathbf{x})=(\mathbf{x}-\mathbf{x}^b)^T.\mathbf{B}^{-1}.(\mathbf{x}-\mathbf{x}^b)+(\mathbf{y}^o-\mathbf{H}.\mathbf{x})^T.\mathbf{R}^{-1}.(\mathbf{y}^o-\mathbf{H}.\mathbf{x})
39 qui est usuellement désignée comme la fonctionnelle "*3D-VAR*" (voir par exemple
42 Commandes requises et optionnelles
43 ++++++++++++++++++++++++++++++++++
45 .. index:: single: Background
46 .. index:: single: BackgroundError
47 .. index:: single: Observation
48 .. index:: single: ObservationError
49 .. index:: single: ObservationOperator
50 .. index:: single: Minimizer
51 .. index:: single: Bounds
52 .. index:: single: MaximumNumberOfSteps
53 .. index:: single: CostDecrementTolerance
54 .. index:: single: ProjectedGradientTolerance
55 .. index:: single: GradientNormTolerance
56 .. index:: single: StoreInternalVariables
57 .. index:: single: StoreSupplementaryCalculations
58 .. index:: single: Quantiles
59 .. index:: single: SetSeed
60 .. index:: single: NumberOfSamplesForQuantiles
61 .. index:: single: SimulationForQuantiles
63 Les commandes requises générales, disponibles dans l'interface en édition, sont
67 *Commande obligatoire*. Elle définit le vecteur d'ébauche ou
68 d'initialisation, noté précédemment :math:`\mathbf{x}^b`. Sa valeur est
69 définie comme un objet de type "*Vector*" ou de type "*VectorSerie*".
72 *Commande obligatoire*. Elle définit la matrice de covariance des erreurs
73 d'ébauche, notée précédemment :math:`\mathbf{B}`. Sa valeur est définie
74 comme un objet de type "*Matrix*", de type "*ScalarSparseMatrix*", ou de
75 type "*DiagonalSparseMatrix*".
78 *Commande obligatoire*. Elle définit le vecteur d'observation utilisé en
79 assimilation de données ou en optimisation, et noté précédemment
80 :math:`\mathbf{y}^o`. Sa valeur est définie comme un objet de type "*Vector*"
81 ou de type "*VectorSerie*".
84 *Commande obligatoire*. Elle définit la matrice de covariance des erreurs
85 d'ébauche, notée précédemment :math:`\mathbf{R}`. Sa valeur est définie
86 comme un objet de type "*Matrix*", de type "*ScalarSparseMatrix*", ou de
87 type "*DiagonalSparseMatrix*".
90 *Commande obligatoire*. Elle indique l'opérateur d'observation, noté
91 précédemment :math:`H`, qui transforme les paramètres d'entrée
92 :math:`\mathbf{x}` en résultats :math:`\mathbf{y}` qui sont à comparer aux
93 observations :math:`\mathbf{y}^o`. Sa valeur est définie comme un objet de
94 type "*Function*" ou de type "*Matrix*". Dans le cas du type "*Function*",
95 différentes formes fonctionnelles peuvent être utilisées, comme décrit dans
96 la section :ref:`section_ref_operator_requirements`. Si un contrôle
97 :math:`U` est inclus dans le modèle d'observation, l'opérateur doit être
98 appliqué à une paire :math:`(X,U)`.
100 Les commandes optionnelles générales, disponibles dans l'interface en édition,
101 sont indiquées dans la :ref:`section_ref_assimilation_keywords`. En particulier,
102 la commande optionnelle "*AlgorithmParameters*" permet d'indiquer les options
103 particulières, décrites ci-après, de l'algorithme. On se reportera à la
104 :ref:`section_ref_options_AlgorithmParameters` pour le bon usage de cette
107 Les options de l'algorithme sont les suivantes:
110 Cette clé permet de changer le minimiseur pour l'optimiseur. Le choix par
111 défaut est "LBFGSB", et les choix possibles sont "LBFGSB" (minimisation non
112 linéaire sous contraintes, voir [Byrd95]_, [Morales11]_ et [Zhu97]_), "TNC"
113 (minimisation non linéaire sous contraintes), "CG" (minimisation non
114 linéaire sans contraintes), "BFGS" (minimisation non linéaire sans
115 contraintes), "NCG" (minimisation de type gradient conjugué de Newton). Il
116 est fortement conseillé de conserver la valeur par défaut.
118 Exemple : ``{"Minimizer":"LBFGSB"}``
121 Cette clé permet de définir des bornes supérieure et inférieure pour chaque
122 variable d'état optimisée. Les bornes doivent être données par une liste de
123 liste de paires de bornes inférieure/supérieure pour chaque variable, avec
124 une valeur ``None`` chaque fois qu'il n'y a pas de borne. Les bornes peuvent
125 toujours être spécifiées, mais seuls les optimiseurs sous contraintes les
128 Exemple : ``{"Bounds":[[2.,5.],[1.e-2,10.],[-30.,None],[None,None]]}``
131 Cette clé indique le nombre maximum d'itérations possibles en optimisation
132 itérative. Le défaut est 15000, qui est très similaire à une absence de
133 limite sur les itérations. Il est ainsi recommandé d'adapter ce paramètre
134 aux besoins pour des problèmes réels. Pour certains optimiseurs, le nombre
135 de pas effectif d'arrêt peut être légèrement différent de la limite à cause
136 d'exigences de contrôle interne de l'algorithme.
138 Exemple : ``{"MaximumNumberOfSteps":100}``
140 CostDecrementTolerance
141 Cette clé indique une valeur limite, conduisant à arrêter le processus
142 itératif d'optimisation lorsque la fonction coût décroît moins que cette
143 tolérance au dernier pas. Le défaut est de 1.e-7, et il est recommandé
144 de l'adapter aux besoins pour des problèmes réels.
146 Exemple : ``{"CostDecrementTolerance":1.e-7}``
148 ProjectedGradientTolerance
149 Cette clé indique une valeur limite, conduisant à arrêter le processus
150 itératif d'optimisation lorsque toutes les composantes du gradient projeté
151 sont en-dessous de cette limite. C'est utilisé uniquement par les
152 optimiseurs sous contraintes. Le défaut est -1, qui désigne le défaut
153 interne de chaque optimiseur (usuellement 1.e-5), et il n'est pas recommandé
156 Exemple : ``{"ProjectedGradientTolerance":-1}``
158 GradientNormTolerance
159 Cette clé indique une valeur limite, conduisant à arrêter le processus
160 itératif d'optimisation lorsque la norme du gradient est en dessous de cette
161 limite. C'est utilisé uniquement par les optimiseurs sans contraintes. Le
162 défaut est 1.e-5 et il n'est pas recommandé de le changer.
164 Exemple : ``{"GradientNormTolerance":1.e-5}``
166 StoreInternalVariables
167 Cette clé booléenne permet de stocker les variables internes par défaut,
168 principalement l'état courant lors d'un processus itératif. Attention, cela
169 peut être un choix numériquement coûteux dans certains cas de calculs. La
170 valeur par défaut est "False".
172 Exemple : ``{"StoreInternalVariables":True}``
174 StoreSupplementaryCalculations
175 Cette liste indique les noms des variables supplémentaires qui peuvent être
176 disponibles à la fin de l'algorithme. Cela implique potentiellement des
177 calculs ou du stockage coûteux. La valeur par défaut est une liste vide,
178 aucune de ces variables n'étant calculée et stockée par défaut. Les noms
179 possibles sont dans la liste suivante : ["APosterioriCovariance", "BMA",
180 "OMA", "OMB", "Innovation", "SigmaObs2", "MahalanobisConsistency",
181 "SimulationQuantiles"].
183 Exemple : ``{"StoreSupplementaryCalculations":["BMA","Innovation"]}``
186 Cette liste indique les valeurs de quantile, entre 0 et 1, à estimer par
187 simulation autour de l'état optimal. L'échantillonnage utilise des tirages
188 aléatoires gaussiens multivariés, dirigés par la matrice de covariance a
189 posteriori. Cette option n'est utile que si le calcul supplémentaire
190 "SimulationQuantiles" a été choisi. La valeur par défaut est une liste vide.
192 Exemple : ``{"Quantiles":[0.1,0.9]}``
195 Cette clé permet de donner un nombre entier pour fixer la graine du
196 générateur aléatoire utilisé pour générer l'ensemble. Un valeur pratique est
197 par exemple 1000. Par défaut, la graine est laissée non initialisée, et elle
198 utilise ainsi l'initialisation par défaut de l'ordinateur.
200 Exemple : ``{"SetSeed":1000}``
202 NumberOfSamplesForQuantiles
203 Cette clé indique le nombre de simulations effectuées pour estimer les
204 quantiles. Cette option n'est utile que si le calcul supplémentaire
205 "SimulationQuantiles" a été choisi. Le défaut est 100, ce qui suffit souvent
206 pour une estimation correcte de quantiles courants à 5%, 10%, 90% ou 95%.
208 Exemple : ``{"NumberOfSamplesForQuantiles":100}``
210 SimulationForQuantiles
211 Cette clé indique le type de simulation, linéaire (avec l'opérateur
212 d'observation tangent appliqué sur des incréments de perturbations autour de
213 l'état optimal) ou non-linéaire (avec l'opérateur d'observation standard
214 appliqué aux états perturbés), que l'on veut faire pour chaque perturbation.
215 Cela change essentiellement le temps de chaque simulation élémentaire,
216 usuellement plus long en non-linéaire qu'en linéaire. Cette option n'est
217 utile que si le calcul supplémentaire "SimulationQuantiles" a été choisi. La
218 valeur par défaut est "Linear", et les choix possibles sont "Linear" et
221 Exemple : ``{"SimulationForQuantiles":"Linear"}``
226 Références vers d'autres sections :
227 - :ref:`section_ref_algorithm_Blue`
228 - :ref:`section_ref_algorithm_ExtendedBlue`
229 - :ref:`section_ref_algorithm_LinearityTest`
231 Références bibliographiques :