doc/source/examples/damagecontinuity.rst

   1 #!/usr/bin/python
   2 from ion.maths.damagecontinuity2 import totaldmg, adjustvsArmor
   3 from numpy import arange
   4 import numpy
   5 import sys
   6 minx, maxx = 0.09, 0.9901
   7 miny, maxy = 0, 100.01
   8
   9 scale = 144.0 # points
  10 iterations = 96
  11 xv, yv, data, color= [],[],[],[]
  12 for x in arange(minx, maxx, 0.05):
  13     sys.stdout.write ('-'); sys.stdout.flush()
  14     for y in arange (miny, maxy, 5):
  15 #        xv.append (x)
  16 #        yv.append (y)
  17         sum = 0.0
  18         minv = 9000000
  19         maxv = -9000000
  20         for i in range(iterations):
  21             tmp = totaldmg(x,y)
  22             sum += tmp
  23             minv = min (tmp, minv)
  24             maxv = max (tmp, maxv)
  25         xv.append (x)
  26         yv.append (y)
  27         data.append (maxv)
  28         color.append ('y')
  29         xv.append(x)
  30         yv.append(y)
  31         data.append (minv)
  32         color.append ('b')
  33         sys.stdout.write ('|'); sys.stdout.flush()
  34 import pylab
  35 maxdata = max(data)
  36 factor = scale / maxdata
  37 import math
  38 data = [(math.sqrt(v) + v) / 4.0 for v in data]
  39 #data = [v * factor for v in data]
  40
  41 pylab.scatter (xv, yv, data, color = color, marker = 'o', linewidth = 0.0, alpha = .5)
  42 pylab.title ('max data value : %.2f' % maxdata)
  43 pylab.show()