2 Copyright (C) 2008-2023 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). Les calculs sont optimisés selon les ressources informatiques
48 disponibles et les options demandées par l'utilisateur. Attention à la taille
49 de l'hypercube (et donc au nombre de calculs) qu'il est possible d'atteindre,
50 elle peut rapidement devenir importante. Lorsqu'un état n'est pas observable,
51 une valeur *"NaN"* est retournée.
53 Il est aussi possible de fournir un ensemble de simulations :math:`\mathbf{y}`
54 déjà établies par ailleurs (donc sans besoin explicite d'un opérateur
55 :math:`\mathcal{H}`), qui sont implicitement associées à un ensemble
56 d'échantillons d'états :math:`\mathbf{x}`. Dans ce cas où l'ensemble de
57 simulations est fourni, il est impératif de fournir aussi l'ensemble des états
58 :math:`\mathbf{x}` par un échantillonnage explicite, dont l'ordre des états
59 correspond à l'ordre des simulations :math:`\mathbf{y}`.
61 Pour accéder aux informations calculées, les résultats de l'échantillonnage ou
62 des simulations doivent être demandés **explicitement** pour éviter les
63 difficultés de stockage (en l'absence de résultats demandés, rien n'est
64 disponible). On utilise pour cela, sur la variable désirée, la sauvegarde
65 finale à l'aide du mot-clé "*UserPostAnalysis*" ou le traitement en cours de
66 calcul à l'aide des "*observer*" adaptés.
68 .. ------------------------------------ ..
69 .. include:: snippets/Header2Algo02.rst
71 .. include:: snippets/CheckingPoint.rst
73 .. include:: snippets/BackgroundError.rst
75 .. include:: snippets/Observation.rst
77 .. include:: snippets/ObservationError.rst
79 .. include:: snippets/ObservationOperator.rst
81 .. ------------------------------------ ..
82 .. include:: snippets/Header2Algo03Chck.rst
84 .. include:: snippets/EnsembleOfSnapshots.rst
86 .. include:: snippets/QualityCriterion.rst
88 .. include:: snippets/SampleAsExplicitHyperCube.rst
90 .. include:: snippets/SampleAsIndependantRandomVariables.rst
92 .. include:: snippets/SampleAsMinMaxLatinHyperCube.rst
94 .. include:: snippets/SampleAsMinMaxStepHyperCube.rst
96 .. include:: snippets/SampleAsnUplet.rst
98 .. include:: snippets/SetDebug.rst
100 .. include:: snippets/SetSeed.rst
102 StoreSupplementaryCalculations
103 .. index:: single: StoreSupplementaryCalculations
105 *Liste de noms*. Cette liste indique les noms des variables supplémentaires,
106 qui peuvent être disponibles au cours du déroulement ou à la fin de
107 l'algorithme, si elles sont initialement demandées par l'utilisateur. Leur
108 disponibilité implique, potentiellement, des calculs ou du stockage coûteux.
109 La valeur par défaut est donc une liste vide, aucune de ces variables n'étant
110 calculée et stockée par défaut (sauf les variables inconditionnelles). Les
111 noms possibles pour les variables supplémentaires sont dans la liste suivante
112 (la description détaillée de chaque variable nommée est donnée dans la suite
113 de cette documentation par algorithme spécifique, dans la sous-partie
114 "*Informations et variables disponibles à la fin de l'algorithme*") : [
119 "EnsembleOfSimulations",
122 "InnovationAtCurrentState",
123 "SimulatedObservationAtCurrentState",
127 ``{"StoreSupplementaryCalculations":["CurrentState", "Residu"]}``
129 .. ------------------------------------ ..
130 .. include:: snippets/Header2Algo04.rst
132 .. include:: snippets/CostFunctionJ.rst
134 .. include:: snippets/CostFunctionJb.rst
136 .. include:: snippets/CostFunctionJo.rst
138 .. ------------------------------------ ..
139 .. include:: snippets/Header2Algo05.rst
141 .. include:: snippets/CostFunctionJ.rst
143 .. include:: snippets/CostFunctionJb.rst
145 .. include:: snippets/CostFunctionJo.rst
147 .. include:: snippets/CurrentState.rst
149 .. include:: snippets/EnsembleOfSimulations.rst
151 .. include:: snippets/EnsembleOfStates.rst
153 .. include:: snippets/Innovation.rst
155 .. include:: snippets/InnovationAtCurrentState.rst
157 .. include:: snippets/SimulatedObservationAtCurrentState.rst
159 .. ------------------------------------ ..
160 .. _section_ref_algorithm_SamplingTest_examples:
162 .. include:: snippets/Header2Algo06.rst
164 - :ref:`section_ref_algorithm_FunctionTest`
165 - :ref:`section_ref_algorithm_LocalSensitivityTest`
167 .. ------------------------------------ ..
168 .. include:: snippets/Header2Algo08.rst
170 - OPENTURNS, voir le *Guide utilisateur du module OPENTURNS* dans le menu principal *Aide* de l'environnement SALOME