diff --git a/bench/Makefile b/bench/Makefile
new file mode 100644
index 00000000..5525d75d
--- /dev/null
+++ b/bench/Makefile
@@ -0,0 +1,17 @@
+CC=cc -std=c99
+CXX=c++
+CFLAGS=-O2 -march=native -Wall
+CPPFLAGS=-I..
+
+all: bench-mlib bench-stl
+
+bench-mlib:
+	$(CC) $(CFLAGS) $(CPPFLAGS) bench-mlib.c -o bench-mlib.exe
+	./bench-mlib.exe
+
+bench-stl:
+	$(CXX) $(CFLAGS) $(CPPFLAGS) bench-stl.cpp -o bench-stl.exe
+	./bench-stl.exe
+
+clean:
+	rm -f *.o bench-mlib.exe bench-stl.exe *~ 
diff --git a/bench/bench-mlib.c b/bench/bench-mlib.c
new file mode 100644
index 00000000..e5c53f72
--- /dev/null
+++ b/bench/bench-mlib.c
@@ -0,0 +1,140 @@
+#include <stdlib.h>
+#include <stdio.h>
+#include <stdbool.h>
+#include <sys/types.h>
+#include <sys/resource.h>
+#include <sys/time.h>
+
+#define NDEBUG
+
+#include "m-array.h"
+#include "m-list.h"
+#include "m-rbtree.h"
+#include "m-dict.h"
+
+static unsigned long long
+cputime (void)
+{
+  struct rusage rus;
+  getrusage (0, &rus);
+  return rus.ru_utime.tv_sec * 1000000ULL + rus.ru_utime.tv_usec;
+}
+
+static unsigned int randValue = 0;
+static void rand_init(void)
+{
+  randValue = 0;
+}
+static unsigned int rand_get(void)
+{
+  randValue = randValue * 31421U + 6927U;
+  return (int) randValue;
+}
+
+static void test_function(const char str[], size_t n, void (*func)(size_t))
+{
+  unsigned long long start, end;
+  //  (*func)(n);
+  start = cputime();
+  (*func)(n);
+  end = cputime();
+  end = (end - start) / 1000U;
+  printf ("%s %Lu ms for n = %lu\n", str, end, n);
+}
+
+/********************************************************************************************/
+
+ARRAY_DEF(uint, unsigned int)
+
+static void test_array(size_t n)
+{
+  rand_init();
+  M_LET(a1, a2, ARRAY_OPLIST(uint)) {
+    for(size_t i = 0; i < n; i++) {
+      array_uint_push_back(a1, rand_get() );
+      array_uint_push_back(a2, rand_get() );
+    }
+    unsigned int s = 0;
+    for(unsigned long i = 0; i < n; i++) {
+      s += *array_uint_cget(a1, i ) * *array_uint_cget(a2, i );
+    }
+  }
+}
+
+/********************************************************************************************/
+
+LIST_DEF(uint, unsigned int)
+
+static void test_list (size_t n)
+{
+  rand_init();
+  M_LET(a1, a2, LIST_OPLIST(uint)) {
+    for(size_t i = 0; i < n; i++) {
+      list_uint_push_back(a1, rand_get() );
+      list_uint_push_back(a2, rand_get() );
+    }
+    unsigned int s = 0;
+    list_it_uint_t it1, it2;
+    for(list_uint_it(it1, a1), list_uint_it(it2,a2); !list_uint_end_p(it1); list_uint_next(it1), list_uint_next(it2)) {
+      s += *list_uint_cref(it1) * *list_uint_cref(it2);
+    }
+  }
+}
+
+/********************************************************************************************/
+
+RBTREE_DEF(ulong, unsigned long)
+
+static void test_rbtree(size_t n)
+{
+  rand_init();
+  M_LET(tree, RBTREE_OPLIST(ulong)) {
+    for (size_t i = 0; i < n; i++) {
+      rbtree_ulong_push(tree, rand_get());
+    }
+    
+    unsigned int s = 0;
+    for (size_t i = 0; i < n; i++) {
+      unsigned long *p = rbtree_ulong_get(tree, rand_get());
+      if (p)
+        s += *p;
+    }
+  }
+}
+
+/********************************************************************************************/
+
+unsigned long ulong_hash(unsigned long t) { return (unsigned long) t; }
+
+#define ULONG_OPLIST M_OPLIST_CAT((HASH(ulong_hash)), M_DEFAULT_OPLIST)
+
+DICT_DEF2(ulong, unsigned long, ULONG_OPLIST, unsigned long, ULONG_OPLIST)
+
+static void
+test_dict(unsigned long  n)
+{
+  rand_init();
+  M_LET(dict, DICT_OPLIST(ulong)) {
+    for (size_t i = 0; i < n; i++) {
+      dict_ulong_set_at(dict, rand_get(), rand_get() );
+    }
+    
+    unsigned int s = 0;
+    for (size_t i = 0; i < n; i++) {
+      unsigned long *p = dict_ulong_get(dict, rand_get());
+      if (p)
+        s += *p;
+    }
+  }
+}
+
+/********************************************************************************************/
+
+int main(int argc, const char *argv[])
+{
+  test_function("List   time",10000000, test_list);
+  test_function("Array  time", 100000000, test_array);
+  test_function("Rbtree time", 1000000, test_rbtree);
+  test_function("Dict   time", 1000000, test_dict);
+}
+
diff --git a/bench/bench-stl.cpp b/bench/bench-stl.cpp
new file mode 100644
index 00000000..91d344b3
--- /dev/null
+++ b/bench/bench-stl.cpp
@@ -0,0 +1,128 @@
+#include <vector>
+#include <list>
+#include <iterator>
+#include <map>
+#include <set>
+#include <unordered_map>
+
+#include <stdlib.h>
+#include <stdio.h>
+#include <stdbool.h>
+#include <sys/types.h>
+#include <sys/resource.h>
+#include <sys/time.h>
+
+#define NDEBUG
+
+using namespace std;
+
+static unsigned long long
+cputime (void)
+{
+  struct rusage rus;
+  getrusage (0, &rus);
+  return rus.ru_utime.tv_sec * 1000000ULL + rus.ru_utime.tv_usec;
+}
+
+static unsigned int randValue = 0;
+static void rand_init(void)
+{
+  randValue = 0;
+}
+static unsigned int rand_get(void)
+{
+  randValue = randValue * 31421U + 6927U;
+  return (int) randValue;
+}
+
+static void test_function(const char str[], size_t n, void (*func)(size_t))
+{
+  unsigned long long start, end;
+  //  (*func)(n);
+  start = cputime();
+  (*func)(n);
+  end = cputime();
+  end = (end - start) / 1000U;
+  printf ("%s %Lu ms for n = %lu\n", str, end, n);
+}
+
+/********************************************************************************************/
+
+static void test_list(size_t n)
+{
+  rand_init();
+  list<unsigned int> a1, a2;
+  for(size_t i = 0; i < n; i++) {
+    a1.push_back(rand_get() );
+    a2.push_back(rand_get() );
+  }
+  unsigned int s = 0;
+  for (list<unsigned int>::const_iterator ci1 = a1.begin(), ci2 = a2.begin (); ci1 != a1.end(); ++ci1, ++ci2) {
+    s += *ci1 * *ci2;
+  }
+}
+
+/********************************************************************************************/
+
+static void test_array(size_t n)
+{
+  rand_init();
+  vector<unsigned int> a1, a2;
+  for(size_t i = 0; i < n; i++) {
+    a1.push_back(rand_get() );
+    a2.push_back(rand_get() );
+  }
+  unsigned int s = 0;
+  for(unsigned long i = 0; i < n; i++) {
+    s += a1[i] * a2[i];
+  }
+}
+
+/********************************************************************************************/
+
+static void test_rbtree(size_t n)
+{
+  rand_init();
+  set<unsigned long> tree;
+
+  for (size_t i = 0; i < n; i++) {
+    tree.insert(rand_get());
+  }
+    
+  unsigned int s = 0;
+  for (size_t i = 0; i < n; i++) {
+    set<unsigned long>::iterator it = tree.find(rand_get());
+    if (it != tree.end())
+      s += *it;
+  }
+}
+
+/********************************************************************************************/
+
+static void
+test_dict(unsigned long  n)
+{
+  rand_init();
+  unordered_map<unsigned long, unsigned long> dict;
+
+  for (size_t i = 0; i < n; i++) {
+    dict[rand_get()] = rand_get();
+  }
+    
+  unsigned int s = 0;
+  for (size_t i = 0; i < n; i++) {
+    unordered_map<unsigned long, unsigned long>::iterator it = dict.find(rand_get());
+    if (it != dict.end())
+      s += it->second;
+  }
+}
+
+/********************************************************************************************/
+
+int main(int argc, const char *argv[])
+{
+  test_function("List   time",10000000, test_list);
+  test_function("Array  time", 100000000, test_array);
+  test_function("Rbtree time", 1000000, test_rbtree);
+  test_function("Dict   time", 1000000, test_dict);
+}