src/tests/daSalome/test017_3DVAR_function_script.py

   1 import numpy
   2 import pickle
   3
   4 print computation["method"]
   5
   6 dimension = 300
   7
   8 H  = numpy.matrix(numpy.core.identity(dimension))
   9
  10 def FunctionH( X ):
  11     return H * X
  12
  13 def AdjointH( (X, Y) ):
  14     return H.T * Y
  15
  16 if computation["method"] == "Direct":
  17   result = FunctionH(computation["data"])
  18
  19 if computation["method"] == "Tangent":
  20   result = FunctionH(computation["data"])
  21
  22 if computation["method"] == "Adjoint":
  23   result = AdjointH(computation["data"])
  24
  25 print "Computation end"