/[escript]/trunk/escript/test/python/UnaryOps.py
ViewVC logotype

Diff of /trunk/escript/test/python/UnaryOps.py

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 559 by jgs, Tue Feb 21 01:59:07 2006 UTC revision 560 by jgs, Thu Feb 23 01:58:35 2006 UTC
# Line 25  tol=1.E-15 Line 25  tol=1.E-15
25  #  #
26  #  list of arguments: a list item has the form [a0,a1,a2]  #  list of arguments: a list item has the form [a0,a1,a2]
27  #  what a0 is the default value and a1 is used for tag Tag1  #  what a0 is the default value and a1 is used for tag Tag1
28  #  and a2 for tag2. a0,a1,a2 are converted into numarrays.  #  and a2 for Tag2.
29  #    #  
30    
31  arglist = [ \  arglist = [ \
# Line 33  arglist = [ \ Line 33  arglist = [ \
33  [ [[1,2],[3,4]], [[5,6],[7,8]], [[-5,-6],[7,8]] ], \  [ [[1,2],[3,4]], [[5,6],[7,8]], [[-5,-6],[7,8]] ], \
34  [ [[15,8],[12,8]], [[-9,9],[13,8]], [[7,34],[19,7]] ], \  [ [[15,8],[12,8]], [[-9,9],[13,8]], [[7,34],[19,7]] ], \
35  [ [[[15,8],[12,8]],[[-9,9],[13,8]]], [[[3,4],[-9,4]],[[1,-9],[7,4]]], [[[5,2],[6,2]],[[-6,4],[7,5]]] ], \  [ [[[15,8],[12,8]],[[-9,9],[13,8]]], [[[3,4],[-9,4]],[[1,-9],[7,4]]], [[[5,2],[6,2]],[[-6,4],[7,5]]] ], \
36  [ 3.0, 6.0, 3 ] \  [ [3.0], [6.0], [3] ] \
37  ]  ]
38    
39  def prepareArg(val,ex,wh):  def prepareArg(val,ex,wh):
# Line 141  for wh in [ContinuousFunction(msh),Funct Line 141  for wh in [ContinuousFunction(msh),Funct
141    
142        # where positive:        # where positive:
143        ref=checkResult("where positive("+ex1+")", \        ref=checkResult("where positive("+ex1+")", \
144                        arg1._wherePositive(), \                        (arg1-3)._wherePositive(), \
145                        numarray.greater(arrays1[0],0.), \                        numarray.greater(arrays1[0],3.), \
146                        numarray.greater(arrays1[1],0.), \                        numarray.greater(arrays1[1],3.), \
147                        numarray.greater(arrays1[2],0.), \                        numarray.greater(arrays1[2],3.), \
148                        wh)                        wh)
149    
150        # where negative:        # where negative:
151        #ref=checkResult("where negative("+ex1+")", \        ref=checkResult("where negative("+ex1+")", \
152        #                (arg1-3)._whereNegative(), \                        (arg1-3)._whereNegative(), \
153        #                numarray.greater(3.,arrays1[0]), \                        numarray.greater(3.,arrays1[0]), \
154        #                numarray.greater(3.,arrays1[1]), \                        numarray.greater(3.,arrays1[1]), \
155        #                numarray.greater(3.,arrays1[2]), \                        numarray.greater(3.,arrays1[2]), \
156        #                wh)                        wh)
157    
158        # where non-negative:        # where non-negative:
159        #ref=checkResult("where nonnegative("+ex1+")", \        ref=checkResult("where nonnegative("+ex1+")", \
160        #                (arg1-3)._whereNonNegative(), \                        (arg1-3)._whereNonNegative(), \
161        #                numarray.greater_equal(arrays1[0],3.), \                        numarray.greater_equal(arrays1[0],3.), \
162        #                numarray.greater_equal(arrays1[1],3.), \                        numarray.greater_equal(arrays1[1],3.), \
163        #                numarray.greater_equal(arrays1[2],3.), \                        numarray.greater_equal(arrays1[2],3.), \
164        #                wh)                        wh)
165    
166        # where non-positive:        # where non-positive:
167        #ref=checkResult("where nonpositive("+ex1+")", \        ref=checkResult("where nonpositive("+ex1+")", \
168        #                (arg1-3)._whereNonPositive(), \                        (arg1-3)._whereNonPositive(), \
169        #                numarray.greater_equal(3.,arrays1[0]), \                        numarray.greater_equal(3.,arrays1[0]), \
170        #                numarray.greater_equal(3.,arrays1[1]), \                        numarray.greater_equal(3.,arrays1[1]), \
171        #                numarray.greater_equal(3.,arrays1[2]), \                        numarray.greater_equal(3.,arrays1[2]), \
172        #                wh)                        wh)
173    
174        # where zero:        # where zero:
175        #ref=checkResult("where zero("+ex1+")", \        ref=checkResult("where zero("+ex1+")", \
176        #                (arg1-3)._whereZero(), \                        (arg1-3)._whereZero(), \
177        #                numarray.less_equal(numarray.abs(arrays1[0]-3.),0.0), \                        numarray.less_equal(numarray.abs(arrays1[0]-3.),0.0), \
178        #                numarray.less_equal(numarray.abs(arrays1[1]-3.),0.0), \                        numarray.less_equal(numarray.abs(arrays1[1]-3.),0.0), \
179        #                numarray.less_equal(numarray.abs(arrays1[2]-3.),0.0), \                        numarray.less_equal(numarray.abs(arrays1[2]-3.),0.0), \
180        #                wh)                        wh)
181    
182        # where non-zero:        # where non-zero:
183        #ref=checkResult("where nonzero("+ex1+")", \        ref=checkResult("where nonzero("+ex1+")", \
184        #                (arg1-3)._whereNonZero(), \                        (arg1-3)._whereNonZero(), \
185        #                numarray.greater(numarray.abs(arrays1[0]-3.),0.0), \                        numarray.greater(numarray.abs(arrays1[0]-3.),0.0), \
186        #                numarray.greater(numarray.abs(arrays1[1]-3.),0.0), \                        numarray.greater(numarray.abs(arrays1[1]-3.),0.0), \
187        #                numarray.greater(numarray.abs(arrays1[2]-3.),0.0), \                        numarray.greater(numarray.abs(arrays1[2]-3.),0.0), \
188        #                wh)                        wh)
189    
190        # exponential function:        # exponential function:
191        ref=checkResult("exp("+ex1+")", \        ref=checkResult("exp("+ex1+")", \
# Line 196  for wh in [ContinuousFunction(msh),Funct Line 196  for wh in [ContinuousFunction(msh),Funct
196                        wh)                        wh)
197    
198        # sqrt        # sqrt
199        #ref=checkResult("sqrt("+ex1+")", \        ref=checkResult("sqrt("+ex1+")", \
200        #                arg1._abs()._sqrt(), \                        arg1._sqrt(), \
201        #                numarray.sqrt(numarray.abs(arrays1[0])), \                        numarray.sqrt(numarray.abs(arrays1[0])), \
202        #                numarray.sqrt(numarray.abs(arrays1[1])), \                        numarray.sqrt(numarray.abs(arrays1[1])), \
203        #                numarray.sqrt(numarray.abs(arrays1[2])), \                        numarray.sqrt(numarray.abs(arrays1[2])), \
204        #                wh)                        wh)
205    
206        # sin:        # sin:
207        ref=checkResult("sin("+ex1+")", \        ref=checkResult("sin("+ex1+")", \
# Line 227  for wh in [ContinuousFunction(msh),Funct Line 227  for wh in [ContinuousFunction(msh),Funct
227                        numarray.tan(arrays1[2]), \                        numarray.tan(arrays1[2]), \
228                        wh)                        wh)
229    
230          # numarray has no asin/acos/atan funcs
231    
232        # asin:        # asin:
233        #ref=checkResult("asin("+ex1+")", \        #ref=checkResult("asin("+ex1+")", \
234        #                arg1.asin(), \        #                arg1._asin(), \
235        #                numarray.asin(arrays1[0]), \        #                numarray.asin(arrays1[0]), \
236        #                numarray.asin(arrays1[1]), \        #                numarray.asin(arrays1[1]), \
237        #                numarray.asin(arrays1[2]), \        #                numarray.asin(arrays1[2]), \
# Line 237  for wh in [ContinuousFunction(msh),Funct Line 239  for wh in [ContinuousFunction(msh),Funct
239    
240        # acos:        # acos:
241        #ref=checkResult("acos("+ex1+")", \        #ref=checkResult("acos("+ex1+")", \
242        #                arg1.acos(), \        #                arg1._acos(), \
243        #                numarray.acos(arrays1[0]), \        #                numarray.acos(arrays1[0]), \
244        #                numarray.acos(arrays1[1]), \        #                numarray.acos(arrays1[1]), \
245        #                numarray.acos(arrays1[2]), \        #                numarray.acos(arrays1[2]), \
# Line 245  for wh in [ContinuousFunction(msh),Funct Line 247  for wh in [ContinuousFunction(msh),Funct
247    
248        # atan:        # atan:
249        #ref=checkResult("atan("+ex1+")", \        #ref=checkResult("atan("+ex1+")", \
250        #                arg1.atan(), \        #                arg1._atan(), \
251        #                numarray.atan(arrays1[0]), \        #                numarray.atan(arrays1[0]), \
252        #                numarray.atan(arrays1[1]), \        #                numarray.atan(arrays1[1]), \
253        #                numarray.atan(arrays1[2]), \        #                numarray.atan(arrays1[2]), \
# Line 275  for wh in [ContinuousFunction(msh),Funct Line 277  for wh in [ContinuousFunction(msh),Funct
277                        numarray.tanh(arrays1[2]), \                        numarray.tanh(arrays1[2]), \
278                        wh)                        wh)
279    
280          # numarray has no asinh/acosh/atanh funcs
281    
282        # asinh:        # asinh:
283        #ref=checkResult("asinh("+ex1+")", \        #ref=checkResult("asinh("+ex1+")", \
284        #                arg1.asinh(), \        #                arg1._asinh(), \
285        #                numarray.asinh(arrays1[0]), \        #                numarray.asinh(arrays1[0]), \
286        #                numarray.asinh(arrays1[1]), \        #                numarray.asinh(arrays1[1]), \
287        #                numarray.asinh(arrays1[2]), \        #                numarray.asinh(arrays1[2]), \
# Line 285  for wh in [ContinuousFunction(msh),Funct Line 289  for wh in [ContinuousFunction(msh),Funct
289    
290        # acosh:        # acosh:
291        #ref=checkResult("acosh("+ex1+")", \        #ref=checkResult("acosh("+ex1+")", \
292        #                arg1.acosh(), \        #                arg1._acosh(), \
293        #                numarray.acosh(arrays1[0]), \        #                numarray.acosh(arrays1[0]), \
294        #                numarray.acosh(arrays1[1]), \        #                numarray.acosh(arrays1[1]), \
295        #                numarray.acosh(arrays1[2]), \        #                numarray.acosh(arrays1[2]), \
# Line 293  for wh in [ContinuousFunction(msh),Funct Line 297  for wh in [ContinuousFunction(msh),Funct
297    
298        # atanh:        # atanh:
299        #ref=checkResult("atanh("+ex1+")", \        #ref=checkResult("atanh("+ex1+")", \
300        #                arg1.atanh(), \        #                arg1._atanh(), \
301        #                numarray.atanh(arrays1[0]), \        #                numarray.atanh(arrays1[0]), \
302        #                numarray.atanh(arrays1[1]), \        #                numarray.atanh(arrays1[1]), \
303        #                numarray.atanh(arrays1[2]), \        #                numarray.atanh(arrays1[2]), \
304        #                wh)        #                wh)
305    
306        # get the maximum value at each data point        # get the maximum value at each data point
307        #ref=checkResult("maxval("+ex1+")", \        ref=checkResult("maxval("+ex1+")", \
308        #                arg1.maxval(), \                        arg1._maxval(), \
309        #                arrays1[0].max(), \                        arrays1[0].max(), \
310        #                arrays1[1].max(), \                        arrays1[1].max(), \
311        #                arrays1[2].max(), \                        arrays1[2].max(), \
312        #                wh)                        wh)
313    
314        # get the minimum value at each data point        # get the minimum value at each data point
315        #ref=checkResult("minval("+ex1+")", \        ref=checkResult("minval("+ex1+")", \
316        #                arg1.minval(), \                        arg1._minval(), \
317        #                arrays1[0].min(), \                        arrays1[0].min(), \
318        #                arrays1[1].min(), \                        arrays1[1].min(), \
319        #                arrays1[2].min(), \                        arrays1[2].min(), \
320        #                wh)                        wh)
321    
322          # length/trace/transpose not yet implemented for Data
323    
324        # get the length at each data point = sqrt(sum_{i,j,k,l} A[i,j,k,l]^2)        # get the length at each data point = sqrt(sum_{i,j,k,l} A[i,j,k,l]^2)
325        #ref=checkResult("length("+ex1+")", \        #ref=checkResult("length("+ex1+")", \
326        #                arg1.length(), \        #                arg1._length(), \
327        #                numarray.sqrt((arrays1[0]**2).sum()), \        #                numarray.sqrt((arrays1[0]**2).sum()), \
328        #                numarray.sqrt((arrays1[1]**2).sum()), \        #                numarray.sqrt((arrays1[1]**2).sum()), \
329        #                numarray.sqrt((arrays1[2]**2).sum()), \        #                numarray.sqrt((arrays1[2]**2).sum()), \
# Line 325  for wh in [ContinuousFunction(msh),Funct Line 331  for wh in [ContinuousFunction(msh),Funct
331    
332        # trace:        # trace:
333        #ref=checkResult("trace("+ex1+")", \        #ref=checkResult("trace("+ex1+")", \
334        #                arg1.trace(), \        #                arg1._trace(), \
335        #                numarray.trace(arrays1[0]), \        #                numarray.trace(arrays1[0]), \
336        #                numarray.trace(arrays1[1]), \        #                numarray.trace(arrays1[1]), \
337        #                numarray.trace(arrays1[2]), \        #                numarray.trace(arrays1[2]), \
# Line 334  for wh in [ContinuousFunction(msh),Funct Line 340  for wh in [ContinuousFunction(msh),Funct
340        # transpose:        # transpose:
341        #axis=arrays1[0]/2        #axis=arrays1[0]/2
342        #ref=checkResult("transpose("+ex1+")", \        #ref=checkResult("transpose("+ex1+")", \
343        #                arg1.transpose(), \        #                arg1._transpose(), \
344        #                numarray.transpose(arrays1[0],axis), \        #                numarray.transpose(arrays1[0],axis), \
345        #                numarray.transpose(arrays1[1],axis), \        #                numarray.transpose(arrays1[1],axis), \
346        #                numarray.transpose(arrays1[2],axis), \        #                numarray.transpose(arrays1[2],axis), \

Legend:
Removed from v.559  
changed lines
  Added in v.560

  ViewVC Help
Powered by ViewVC 1.1.26