examples/daSalome/test006_Observers.py.in

   1 #-*-coding:iso-8859-1-*-
   2 study_config = {}
   3 study_config['StudyType'] = 'ASSIMILATION_STUDY'
   4 study_config['Name'] = 'test_observers'
   5 study_config['Debug'] = '0'
   6 study_config['Algorithm'] = '3DVAR'
   7 AlgorithmParameters_config = {}
   8 AlgorithmParameters_config['Type'] = 'Dict'
   9 AlgorithmParameters_config['From'] = 'Script'
  10 AlgorithmParameters_config['Data'] = 'test006_Observers_var.py'
  11 study_config['AlgorithmParameters'] = AlgorithmParameters_config
  12 Background_config = {}
  13 Background_config['Type'] = 'Vector'
  14 Background_config['From'] = 'Script'
  15 Background_config['Data'] = 'test006_Observers_var.py'
  16 Background_config['Stored'] = '0'
  17 study_config['Background'] = Background_config
  18 BackgroundError_config = {}
  19 BackgroundError_config['Type'] = 'Matrix'
  20 BackgroundError_config['From'] = 'Script'
  21 BackgroundError_config['Data'] = 'test006_Observers_var.py'
  22 BackgroundError_config['Stored'] = '0'
  23 study_config['BackgroundError'] = BackgroundError_config
  24 Observation_config = {}
  25 Observation_config['Type'] = 'Vector'
  26 Observation_config['From'] = 'Script'
  27 Observation_config['Data'] = 'test006_Observers_var.py'
  28 Observation_config['Stored'] = '0'
  29 study_config['Observation'] = Observation_config
  30 ObservationError_config = {}
  31 ObservationError_config['Type'] = 'Matrix'
  32 ObservationError_config['From'] = 'Script'
  33 ObservationError_config['Data'] = 'test006_Observers_var.py'
  34 ObservationError_config['Stored'] = '0'
  35 study_config['ObservationError'] = ObservationError_config
  36 ObservationOperator_ScriptWithSwitch = {}
  37 ObservationOperator_ScriptWithSwitch['Function'] = ['Direct', 'Tangent', 'Adjoint']
  38 ObservationOperator_ScriptWithSwitch['Script'] = {}
  39 ObservationOperator_ScriptWithSwitch['Script']['Direct'] = 'test006_Observers_Observation_Operator.py'
  40 ObservationOperator_ScriptWithSwitch['Script']['Tangent'] = 'test006_Observers_Observation_Operator.py'
  41 ObservationOperator_ScriptWithSwitch['Script']['Adjoint'] = 'test006_Observers_Observation_Operator.py'
  42 ObservationOperator_config = {}
  43 ObservationOperator_config['Type'] = 'Function'
  44 ObservationOperator_config['From'] = 'ScriptWithSwitch'
  45 ObservationOperator_config['Data'] = ObservationOperator_ScriptWithSwitch
  46 study_config['ObservationOperator'] = ObservationOperator_config
  47 inputvariables_config = {}
  48 inputvariables_config['Order'] =['adao_default']
  49 inputvariables_config['adao_default'] = -1
  50 study_config['InputVariables'] = inputvariables_config
  51 outputvariables_config = {}
  52 outputvariables_config['Order'] = ['adao_default']
  53 outputvariables_config['adao_default'] = -1
  54 study_config['OutputVariables'] = outputvariables_config
  55 study_config['Repertory'] = '@prefix@/share/salome/adao_examples/daSalome'
  56 Analysis_config = {}
  57 Analysis_config['From'] = 'String'
  58 Analysis_config['Data'] = """import numpy
  59 xa=numpy.ravel(ADD.get('Analysis')[-1])
  60 print 'Analysis:',xa"""
  61 study_config['UserPostAnalysis'] = Analysis_config
  62 Init_config = {}
  63 Init_config['Type'] = 'Dict'
  64 Init_config['From'] = 'Script'
  65 Init_config['Data'] = 'test006_Observers_init.py'
  66 Init_config['Target'] = ['Background','BackgroundError','Observation','ObservationError','AlgorithmParameters',]
  67 study_config['UserDataInit'] = Init_config
  68 observers = {}
  69 observers["CostFunctionJ"] = {}
  70 observers["CostFunctionJ"]["number"] = 3
  71 observers["CostFunctionJ"]["nodetype"] = "String"
  72 observers["CostFunctionJ"]["String"] = """print "  ---> observerCost"
  73 print "       var  =",var[:]
  74 print "       info =",info
  75 #
  76 import Gnuplot
  77 import os
  78 try:
  79   numero
  80 except NameError:
  81   numero = 0
  82 gp = Gnuplot.Gnuplot()
  83 gp('set style data lines')
  84 gp('set title  "'+str(info)+'"')
  85 gp.plot( Gnuplot.Data( var[:] ) )
  86 filename = os.path.join("/tmp", "imageCost_%02i.ps"%numero)
  87 print "       imageCost %s"%filename
  88 gp.hardcopy(filename=filename, color=1)
  89 numero += 1
  90 """
  91 observers["CostFunctionJ"]["info"] = """CostFunctionJ"""
  92 observers["CurrentState"] = {}
  93 observers["CurrentState"]["number"] = 4
  94 observers["CurrentState"]["nodetype"] = "Script"
  95 observers["CurrentState"]["Script"] = "test006_Observers_observer_with_file.py"
  96 observers["CurrentState"]["info"] = """CurrentState"""
  97 study_config['Observers'] = observers