resources/ADAOSchemaCatalog.xml

   1 <?xml version='1.0' encoding='iso-8859-1' ?>
   2 <!--
   3   Copyright (C) 2010-2011 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: Andre Ribes, andre.ribes@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 # Data
  80 # print "Data entered are:"
  81
  82 # Background
  83 try:
  84   Background
  85 except NameError:
  86   pass
  87 else:
  88   logging.debug("CREATE Background is %s"%Background)
  89   logging.debug("CREATE BackgroundType is %s"%BackgroundType)
  90   logging.debug("CREATE BackgroundStored is %s"%BackgroundStored)
  91   assim_study.setBackgroundType(BackgroundType)
  92   assim_study.setBackgroundStored(BackgroundStored)
  93   assim_study.setBackground(Background)
  94
  95 # CheckingPoint
  96 try:
  97   CheckingPoint
  98 except NameError:
  99   pass
 100 else:
 101   logging.debug("CREATE CheckingPoint is %s"%CheckingPoint)
 102   logging.debug("CREATE CheckingPointType is %s"%CheckingPointType)
 103   logging.debug("CREATE CheckingPointStored is %s"%CheckingPointStored)
 104   assim_study.setCheckingPointType(CheckingPointType)
 105   assim_study.setCheckingPointStored(CheckingPointStored)
 106   assim_study.setCheckingPoint(CheckingPoint)
 107
 108 # BackgroundError
 109 try:
 110   BackgroundError
 111 except NameError:
 112   pass
 113 else:
 114   logging.debug("CREATE BackgroundError is %s"%BackgroundError)
 115   logging.debug("CREATE BackgroundErrorStored is %s"%BackgroundErrorStored)
 116   assim_study.setBackgroundErrorStored(BackgroundErrorStored)
 117   assim_study.setBackgroundError(BackgroundError)
 118
 119 # Observation
 120 try:
 121   Observation
 122 except NameError:
 123   pass
 124 else:
 125   logging.debug("CREATE Observation is %s"%Observation)
 126   logging.debug("CREATE ObservationType is %s"%ObservationType)
 127   logging.debug("CREATE ObservationStored is %s"%ObservationStored)
 128   assim_study.setObservationType(ObservationType)
 129   assim_study.setObservationStored(ObservationStored)
 130   assim_study.setObservation(Observation)
 131
 132 # ObservationError
 133 try:
 134   ObservationError
 135 except NameError:
 136   pass
 137 else:
 138   logging.debug("CREATE ObservationError is %s"%ObservationError)
 139   logging.debug("CREATE ObservationErrorStored is %s"%ObservationErrorStored)
 140   assim_study.setObservationErrorStored(ObservationErrorStored)
 141   assim_study.setObservationError(ObservationError)
 142
 143 # ObservationOperator
 144 ObservationOperatorOk = 0
 145 try:
 146   ObservationOperator
 147 except NameError:
 148   pass
 149 else:
 150   logging.debug("CREATE ObservationOperator is %s"%ObservationOperator)
 151   logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
 152   assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
 153   assim_study.setObservationOperator("Matrix", ObservationOperator)
 154   ObservationOperatorOk = 1
 155
 156 if ObservationOperatorOk == 0:
 157   try:
 158     ObservationOperatorDirect
 159   except NameError:
 160     pass
 161   else:
 162     logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
 163     assim_study.setObservationOperatorType("Direct", "Function")
 164     assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
 165   try:
 166     ObservationOperatorTangent
 167   except NameError:
 168     pass
 169   else:
 170     logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
 171     assim_study.setObservationOperatorType("Tangent", "Function")
 172     assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
 173   try:
 174     ObservationOperatorAdjoint
 175   except NameError:
 176     pass
 177   else:
 178     logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
 179     assim_study.setObservationOperatorType("Adjoint", "Function")
 180     assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
 181
 182 # Variables
 183 for name, size in zip(InputVariablesNames, InputVariablesSizes):
 184   assim_study.setInputVariable(name, size)
 185 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
 186   assim_study.setOutputVariable(name, size)
 187
 188 if has_observers:
 189   logging.debug("CREATE Observers is %s"%observers.keys())
 190   # Adding observers to the study
 191   for observer_name in observers.keys():
 192     scheduler = ""
 193     info = ""
 194     number = str(observers[observer_name]["number"])
 195     if "scheduler" in observers[observer_name].keys():
 196       scheduler = observers[observer_name]["scheduler"]
 197     if "info" in observers[observer_name].keys():
 198       info = observers[observer_name]["info"]
 199     assim_study.addObserver(observer_name, scheduler, info, number)
 200 Study = assim_study
 201 ]]>
 202
 203 </code></script>
 204     <inport name="Name" type="string"/>
 205     <inport name="Algorithm" type="string"/>
 206     <inport name="Debug" type="bool"/>
 207     <inport name="InputVariablesNames" type="stringvec"/>
 208     <inport name="InputVariablesSizes" type="intvec"/>
 209     <inport name="OutputVariablesNames" type="stringvec"/>
 210     <inport name="OutputVariablesSizes" type="intvec"/>
 211     <inport name="has_observers" type="bool"/>
 212     <inport name="observers" type="pyobj"/>
 213     <outport name="Study" type="pyobj"/>
 214   </inline>
 215
 216   <inline name="CreateNumpyMatrixFromString">
 217     <script><code><![CDATA[
 218 import numpy, logging
 219 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
 220 matrix = numpy.matrix(matrix_in_string)
 221 type = "Matrix"
 222 logging.debug("CREATE Matrix is %s"%matrix)
 223 ]]></code></script>
 224     <inport name="matrix_in_string" type="string"/>
 225     <outport name="matrix" type="pyobj"/>
 226     <outport name="type" type="string"/>
 227     <outport name="stored" type="bool"/>
 228   </inline>
 229
 230   <inline name="CreateNumpyMatrixFromScript">
 231     <script><code><![CDATA[
 232 import logging
 233 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
 234 type = "Matrix"
 235
 236 # Get file path and filename
 237 import sys
 238 import os
 239 filepath = os.path.dirname(script)
 240 filename = os.path.basename(script)
 241 module_name = os.path.splitext(filename)[0]
 242 sys.path.insert(0,filepath)
 243
 244 # Import script
 245 __import__(module_name)
 246 user_script_module = sys.modules[module_name]
 247
 248 # Get Data from script
 249 ]]></code></script>
 250     <inport name="script" type="string"/>
 251     <outport name="type" type="string"/>
 252     <outport name="stored" type="bool"/>
 253   </inline>
 254
 255   <inline name="CreateNumpyVectorFromString">
 256     <script><code><![CDATA[
 257 import numpy, logging
 258 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
 259 vector = numpy.matrix(vector_in_string)
 260 type = "Vector"
 261 logging.debug("Vector is %s"%vector)
 262 ]]></code></script>
 263     <inport name="vector_in_string" type="string"/>
 264     <outport name="vector" type="pyobj"/>
 265     <outport name="type" type="string"/>
 266     <outport name="stored" type="bool"/>
 267   </inline>
 268
 269   <inline name="CreateNumpyVectorFromScript">
 270     <script><code><![CDATA[
 271 import logging
 272 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
 273 type = "Vector"
 274
 275 # Get file path and filename
 276 import sys
 277 import os
 278 filepath = os.path.dirname(script)
 279 filename = os.path.basename(script)
 280 module_name = os.path.splitext(filename)[0]
 281 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
 282     sys.path.insert(0,filepath)
 283
 284 # Import script
 285 __import__(module_name)
 286 user_script_module = sys.modules[module_name]
 287
 288 # Get Data from script
 289 ]]></code></script>
 290     <inport name="script" type="string"/>
 291     <outport name="type" type="string"/>
 292     <outport name="stored" type="bool"/>
 293   </inline>
 294
 295   <inline name="SimpleExecuteDirectAlgorithm">
 296     <script><code><![CDATA[
 297 import logging
 298 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
 299 from daYacsIntegration.daStudy import *
 300 ADD = Study.getAssimilationStudy()
 301 ADD.analyze()
 302 ]]></code></script>
 303     <inport name="Study" type="pyobj"/>
 304     <outport name="Study" type="pyobj"/>
 305   </inline>
 306
 307   <inline name="SimpleUserAnalysis">
 308     <script><code><![CDATA[
 309 #-*-coding:iso-8859-1-*-
 310 import logging
 311 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
 312 from daYacsIntegration.daStudy import *
 313 ADD = Study.getAssimilationStudy()
 314 # User code is below
 315
 316 ]]></code></script>
 317     <inport name="Study" type="pyobj"/>
 318   </inline>
 319
 320   <inline name="FakeOptimizerLoopNode">
 321     <script><code><![CDATA[
 322 import logging
 323 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
 324 result = None
 325 ]]></code></script>
 326     <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
 327     <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
 328   </inline>
 329
 330   <inline name="CreateDictFromScript">
 331     <script><code><![CDATA[
 332 import logging
 333 logging.debug("CREATE Entering in CreateDictFromScript")
 334
 335 # Get file path and filename
 336 import sys
 337 import os
 338 filepath = os.path.dirname(script)
 339 filename = os.path.basename(script)
 340 module_name = os.path.splitext(filename)[0]
 341 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
 342     sys.path.insert(0,filepath)
 343
 344 # Import script
 345 __import__(module_name)
 346 user_script_module = sys.modules[module_name]
 347
 348 # Get Data from script
 349 ]]></code></script>
 350     <inport name="script" type="string"/>
 351   </inline>
 352
 353   <inline name="UserDataInitFromScript">
 354     <script><code><![CDATA[
 355 import logging
 356 logging.debug("CREATE Entering in UserDataInitFromScript")
 357
 358 # Get file path and filename
 359 import sys
 360 import os
 361 filepath = os.path.dirname(script)
 362 filename = os.path.basename(script)
 363 module_name = os.path.splitext(filename)[0]
 364 sys.path.insert(0,filepath)
 365
 366 # Import script
 367 __import__(module_name)
 368 user_script_module = sys.modules[module_name]
 369
 370 # Get Data from script
 371 ]]></code></script>
 372     <inport name="script" type="string"/>
 373     <outport name="init_data" type="pyobj"/>
 374   </inline>
 375
 376   <inline name="ReadForSwitchNode">
 377     <script><code><![CDATA[
 378 import logging
 379 logging.debug("CREATE Entering in ReadForSwitchNode")
 380 switch_value = -1
 381 for param in data["specificParameters"]:
 382   if param["name"] == "switch_value":
 383     switch_value = int(param["value"])
 384 ]]></code></script>
 385     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 386     <outport name="data" type="SALOME_TYPES/ParametricInput"/>
 387     <outport name="switch_value" type="int"/>
 388   </inline>
 389
 390   <inline name="ExtractDataNode">
 391     <script><code><![CDATA[
 392 import logging
 393 logging.debug("TERMINATE Entering in ExtractDataNode")
 394 import pickle
 395 from daCore.AssimilationStudy import AssimilationStudy
 396 var = None
 397 info = None
 398 for param in data["specificParameters"]:
 399   if param["name"] == "var":
 400     var = pickle.loads(param["value"])
 401   if param["name"] == "info":
 402     info = param["value"]
 403 ]]></code></script>
 404     <inport name="data" type="SALOME_TYPES/ParametricInput"/>
 405     <outport name="var" type="pyobj"/>
 406     <outport name="info" type="pyobj"/>
 407   </inline>
 408
 409   <inline name="ObservationNodeString">
 410     <script><code><![CDATA[
 411 #print "Entering in Observation"
 412
 413 ]]></code></script>
 414     <inport name="var" type="pyobj"/>
 415     <inport name="info" type="pyobj"/>
 416   </inline>
 417
 418   <inline name="ObservationNodeFile">
 419     <script><code><![CDATA[
 420 #print "Entering in Observation"
 421 execfile(script)
 422
 423 ]]></code></script>
 424     <inport name="var"    type="pyobj"/>
 425     <inport name="info"   type="pyobj"/>
 426     <inport name="script" type="string"/>
 427   </inline>
 428
 429   <inline name="EndObservationNode">
 430     <script><code><![CDATA[
 431 # Create a fake output object.
 432 # An observer is always successful.
 433 output = {}
 434 output["outputValues"]        = [[[[]]]]
 435 output["specificOutputInfos"] = []
 436 output["returnCode"]          = 0
 437 output["errorMessage"]        = ""
 438 ]]></code></script>
 439     <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
 440   </inline>
 441
 442   <inline name="SetObserversNode">
 443     <script><code><![CDATA[
 444 #print "Setting observers"
 445 ]]></code></script>
 446     <outport name="has_observers" type="bool"/>
 447     <outport name="observers" type="pyobj"/>
 448   </inline>
 449 </proc>