summaryrefslogtreecommitdiff
path: root/source3/lua-5.1.4/test/bisect.lua
diff options
context:
space:
mode:
authorVolker Lendecke <vl@samba.org>2008-09-19 08:01:23 +0200
committerVolker Lendecke <vl@samba.org>2008-12-17 12:57:18 +0100
commit390f19b649b85a8f8ab762c6a4e9d1a3d23d3f47 (patch)
tree9f9506e36222640d23b086d62b32de940c526a98 /source3/lua-5.1.4/test/bisect.lua
parent081f8883bafc11ca7c08f868ec19f1ea32071837 (diff)
downloadsamba-390f19b649b85a8f8ab762c6a4e9d1a3d23d3f47.tar.gz
samba-390f19b649b85a8f8ab762c6a4e9d1a3d23d3f47.tar.bz2
samba-390f19b649b85a8f8ab762c6a4e9d1a3d23d3f47.zip
Add the Lua distibution from http://www.lua.org/ftp/lua-5.1.4.tar.gz
Available under the MIT license. Adding it to see how the build farm likes it. They claim to be 100% pure ANSI C and compile everywhere. Lets see. If it breaks badly, we can remove it again.
Diffstat (limited to 'source3/lua-5.1.4/test/bisect.lua')
-rw-r--r--source3/lua-5.1.4/test/bisect.lua27
1 files changed, 27 insertions, 0 deletions
diff --git a/source3/lua-5.1.4/test/bisect.lua b/source3/lua-5.1.4/test/bisect.lua
new file mode 100644
index 0000000000..f91e69bfba
--- /dev/null
+++ b/source3/lua-5.1.4/test/bisect.lua
@@ -0,0 +1,27 @@
+-- bisection method for solving non-linear equations
+
+delta=1e-6 -- tolerance
+
+function bisect(f,a,b,fa,fb)
+ local c=(a+b)/2
+ io.write(n," c=",c," a=",a," b=",b,"\n")
+ if c==a or c==b or math.abs(a-b)<delta then return c,b-a end
+ n=n+1
+ local fc=f(c)
+ if fa*fc<0 then return bisect(f,a,c,fa,fc) else return bisect(f,c,b,fc,fb) end
+end
+
+-- find root of f in the inverval [a,b]. needs f(a)*f(b)<0
+function solve(f,a,b)
+ n=0
+ local z,e=bisect(f,a,b,f(a),f(b))
+ io.write(string.format("after %d steps, root is %.17g with error %.1e, f=%.1e\n",n,z,e,f(z)))
+end
+
+-- our function
+function f(x)
+ return x*x*x-x-1
+end
+
+-- find zero in [1,2]
+solve(f,1,2)