1 <?xml version='1.0' encoding='iso-8859-1' ?>
3 Copyright (C) 2010-2013 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: Jean-Philippe Argaud, jean-philippe.argaud@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
66 # Create Assimilation study
67 from daYacsIntegration.daStudy import *
68 assim_study = daStudy(Name, Algorithm, Debug)
70 logging.debug("CREATE Data entered are:")
72 # Algorithm parameters
78 logging.debug("CREATE AlgorithmParameters is %s"%AlgorithmParameters)
79 assim_study.setAlgorithmParameters(AlgorithmParameters)
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)
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)
113 logging.debug("CREATE ControlInput is set")
114 logging.debug("CREATE ControlInputType is %s"%ControlInputType)
115 logging.debug("CREATE ControlInputStored is %s"%ControlInputStored)
116 assim_study.setControlInputType(ControlInputType)
117 assim_study.setControlInputStored(ControlInputStored)
118 assim_study.setControlInput(ControlInput)
126 logging.debug("CREATE Observation is set")
127 logging.debug("CREATE ObservationType is %s"%ObservationType)
128 logging.debug("CREATE ObservationStored is %s"%ObservationStored)
129 assim_study.setObservationType(ObservationType)
130 assim_study.setObservationStored(ObservationStored)
131 assim_study.setObservation(Observation)
139 logging.debug("CREATE BackgroundError is set")
140 logging.debug("CREATE BackgroundErrorStored is %s"%BackgroundErrorStored)
141 assim_study.setBackgroundErrorStored(BackgroundErrorStored)
142 assim_study.setBackgroundError(BackgroundError)
150 logging.debug("CREATE ObservationError is set")
151 logging.debug("CREATE ObservationErrorStored is %s"%ObservationErrorStored)
152 assim_study.setObservationErrorStored(ObservationErrorStored)
153 assim_study.setObservationError(ObservationError)
161 logging.debug("CREATE EvolutionError is set")
162 logging.debug("CREATE EvolutionErrorStored is %s"%EvolutionErrorStored)
163 assim_study.setEvolutionErrorStored(EvolutionErrorStored)
164 assim_study.setEvolutionError(EvolutionError)
166 # ObservationOperator
167 ObservationOperatorOk = 0
173 logging.debug("CREATE ObservationOperator is set")
174 logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
175 assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
176 assim_study.setObservationOperator("Matrix", ObservationOperator)
177 ObservationOperatorOk = 1
179 if ObservationOperatorOk == 0:
181 ObservationOperatorDirect
185 logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
186 assim_study.setObservationOperatorType("Direct", "Function")
187 assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
189 ObservationOperatorTangent
193 logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
194 assim_study.setObservationOperatorType("Tangent", "Function")
195 assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
197 ObservationOperatorAdjoint
201 logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
202 assim_study.setObservationOperatorType("Adjoint", "Function")
203 assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
212 logging.debug("CREATE EvolutionModel is set")
213 logging.debug("CREATE EvolutionModelType is %s"%EvolutionModelType)
214 assim_study.setEvolutionModelType("Matrix", EvolutionModelType)
215 assim_study.setEvolutionModel("Matrix", EvolutionModel)
218 if EvolutionModelOk == 0:
224 logging.debug("CREATE EvolutionModelDirect is %s"%EvolutionModelDirect)
225 assim_study.setEvolutionModelType("Direct", "Function")
226 assim_study.setEvolutionModel("Direct", EvolutionModelDirect)
228 EvolutionModelTangent
232 logging.debug("CREATE EvolutionModelTangent is %s"%EvolutionModelTangent)
233 assim_study.setEvolutionModelType("Tangent", "Function")
234 assim_study.setEvolutionModel("Tangent", EvolutionModelTangent)
236 EvolutionModelAdjoint
240 logging.debug("CREATE EvolutionModelAdjoint is %s"%EvolutionModelAdjoint)
241 assim_study.setEvolutionModelType("Adjoint", "Function")
242 assim_study.setEvolutionModel("Adjoint", EvolutionModelAdjoint)
245 for name, size in zip(InputVariablesNames, InputVariablesSizes):
246 assim_study.setInputVariable(name, size)
247 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
248 assim_study.setOutputVariable(name, size)
251 logging.debug("CREATE Observers keys are %s"%observers.keys())
252 # Adding observers to the study
253 for observer_name in observers.keys():
256 number = str(observers[observer_name]["number"])
257 if "scheduler" in observers[observer_name].keys():
258 scheduler = observers[observer_name]["scheduler"]
259 if "info" in observers[observer_name].keys():
260 info = observers[observer_name]["info"]
261 assim_study.addObserver(observer_name, scheduler, info, number)
266 <inport name="Name" type="string"/>
267 <inport name="Algorithm" type="string"/>
268 <inport name="Debug" type="bool"/>
269 <inport name="InputVariablesNames" type="stringvec"/>
270 <inport name="InputVariablesSizes" type="intvec"/>
271 <inport name="OutputVariablesNames" type="stringvec"/>
272 <inport name="OutputVariablesSizes" type="intvec"/>
273 <inport name="has_observers" type="bool"/>
274 <inport name="observers" type="pyobj"/>
275 <outport name="Study" type="pyobj"/>
278 <inline name="CreateNumpyMatrixFromString">
279 <script><code><![CDATA[
280 import numpy, logging
281 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
282 matrix = numpy.matrix(matrix_in_string)
284 logging.debug("CREATE Matrix is %s"%matrix)
286 <inport name="matrix_in_string" type="string"/>
287 <outport name="matrix" type="pyobj"/>
288 <outport name="type" type="string"/>
289 <outport name="stored" type="bool"/>
292 <inline name="CreateNumpyMatrixFromScript">
293 <script><code><![CDATA[
295 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
298 # Get file path and filename
301 filepath = os.path.dirname(script)
302 filename = os.path.basename(script)
303 module_name = os.path.splitext(filename)[0]
304 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
305 sys.path.insert(0,filepath)
308 __import__(module_name)
309 user_script_module = sys.modules[module_name]
311 # Get Data from script
313 <inport name="script" type="string"/>
314 <outport name="type" type="string"/>
315 <outport name="stored" type="bool"/>
318 <inline name="CreateNumpyVectorFromString">
319 <script><code><![CDATA[
320 import numpy, logging
321 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
322 vector = numpy.matrix(vector_in_string)
324 logging.debug("Vector is %s"%vector)
326 <inport name="vector_in_string" type="string"/>
327 <outport name="vector" type="pyobj"/>
328 <outport name="type" type="string"/>
329 <outport name="stored" type="bool"/>
332 <inline name="CreateNumpyVectorFromScript">
333 <script><code><![CDATA[
335 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
338 # Get file path and filename
341 filepath = os.path.dirname(script)
342 filename = os.path.basename(script)
343 module_name = os.path.splitext(filename)[0]
344 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
345 sys.path.insert(0,filepath)
348 __import__(module_name)
349 user_script_module = sys.modules[module_name]
351 # Get Data from script
353 <inport name="script" type="string"/>
354 <outport name="type" type="string"/>
355 <outport name="stored" type="bool"/>
358 <inline name="CreateNumpyVectorSerieFromString">
359 <script><code><![CDATA[
360 import numpy, logging
361 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromString")
362 vector_in_list = eval(str(vector_in_string),{},{})
363 vector = numpy.matrix(vector_in_list)
365 logging.debug("VectorSerie is %s"%vector)
367 <inport name="vector_in_string" type="string"/>
368 <outport name="vector" type="pyobj"/>
369 <outport name="type" type="string"/>
370 <outport name="stored" type="bool"/>
373 <inline name="CreateNumpyVectorSerieFromScript">
374 <script><code><![CDATA[
376 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromScript")
379 # Get file path and filename
382 filepath = os.path.dirname(script)
383 filename = os.path.basename(script)
384 module_name = os.path.splitext(filename)[0]
385 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
386 sys.path.insert(0,filepath)
389 __import__(module_name)
390 user_script_module = sys.modules[module_name]
392 # Get Data from script
394 <inport name="script" type="string"/>
395 <outport name="type" type="string"/>
396 <outport name="stored" type="bool"/>
399 <inline name="SimpleExecuteDirectAlgorithm">
400 <script><code><![CDATA[
402 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
403 from daYacsIntegration.daStudy import *
404 ADD = Study.getAssimilationStudy()
407 <inport name="Study" type="pyobj"/>
408 <outport name="Study" type="pyobj"/>
411 <inline name="SimpleUserAnalysis">
412 <script><code><![CDATA[
413 #-*-coding:iso-8859-1-*-
415 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
416 from daYacsIntegration.daStudy import *
417 ADD = Study.getAssimilationStudy()
421 <inport name="Study" type="pyobj"/>
424 <inline name="FakeOptimizerLoopNode">
425 <script><code><![CDATA[
427 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
430 <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
431 <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
434 <inline name="CreateDictFromScript">
435 <script><code><![CDATA[
437 logging.debug("CREATE Entering in CreateDictFromScript")
439 # Get file path and filename
442 filepath = os.path.dirname(script)
443 filename = os.path.basename(script)
444 module_name = os.path.splitext(filename)[0]
445 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
446 sys.path.insert(0,filepath)
449 __import__(module_name)
450 user_script_module = sys.modules[module_name]
452 # Get Data from script
454 <inport name="script" type="string"/>
457 <inline name="UserDataInitFromScript">
458 <script><code><![CDATA[
460 logging.debug("CREATE Entering in UserDataInitFromScript")
462 # Get file path and filename
465 filepath = os.path.dirname(script)
466 filename = os.path.basename(script)
467 module_name = os.path.splitext(filename)[0]
468 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
469 sys.path.insert(0,filepath)
472 __import__(module_name)
473 user_script_module = sys.modules[module_name]
475 # Get Data from script
477 <inport name="script" type="string"/>
478 <outport name="init_data" type="pyobj"/>
481 <inline name="ReadForSwitchNode">
482 <script><code><![CDATA[
484 logging.debug("CREATE Entering in ReadForSwitchNode")
485 logging.debug(" with input data : "+str(data["specificParameters"]))
487 for param in data["specificParameters"]:
488 if param["name"] == "switch_value":
489 switch_value = int(param["value"])
490 logging.debug(" switching to value : "+str(switch_value))
492 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
493 <outport name="data" type="SALOME_TYPES/ParametricInput"/>
494 <outport name="switch_value" type="int"/>
497 <inline name="ExtractDataNode">
498 <script><code><![CDATA[
500 logging.debug("TERMINATE Entering in ExtractDataNode")
502 from daCore.AssimilationStudy import AssimilationStudy
505 for param in data["specificParameters"]:
506 if param["name"] == "var":
507 var = cPickle.loads(param["value"])
508 if param["name"] == "info":
509 info = param["value"]
511 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
512 <outport name="var" type="pyobj"/>
513 <outport name="info" type="pyobj"/>
516 <inline name="ObservationNodeString">
517 <script><code><![CDATA[
518 #print "Entering in Observation"
521 <inport name="var" type="pyobj"/>
522 <inport name="info" type="pyobj"/>
525 <inline name="ObservationNodeFile">
526 <script><code><![CDATA[
527 #print "Entering in Observation"
531 <inport name="var" type="pyobj"/>
532 <inport name="info" type="pyobj"/>
533 <inport name="script" type="string"/>
536 <inline name="EndObservationNode">
537 <script><code><![CDATA[
538 # Create a fake output object.
539 # An observer is always successful.
541 output["outputValues"] = [[[[]]]]
542 output["specificOutputInfos"] = []
543 output["returnCode"] = 0
544 output["errorMessage"] = ""
546 <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
549 <inline name="SetObserversNode">
550 <script><code><![CDATA[
551 #print "Setting observers"
553 <outport name="has_observers" type="bool"/>
554 <outport name="observers" type="pyobj"/>