Benchmark: bigarrays

Author: vouillon

benchmarks/benchmark-others/bigarrays/Makefile

@@ -0,0 +1,18 @@
+.PHONY: bench perform
+
+export NAME=Others/bigarrays
+
+SHELL=/bin/bash -o pipefail
+
+bench:
+	@date -u +"%FT%TZ - $(NAME): starting"
+	ocamlc bench.ml -o bench
+	$(MAKE) perform COMPILER=js_of_ocaml SCRIPT=bench.js KIND=js
+	$(MAKE) perform COMPILER=wasm_of_ocaml SCRIPT=bench.wasm.js KIND=wasm
+	@date -u +"%FT%TZ - $(NAME): done"
+
+perform:
+	$(COMPILER) --opt 2 --pretty bench -o $(SCRIPT)
+	/usr/bin/time -f '{"compiler": "$(COMPILER)", "time":"%E"}' node $(SCRIPT) 2>&1 | \
+	sh ../../utils/format_metrics.sh exec | \
+	sh ../../utils/aggregate.sh $(KIND)

benchmarks/benchmark-others/bigarrays/bench.ml

@@ -0,0 +1,36 @@
+let pi = 4. *. atan 1.
+
+let deltay = 40_000. /. 360. /. 3600. *. 1000.
+
+let deltax = deltay *. cos (44. *. pi /. 180.)
+
+let precompute tile_height tile_width tile =
+  let normals =
+    Bigarray.(Array3.create Int8_signed C_layout) (tile_height - 2) (tile_width - 2) 3
+  in
+  let heights =
+    Bigarray.(Array2.create Float32 C_layout) (tile_height - 2) (tile_width - 2)
+  in
+  for y = 1 to tile_height - 2 do
+    for x = 1 to tile_width - 2 do
+      let nx = (tile.{y, x - 1} -. tile.{y, x + 1}) *. deltay in
+      let ny = (tile.{y - 1, x} -. tile.{y + 1, x}) *. deltax in
+      let nz = 2. *. deltax *. deltay in
+      let n = 127. /. sqrt ((nx *. nx) +. (ny *. ny) +. (nz *. nz)) in
+      normals.{tile_height - 2 - y, x - 1, 0} <- truncate (nx *. n);
+      normals.{tile_height - 2 - y, x - 1, 1} <- truncate (ny *. n);
+      normals.{tile_height - 2 - y, x - 1, 2} <- truncate (nz *. n);
+      heights.{tile_height - 2 - y, x - 1} <- tile.{y, x}
+    done
+  done
+
+let tile_height = 1024
+
+let tile_width = 1024
+
+let tile = Bigarray.(Array2.create Float32 C_layout) tile_height tile_width
+
+let () =
+  for _ = 1 to 10 do
+    precompute tile_height tile_width tile
+  done