resources/ADAOSchemaCatalog.xml

   1 <?xml version='1.0' encoding='iso-8859-1' ?>
   2 <!--
   3   Copyright (C) 2010-2012 EDF R&D
   4
   5   This library is free software; you can redistribute it and/or
   6   modify it under the terms of the GNU Lesser General Public
   7   License as published by the Free Software Foundation; either
   8   version 2.1 of the License.
   9
  10   This library is distributed in the hope that it will be useful,
  11   but WITHOUT ANY WARRANTY; without even the implied warranty of
  12   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13   Lesser General Public License for more details.
  14
  15   You should have received a copy of the GNU Lesser General Public
  16   License along with this library; if not, write to the Free Software
  17   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA
  18
  19   See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  20
  21   Author: Jean-Philippe Argaud, jean-philippe.argaud@edf.fr, EDF R&D
  22 -->
  23 <proc>
  24
  25   <objref name="computeAD" id="python:computeAD:1.0">
  26     <base>pyobj</base>
  27   </objref>
  28
  29   <!-- Types for parametric computations -->
  30   <!-- TODO On devrait pouvoir le lire depuis le KERNEL !!!-->
  31   <type name="long" kind="int"/>
  32   <struct name="SALOME_TYPES/Parameter">
  33     <member type="string" name="name"></member>
  34     <member type="string" name="value"></member>
  35   </struct>
  36   <sequence content="SALOME_TYPES/Parameter" name="SALOME_TYPES/ParameterList"></sequence>
  37   <sequence content="double" name="SALOME_TYPES/Variable"></sequence>
  38   <sequence content="SALOME_TYPES/Variable" name="SALOME_TYPES/VariableSequence"></sequence>
  39   <sequence content="SALOME_TYPES/VariableSequence" name="SALOME_TYPES/StateSequence"></sequence>
  40   <sequence content="SALOME_TYPES/StateSequence" name="SALOME_TYPES/TimeSequence"></sequence>
  41   <sequence content="string" name="SALOME_TYPES/VarList"></sequence>
  42   <struct name="SALOME_TYPES/ParametricInput">
  43     <member type="SALOME_TYPES/VarList" name="inputVarList"></member>
  44     <member type="SALOME_TYPES/VarList" name="outputVarList"></member>
  45     <member type="SALOME_TYPES/TimeSequence" name="inputValues"></member>
  46     <member type="SALOME_TYPES/ParameterList" name="specificParameters"></member>
  47   </struct>
  48   <struct name="SALOME_TYPES/ParametricOutput">
  49     <member type="SALOME_TYPES/TimeSequence" name="outputValues"></member>
  50     <member type="SALOME_TYPES/ParameterList" name="specificOutputInfos"></member>
  51     <member type="long" name="returnCode"></member>
  52     <member type="string" name="errorMessage"></member>
  53   </struct>
  54
  55
  56   <inline name="CreateAssimilationStudy">
  57     <script><code>
  58
  59 <![CDATA[
  60 import numpy, logging
  61 logging.debug("CREATE Entering in CreateAssimilationStudy")
  62 print "Entering in the assimilation study"
  63 print "Name is set to........:", Name
  64 print "Algorithm is set to...:", Algorithm
  65 print "Debug is set to.......:", Debug
  66
  67 # Create Assimilation study
  68 from daYacsIntegration.daStudy import *
  69 assim_study = daStudy(Name, Algorithm, Debug)
  70
  71 # Algorithm parameters
  72 try:
  73   AlgorithmParameters
  74 except NameError:
  75   pass
  76 else:
  77   assim_study.setAlgorithmParameters(AlgorithmParameters)
  78
  79 logging.debug("CREATE Data entered are:")
  80
  81 # Background
  82 try:
  83   Background
  84 except NameError:
  85   pass
  86 else:
  87   logging.debug("CREATE Background is set")
  88   logging.debug("CREATE BackgroundType is %s"%BackgroundType)
  89   logging.debug("CREATE BackgroundStored is %s"%BackgroundStored)
  90   assim_study.setBackgroundType(BackgroundType)
  91   assim_study.setBackgroundStored(BackgroundStored)
  92   assim_study.setBackground(Background)
  93
  94 # CheckingPoint
  95 try:
  96   CheckingPoint
  97 except NameError:
  98   pass
  99 else:
 100   logging.debug("CREATE CheckingPoint is set")
 101   logging.debug("CREATE CheckingPointType is %s"%CheckingPointType)
 102   logging.debug("CREATE CheckingPointStored is %s"%CheckingPointStored)
 103   assim_study.setCheckingPointType(CheckingPointType)
 104   assim_study.setCheckingPointStored(CheckingPointStored)
 105   assim_study.setCheckingPoint(CheckingPoint)
 106
 107 # BackgroundError
 108 try:
 109   BackgroundError
 110 except NameError:
 111   pass
 112 else:
 113   logging.debug("CREATE BackgroundError is set")
 114   logging.debug("CREATE BackgroundErrorStored is %s"%BackgroundErrorStored)
 115   assim_study.setBackgroundErrorStored(BackgroundErrorStored)
 116   assim_study.setBackgroundError(BackgroundError)
 117
 118 # Observation
 119 try:
 120   Observation
 121 except NameError:
 122   pass
 123 else:
 124   logging.debug("CREATE Observation is set")
 125   logging.debug("CREATE ObservationType is %s"%ObservationType)
 126   logging.debug("CREATE ObservationStored is %s"%ObservationStored)
 127   assim_study.setObservationType(ObservationType)
 128   assim_study.setObservationStored(ObservationStored)
 129   assim_study.setObservation(Observation)
 130
 131 # ObservationError
 132 try:
 133   ObservationError
 134 except NameError:
 135   pass
 136 else:
 137   logging.debug("CREATE ObservationError is set")
 138   logging.debug("CREATE ObservationErrorStored is %s"%ObservationErrorStored)
 139   assim_study.setObservationErrorStored(ObservationErrorStored)
 140   assim_study.setObservationError(ObservationError)
 141
 142 # EvolutionError
 143 try:
 144   EvolutionError
 145 except NameError:
 146   pass
 147 else:
 148   logging.debug("CREATE EvolutionError is set")
 149   logging.debug("CREATE EvolutionErrorStored is %s"%EvolutionErrorStored)
 150   assim_study.setEvolutionErrorStored(EvolutionErrorStored)
 151   assim_study.setEvolutionError(EvolutionError)
 152
 153 # ObservationOperator
 154 ObservationOperatorOk = 0
 155 try:
 156   ObservationOperator
 157 except NameError:
 158   pass
 159 else:
 160   logging.debug("CREATE ObservationOperator is set")
 161   logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
 162   assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
 163   assim_study.setObservationOperator("Matrix", ObservationOperator)
 164   ObservationOperatorOk = 1
 165
 166 if ObservationOperatorOk == 0:
 167   try:
 168     ObservationOperatorDirect
 169   except NameError:
 170     pass
 171   else:
 172     logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
 173     assim_study.setObservationOperatorType("Direct", "Function")
 174     assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
 175   try:
 176     ObservationOperatorTangent
 177   except NameError:
 178     pass
 179   else:
 180     logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
 181     assim_study.setObservationOperatorType("Tangent", "Function")
 182     assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
 183   try:
 184     ObservationOperatorAdjoint
 185   except NameError:
 186     pass
 187   else:
 188     logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
 189     assim_study.setObservationOperatorType("Adjoint", "Function")
 190     assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
 191
 192 # EvolutionModel
 193 EvolutionModelOk = 0
 194 try:
 195   EvolutionModel
 196 except NameError:
 197   pass
 198 else:
 199   logging.debug("CREATE EvolutionModel is set")
 200   logging.debug("CREATE EvolutionModelType is %s"%EvolutionModelType)
 201   assim_study.setEvolutionModelType("Matrix", EvolutionModelType)
 202   assim_study.setEvolutionModel("Matrix", EvolutionModel)
 203   EvolutionModelOk = 1
 204
 205 if EvolutionModelOk == 0:
 206   try:
 207     EvolutionModelDirect
 208   except NameError:
 209     pass
 210   else:
 211     logging.debug("CREATE EvolutionModelDirect is %s"%EvolutionModelDirect)
 212     assim_study.setEvolutionModelType("Direct", "Function")
 213     assim_study.setEvolutionModel("Direct", EvolutionModelDirect)
 214   try:
 215     EvolutionModelTangent
 216   except NameError:
 217     pass
 218   else:
 219     logging.debug("CREATE EvolutionModelTangent is %s"%EvolutionModelTangent)
 220     assim_study.setEvolutionModelType("Tangent", "Function")
 221     assim_study.setEvolutionModel("Tangent", EvolutionModelTangent)
 222   try:
 223     EvolutionModelAdjoint
 224   except NameError:
 225     pass
 226   else:
 227     logging.debug("CREATE EvolutionModelAdjoint is %s"%EvolutionModelAdjoint)
 228     assim_study.setEvolutionModelType("Adjoint", "Function")
 229     assim_study.setEvolutionModel("Adjoint", EvolutionModelAdjoint)
 230
 231 # Variables
 232 for name, size in zip(InputVariablesNames, InputVariablesSizes):
 233   assim_study.setInputVariable(name, size)
 234 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
 235   assim_study.setOutputVariable(name, size)
 236
 237 if has_observers:
 238   logging.debug("CREATE Observers keys are %s"%observers.keys())
 239   # Adding observers to the study
 240   for observer_name in observers.keys():
 241     scheduler = ""
 242     info = ""
 243     number = str(observers[observer_name]["number"])
 244     if "scheduler" in observers[observer_name].keys():
 245       scheduler = observers[observer_name]["scheduler"]
 246     if "info" in observers[observer_name].keys():
 247       info = observers[observer_name]["info"]
 248     assim_study.addObserver(observer_name, scheduler, info, number)
 249 Study = assim_study
 250 ]]>
 251
 252 </code></script>
 253     <inport name="Name" type="string"/>
 254     <inport name="Algorithm" type="string"/>
 255     <inport name="Debug" type="bool"/>
 256     <inport name="InputVariablesNames" type="stringvec"/>
 257     <inport name="InputVariablesSizes" type="intvec"/>
 258     <inport name="OutputVariablesNames" type="stringvec"/>
 259     <inport name="OutputVariablesSizes" type="intvec"/>
 260     <inport name="has_observers" type="bool"/>
 261     <inport name="observers" type="pyobj"/>
 262     <outport name="Study" type="pyobj"/>
 263   </inline>
 264
 265   <inline name="CreateNumpyMatrixFromString">
 266     <script><code><![CDATA[
 267 import numpy, logging
 268 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
 269 matrix = numpy.matrix(matrix_in_string)
 270 type = "Matrix"
 271 logging.debug("CREATE Matrix is %s"%matrix)
 272 ]]></code></script>
 273     <inport name="matrix_in_string" type="string"/>
 274     <outport name="matrix" type="pyobj"/>
 275     <outport name="type" type="string"/>
 276     <outport name="stored" type="bool"/>
 277   </inline>
 278
 279   <inline name="CreateNumpyMatrixFromScript">
 280     <script><code><![CDATA[
 281 import logging
 282 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
 283 type = "Matrix"
 284
 285 # Get file path and filename
 286 import sys
 287 import os
 288 filepath = os.path.dirname(script)
 289 filename = os.path.basename(script)
 290 module_name = os.path.splitext(filename)[0]
 291 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
 292     sys.path.insert(0,filepath)
 293
 294 # Import script
 295 __import__(module_name)
 296 user_script_module = sys.modules[module_name]
 297
 298 # Get Data from script
 299 ]]></code></script>
 300     <inport name="script" type="string"/>
 301     <outport name="type" type="string"/>
 302     <outport name="stored" type="bool"/>
 303   </inline>
 304
 305   <inline name="CreateNumpyVectorFromString">
 306     <script><code><![CDATA[
 307 import numpy, logging
 308 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
 309 vector = numpy.matrix(vector_in_string)
 310 type = "Vector"
 311 logging.debug("Vector is %s"%vector)
 312 ]]></code></script>
 313     <inport name="vector_in_string" type="string"/>
 314     <outport name="vector" type="pyobj"/>
 315     <outport name="type" type="string"/>
 316     <outport name="stored" type="bool"/>
 317   </inline>
 318
 319   <inline name="CreateNumpyVectorFromScript">
 320     <script><code><![CDATA[
 321 import logging
 322 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
 323 type = "Vector"
 324
 325 # Get file path and filename
 326 import sys
 327 import os
 328 filepath = os.path.dirname(script)
 329 filename = os.path.basename(script)
 330 module_name = os.path.splitext(filename)[0]
 331 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
 332     sys.path.insert(0,filepath)
 333
 334 # Import script
 335 __import__(module_name)
 336 user_script_module = sys.modules[module_name]
 337
 338 # Get Data from script
 339 ]]></code></script>
 340     <inport name="script" type="string"/>
 341     <outport name="type" type="string"/>
 342     <outport name="stored" type="bool"/>
 343   </inline>
 344
 345   <inline name="CreateNumpyVectorSerieFromString">
 346     <script><code><![CDATA[
 347 import numpy, logging
 348 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromString")
 349 vector = numpy.matrix(vector_in_string)
 350 type = "VectorSerie"
 351 logging.debug("VectorSerie is %s"%vector)
 352 ]]></code></script>
 353     <inport name="vector_in_string" type="string"/>
 354     <outport name="vector" type="pyobj"/>
 355     <outport name="type" type="string"/>
 356     <outport name="stored" type="bool"/>
 357   </inline>
 358
 359   <inline name="CreateNumpyVectorSerieFromScript">
 360     <script><code><![CDATA[
 361 import logging
 362 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromScript")
 363 type = "VectorSerie"
 364
 365 # Get file path and filename
 366 import sys
 367 import os
 368 filepath = os.path.dirname(script)
 369 filename = os.path.basename(script)
 370 module_name = os.path.splitext(filename)[0]
 371 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
 372     sys.path.insert(0,filepath)
 373
 374 # Import script
 375 __import__(module_name)
 376 user_script_module = sys.modules[module_name]
 377
 378 # Get Data from script
 379 ]]></code></script>
 380     <inport name="script" type="string"/>
 381     <outport name="type" type="string"/>
 382     <outport name="stored" type="bool"/>
 383   </inline>
 384
 385   <inline name="SimpleExecuteDirectAlgorithm">
 386     <script><code><![CDATA[
 387 import logging
 388 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
 389 from daYacsIntegration.daStudy import *
 390 ADD = Study.getAssimilationStudy()
 391 ADD.analyze()
 392 ]]></code></script>
 393     <inport name="Study" type="pyobj"/>
 394     <outport name="Study" type="pyobj"/>
 395   </inline>
 396
 397   <inline name="SimpleUserAnalysis">
 398     <script><code><![CDATA[
 399 #-*-coding:iso-8859-1-*-
 400 import logging
 401 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
 402 from daYacsIntegration.daStudy import *
 403 ADD = Study.getAssimilationStudy()
 404 # User code is below
 405
 406 ]]></code></script>
 407     <inport name="Study" type="pyobj"/>
 408   </inline>
 409
 410   <inline name="FakeOptimizerLoopNode">
 411     <script><code><![CDATA[
 412 import logging
 413 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
 414 result = None
 415 ]]></code></script>
 416     <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
 417     <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
 418   </inline>
 419
 420   <inline name="CreateDictFromScript">
 421     <script><code><![CDATA[
 422 import logging
 423 logging.debug("CREATE Entering in CreateDictFromScript")
 424
 425 # Get file path and filename
 426 import sys
 427 import os
 428 filepath = os.path.dirname(script)
 429 filename = os.path.basename(script)
 430 module_name = os.path.splitext(filename)[0]
 431 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
 432     sys.path.insert(0,filepath)
 433
 434 # Import script
 435 __import__(module_name)
 436 user_script_module = sys.modules[module_name]
 437
 438 # Get Data from script
 439 ]]></code></script>
 440     <inport name="script" type="string"/>
 441   </inline>
 442
 443   <inline name="UserDataInitFromScript">
 444     <script><code><![CDATA[
 445 import logging
 446 logging.debug("CREATE Entering in UserDataInitFromScript")
 447
 448 # Get file path and filename
 449 import sys
 450 import os
 451 filepath = os.path.dirname(script)
 452 filename = os.path.basename(script)
 453 module_name = os.path.splitext(filename)[0]
 454 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
 455     sys.path.insert(0,filepath)
 456
 457 # Import script
 458 __import__(module_name)
 459 user_script_module = sys.modules[module_name]
 460
 461 # Get Data from script
 462 ]]></code></script>
 463     <inport name="script" type="string"/>
 464     <outport name="init_data" type="pyobj"/>
 465   </inline>
 466
 467   <inline name="ReadForSwitchNode">
 468     <script><code><![CDATA[
 469 import logging
 470 logging.debug("CREATE Entering in ReadForSwitchNode")
 471 logging.debug("       with input data : "+str(data["specificParameters"]))
 472 switch_value = -1
 473 for param in data["specificParameters"]:
 474   if param["name"] == "switch_value":
 475     switch_value = int(param["value"])
 476 logging.debug("       switching to value : "+str(switch_value))
 477 ]]></code></script>
 478     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 479     <outport name="data" type="SALOME_TYPES/ParametricInput"/>
 480     <outport name="switch_value" type="int"/>
 481   </inline>
 482
 483   <inline name="ExtractDataNode">
 484     <script><code><![CDATA[
 485 import logging
 486 logging.debug("TERMINATE Entering in ExtractDataNode")
 487 import cPickle
 488 from daCore.AssimilationStudy import AssimilationStudy
 489 var = None
 490 info = None
 491 for param in data["specificParameters"]:
 492   if param["name"] == "var":
 493     var = cPickle.loads(param["value"])
 494   if param["name"] == "info":
 495     info = param["value"]
 496 ]]></code></script>
 497     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 498     <outport name="var" type="pyobj"/>
 499     <outport name="info" type="pyobj"/>
 500   </inline>
 501
 502   <inline name="ObservationNodeString">
 503     <script><code><![CDATA[
 504 #print "Entering in Observation"
 505
 506 ]]></code></script>
 507     <inport name="var" type="pyobj"/>
 508     <inport name="info" type="pyobj"/>
 509   </inline>
 510
 511   <inline name="ObservationNodeFile">
 512     <script><code><![CDATA[
 513 #print "Entering in Observation"
 514 execfile(script)
 515
 516 ]]></code></script>
 517     <inport name="var"    type="pyobj"/>
 518     <inport name="info"   type="pyobj"/>
 519     <inport name="script" type="string"/>
 520   </inline>
 521
 522   <inline name="EndObservationNode">
 523     <script><code><![CDATA[
 524 # Create a fake output object.
 525 # An observer is always successful.
 526 output = {}
 527 output["outputValues"]        = [[[[]]]]
 528 output["specificOutputInfos"] = []
 529 output["returnCode"]          = 0
 530 output["errorMessage"]        = ""
 531 ]]></code></script>
 532     <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
 533   </inline>
 534
 535   <inline name="SetObserversNode">
 536     <script><code><![CDATA[
 537 #print "Setting observers"
 538 ]]></code></script>
 539     <outport name="has_observers" type="bool"/>
 540     <outport name="observers" type="pyobj"/>
 541   </inline>
 542 </proc>