1 <?xml version='1.0' encoding='iso-8859-1' ?>
3 Copyright (C) 2010-2014 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 BackgroundErrorType is %s"%BackgroundErrorType)
141 logging.debug("CREATE BackgroundErrorStored is %s"%BackgroundErrorStored)
142 assim_study.setBackgroundErrorType(BackgroundErrorType)
143 assim_study.setBackgroundErrorStored(BackgroundErrorStored)
144 assim_study.setBackgroundError(BackgroundError)
152 logging.debug("CREATE ObservationError is set")
153 logging.debug("CREATE ObservationErrorType is %s"%ObservationErrorType)
154 logging.debug("CREATE ObservationErrorStored is %s"%ObservationErrorStored)
155 assim_study.setObservationErrorType(ObservationErrorType)
156 assim_study.setObservationErrorStored(ObservationErrorStored)
157 assim_study.setObservationError(ObservationError)
165 logging.debug("CREATE EvolutionError is set")
166 logging.debug("CREATE EvolutionErrorType is %s"%EvolutionErrorType)
167 logging.debug("CREATE EvolutionErrorStored is %s"%EvolutionErrorStored)
168 assim_study.setEvolutionErrorType(EvolutionErrorType)
169 assim_study.setEvolutionErrorStored(EvolutionErrorStored)
170 assim_study.setEvolutionError(EvolutionError)
172 # ObservationOperator
173 ObservationOperatorOk = 0
179 logging.debug("CREATE ObservationOperator is set")
180 logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
181 assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
182 assim_study.setObservationOperator("Matrix", ObservationOperator)
183 ObservationOperatorOk = 1
185 if ObservationOperatorOk == 0:
187 ObservationOperatorDirect
191 logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
192 assim_study.setObservationOperatorType("Direct", "Function")
193 assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
195 ObservationOperatorTangent
199 logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
200 assim_study.setObservationOperatorType("Tangent", "Function")
201 assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
203 ObservationOperatorAdjoint
207 logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
208 assim_study.setObservationOperatorType("Adjoint", "Function")
209 assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
218 logging.debug("CREATE EvolutionModel is set")
219 logging.debug("CREATE EvolutionModelType is %s"%EvolutionModelType)
220 assim_study.setEvolutionModelType("Matrix", EvolutionModelType)
221 assim_study.setEvolutionModel("Matrix", EvolutionModel)
224 if EvolutionModelOk == 0:
230 logging.debug("CREATE EvolutionModelDirect is %s"%EvolutionModelDirect)
231 assim_study.setEvolutionModelType("Direct", "Function")
232 assim_study.setEvolutionModel("Direct", EvolutionModelDirect)
234 EvolutionModelTangent
238 logging.debug("CREATE EvolutionModelTangent is %s"%EvolutionModelTangent)
239 assim_study.setEvolutionModelType("Tangent", "Function")
240 assim_study.setEvolutionModel("Tangent", EvolutionModelTangent)
242 EvolutionModelAdjoint
246 logging.debug("CREATE EvolutionModelAdjoint is %s"%EvolutionModelAdjoint)
247 assim_study.setEvolutionModelType("Adjoint", "Function")
248 assim_study.setEvolutionModel("Adjoint", EvolutionModelAdjoint)
251 for name, size in zip(InputVariablesNames, InputVariablesSizes):
252 assim_study.setInputVariable(name, size)
253 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
254 assim_study.setOutputVariable(name, size)
257 logging.debug("CREATE Observers keys are %s"%observers.keys())
258 # Adding observers to the study
259 for observer_name in observers.keys():
262 number = str(observers[observer_name]["number"])
263 if "scheduler" in observers[observer_name].keys():
264 scheduler = observers[observer_name]["scheduler"]
265 if "info" in observers[observer_name].keys():
266 info = observers[observer_name]["info"]
267 assim_study.addObserver(observer_name, scheduler, info, number)
272 <inport name="Name" type="string"/>
273 <inport name="Algorithm" type="string"/>
274 <inport name="Debug" type="bool"/>
275 <inport name="InputVariablesNames" type="stringvec"/>
276 <inport name="InputVariablesSizes" type="intvec"/>
277 <inport name="OutputVariablesNames" type="stringvec"/>
278 <inport name="OutputVariablesSizes" type="intvec"/>
279 <inport name="has_observers" type="bool"/>
280 <inport name="observers" type="pyobj"/>
281 <outport name="Study" type="pyobj"/>
284 <inline name="CreateNumpyMatrixFromString">
285 <script><code><![CDATA[
286 import numpy, logging
287 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
288 matrix = numpy.matrix(matrix_in_string)
290 logging.debug("CREATE Matrix is %s"%matrix)
292 <inport name="matrix_in_string" type="string"/>
293 <outport name="matrix" type="pyobj"/>
294 <outport name="type" type="string"/>
295 <outport name="stored" type="bool"/>
298 <inline name="CreateNumpyScalarSparseMatrixFromString">
299 <script><code><![CDATA[
300 import numpy, logging
301 logging.debug("CREATE Entering in CreateNumpyScalarSparseMatrixFromString")
302 matrix = numpy.matrix(matrix_in_string)
303 type = "ScalarSparseMatrix"
304 logging.debug("CREATE ScalarSparseMatrix is %s"%matrix)
306 <inport name="matrix_in_string" type="string"/>
307 <outport name="matrix" type="pyobj"/>
308 <outport name="type" type="string"/>
309 <outport name="stored" type="bool"/>
312 <inline name="CreateNumpyDiagonalSparseMatrixFromString">
313 <script><code><![CDATA[
314 import numpy, logging
315 logging.debug("CREATE Entering in CreateNumpyDiagonalSparseMatrixFromString")
316 matrix = numpy.matrix(matrix_in_string)
317 type = "DiagonalSparseMatrix"
318 logging.debug("CREATE DiagonalSparseMatrix is %s"%matrix)
320 <inport name="matrix_in_string" type="string"/>
321 <outport name="matrix" type="pyobj"/>
322 <outport name="type" type="string"/>
323 <outport name="stored" type="bool"/>
326 <inline name="CreateNumpyMatrixFromScript">
327 <script><code><![CDATA[
328 import logging, sys, os
329 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
332 # Get file path and filename
334 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
335 sys.path.insert(0,studydir)
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):
342 sys.path.insert(0,filepath)
345 <inport name="script" type="string"/>
346 <outport name="type" type="string"/>
347 <outport name="stored" type="bool"/>
350 <inline name="CreateNumpyScalarSparseMatrixFromScript">
351 <script><code><![CDATA[
352 import logging, sys, os
353 logging.debug("CREATE Entering in CreateNumpyScalarSparseMatrixFromScript")
354 type = "ScalarSparseMatrix"
356 # Get file path and filename
358 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
359 sys.path.insert(0,studydir)
362 filepath = os.path.dirname(script)
363 filename = os.path.basename(script)
364 module_name = os.path.splitext(filename)[0]
365 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
366 sys.path.insert(0,filepath)
369 <inport name="script" type="string"/>
370 <outport name="type" type="string"/>
371 <outport name="stored" type="bool"/>
374 <inline name="CreateNumpyDiagonalSparseMatrixFromScript">
375 <script><code><![CDATA[
376 import logging, sys, os
377 logging.debug("CREATE Entering in CreateNumpyDiagonalSparseMatrixFromScript")
378 type = "DiagonalSparseMatrix"
380 # Get file path and filename
382 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
383 sys.path.insert(0,studydir)
386 filepath = os.path.dirname(script)
387 filename = os.path.basename(script)
388 module_name = os.path.splitext(filename)[0]
389 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
390 sys.path.insert(0,filepath)
393 <inport name="script" type="string"/>
394 <outport name="type" type="string"/>
395 <outport name="stored" type="bool"/>
398 <inline name="CreateNumpyVectorFromString">
399 <script><code><![CDATA[
400 import numpy, logging
401 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
402 vector = numpy.matrix(vector_in_string)
404 logging.debug("Vector is %s"%vector)
406 <inport name="vector_in_string" type="string"/>
407 <outport name="vector" type="pyobj"/>
408 <outport name="type" type="string"/>
409 <outport name="stored" type="bool"/>
412 <inline name="CreateNumpyVectorFromScript">
413 <script><code><![CDATA[
414 import logging, sys, os
415 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
418 # Get file path and filename
420 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
421 sys.path.insert(0,studydir)
424 filepath = os.path.dirname(script)
425 filename = os.path.basename(script)
426 module_name = os.path.splitext(filename)[0]
427 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
428 sys.path.insert(0,filepath)
431 <inport name="script" type="string"/>
432 <outport name="type" type="string"/>
433 <outport name="stored" type="bool"/>
436 <inline name="CreateNumpyVectorSerieFromString">
437 <script><code><![CDATA[
438 import numpy, logging
439 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromString")
440 vector_in_list = eval(str(vector_in_string),{},{})
441 vector = numpy.matrix(vector_in_list)
443 logging.debug("VectorSerie is %s"%vector)
445 <inport name="vector_in_string" type="string"/>
446 <outport name="vector" type="pyobj"/>
447 <outport name="type" type="string"/>
448 <outport name="stored" type="bool"/>
451 <inline name="CreateNumpyVectorSerieFromScript">
452 <script><code><![CDATA[
453 import logging, sys, os
454 logging.debug("CREATE Entering in CreateNumpyVectorSerieFromScript")
457 # Get file path and filename
459 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
460 sys.path.insert(0,studydir)
463 filepath = os.path.dirname(script)
464 filename = os.path.basename(script)
465 module_name = os.path.splitext(filename)[0]
466 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
467 sys.path.insert(0,filepath)
470 <inport name="script" type="string"/>
471 <outport name="type" type="string"/>
472 <outport name="stored" type="bool"/>
475 <inline name="SimpleExecuteDirectAlgorithm">
476 <script><code><![CDATA[
478 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
479 from daYacsIntegration.daStudy import *
480 ADD = Study.getAssimilationStudy()
483 <inport name="Study" type="pyobj"/>
484 <outport name="Study" type="pyobj"/>
487 <inline name="SimpleUserAnalysis">
488 <script><code><![CDATA[
489 #-*-coding:iso-8859-1-*-
491 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
492 ADD = Study.getAssimilationStudy()
496 <inport name="Study" type="pyobj"/>
499 <inline name="FakeOptimizerLoopNode">
500 <script><code><![CDATA[
502 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
505 <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
506 <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
509 <inline name="CreateDictFromString">
510 <script><code><![CDATA[
511 import logging, sys, os
512 logging.debug("CREATE Entering in CreateDictFromString")
513 dico = eval(dict_in_string)
515 <inport name="dict_in_string" type="string"/>
518 <inline name="CreateDictFromScript">
519 <script><code><![CDATA[
520 import logging, sys, os
521 logging.debug("CREATE Entering in CreateDictFromScript")
523 # Get file path and filename
525 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
526 sys.path.insert(0,studydir)
529 filepath = os.path.dirname(script)
530 filename = os.path.basename(script)
531 module_name = os.path.splitext(filename)[0]
532 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
533 sys.path.insert(0,filepath)
536 <inport name="script" type="string"/>
539 <inline name="UserDataInitFromScript">
540 <script><code><![CDATA[
541 import logging, sys, os
542 logging.debug("CREATE Entering in UserDataInitFromScript")
544 # Get file path and filename
546 if sys.path.count(studydir)==0 or (sys.path.count(studydir)>0 and sys.path.index(studydir)>0):
547 sys.path.insert(0,studydir)
550 filepath = os.path.dirname(script)
551 filename = os.path.basename(script)
552 module_name = os.path.splitext(filename)[0]
553 if sys.path.count(filepath)==0 or (sys.path.count(filepath)>0 and sys.path.index(filepath)>0):
554 sys.path.insert(0,filepath)
557 <inport name="script" type="string"/>
558 <outport name="init_data" type="pyobj"/>
561 <inline name="ReadForSwitchNode">
562 <script><code><![CDATA[
564 logging.debug("CREATE Entering in ReadForSwitchNode")
565 logging.debug(" with input data : "+str(data["specificParameters"]))
567 for param in data["specificParameters"]:
568 if param["name"] == "switch_value":
569 switch_value = int(param["value"])
570 logging.debug(" switching to value : "+str(switch_value))
572 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
573 <outport name="data" type="SALOME_TYPES/ParametricInput"/>
574 <outport name="switch_value" type="int"/>
577 <inline name="ExtractDataNode">
578 <script><code><![CDATA[
580 logging.debug("TERMINATE Entering in ExtractDataNode")
582 from daCore.AssimilationStudy import AssimilationStudy
585 for param in data["specificParameters"]:
586 if param["name"] == "var":
587 var = cPickle.loads(param["value"])
588 if param["name"] == "info":
589 info = param["value"]
591 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
592 <outport name="var" type="pyobj"/>
593 <outport name="info" type="pyobj"/>
596 <inline name="ObservationNodeString">
597 <script><code><![CDATA[]]></code></script>
598 <inport name="var" type="pyobj"/>
599 <inport name="info" type="pyobj"/>
602 <inline name="ObservationNodeFile">
603 <script><code><![CDATA[execfile(script)]]></code></script>
604 <inport name="var" type="pyobj"/>
605 <inport name="info" type="pyobj"/>
606 <inport name="script" type="string"/>
609 <inline name="EndObservationNode">
610 <script><code><![CDATA[
611 # Create a fake output object.
612 # An observer is always successful.
614 output["outputValues"] = [[[[]]]]
615 output["specificOutputInfos"] = []
616 output["returnCode"] = 0
617 output["errorMessage"] = ""
619 <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
622 <inline name="SetObserversNode">
623 <script><code><![CDATA[
624 #print "Setting observers"
626 <outport name="has_observers" type="bool"/>
627 <outport name="observers" type="pyobj"/>