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
  11 dimension = 300
  12 H  = numpy.matrix(numpy.core.identity(dimension))
  13
  14 def FunctionH( X ):
  15     return H * X
  16
  17 def AdjointH( (X, Y) ):
  18     return H.T * Y
  19
  20 if method == "Direct":
  21   result = FunctionH(numpy.matrix(computation["inputValues"][0][0]).T)
  22
  23 if method == "Tangent":
  24   result = FunctionH(numpy.matrix(computation["inputValues"][0][0]).T)
  25
  26 if method == "Adjoint":
  27   result = AdjointH((numpy.matrix(computation["inputValues"][0][0]).T, numpy.matrix(computation["inputValues"][0][1]).T))
  28
  29 print "Computation end"