-================================================================================
-ADAO, a module for Data Assimilation and Optimization
-================================================================================
+=====================================================
+ADAO: A module for Data Assimilation and Optimization
+=====================================================
About
-=====
+-----
-The ADAO module provides data assimilation and optimization features in SALOME
-context or alone. Briefly stated, Data Assimilation is a methodological
-framework to compute the optimal estimate of the inaccessible true value of a
-system state over time. It uses information coming from experimental
-measurements or observations, and from numerical *a priori* models, including
-information about their errors. Parts of the framework are also known under the
-names of *parameter estimation*, *inverse problems*, *Bayesian estimation*,
-*optimal interpolation*, etc.
+**The ADAO module provides data assimilation and optimization**
+features in Python or SALOME context (see
+http://www.salome-platform.org/). Briefly stated, Data Assimilation is
+a methodological framework to compute the optimal estimate of the
+inaccessible true value of a system state, eventually over time. It
+uses information coming from experimental measurements or observations,
+and from numerical *a priori* models, including information about their
+errors. Parts of the framework are also known under the names of
+*parameter estimation*, *inverse problems*, *Bayesian estimation*,
+*optimal interpolation*, etc. More details can be found in the full
+ADAO documentation (see http://www.salome-platform.org/).
+
+Only the use of ADAO text programming interface (API/TUI) is introduced
+here. This interface gives ability to create a calculation object in a
+similar way than the case building obtained through the graphical
+interface (GUI). When one wants to elaborate "by hand" the TUI
+calculation case, it is recommended to extensively use all the ADAO
+module documentation, and to go back if necessary to the graphical
+interface (GUI), to get all the elements allowing to correctly set the
+commands.
+
+A simple setup example of an ADAO TUI calculation case
+------------------------------------------------------
+
+To introduce the TUI interface, lets begin by a simple but complete
+example of ADAO calculation case. All the data are explicitly defined
+inside the script in order to make the reading easier. The whole set of
+commands is the following one::
+
+ from numpy import array, matrix
+ from adao import adaoBuilder
+ case = adaoBuilder.New()
+ case.set( 'AlgorithmParameters', Algorithm='3DVAR' )
+ case.set( 'Background', Vector=[0, 1, 2] )
+ case.set( 'BackgroundError', ScalarSparseMatrix=1.0 )
+ case.set( 'Observation', Vector=array([0.5, 1.5, 2.5]) )
+ case.set( 'ObservationError', DiagonalSparseMatrix='1 1 1' )
+ case.set( 'ObservationOperator', Matrix='1 0 0;0 2 0;0 0 3' )
+ case.set( 'Observer', Variable="Analysis", Template="ValuePrinter" )
+ case.execute()
+
+The result of running these commands in SALOME (either as a SALOME
+"*shell*" command, in the Python command window of the interface, or by
+the script execution entry of the menu) is the following::
+
+ Analysis [ 0.25000264 0.79999797 0.94999939]
+
+More advanced examples of ADAO TUI calculation case
+---------------------------------------------------
+
+Real cases involves observations loaded from files, operators
+explicitly defined as generic functions including physical simulators,
+time dependant information in order to deal with forecast analysis in
+addition to calibration or re-analysis. More details can be found in
+the full ADAO documentation (see http://www.salome-platform.org/).
License and requirements
-========================
+------------------------
-The license for this module is the GNU Lesser General Public License (Lesser
-GPL), as stated here and in the source files::
+The license for this module is the GNU Lesser General Public License
+(Lesser GPL), as stated here and in the source files::
<ADAO, a module for Data Assimilation and Optimization>
This library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
- MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with this library; if not, write to the Free Software
- Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+ Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
- See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
+ See http://www.salome-platform.org/
-In addition, we require that all publication or presentation describing work
-using this software, or all commercial or not products using it, quote at least
-one of the references given below:
+In addition, we require that all publication or presentation describing
+work using this software, or all commercial or not products using it,
+quote at least one of the references given below:
* *ADAO, a module for Data Assimilation and Optimization*,
http://www.salome-platform.org/
salomeStudyItem = adaoGuiHelper.getSelectedItem()
if salomeStudyItem is not None:
for case_editor, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
self.eficas_manager.selectCase(adao_case.eficas_editor)
break
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
if not adao_case.isOk():
adaoLogger.debug("Cas invalide, donc il est sauvegarde, mais il ne peut pas etre exporte vers YACS ensuite")
self.eficas_manager.adaoFileSave(adao_case)
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
if not adao_case.isOk():
adaoLogger.debug("Cas invalide, donc il est sauvegarde, mais il ne peut pas etre exporte vers YACS ensuite")
self.eficas_manager.adaoFileSaveAs(adao_case)
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
self.eficas_manager.adaoFileClose(adao_case)
break
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
msg = adao_case.validationReportforJDC()
adaoGuiHelper.gui_information(SalomePyQt.SalomePyQt().getDesktop(), msg)
break
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
msg = adao_case.showTreeAdaoCase()
break
self.harmonizeSelectionFromEficas()
salomeStudyItem = adaoGuiHelper.getSelectedItem()
for case_name, adao_case in self.cases.items():
- if adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
+ if hasattr(salomeStudyItem,"GetID") and adao_case.salome_study_item.GetID() == salomeStudyItem.GetID():
if adao_case.isOk():
msg = adao_case.exportCaseToYACS()
# If msg is not empty -> error found
