src/tests/daSalome/test017_3DVAR_init_data.py

   1 import numpy
   2
   3 numpy.random.seed(1000)
   4 dimension = 300
   5
   6 xt = numpy.matrix(numpy.arange(dimension)).T
   7 Eo = numpy.matrix(numpy.zeros((dimension,))).T
   8 Eb = numpy.matrix(numpy.random.normal(0.,1.,size=(dimension,))).T
   9 H  = numpy.matrix(numpy.core.identity(dimension))
  10 B = numpy.matrix(numpy.core.identity(dimension)).T
  11 R = numpy.matrix(numpy.core.identity(dimension)).T
  12
  13 def FunctionH( X ):
  14     return H * X
  15
  16 xb = xt + Eb
  17 xb = xb.A1
  18 yo = FunctionH( xt ) + Eo
  19 yo = yo.A1
  20
  21 Background = xb
  22 BackgroundError = B
  23 Observation = yo
  24 ObservationError = R