1 <?xml version='1.0' encoding='iso-8859-1' ?>
3 Copyright (C) 2010-2011 EDF R&D
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.
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.
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
19 See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
21 Author: Andre Ribes, andre.ribes@edf.fr, EDF R&D
25 <objref name="computeAD" id="python:computeAD:1.0">
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>
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>
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>
56 <inline name="CreateAssimilationStudy">
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
67 # Create Assimilation study
68 from daYacsIntegration.daStudy import *
69 assim_study = daStudy(Name, Algorithm, Debug)
71 # Algorithm parameters
77 assim_study.setAlgorithmParameters(AlgorithmParameters)
80 # print "Data entered are:"
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)
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)
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)
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)
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)
149 logging.debug("CREATE EvolutionError is %s"%EvolutionError)
150 logging.debug("CREATE EvolutionErrorStored is %s"%EvolutionErrorStored)
151 assim_study.setEvolutionErrorStored(EvolutionErrorStored)
152 assim_study.setEvolutionError(EvolutionError)
154 # ObservationOperator
155 ObservationOperatorOk = 0
161 logging.debug("CREATE ObservationOperator is %s"%ObservationOperator)
162 logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
163 assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
164 assim_study.setObservationOperator("Matrix", ObservationOperator)
165 ObservationOperatorOk = 1
167 if ObservationOperatorOk == 0:
169 ObservationOperatorDirect
173 logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
174 assim_study.setObservationOperatorType("Direct", "Function")
175 assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
177 ObservationOperatorTangent
181 logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
182 assim_study.setObservationOperatorType("Tangent", "Function")
183 assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
185 ObservationOperatorAdjoint
189 logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
190 assim_study.setObservationOperatorType("Adjoint", "Function")
191 assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
200 logging.debug("CREATE EvolutionModel is %s"%EvolutionModel)
201 logging.debug("CREATE EvolutionModelType is %s"%EvolutionModelType)
202 assim_study.setEvolutionModelType("Matrix", EvolutionModelType)
203 assim_study.setEvolutionModel("Matrix", EvolutionModel)
206 if EvolutionModelOk == 0:
212 logging.debug("CREATE EvolutionModelDirect is %s"%EvolutionModelDirect)
213 assim_study.setEvolutionModelType("Direct", "Function")
214 assim_study.setEvolutionModel("Direct", EvolutionModelDirect)
216 EvolutionModelTangent
220 logging.debug("CREATE EvolutionModelTangent is %s"%EvolutionModelTangent)
221 assim_study.setEvolutionModelType("Tangent", "Function")
222 assim_study.setEvolutionModel("Tangent", EvolutionModelTangent)
224 EvolutionModelAdjoint
228 logging.debug("CREATE EvolutionModelAdjoint is %s"%EvolutionModelAdjoint)
229 assim_study.setEvolutionModelType("Adjoint", "Function")
230 assim_study.setEvolutionModel("Adjoint", EvolutionModelAdjoint)
233 for name, size in zip(InputVariablesNames, InputVariablesSizes):
234 assim_study.setInputVariable(name, size)
235 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
236 assim_study.setOutputVariable(name, size)
239 logging.debug("CREATE Observers is %s"%observers.keys())
240 # Adding observers to the study
241 for observer_name in observers.keys():
244 number = str(observers[observer_name]["number"])
245 if "scheduler" in observers[observer_name].keys():
246 scheduler = observers[observer_name]["scheduler"]
247 if "info" in observers[observer_name].keys():
248 info = observers[observer_name]["info"]
249 assim_study.addObserver(observer_name, scheduler, info, number)
254 <inport name="Name" type="string"/>
255 <inport name="Algorithm" type="string"/>
256 <inport name="Debug" type="bool"/>
257 <inport name="InputVariablesNames" type="stringvec"/>
258 <inport name="InputVariablesSizes" type="intvec"/>
259 <inport name="OutputVariablesNames" type="stringvec"/>
260 <inport name="OutputVariablesSizes" type="intvec"/>
261 <inport name="has_observers" type="bool"/>
262 <inport name="observers" type="pyobj"/>
263 <outport name="Study" type="pyobj"/>
266 <inline name="CreateNumpyMatrixFromString">
267 <script><code><![CDATA[
268 import numpy, logging
269 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
270 matrix = numpy.matrix(matrix_in_string)
272 logging.debug("CREATE Matrix is %s"%matrix)
274 <inport name="matrix_in_string" type="string"/>
275 <outport name="matrix" type="pyobj"/>
276 <outport name="type" type="string"/>
277 <outport name="stored" type="bool"/>
280 <inline name="CreateNumpyMatrixFromScript">
281 <script><code><![CDATA[
283 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
286 # Get file path and filename
289 filepath = os.path.dirname(script)
290 filename = os.path.basename(script)
291 module_name = os.path.splitext(filename)[0]
292 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
293 sys.path.insert(0,filepath)
296 __import__(module_name)
297 user_script_module = sys.modules[module_name]
299 # Get Data from script
301 <inport name="script" type="string"/>
302 <outport name="type" type="string"/>
303 <outport name="stored" type="bool"/>
306 <inline name="CreateNumpyVectorFromString">
307 <script><code><![CDATA[
308 import numpy, logging
309 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
310 vector = numpy.matrix(vector_in_string)
312 logging.debug("Vector is %s"%vector)
314 <inport name="vector_in_string" type="string"/>
315 <outport name="vector" type="pyobj"/>
316 <outport name="type" type="string"/>
317 <outport name="stored" type="bool"/>
320 <inline name="CreateNumpyVectorFromScript">
321 <script><code><![CDATA[
323 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
326 # Get file path and filename
329 filepath = os.path.dirname(script)
330 filename = os.path.basename(script)
331 module_name = os.path.splitext(filename)[0]
332 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
333 sys.path.insert(0,filepath)
336 __import__(module_name)
337 user_script_module = sys.modules[module_name]
339 # Get Data from script
341 <inport name="script" type="string"/>
342 <outport name="type" type="string"/>
343 <outport name="stored" type="bool"/>
346 <inline name="CreateNumpyVectorSerieFromString">
347 <script><code><![CDATA[
348 import numpy, logging
349 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromString")
350 vector = numpy.matrix(vector_in_string)
352 logging.debug("VectorSerie is %s"%vector)
354 <inport name="vector_in_string" type="string"/>
355 <outport name="vector" type="pyobj"/>
356 <outport name="type" type="string"/>
357 <outport name="stored" type="bool"/>
360 <inline name="CreateNumpyVectorSerieFromScript">
361 <script><code><![CDATA[
363 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromScript")
366 # Get file path and filename
369 filepath = os.path.dirname(script)
370 filename = os.path.basename(script)
371 module_name = os.path.splitext(filename)[0]
372 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
373 sys.path.insert(0,filepath)
376 __import__(module_name)
377 user_script_module = sys.modules[module_name]
379 # Get Data from script
381 <inport name="script" type="string"/>
382 <outport name="type" type="string"/>
383 <outport name="stored" type="bool"/>
386 <inline name="SimpleExecuteDirectAlgorithm">
387 <script><code><![CDATA[
389 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
390 from daYacsIntegration.daStudy import *
391 ADD = Study.getAssimilationStudy()
394 <inport name="Study" type="pyobj"/>
395 <outport name="Study" type="pyobj"/>
398 <inline name="SimpleUserAnalysis">
399 <script><code><![CDATA[
400 #-*-coding:iso-8859-1-*-
402 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
403 from daYacsIntegration.daStudy import *
404 ADD = Study.getAssimilationStudy()
408 <inport name="Study" type="pyobj"/>
411 <inline name="FakeOptimizerLoopNode">
412 <script><code><![CDATA[
414 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
417 <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
418 <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
421 <inline name="CreateDictFromScript">
422 <script><code><![CDATA[
424 logging.debug("CREATE Entering in CreateDictFromScript")
426 # Get file path and filename
429 filepath = os.path.dirname(script)
430 filename = os.path.basename(script)
431 module_name = os.path.splitext(filename)[0]
432 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
433 sys.path.insert(0,filepath)
436 __import__(module_name)
437 user_script_module = sys.modules[module_name]
439 # Get Data from script
441 <inport name="script" type="string"/>
444 <inline name="UserDataInitFromScript">
445 <script><code><![CDATA[
447 logging.debug("CREATE Entering in UserDataInitFromScript")
449 # Get file path and filename
452 filepath = os.path.dirname(script)
453 filename = os.path.basename(script)
454 module_name = os.path.splitext(filename)[0]
455 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath,0,1)>0):
456 sys.path.insert(0,filepath)
459 __import__(module_name)
460 user_script_module = sys.modules[module_name]
462 # Get Data from script
464 <inport name="script" type="string"/>
465 <outport name="init_data" type="pyobj"/>
468 <inline name="ReadForSwitchNode">
469 <script><code><![CDATA[
471 logging.debug("CREATE Entering in ReadForSwitchNode")
472 logging.debug(" with input data : "+str(data["specificParameters"]))
474 for param in data["specificParameters"]:
475 if param["name"] == "switch_value":
476 switch_value = int(param["value"])
477 logging.debug(" switching to value : "+str(switch_value))
479 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
480 <outport name="data" type="SALOME_TYPES/ParametricInput"/>
481 <outport name="switch_value" type="int"/>
484 <inline name="ExtractDataNode">
485 <script><code><![CDATA[
487 logging.debug("TERMINATE Entering in ExtractDataNode")
489 from daCore.AssimilationStudy import AssimilationStudy
492 for param in data["specificParameters"]:
493 if param["name"] == "var":
494 var = pickle.loads(param["value"])
495 if param["name"] == "info":
496 info = param["value"]
498 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
499 <outport name="var" type="pyobj"/>
500 <outport name="info" type="pyobj"/>
503 <inline name="ObservationNodeString">
504 <script><code><![CDATA[
505 #print "Entering in Observation"
508 <inport name="var" type="pyobj"/>
509 <inport name="info" type="pyobj"/>
512 <inline name="ObservationNodeFile">
513 <script><code><![CDATA[
514 #print "Entering in Observation"
518 <inport name="var" type="pyobj"/>
519 <inport name="info" type="pyobj"/>
520 <inport name="script" type="string"/>
523 <inline name="EndObservationNode">
524 <script><code><![CDATA[
525 # Create a fake output object.
526 # An observer is always successful.
528 output["outputValues"] = [[[[]]]]
529 output["specificOutputInfos"] = []
530 output["returnCode"] = 0
531 output["errorMessage"] = ""
533 <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
536 <inline name="SetObserversNode">
537 <script><code><![CDATA[
538 #print "Setting observers"
540 <outport name="has_observers" type="bool"/>
541 <outport name="observers" type="pyobj"/>