src/tests/daSalome/test017_3DVAR_function_script.py

   1 import numpy
   2 import pickle
   3
   4 print computation
   5 method = ""
   6 for param in computation["specificParameters"]:
   7   if param["name"] == "method":
   8     method = param["value"]
   9 print "Method found is", method
  10 print computation["inputValues"]
  11 print computation["inputValues"][0][0]
  12
  13 dimension = 300
  14 H  = numpy.matrix(numpy.core.identity(dimension))
  15
  16 def FunctionH( X ):
  17     return H * X
  18
  19 def AdjointH( (X, Y) ):
  20     return H.T * Y
  21
  22 if method == "Direct":
  23   result = FunctionH(numpy.matrix(computation["inputValues"][0][0]).T)
  24
  25 if method == "Tangent":
  26   result = FunctionH(numpy.matrix(computation["inputValues"][0][0]).T)
  27
  28 if method == "Adjoint":
  29   result = AdjointH(numpy.matrix(computation["inputValues"][0][0]).T)
  30
  31 print "Computation end"