return val
black = [0,0,0]
red = [255,0,0]
green = [0,255,0]
blue = [0,0,255]
w = 600
h = 600
# display list of line segments from contouring routine
gr = racg.Graphics(w=600,h=600)
gr.Clear()
s = 0.5*10**8
k = 4
x0,y0 = 200,100
gr.Clear()
G = lambda s: lambda x,y: F(x,y)-s
ix = 100
iy = 100
for y0 in range(1,600,iy):
for x0 in range(1,600,ix):
ctr.ShowLevelSet(gr,G(s), x0, y0, blue, tmax = 50)
gr.Point([x0,y0],red)
ch = gr.Show("result",15)
if ch == ord('e'):
break
gr.Show("result",-1)
fn = "ImplicitFunction-ContourVelocity-%02d.jpg" % k
#gr.Save(fn)
gr.Close()
black = [0, 0, 0]
red = [255, 0, 0]
green = [0, 255, 0]
blue = [0, 0, 255]
w = 600
h = 600
# display list of line segments from contouring routine
gr = racg.Graphics(w=600, h=600)
gr.Clear()
s = 0.5 * 10**8
k = 3
x0, y0 = 200, 100
gr.Clear()
G = lambda s: lambda x, y: F(x, y) - s
ix = 50
iy = 50
for y0 in range(1, 600, iy):
for x0 in range(1, 600, ix):
ctr.ShowLevelSet(gr, G(s), x0, y0, blue)
gr.Point([x0, y0], red)
ch = gr.Show("result", 15)
if ch == ord('e'):
break
gr.Show("result", -1)
fn = "ImplicitFunction-ContourVelocity-%02d.jpg" % k
#gr.Save(fn)
gr.Close()
