1 <?xml version='1.0' encoding='utf-8' ?>
3 Copyright (C) 2008-2018 EDF R&D
5 This file is part of SALOME ADAO module
7 This library is free software; you can redistribute it and/or
8 modify it under the terms of the GNU Lesser General Public
9 License as published by the Free Software Foundation; either
10 version 2.1 of the License.
12 This library is distributed in the hope that it will be useful,
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
15 Lesser General Public License for more details.
17 You should have received a copy of the GNU Lesser General Public
18 License along with this library; if not, write to the Free Software
19 Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
21 See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
23 Author: Jean-Philippe Argaud, jean-philippe.argaud@edf.fr, EDF R&D
27 <objref name="computeAD" id="python:computeAD:1.0">
31 <!-- Types for parametric computations -->
32 <!-- TODO On devrait pouvoir le lire depuis le KERNEL !!!-->
33 <type name="long" kind="int"/>
34 <struct name="SALOME_TYPES/Parameter">
35 <member type="string" name="name"></member>
36 <member type="string" name="value"></member>
38 <sequence content="SALOME_TYPES/Parameter" name="SALOME_TYPES/ParameterList"></sequence>
39 <sequence content="double" name="SALOME_TYPES/Variable"></sequence>
40 <sequence content="SALOME_TYPES/Variable" name="SALOME_TYPES/VariableSequence"></sequence>
41 <sequence content="SALOME_TYPES/VariableSequence" name="SALOME_TYPES/StateSequence"></sequence>
42 <sequence content="SALOME_TYPES/StateSequence" name="SALOME_TYPES/TimeSequence"></sequence>
43 <sequence content="string" name="SALOME_TYPES/VarList"></sequence>
44 <struct name="SALOME_TYPES/ParametricInput">
45 <member type="SALOME_TYPES/VarList" name="inputVarList"></member>
46 <member type="SALOME_TYPES/VarList" name="outputVarList"></member>
47 <member type="SALOME_TYPES/TimeSequence" name="inputValues"></member>
48 <member type="SALOME_TYPES/ParameterList" name="specificParameters"></member>
50 <struct name="SALOME_TYPES/ParametricOutput">
51 <member type="SALOME_TYPES/TimeSequence" name="outputValues"></member>
52 <member type="SALOME_TYPES/ParameterList" name="specificOutputInfos"></member>
53 <member type="long" name="returnCode"></member>
54 <member type="string" name="errorMessage"></member>
58 <inline name="CreateAssimilationStudy">
62 import sys, numpy, logging
63 logging.debug("CREATE YI Entering in CreateAssimilationStudy")
64 print("Entering in the assimilation study")
65 print("Name is set to........:", Name)
66 print("Algorithm is set to...:", Algorithm)
69 # Create Assimilation study
70 if sys.version_info.major > 2: import adao
71 from daYacsIntegration.daStudy import *
72 assim_study = daStudy(Name, Algorithm, Debug)
74 logging.debug("CREATE YI Data entered are:")
76 # Algorithm parameters
77 if "AlgorithmParameters" in locals():
78 logging.debug("CREATE YI AlgorithmParameters is %s"%AlgorithmParameters)
79 assim_study.setYIAlgorithmParameters(AlgorithmParameters)
82 if "Background" in locals():
83 logging.debug("CREATE YI Background is set")
84 logging.debug("CREATE YI BackgroundType is %s"%BackgroundType)
85 logging.debug("CREATE YI BackgroundStored is %s"%BackgroundStored)
86 assim_study.setYIBackgroundType(BackgroundType)
87 assim_study.setYIBackgroundStored(BackgroundStored)
88 assim_study.setYIBackground(Background)
91 if "CheckingPoint" in locals():
92 logging.debug("CREATE YI CheckingPoint is set")
93 logging.debug("CREATE YI CheckingPointType is %s"%CheckingPointType)
94 logging.debug("CREATE YI CheckingPointStored is %s"%CheckingPointStored)
95 assim_study.setYICheckingPointType(CheckingPointType)
96 assim_study.setYICheckingPointStored(CheckingPointStored)
97 assim_study.setYICheckingPoint(CheckingPoint)
100 if "ControlInput" in locals():
101 logging.debug("CREATE YI ControlInput is set")
102 logging.debug("CREATE YI ControlInputType is %s"%ControlInputType)
103 logging.debug("CREATE YI ControlInputStored is %s"%ControlInputStored)
104 assim_study.setYIControlInputType(ControlInputType)
105 assim_study.setYIControlInputStored(ControlInputStored)
106 assim_study.setYIControlInput(ControlInput)
109 if "Observation" in locals():
110 logging.debug("CREATE YI Observation is set")
111 logging.debug("CREATE YI ObservationType is %s"%ObservationType)
112 logging.debug("CREATE YI ObservationStored is %s"%ObservationStored)
113 assim_study.setYIObservationType(ObservationType)
114 assim_study.setYIObservationStored(ObservationStored)
115 assim_study.setYIObservation(Observation)
118 if "BackgroundError" in locals():
119 logging.debug("CREATE YI BackgroundError is set")
120 logging.debug("CREATE YI BackgroundErrorType is %s"%BackgroundErrorType)
121 logging.debug("CREATE YI BackgroundErrorStored is %s"%BackgroundErrorStored)
122 assim_study.setYIBackgroundErrorType(BackgroundErrorType)
123 assim_study.setYIBackgroundErrorStored(BackgroundErrorStored)
124 assim_study.setYIBackgroundError(BackgroundError)
127 if "ObservationError" in locals():
128 logging.debug("CREATE YI ObservationError is set")
129 logging.debug("CREATE YI ObservationErrorType is %s"%ObservationErrorType)
130 logging.debug("CREATE YI ObservationErrorStored is %s"%ObservationErrorStored)
131 assim_study.setYIObservationErrorType(ObservationErrorType)
132 assim_study.setYIObservationErrorStored(ObservationErrorStored)
133 assim_study.setYIObservationError(ObservationError)
136 if "EvolutionError" in locals():
137 logging.debug("CREATE YI EvolutionError is set")
138 logging.debug("CREATE YI EvolutionErrorType is %s"%EvolutionErrorType)
139 logging.debug("CREATE YI EvolutionErrorStored is %s"%EvolutionErrorStored)
140 assim_study.setYIEvolutionErrorType(EvolutionErrorType)
141 assim_study.setYIEvolutionErrorStored(EvolutionErrorStored)
142 assim_study.setYIEvolutionError(EvolutionError)
144 # ObservationOperator
145 ObservationOperatorOk = 0
146 if "ObservationOperator" in locals():
147 logging.debug("CREATE YI ObservationOperator is set")
148 logging.debug("CREATE YI ObservationOperatorType is %s"%ObservationOperatorType)
149 assim_study.setYIObservationOperatorType("Matrix", ObservationOperatorType)
150 assim_study.setYIObservationOperator("Matrix", ObservationOperator)
151 ObservationOperatorOk = 1
153 if ObservationOperatorOk == 0:
154 if "ObservationOperatorDirect" in locals():
155 logging.debug("CREATE YI ObservationOperatorDirect is %s"%ObservationOperatorDirect)
156 assim_study.setYIObservationOperatorType("Direct", "Function")
157 assim_study.setYIObservationOperator("Direct", ObservationOperatorDirect)
158 if "ObservationOperatorTangent" in locals():
159 logging.debug("CREATE YI ObservationOperatorTangent is %s"%ObservationOperatorTangent)
160 assim_study.setYIObservationOperatorType("Tangent", "Function")
161 assim_study.setYIObservationOperator("Tangent", ObservationOperatorTangent)
162 if "ObservationOperatorAdjoint" in locals():
163 logging.debug("CREATE YI ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
164 assim_study.setYIObservationOperatorType("Adjoint", "Function")
165 assim_study.setYIObservationOperator("Adjoint", ObservationOperatorAdjoint)
169 if "EvolutionModel" in locals():
170 logging.debug("CREATE YI EvolutionModel is set")
171 logging.debug("CREATE YI EvolutionModelType is %s"%EvolutionModelType)
172 assim_study.setYIEvolutionModelType("Matrix", EvolutionModelType)
173 assim_study.setYIEvolutionModel("Matrix", EvolutionModel)
176 if EvolutionModelOk == 0:
177 if "EvolutionModelDirect" in locals():
178 logging.debug("CREATE YI EvolutionModelDirect is %s"%EvolutionModelDirect)
179 assim_study.setYIEvolutionModelType("Direct", "Function")
180 assim_study.setYIEvolutionModel("Direct", EvolutionModelDirect)
181 if "EvolutionModelTangent" in locals():
182 logging.debug("CREATE YI EvolutionModelTangent is %s"%EvolutionModelTangent)
183 assim_study.setYIEvolutionModelType("Tangent", "Function")
184 assim_study.setYIEvolutionModel("Tangent", EvolutionModelTangent)
185 if "EvolutionModelAdjoint" in locals():
186 logging.debug("CREATE YI EvolutionModelAdjoint is %s"%EvolutionModelAdjoint)
187 assim_study.setYIEvolutionModelType("Adjoint", "Function")
188 assim_study.setYIEvolutionModel("Adjoint", EvolutionModelAdjoint)
191 for name, size in zip(InputVariablesNames, InputVariablesSizes):
192 assim_study.setYIInputVariable(name, size)
193 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
194 assim_study.setYIOutputVariable(name, size)
197 logging.debug("CREATE YI Observers keys are %s"%observers.keys())
198 # Adding observers to the study
199 for observer_name in observers.keys():
202 number = str(observers[observer_name]["number"])
203 if "scheduler" in observers[observer_name].keys():
204 scheduler = observers[observer_name]["scheduler"]
205 if "info" in observers[observer_name].keys():
206 info = observers[observer_name]["info"]
207 assim_study.addYIObserver(observer_name, scheduler, info, number)
209 except Exception as e:
210 if isinstance(e, SyntaxError): msg = "at %s: %s"%(e.offset, e.text)
212 raise ValueError("during execution, the following error occurs:\n\n%s %s\n\nSee also the potential messages, which can show the origin of the above error, in the launching terminal."%(str(e),msg))
218 <inport name="Name" type="string"/>
219 <inport name="Algorithm" type="string"/>
220 <inport name="Debug" type="bool"/>
221 <inport name="InputVariablesNames" type="stringvec"/>
222 <inport name="InputVariablesSizes" type="intvec"/>
223 <inport name="OutputVariablesNames" type="stringvec"/>
224 <inport name="OutputVariablesSizes" type="intvec"/>
225 <inport name="has_observers" type="bool"/>
226 <inport name="observers" type="pyobj"/>
227 <outport name="Study" type="pyobj"/>
230 <inline name="CreateNumpyMatrixFromString">
231 <script><code><![CDATA[
232 import numpy, logging
233 logging.debug("CREATE YI Entering in CreateNumpyMatrixFromString")
234 if (matrix_in_string.find("array")>-1) or (matrix_in_string.find("matrix")>-1):
235 matrix_in_string=matrix_in_string.replace("array","numpy.array")
236 matrix_in_string=matrix_in_string.replace("matrix","numpy.matrix")
237 exec("matrix="+matrix_in_string)
238 matrix = numpy.matrix(matrix)
240 matrix_in_string=matrix_in_string.replace("'","").replace('"','')
241 matrix = numpy.matrix(matrix_in_string)
243 logging.debug("CREATE YI Matrix is %s"%matrix)
245 <inport name="matrix_in_string" type="string"/>
246 <outport name="matrix" type="pyobj"/>
247 <outport name="type" type="string"/>
248 <outport name="stored" type="bool"/>
251 <inline name="CreateNumpyScalarSparseMatrixFromString">
252 <script><code><![CDATA[
253 import numpy, logging
254 logging.debug("CREATE YI Entering in CreateNumpyScalarSparseMatrixFromString")
255 if (matrix_in_string.find("array")>-1) or (matrix_in_string.find("matrix")>-1):
256 matrix_in_string=matrix_in_string.replace("array","numpy.array")
257 matrix_in_string=matrix_in_string.replace("matrix","numpy.matrix")
258 exec("matrix="+matrix_in_string)
259 matrix = numpy.matrix(matrix)
261 matrix_in_string=matrix_in_string.replace("'","").replace('"','')
262 matrix = numpy.matrix(matrix_in_string)
263 type = "ScalarSparseMatrix"
264 logging.debug("CREATE YI ScalarSparseMatrix is %s"%matrix)
266 <inport name="matrix_in_string" type="string"/>
267 <outport name="matrix" type="pyobj"/>
268 <outport name="type" type="string"/>
269 <outport name="stored" type="bool"/>
272 <inline name="CreateNumpyDiagonalSparseMatrixFromString">
273 <script><code><![CDATA[
274 import numpy, logging
275 logging.debug("CREATE YI Entering in CreateNumpyDiagonalSparseMatrixFromString")
276 if (matrix_in_string.find("array")>-1) or (matrix_in_string.find("matrix")>-1):
277 matrix_in_string=matrix_in_string.replace("array","numpy.array")
278 matrix_in_string=matrix_in_string.replace("matrix","numpy.matrix")
279 exec("matrix="+matrix_in_string)
280 matrix = numpy.matrix(matrix)
282 matrix_in_string=matrix_in_string.replace("'","").replace('"','')
283 matrix = numpy.matrix(matrix_in_string)
284 type = "DiagonalSparseMatrix"
285 logging.debug("CREATE YI DiagonalSparseMatrix is %s"%matrix)
287 <inport name="matrix_in_string" type="string"/>
288 <outport name="matrix" type="pyobj"/>
289 <outport name="type" type="string"/>
290 <outport name="stored" type="bool"/>
293 <inline name="CreateNumpyMatrixFromScript">
294 <script><code><![CDATA[
295 import logging, sys, os
296 logging.debug("CREATE YI Entering in CreateNumpyMatrixFromScript")
299 # Get file path and filename
301 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
302 sys.path.insert(0,studydir)
305 filepath = os.path.dirname(script)
306 filename = os.path.basename(script)
307 module_name = os.path.splitext(filename)[0]
308 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
309 sys.path.insert(0,filepath)
312 <inport name="script" type="string"/>
313 <outport name="type" type="string"/>
314 <outport name="stored" type="bool"/>
317 <inline name="CreateNumpyScalarSparseMatrixFromScript">
318 <script><code><![CDATA[
319 import logging, sys, os
320 logging.debug("CREATE YI Entering in CreateNumpyScalarSparseMatrixFromScript")
321 type = "ScalarSparseMatrix"
323 # Get file path and filename
325 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
326 sys.path.insert(0,studydir)
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):
333 sys.path.insert(0,filepath)
336 <inport name="script" type="string"/>
337 <outport name="type" type="string"/>
338 <outport name="stored" type="bool"/>
341 <inline name="CreateNumpyDiagonalSparseMatrixFromScript">
342 <script><code><![CDATA[
343 import logging, sys, os
344 logging.debug("CREATE YI Entering in CreateNumpyDiagonalSparseMatrixFromScript")
345 type = "DiagonalSparseMatrix"
347 # Get file path and filename
349 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
350 sys.path.insert(0,studydir)
353 filepath = os.path.dirname(script)
354 filename = os.path.basename(script)
355 module_name = os.path.splitext(filename)[0]
356 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
357 sys.path.insert(0,filepath)
360 <inport name="script" type="string"/>
361 <outport name="type" type="string"/>
362 <outport name="stored" type="bool"/>
365 <inline name="CreateNumpyVectorFromString">
366 <script><code><![CDATA[
367 import numpy, logging
368 logging.debug("CREATE YI Entering in CreateNumpyVectorFromString")
369 if (vector_in_string.find("array")>-1) or (vector_in_string.find("matrix")>-1):
370 vector_in_string=vector_in_string.replace("array","numpy.array")
371 vector_in_string=vector_in_string.replace("matrix","numpy.matrix")
372 exec("vector="+vector_in_string)
373 vector = numpy.matrix(vector)
375 vector_in_string=vector_in_string.replace("'","").replace('"','')
376 vector = numpy.matrix(vector_in_string)
378 logging.debug("Vector is %s"%vector)
380 <inport name="vector_in_string" type="string"/>
381 <outport name="vector" type="pyobj"/>
382 <outport name="type" type="string"/>
383 <outport name="stored" type="bool"/>
386 <inline name="CreateNumpyVectorFromScript">
387 <script><code><![CDATA[
388 import logging, sys, os
389 logging.debug("CREATE YI Entering in CreateNumpyVectorFromScript")
392 # Get file path and filename
394 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
395 sys.path.insert(0,studydir)
398 filepath = os.path.dirname(script)
399 filename = os.path.basename(script)
400 module_name = os.path.splitext(filename)[0]
401 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
402 sys.path.insert(0,filepath)
405 <inport name="script" type="string"/>
406 <outport name="type" type="string"/>
407 <outport name="stored" type="bool"/>
410 <inline name="CreateNumpyVectorSerieFromString">
411 <script><code><![CDATA[
412 import numpy, logging
413 logging.debug("CREATE YI Entering in CreateNumpyVectorSerieFromString")
414 vector_in_list = eval(str(vector_in_string),{},{})
415 vector = numpy.matrix(vector_in_list)
417 logging.debug("VectorSerie is %s"%vector)
419 <inport name="vector_in_string" type="string"/>
420 <outport name="vector" type="pyobj"/>
421 <outport name="type" type="string"/>
422 <outport name="stored" type="bool"/>
425 <inline name="CreateNumpyVectorSerieFromScript">
426 <script><code><![CDATA[
427 import logging, sys, os
428 logging.debug("CREATE YI Entering in CreateNumpyVectorSerieFromScript")
431 # Get file path and filename
433 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
434 sys.path.insert(0,studydir)
437 filepath = os.path.dirname(script)
438 filename = os.path.basename(script)
439 module_name = os.path.splitext(filename)[0]
440 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
441 sys.path.insert(0,filepath)
444 <inport name="script" type="string"/>
445 <outport name="type" type="string"/>
446 <outport name="stored" type="bool"/>
449 <inline name="SimpleExecuteDirectAlgorithm">
450 <script><code><![CDATA[
452 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
453 from daYacsIntegration.daStudy import *
454 ADD = Study.getResults()
457 <inport name="Study" type="pyobj"/>
458 <outport name="Study" type="pyobj"/>
461 <inline name="SimpleUserAnalysis">
462 <script><code><![CDATA[
463 #-*- coding: utf-8 -*-
465 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
466 if sys.version_info.major < 3:
467 ADD = Study.getResults()
473 <inport name="Study" type="pyobj"/>
476 <inline name="FakeOptimizerLoopNode">
477 <script><code><![CDATA[
479 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
482 <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
483 <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
486 <inline name="CreateDictFromString">
487 <script><code><![CDATA[
488 import logging, sys, os
489 logging.debug("CREATE YI Entering in CreateDictFromString")
490 dico = eval(dict_in_string)
492 <inport name="dict_in_string" type="string"/>
495 <inline name="CreateDictFromScript">
496 <script><code><![CDATA[
497 import logging, sys, os
498 logging.debug("CREATE YI Entering in CreateDictFromScript")
500 # Get file path and filename
502 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
503 sys.path.insert(0,studydir)
506 filepath = os.path.dirname(script)
507 filename = os.path.basename(script)
508 module_name = os.path.splitext(filename)[0]
509 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
510 sys.path.insert(0,filepath)
513 <inport name="script" type="string"/>
516 <inline name="UserDataInitFromScript">
517 <script><code><![CDATA[
518 import logging, sys, os
519 logging.debug("CREATE YI Entering in UserDataInitFromScript")
521 # Get file path and filename
523 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
524 sys.path.insert(0,studydir)
527 filepath = os.path.dirname(script)
528 filename = os.path.basename(script)
529 module_name = os.path.splitext(filename)[0]
530 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
531 sys.path.insert(0,filepath)
534 <inport name="script" type="string"/>
535 <outport name="init_data" type="pyobj"/>
538 <inline name="ReadForSwitchNode">
539 <script><code><![CDATA[import sys, logging, codecs, pickle
540 if sys.version_info.major > 2: import adao
542 return pickle.loads(codecs.decode(data.encode(), "base64"))
543 logging.debug("CREATE YI Entering in ReadForSwitchNode")
544 # logging.debug(" with input data : ",list(data["specificParameters"]))
546 for param in data["specificParameters"]:
547 # logging.debug(" ReadForSwitchNode specificParameters : ",param)
548 if param["name"] == "switch_value":
549 switch_value = int(loads(param["value"]))
550 logging.debug(" switching to value : "+str(switch_value))
552 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
553 <outport name="data" type="SALOME_TYPES/ParametricInput"/>
554 <outport name="switch_value" type="int"/>
557 <inline name="ExtractDataNode">
558 <script><code><![CDATA[import sys, logging, codecs, pickle
559 if sys.version_info.major > 2: import adao
561 return pickle.loads(codecs.decode(data.encode(), "base64"))
562 logging.debug("TERMINATE Entering in ExtractData/Node")
563 from daCore.AssimilationStudy import AssimilationStudy
566 for param in data["specificParameters"]:
567 if param["name"] == "var":
568 var = loads(param["value"])
569 if param["name"] == "info":
570 info = loads(param["value"])
572 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
573 <outport name="var" type="pyobj"/>
574 <outport name="info" type="pyobj"/>
577 <inline name="ObservationNodeString">
578 <script><code><![CDATA[]]></code></script>
579 <inport name="var" type="pyobj"/>
580 <inport name="info" type="pyobj"/>
583 <inline name="ObservationNodeFile">
584 <script><code><![CDATA[execfile(script)]]></code></script>
585 <inport name="var" type="pyobj"/>
586 <inport name="info" type="pyobj"/>
587 <inport name="script" type="string"/>
590 <inline name="EndObservationNode">
591 <script><code><![CDATA[
592 # Create a fake output object.
593 # An observer is always successful.
595 output["outputValues"] = [[[[]]]]
596 output["specificOutputInfos"] = []
597 output["returnCode"] = 0
598 output["errorMessage"] = ""
600 <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
603 <inline name="SetObserversNode">
604 <script><code><![CDATA[
605 #print("Setting observers")
607 <outport name="has_observers" type="bool"/>
608 <outport name="observers" type="pyobj"/>