Changeset 3190 for production/sydney_2006/find_max.py

Timestamp:

Jun 21, 2006, 9:31:57 AM (18 years ago)

Author:

sexton

Message:

MOST and ANUGA comparisons and updates

File:

• production/sydney_2006/find_max.py (modified) (1 diff)

production/sydney_2006/find_max.py

@@ -35 +35 @@
 c = 0
 direction = 'negative'
+xstar = x0
+print 'start of loop', deriv
 while c < 100000 and deriv > 0:
     deriv = dfuncdx(x,0.83)
-    if deriv < 0: xstar = x
+    if deriv < 0:
+        xstar = x
+        print 'hello', deriv, xstar/lam0, c
     if direction == 'positive': x += step
     if direction == 'negative': x -= step
     c += 1
+    
 
 print 'location of maximum of surface elevation function: xstar = %f' % (xstar/lam0)
 const = 1.0 #a3D = ? #sydney = 86? and kappad=1
 
-x = arange(-20,25,0.001)
-y = arange(-10,10,0.1)
-X,Y = meshgrid(x,y)
+#x = arange(-20,25,0.001)
+#y = arange(-10,10,0.1)
+#X,Y = meshgrid(x,y)
 
-test = 0.0
-for xi in x:
-    testi = func(xi,0.0,0.83)
-    if direction == 'positive':
-        if testi > test:
-            test = testi
-            xstar = xi
-    else:
-        if testi < test:
-            test = testi
-            xstar = xi
+#test = 0.0
+#for xi in x:
+#    testi = func(xi,0.0,0.83)
+#    if direction == 'positive':
+#        if testi > test:
+#            test = testi
+#            xstar = xi
+#    else:
+#        if testi < test:
+#            test = testi
+#            xstar = xi
     
-print 'check: xstar = %f and eta(xstar) = %f' %(xstar/lam0, test)
+#print 'check: xstar = %f and eta(xstar) = %f' %(xstar/lam0, test)