Mandelbrot example from benchmarks game

Signed-off-by: BlackMinou <romain.chanoir@viacesi.fr>

diff --git a/lib/actors/examples/mandelbrot/makefile b/lib/actors/examples/mandelbrot/makefile
new file mode 100644 (file)
index 0000000..cd626b7
--- /dev/null
+++ b/lib/actors/examples/mandelbrot/makefile
@@ -0,0 +1,16 @@
+file= mandelbrot
+
+default: threaded
+
+threaded:
+       nitc $(file).nit
+
+test:
+       ./$(file) 200
+
+bm:
+       time ./$(file) 16000
+
+clean:
+       rm $(file)
+       rm actors_$(file).nit

diff --git a/lib/actors/examples/mandelbrot/mandelbrot.nit b/lib/actors/examples/mandelbrot/mandelbrot.nit
new file mode 100644 (file)
index 0000000..316b54a
--- /dev/null
+++ b/lib/actors/examples/mandelbrot/mandelbrot.nit
@@ -0,0 +1,121 @@
+# This file is part of NIT (http://www.nitlanguage.org).
+#
+# Licensed under the Apache License, Version 2.0 (the "License");
+# you may not use this file except in compliance with the License.
+# You may obtain a copy of the License at
+#
+#     http://www.apache.org/licenses/LICENSE-2.0
+#
+# Unless required by applicable law or agreed to in writing, software
+# distributed under the License is distributed on an "AS IS" BASIS,
+# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
+# See the License for the specific language governing permissions and
+# limitations under the License.
+
+# Example implemented from "The computer Language Benchmarks Game" - Mandelbrot
+# http://benchmarksgame.alioth.debian.org/
+#
+# Complete description of mandelbrot :
+# https://benchmarksgame.alioth.debian.org/u64q/mandelbrot-description.html#mandelbrot
+module mandelbrot is no_warning("missing-doc")
+
+import actors
+
+class Worker
+       actor
+
+       fun get_byte(x, y: Int): Int do
+               var res = 0
+               for i in [0..8[.step(2) do
+                       var zr1 = crb[x + i]
+                       var zi1 = cib[y]
+
+                       var zr2 = crb [x + i + 1]
+                       var zi2 = cib[y]
+
+                       var b = 0
+                       for j in [0..49[ do
+                               var nzr1 = zr1 * zr1 - zi1 * zi1 + crb[x+i]
+                               var nzi1 = zr1 * zi1 + zr1 * zi1 + cib[y]
+                               zr1 = nzr1
+                               zi1 = nzi1
+
+                               var nzr2 = zr2 * zr2 - zi2 * zi2 + crb[x + i + 1]
+                               var nzi2 = zr2 * zi2 + zr2 * zi2 + cib[y]
+                               zr2 = nzr2
+                               zi2 = nzi2
+
+                               if zr1 * zr1 + zi1 * zi1 > 4.0 then b = b | 2
+                               if zr2 * zr2 + zi2 * zi2 > 4.0 then b = b | 1
+                               if b == 3 then break
+                       end
+                       res = (res << 2) + b
+               end
+               return res ^-1
+       end
+
+       fun put_line(y: Int, line: Array[Byte]) do
+               for i in [0..line.length[ do line[i] = get_byte(i * 8, y).to_b
+       end
+
+       fun work do
+               var line = 0
+               loop
+                       line = atomic.get_and_increment
+                       if line < n then
+                               put_line(line, data[line])
+                       else
+                               break
+                       end
+               end
+       end
+end
+
+redef class Sys
+       var n = 0
+       var inv_n: Float is noautoinit
+       var data: Array[Array[Byte]] is noautoinit
+       var crb: Array[Float] is noautoinit
+       var cib: Array[Float] is noautoinit
+       var atomic = new AtomicInt(0)
+       var nb_threads = 8
+end
+
+if args.is_empty then
+       n = 200
+else
+       n = args[0].to_i
+end
+
+sys.crb = new Array[Float].with_capacity(n + 7)
+sys.cib = new Array[Float].with_capacity(n + 7)
+sys.inv_n = 2.0 / n.to_f
+for i in [0..n[ do
+       sys.cib[i] = i.to_f * inv_n - 1.0
+       sys.crb[i] = i.to_f * inv_n - 1.5
+end
+sys.data = new Array[Array[Byte]].with_capacity(n)
+for i in [0..n[ do sys.data[i] = new Array[Byte].filled_with(0.to_b, (n + 7) / 8)
+
+# Parallel Approach
+var actors = new Array[Worker]
+for i in [0..nb_threads[ do
+       var a = new Worker
+       actors.add(a)
+       a.async.work
+end
+
+for a in actors do
+       a.async.terminate
+       a.async.wait_termination
+end
+
+var filename = "WRITE".environ
+if filename == "" then filename = "out"
+var output = new FileWriter.open(filename)
+output.write_bytes("P4\n{n} {n}\n".to_bytes)
+for i in [0..n[ do
+       var length = data[i].length
+       for j in [0..length[ do output.write_byte(data[i][j])
+end
+output.close

diff --git a/tests/sav/mandelbrot.res b/tests/sav/mandelbrot.res
new file mode 100644 (file)
index 0000000..2f7bbbc
Binary files /dev/null and b/tests/sav/mandelbrot.res differ