2 Copyright (C) 2008-2024 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: SamplingTest
25 .. _section_ref_algorithm_SamplingTest:
27 Algorithme de vérification "*SamplingTest*"
28 -------------------------------------------
30 .. ------------------------------------ ..
31 .. include:: snippets/Header2Algo01.rst
33 Cet algorithme de test permet d'établir la collection des valeurs d'une
34 fonctionnelle d'erreur :math:`J` de type :math:`L^1`, :math:`L^2` ou
35 :math:`L^{\infty}`, avec ou sans pondérations, à l'aide de l'opérateur
36 d'observation :math:`\mathcal{H}`, pour un échantillon donné a priori d'états
37 :math:`\mathbf{x}`. La fonctionnelle d'erreur par défaut est celle de moindres
38 carrés pondérés augmentés, classiquement utilisée en assimilation de données,
39 utilisant en plus des observations :math:`\mathbf{y}^o`.
41 Ce test est utile pour analyser la sensibilité de la fonctionnelle :math:`J`
42 aux variations de l'état :math:`\mathbf{x}` en particulier.
44 L'échantillonnage des états :math:`\mathbf{x}` peut être fourni explicitement
45 ou sous la forme d'hypercubes, explicites ou échantillonnés selon des
46 distributions courantes, ou à l'aide d'un échantillonnage par hypercube latin
47 (LHS) ou par séquence de Sobol. Les calculs sont optimisés selon les ressources
48 informatiques disponibles et les options demandées par l'utilisateur. On pourra
49 se reporter aux :ref:`section_ref_sampling_requirements` pour une illustration
50 de l'échantillonnage. Attention à la taille de l'hypercube (et donc au nombre
51 de calculs) qu'il est possible d'atteindre, elle peut rapidement devenir
52 importante. Lorsqu'un état n'est pas observable, une valeur *"NaN"* est
55 Il est aussi possible de fournir un ensemble de simulations :math:`\mathbf{y}`
56 déjà établies par ailleurs (donc sans besoin explicite d'un opérateur
57 :math:`\mathcal{H}`), qui sont implicitement associées à un ensemble
58 d'échantillons d'états :math:`\mathbf{x}`. Dans ce cas où l'ensemble de
59 simulations est fourni, il est impératif de fournir aussi l'ensemble des états
60 :math:`\mathbf{x}` par un échantillonnage explicite, dont l'ordre des états
61 correspond à l'ordre des simulations :math:`\mathbf{y}`.
63 Pour accéder aux informations calculées, les résultats de l'échantillonnage ou
64 des simulations doivent être demandés **explicitement** pour éviter les
65 difficultés de stockage (en l'absence de résultats demandés, rien n'est
66 disponible). On utilise pour cela, sur la variable désirée, la sauvegarde
67 finale à l'aide du mot-clé "*UserPostAnalysis*" ou le traitement en cours de
68 calcul à l'aide des "*observer*" adaptés.
70 .. ------------------------------------ ..
71 .. include:: snippets/Header2Algo12.rst
73 .. include:: snippets/FeaturePropDerivativeFree.rst
75 .. include:: snippets/FeaturePropParallelAlgorithm.rst
77 .. ------------------------------------ ..
78 .. include:: snippets/Header2Algo02.rst
80 .. include:: snippets/CheckingPoint.rst
82 .. include:: snippets/BackgroundError.rst
84 .. include:: snippets/Observation.rst
86 .. include:: snippets/ObservationError.rst
88 .. include:: snippets/ObservationOperator.rst
90 .. ------------------------------------ ..
91 .. include:: snippets/Header2Algo03Chck.rst
93 .. include:: snippets/EnsembleOfSnapshots.rst
95 .. include:: snippets/QualityCriterion.rst
97 .. include:: snippets/SampleAsExplicitHyperCube.rst
99 .. include:: snippets/SampleAsIndependantRandomVariables.rst
101 .. include:: snippets/SampleAsMinMaxLatinHyperCube.rst
103 .. include:: snippets/SampleAsMinMaxSobolSequence.rst
105 .. include:: snippets/SampleAsMinMaxStepHyperCube.rst
107 .. include:: snippets/SampleAsnUplet.rst
109 .. include:: snippets/SetDebug.rst
111 .. include:: snippets/SetSeed.rst
113 StoreSupplementaryCalculations
114 .. index:: single: StoreSupplementaryCalculations
116 *Liste de noms*. Cette liste indique les noms des variables supplémentaires,
117 qui peuvent être disponibles au cours du déroulement ou à la fin de
118 l'algorithme, si elles sont initialement demandées par l'utilisateur. Leur
119 disponibilité implique, potentiellement, des calculs ou du stockage coûteux.
120 La valeur par défaut est donc une liste vide, aucune de ces variables n'étant
121 calculée et stockée par défaut (sauf les variables inconditionnelles). Les
122 noms possibles pour les variables supplémentaires sont dans la liste suivante
123 (la description détaillée de chaque variable nommée est donnée dans la suite
124 de cette documentation par algorithme spécifique, dans la sous-partie
125 "*Informations et variables disponibles à la fin de l'algorithme*") : [
130 "EnsembleOfSimulations",
133 "InnovationAtCurrentState",
134 "SimulatedObservationAtCurrentState",
138 ``{"StoreSupplementaryCalculations":["CurrentState", "Residu"]}``
140 .. ------------------------------------ ..
141 .. include:: snippets/Header2Algo04.rst
143 .. include:: snippets/CostFunctionJ.rst
145 .. include:: snippets/CostFunctionJb.rst
147 .. include:: snippets/CostFunctionJo.rst
149 .. ------------------------------------ ..
150 .. include:: snippets/Header2Algo05.rst
152 .. include:: snippets/CostFunctionJ.rst
154 .. include:: snippets/CostFunctionJb.rst
156 .. include:: snippets/CostFunctionJo.rst
158 .. include:: snippets/CurrentState.rst
160 .. include:: snippets/EnsembleOfSimulations.rst
162 .. include:: snippets/EnsembleOfStates.rst
164 .. include:: snippets/Innovation.rst
166 .. include:: snippets/InnovationAtCurrentState.rst
168 .. include:: snippets/SimulatedObservationAtCurrentState.rst
170 .. ------------------------------------ ..
171 .. _section_ref_algorithm_SamplingTest_examples:
173 .. include:: snippets/Header2Algo06.rst
175 - :ref:`section_ref_algorithm_FunctionTest`
176 - :ref:`section_ref_algorithm_LocalSensitivityTest`
178 .. ------------------------------------ ..
179 .. include:: snippets/Header2Algo08.rst
181 - OPENTURNS, voir le *Guide utilisateur du module OPENTURNS* dans le menu principal *Aide* de l'environnement SALOME