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 assim_study.setBackgroundType(BackgroundType)
91 assim_study.setBackground(Background)
99 logging.debug("CREATE CheckingPoint is %s"%CheckingPoint)
100 logging.debug("CREATE CheckingPointType is %s"%CheckingPointType)
101 assim_study.setCheckingPointType(CheckingPointType)
102 assim_study.setCheckingPoint(CheckingPoint)
110 logging.debug("CREATE BackgroundError is %s"%BackgroundError)
111 logging.debug("CREATE BackgroundErrorType is %s"%BackgroundErrorType)
112 assim_study.setBackgroundError(BackgroundError)
120 logging.debug("CREATE Observation is %s"%Observation)
121 logging.debug("CREATE ObservationType is %s"%ObservationType)
122 assim_study.setObservationType(ObservationType)
123 assim_study.setObservation(Observation)
131 logging.debug("CREATE ObservationError is %s"%ObservationError)
132 logging.debug("CREATE ObservationErrorType is %s"%ObservationErrorType)
133 assim_study.setObservationError(ObservationError)
135 # ObservationOperator
136 ObservationOperatorOk = 0
142 logging.debug("CREATE ObservationOperator is %s"%ObservationOperator)
143 logging.debug("CREATE ObservationOperatorType is %s"%ObservationOperatorType)
144 assim_study.setObservationOperatorType("Matrix", ObservationOperatorType)
145 assim_study.setObservationOperator("Matrix", ObservationOperator)
146 ObservationOperatorOk = 1
148 if ObservationOperatorOk == 0:
150 ObservationOperatorDirect
154 logging.debug("CREATE ObservationOperatorDirect is %s"%ObservationOperatorDirect)
155 assim_study.setObservationOperatorType("Direct", "Function")
156 assim_study.setObservationOperator("Direct", ObservationOperatorDirect)
158 ObservationOperatorTangent
162 logging.debug("CREATE ObservationOperatorTangent is %s"%ObservationOperatorTangent)
163 assim_study.setObservationOperatorType("Tangent", "Function")
164 assim_study.setObservationOperator("Tangent", ObservationOperatorTangent)
166 ObservationOperatorAdjoint
170 logging.debug("CREATE ObservationOperatorAdjoint is %s"%ObservationOperatorAdjoint)
171 assim_study.setObservationOperatorType("Adjoint", "Function")
172 assim_study.setObservationOperator("Adjoint", ObservationOperatorAdjoint)
175 for name, size in zip(InputVariablesNames, InputVariablesSizes):
176 assim_study.setInputVariable(name, size)
177 for name, size in zip(OutputVariablesNames, OutputVariablesSizes):
178 assim_study.setOutputVariable(name, size)
181 logging.debug("CREATE Observers is %s"%observers.keys())
182 # Adding observers to the study
183 for observer_name in observers.keys():
186 number = str(observers[observer_name]["number"])
187 if "scheduler" in observers[observer_name].keys():
188 scheduler = observers[observer_name]["scheduler"]
189 if "info" in observers[observer_name].keys():
190 info = observers[observer_name]["info"]
191 assim_study.addObserver(observer_name, scheduler, info, number)
196 <inport name="Name" type="string"/>
197 <inport name="Algorithm" type="string"/>
198 <inport name="Debug" type="bool"/>
199 <inport name="InputVariablesNames" type="stringvec"/>
200 <inport name="InputVariablesSizes" type="intvec"/>
201 <inport name="OutputVariablesNames" type="stringvec"/>
202 <inport name="OutputVariablesSizes" type="intvec"/>
203 <inport name="has_observers" type="bool"/>
204 <inport name="observers" type="pyobj"/>
205 <outport name="Study" type="pyobj"/>
208 <inline name="CreateNumpyMatrixFromString">
209 <script><code><![CDATA[
210 import numpy, logging
211 logging.debug("CREATE Entering in CreateNumpyMatrixFromString")
212 matrix = numpy.matrix(matrix_in_string)
214 logging.debug("CREATE Matrix is %s"%matrix)
216 <inport name="matrix_in_string" type="string"/>
217 <outport name="matrix" type="pyobj"/>
218 <outport name="type" type="string"/>
221 <inline name="CreateNumpyMatrixFromScript">
222 <script><code><![CDATA[
224 logging.debug("CREATE Entering in CreateNumpyMatrixFromScript")
227 # Get file path and filename
230 filepath = os.path.dirname(script)
231 filename = os.path.basename(script)
232 module_name = os.path.splitext(filename)[0]
233 sys.path.insert(0,filepath)
236 __import__(module_name)
237 user_script_module = sys.modules[module_name]
239 # Get Data from script
241 <inport name="script" type="string"/>
242 <outport name="type" type="string"/>
245 <inline name="CreateNumpyVectorFromString">
246 <script><code><![CDATA[
247 import numpy, logging
248 logging.debug("CREATE Entering in CreateNumpyVectorFromString")
249 vector = numpy.matrix(vector_in_string)
251 logging.debug("Vector is %s"%vector)
253 <inport name="vector_in_string" type="string"/>
254 <outport name="vector" type="pyobj"/>
255 <outport name="type" type="string"/>
258 <inline name="CreateNumpyVectorFromScript">
259 <script><code><![CDATA[
261 logging.debug("CREATE Entering in CreateNumpyVectorFromScript")
264 # Get file path and filename
267 filepath = os.path.dirname(script)
268 filename = os.path.basename(script)
269 module_name = os.path.splitext(filename)[0]
270 sys.path.insert(0,filepath)
273 __import__(module_name)
274 user_script_module = sys.modules[module_name]
276 # Get Data from script
278 <inport name="script" type="string"/>
279 <outport name="type" type="string"/>
282 <inline name="SimpleExecuteDirectAlgorithm">
283 <script><code><![CDATA[
285 logging.debug("EXECUTE Entering in SimpleExecuteDirectAlgorithm")
286 from daYacsIntegration.daStudy import *
287 ADD = Study.getAssimilationStudy()
290 <inport name="Study" type="pyobj"/>
291 <outport name="Study" type="pyobj"/>
294 <inline name="SimpleUserAnalysis">
295 <script><code><![CDATA[
296 #-*-coding:iso-8859-1-*-
298 logging.debug("TERMINATE Entering in SimpleUserAnalysis")
299 from daYacsIntegration.daStudy import *
300 ADD = Study.getAssimilationStudy()
304 <inport name="Study" type="pyobj"/>
307 <inline name="FakeOptimizerLoopNode">
308 <script><code><![CDATA[
310 logging.debug("EXECUTE Entering in FakeOptimizerLoopNode")
313 <inport name="computation" type="SALOME_TYPES/ParametricInput"/>
314 <outport name="result" type="SALOME_TYPES/ParametricOutput"/>
317 <inline name="CreateDictFromScript">
318 <script><code><![CDATA[
320 logging.debug("CREATE Entering in CreateDictFromScript")
322 # Get file path and filename
325 filepath = os.path.dirname(script)
326 filename = os.path.basename(script)
327 module_name = os.path.splitext(filename)[0]
328 sys.path.insert(0,filepath)
331 __import__(module_name)
332 user_script_module = sys.modules[module_name]
334 # Get Data from script
336 <inport name="script" type="string"/>
339 <inline name="UserDataInitFromScript">
340 <script><code><![CDATA[
342 logging.debug("CREATE Entering in UserDataInitFromScript")
344 # Get file path and filename
347 filepath = os.path.dirname(script)
348 filename = os.path.basename(script)
349 module_name = os.path.splitext(filename)[0]
350 sys.path.insert(0,filepath)
353 __import__(module_name)
354 user_script_module = sys.modules[module_name]
356 # Get Data from script
358 <inport name="script" type="string"/>
359 <outport name="init_data" type="pyobj"/>
362 <inline name="ReadForSwitchNode">
363 <script><code><![CDATA[
365 logging.debug("CREATE Entering in ReadForSwitchNode")
367 for param in data["specificParameters"]:
368 if param["name"] == "switch_value":
369 switch_value = int(param["value"])
371 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
372 <outport name="data" type="SALOME_TYPES/ParametricInput"/>
373 <outport name="switch_value" type="int"/>
376 <inline name="ExtractDataNode">
377 <script><code><![CDATA[
379 logging.debug("TERMINATE Entering in ExtractDataNode")
381 from daCore.AssimilationStudy import AssimilationStudy
384 for param in data["specificParameters"]:
385 if param["name"] == "var":
386 var = pickle.loads(param["value"])
387 if param["name"] == "info":
388 info = param["value"]
390 <inport name="data" type="SALOME_TYPES/ParametricInput"/>
391 <outport name="var" type="pyobj"/>
392 <outport name="info" type="pyobj"/>
395 <inline name="ObservationNodeString">
396 <script><code><![CDATA[
397 #print "Entering in Observation"
400 <inport name="var" type="pyobj"/>
401 <inport name="info" type="pyobj"/>
404 <inline name="ObservationNodeFile">
405 <script><code><![CDATA[
406 #print "Entering in Observation"
410 <inport name="var" type="pyobj"/>
411 <inport name="info" type="pyobj"/>
412 <inport name="script" type="string"/>
415 <inline name="EndObservationNode">
416 <script><code><![CDATA[
417 # Create a fake output object.
418 # An observer is always successful.
420 output["outputValues"] = [[[[]]]]
421 output["specificOutputInfos"] = []
422 output["returnCode"] = 0
423 output["errorMessage"] = ""
425 <outport name="output" type="SALOME_TYPES/ParametricOutput"/>
428 <inline name="SetObserversNode">
429 <script><code><![CDATA[
430 #print "Setting observers"
432 <outport name="has_observers" type="bool"/>
433 <outport name="observers" type="pyobj"/>