/[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 102 by jgs, Wed Dec 15 07:08:39 2004 UTC revision 150 by jgs, Thu Sep 15 03:44:45 2005 UTC
# Line 81  def checkResult(text,res,val0,val1,val2, Line 81  def checkResult(text,res,val0,val1,val2,
81       error=Lsup(ref-res)/norm       error=Lsup(ref-res)/norm
82       print "@@ %s, shape %s: error = %e"%(text,ref.getShape(),error)       print "@@ %s, shape %s: error = %e"%(text,ref.getShape(),error)
83       if error>tol:       if error>tol:
        #raise SystemError,"@@ %s at %s: error is too large"%(text,wh)  
84         print "**** %s: error is too large"%(text)         print "**** %s: error is too large"%(text)
85           raise SystemError,"@@ %s at %s: error is too large"%(text,wh)
86           sys.exit(1)
87    
88  def getRank(arg):  def getRank(arg):
89      if isinstance(arg,Data):      if isinstance(arg,Data):
# Line 186  for wh in [ContinuousFunction(msh),Funct Line 187  for wh in [ContinuousFunction(msh),Funct
187        # where zero:        # where zero:
188        ref=checkResult("where zero("+ex1+")", \        ref=checkResult("where zero("+ex1+")", \
189                        (arg1-3).whereZero(), \                        (arg1-3).whereZero(), \
190                        numarray.less_equal(numarray.abs(arrays1[0]-3.),EPSILON), \                        numarray.less_equal(numarray.abs(arrays1[0]-3.),0.0), \
191                        numarray.less_equal(numarray.abs(arrays1[1]-3.),EPSILON), \                        numarray.less_equal(numarray.abs(arrays1[1]-3.),0.0), \
192                        numarray.less_equal(numarray.abs(arrays1[2]-3.),EPSILON), \                        numarray.less_equal(numarray.abs(arrays1[2]-3.),0.0), \
193                        wh)                        wh)
194    
195        # where non-zero:        # where non-zero:
196        ref=checkResult("where nonzero("+ex1+")", \        ref=checkResult("where nonzero("+ex1+")", \
197                        (arg1-3).whereNonZero(), \                        (arg1-3).whereNonZero(), \
198                        numarray.greater(numarray.abs(arrays1[0]-3.),EPSILON), \                        numarray.greater(numarray.abs(arrays1[0]-3.),0.0), \
199                        numarray.greater(numarray.abs(arrays1[1]-3.),EPSILON), \                        numarray.greater(numarray.abs(arrays1[1]-3.),0.0), \
200                        numarray.greater(numarray.abs(arrays1[2]-3.),EPSILON), \                        numarray.greater(numarray.abs(arrays1[2]-3.),0.0), \
201                        wh)                        wh)
202    
       # transpose:  
       #axis=arrays1[0]/2  
       #ref=checkResult("transpose("+ex1+")", \  
       #                arg1.transpose(), \  
       #                numarray.transpose(arrays1[0],axis), \  
       #                numarray.transpose(arrays1[1],axis), \  
       #                numarray.transpose(arrays1[2],axis), \  
       #                wh)  
   
       # trace:  
       #ref=checkResult("trace("+ex1+")", \  
       #                arg1.trace(), \  
       #                numarray.trace(arrays1[0]), \  
       #                numarray.trace(arrays1[1]), \  
       #                numarray.trace(arrays1[2]), \  
       #                wh)  
   
203        # exponential function:        # exponential function:
204        ref=checkResult("exp("+ex1+")", \        ref=checkResult("exp("+ex1+")", \
205                        arg1.exp(), \                        arg1.exp(), \
# Line 248  for wh in [ContinuousFunction(msh),Funct Line 232  for wh in [ContinuousFunction(msh),Funct
232                        numarray.cos(arrays1[2]), \                        numarray.cos(arrays1[2]), \
233                        wh)                        wh)
234    
235          # tan:
236          ref=checkResult("tan("+ex1+")", \
237                          arg1.tan(), \
238                          numarray.tan(arrays1[0]), \
239                          numarray.tan(arrays1[1]), \
240                          numarray.tan(arrays1[2]), \
241                          wh)
242    
243          # asin:
244          #ref=checkResult("asin("+ex1+")", \
245          #                arg1.asin(), \
246          #                numarray.asin(arrays1[0]), \
247          #                numarray.asin(arrays1[1]), \
248          #                numarray.asin(arrays1[2]), \
249          #                wh)
250    
251          # acos:
252          #ref=checkResult("acos("+ex1+")", \
253          #                arg1.acos(), \
254          #                numarray.acos(arrays1[0]), \
255          #                numarray.acos(arrays1[1]), \
256          #                numarray.acos(arrays1[2]), \
257          #                wh)
258    
259          # atan:
260          #ref=checkResult("atan("+ex1+")", \
261          #                arg1.atan(), \
262          #                numarray.atan(arrays1[0]), \
263          #                numarray.atan(arrays1[1]), \
264          #                numarray.atan(arrays1[2]), \
265          #                wh)
266    
267          # sinh:
268          ref=checkResult("sinh("+ex1+")", \
269                          arg1.sinh(), \
270                          numarray.sinh(arrays1[0]), \
271                          numarray.sinh(arrays1[1]), \
272                          numarray.sinh(arrays1[2]), \
273                          wh)
274    
275          # cosh:
276          ref=checkResult("cosh("+ex1+")", \
277                          arg1.cosh(), \
278                          numarray.cosh(arrays1[0]), \
279                          numarray.cosh(arrays1[1]), \
280                          numarray.cosh(arrays1[2]), \
281                          wh)
282    
283          # tanh:
284          ref=checkResult("tanh("+ex1+")", \
285                          arg1.tanh(), \
286                          numarray.tanh(arrays1[0]), \
287                          numarray.tanh(arrays1[1]), \
288                          numarray.tanh(arrays1[2]), \
289                          wh)
290    
291          # asinh:
292          #ref=checkResult("asinh("+ex1+")", \
293          #                arg1.asinh(), \
294          #                numarray.asinh(arrays1[0]), \
295          #                numarray.asinh(arrays1[1]), \
296          #                numarray.asinh(arrays1[2]), \
297          #                wh)
298    
299          # acosh:
300          #ref=checkResult("acosh("+ex1+")", \
301          #                arg1.acosh(), \
302          #                numarray.acosh(arrays1[0]), \
303          #                numarray.acosh(arrays1[1]), \
304          #                numarray.acosh(arrays1[2]), \
305          #                wh)
306    
307          # atanh:
308          #ref=checkResult("atanh("+ex1+")", \
309          #                arg1.atanh(), \
310          #                numarray.atanh(arrays1[0]), \
311          #                numarray.atanh(arrays1[1]), \
312          #                numarray.atanh(arrays1[2]), \
313          #                wh)
314    
315        # get the maximum value at each data point        # get the maximum value at each data point
316        #ref=checkResult("maxval("+ex1+")", \        #ref=checkResult("maxval("+ex1+")", \
317        #                arg1.maxval(), \        #                arg1.maxval(), \
# Line 272  for wh in [ContinuousFunction(msh),Funct Line 336  for wh in [ContinuousFunction(msh),Funct
336        #                numarray.sqrt((arrays1[2]**2).sum()), \        #                numarray.sqrt((arrays1[2]**2).sum()), \
337        #                wh)        #                wh)
338    
339          # trace:
340          #ref=checkResult("trace("+ex1+")", \
341          #                arg1.trace(), \
342          #                numarray.trace(arrays1[0]), \
343          #                numarray.trace(arrays1[1]), \
344          #                numarray.trace(arrays1[2]), \
345          #                wh)
346    
347          # transpose:
348          #axis=arrays1[0]/2
349          #ref=checkResult("transpose("+ex1+")", \
350          #                arg1.transpose(), \
351          #                numarray.transpose(arrays1[0],axis), \
352          #                numarray.transpose(arrays1[1],axis), \
353          #                numarray.transpose(arrays1[2],axis), \
354          #                wh)
355    
356        # get the signs of the values:        # get the signs of the values:
357        ref=checkResult("sign("+ex1+")", \        ref=checkResult("sign("+ex1+")", \
358                        arg1.sign(), \                        arg1.sign(), \
# Line 283  for wh in [ContinuousFunction(msh),Funct Line 364  for wh in [ContinuousFunction(msh),Funct
364                          -numarray.less(arrays1[2],numarray.zeros(arrays1[2].shape)),\                          -numarray.less(arrays1[2],numarray.zeros(arrays1[2].shape)),\
365                        wh)                        wh)
366    
367    sys.exit(0)
368  # end  # end

Legend:
Removed from v.102  
changed lines
  Added in v.150

  ViewVC Help
Powered by ViewVC 1.1.26