2 ASSIMILATION_STUDY(Study_name='test017',
3 Study_repertory='@prefix@/tests/daSalome',
6 Background=_F(INPUT_TYPE='Vector',
9 'test017_3DVAR_init_data.py',),),
10 BackgroundError=_F(INPUT_TYPE='Matrix',
11 data=_F(FROM='Script',
13 'test017_3DVAR_init_data.py',),),
14 Observation=_F(INPUT_TYPE='Vector',
15 data=_F(FROM='Script',
17 'test017_3DVAR_init_data.py',),),
18 ObservationError=_F(INPUT_TYPE='Matrix',
19 data=_F(FROM='Script',
21 'test017_3DVAR_init_data.py',),),
22 ObservationOperator=_F(INPUT_TYPE='Function',
23 data=_F(FROM='FunctionDict',
25 'test017_3DVAR_function_script.py',),),
26 UserPostAnalysis=_F(FROM='String',
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])
36 # Verification du resultat
37 # ------------------------
38 if max(abs(xa - (xb+xt.A1)/2)) > precision:
39 raise ValueError("Resultat du test errone (1)")
41 print " Test correct, erreur maximale inferieure à %s"%precision