Changeset 7926

Timestamp:

Aug 6, 2010, 1:21:08 PM (15 years ago)

Author:

mungkasi

Message:

 

File:

• trunk/anuga_work/development/sudi/sw_1d/periodic_waves/johns/bisect.py (modified) (1 diff)

trunk/anuga_work/development/sudi/sw_1d/periodic_waves/johns/bisect.py

@@ -16 +16 @@
     if f1*f2 > 0.0: error.err('Root is not bracketed')
     n = ceil(log(abs(x2 - x1)/epsilon)/log(2.0))
+    n = int(n)
     for i in range(n):
         x3 = 0.5*(x1 + x2); f3 = f(x3)