src/tests/daSalome/test017.comm.in

   1
   2 ASSIMILATION_STUDY(Study_name='test017',
   3                    Study_repertory='@prefix@/tests/daSalome',
   4                    Debug=0,
   5                    Algorithm='3DVAR',
   6                    Background=_F(INPUT_TYPE='Vector',
   7                                  data=_F(FROM='Script',
   8                                          SCRIPT_FILE=
   9                            'test017_3DVAR_init_data.py',),),
  10                    BackgroundError=_F(INPUT_TYPE='Matrix',
  11                                       data=_F(FROM='Script',
  12                                               SCRIPT_FILE=
  13                            'test017_3DVAR_init_data.py',),),
  14                    Observation=_F(INPUT_TYPE='Vector',
  15                                   data=_F(FROM='Script',
  16                                           SCRIPT_FILE=
  17                            'test017_3DVAR_init_data.py',),),
  18                    ObservationError=_F(INPUT_TYPE='Matrix',
  19                                        data=_F(FROM='Script',
  20                                                SCRIPT_FILE=
  21                            'test017_3DVAR_init_data.py',),),
  22                    ObservationOperator=_F(INPUT_TYPE='Function',
  23                                           data=_F(FROM='FunctionDict',
  24                                                   FUNCTIONDICT_FILE=
  25                      'test017_3DVAR_function_script.py',),),
  26                    UserPostAnalysis=_F(FROM='String',
  27                                        STRING=
  28 """import numpy
  29 dimension = 300
  30 precision = 1.e-10
  31 xt = numpy.matrix(numpy.arange(dimension)).T
  32 xb = Study.getBackground()
  33 xa = numpy.array(ADD.get("Analysis")[0])
  34 d  = numpy.array(ADD.get("Innovation")[0])
  35 #
  36 # Verification du resultat
  37 # ------------------------
  38 if max(abs(xa - (xb+xt.A1)/2)) > precision:
  39     raise ValueError("Resultat du test errone (1)")
  40 else:
  41     print "    Test correct, erreur maximale inferieure à %s"%precision
  42 """,),);