# test_spiral.py find min and max for various ranges 
#                spiral(x,y) function included below
import math
from spiral_trough import spiral

print "test_spiral.py running"
xmin = -1.0
xmax = 1.0
ymin = -1.0
ymax = 1.0

def run_range():
  n=101

  xh = (xmax-xmin)/n
  yh = (ymax-ymin)/n

  z = spiral(xmin,ymin)
  zmin = z
  zmax = z
  atx = xmin
  aty = ymin

  for i in range(n):
    x = xmin+i*xh
    for j in range(n):
      y = ymin+j*yh
      z = spiral(x,y)
      if z<zmin:
        zmin = z
        atx = x
        aty = y

      if z>zmax:
        zmax = z
  return zmin, zmax, atx, aty

  # end run_range


zmin, zmax, atx, aty = run_range()
print "xmin=",
print  xmin,
print ", xmax=",
print  xmax,
print ", ymin=",
print  ymin,
print ", ymax=",
print  +ymax
print "zmin=",
print  zmin,
print ", zmax=",
print  zmax
print " atx=",
print  +atx,
print ",  aty=",
print  aty

xmin = -2.0
xmax = 2.0
ymin = -2.0
ymax = 2.0
zmin, zmax, atx, aty = run_range()
print "xmin=",
print  xmin,
print ", xmax=",
print  xmax,
print ", ymin=",
print  ymin,
print ", ymax=",
print  +ymax
print "zmin=",
print  zmin,
print ", zmax=",
print  zmax
print " atx=",
print  +atx,
print ",  aty=",
print  aty

xmin = -0.01
xmax = -1.0
ymin = 0.01
ymax = 1.0
zmin, zmax, atx, aty = run_range()
print "xmin=",
print  xmin,
print ", xmax=",
print  xmax,
print ", ymin=",
print  ymin,
print ", ymax=",
print  +ymax
print "zmin=",
print  zmin,
print ", zmax=",
print  zmax
print " atx=",
print  +atx,
print ",  aty=",
print  aty

xmin = 0.01
xmax = 1.0
ymin = 0.01
ymax = 1.0
zmin, zmax, atx, aty = run_range()
print "xmin=",
print  xmin,
print ", xmax=",
print  xmax,
print ", ymin=",
print  ymin,
print ", ymax=",
print  +ymax
print "zmin=",
print  zmin,
print ", zmax=",
print  zmax
print " atx=",
print  +atx,
print ",  aty=",
print  aty

xmin = -0.01
xmax = -1.0
ymin = -0.01
ymax = -1.0
zmin, zmax, atx, aty = run_range()
print "xmin=",
print  xmin,
print ", xmax=",
print  xmax,
print ", ymin=",
print  ymin,
print ", ymax=",
print  +ymax
print "zmin=",
print  zmin,
print ", zmax=",
print  zmax
print " atx=",
print  +atx,
print ",  aty=",
print  aty

# end test_spiral