2 Copyright (C) 2008-2021 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: AdjointTest
25 .. _section_ref_algorithm_AdjointTest:
27 Algorithme de vérification "*AdjointTest*"
28 ------------------------------------------
30 .. ------------------------------------ ..
31 .. include:: snippets/Header2Algo01.rst
33 Cet algorithme permet de vérifier la qualité de l'opérateur adjoint, en
34 calculant un résidu dont les propriétés théoriques sont connues.
36 On observe le résidu suivant, qui est la différence de deux produits scalaires :
38 .. math:: R(\alpha) = | < TangentF_x(\mathbf{dx}) , \mathbf{y} > - < \mathbf{dx} , AdjointF_x(\mathbf{y}) > |
40 qui doit rester constamment égal à zéro à la précision du calcul. On prend
41 :math:`\mathbf{dx}_0=Normal(0,\mathbf{x})` et
42 :math:`\mathbf{dx}=\alpha*\mathbf{dx}_0`. :math:`F` est le code de calcul.
43 :math:`\mathbf{y}` doit être dans l'image de :math:`F`. S'il n'est pas donné,
44 on prend :math:`\mathbf{y} = F(\mathbf{x})`.
46 .. ------------------------------------ ..
47 .. include:: snippets/Header2Algo02.rst
49 .. include:: snippets/CheckingPoint.rst
51 .. include:: snippets/ObservationOperator.rst
53 .. ------------------------------------ ..
54 .. include:: snippets/Header2Algo03Chck.rst
56 .. include:: snippets/AmplitudeOfInitialDirection.rst
58 .. include:: snippets/EpsilonMinimumExponent.rst
60 .. include:: snippets/InitialDirection.rst
62 .. include:: snippets/SetSeed.rst
64 StoreSupplementaryCalculations
65 .. index:: single: StoreSupplementaryCalculations
67 *Liste de noms*. Cette liste indique les noms des variables supplémentaires
68 qui peuvent être disponibles au cours du déroulement ou à la fin de
69 l'algorithme, si elles sont initialement demandées par l'utilisateur. Cela
70 implique potentiellement des calculs ou du stockage coûteux. La valeur par
71 défaut est une liste vide, aucune de ces variables n'étant calculée et
72 stockée par défaut sauf les variables inconditionnelles. Les noms possibles
73 sont dans la liste suivante : [
76 "SimulatedObservationAtCurrentState",
80 ``{"StoreSupplementaryCalculations":["BMA", "CurrentState"]}``
82 .. ------------------------------------ ..
83 .. include:: snippets/Header2Algo04.rst
85 .. include:: snippets/Residu.rst
87 .. ------------------------------------ ..
88 .. include:: snippets/Header2Algo05.rst
90 .. include:: snippets/CurrentState.rst
92 .. include:: snippets/Residu.rst
94 .. include:: snippets/SimulatedObservationAtCurrentState.rst
96 .. ------------------------------------ ..
97 .. include:: snippets/Header2Algo06.rst
99 - :ref:`section_ref_algorithm_FunctionTest`
100 - :ref:`section_ref_algorithm_LinearityTest`
101 - :ref:`section_ref_algorithm_TangentTest`
102 - :ref:`section_ref_algorithm_GradientTest`
103 - :ref:`section_ref_algorithm_LocalSensitivityTest`