diff --git a/library/mpi/examples/demo.sh b/library/mpi/examples/demo.sh
new file mode 100755
index 000000000..fef079486
--- /dev/null
+++ b/library/mpi/examples/demo.sh
@@ -0,0 +1,53 @@
+#!/bin/sh
+# demo.sh
+#
+# This file was originally created on 3/7/2003
+# by Stasinos Konstantopoulos konstant@let.rug.nl
+# as part of the YAP Prolog distribution.
+#
+# This file is in the Public Domain.
+
+# Arguments: $1 will be the Prolog programme to load.
+
+PLFILE=$1
+shift 1
+
+YAP=${HOME}/opt/yap-cur/mpich-gm-gcc-nodebug/bin/yap
+
+# accumulator for stuff before the --
+before=""
+# accumulator for stuff after the --
+after=""
+# the previous parameter is expecting an argument
+argument=""
+# i'm still busy with the "before" stuff
+flag=yes
+
+while test -n "$1"; do
+    #echo "DEBUG: $*"
+    #echo "DEBUG: flag = $flag"
+    if test -n "$flag"; then
+	#echo "DEBUG: argument = $argument"
+	if test -n "$argument"; then
+	    argument=""
+	    before="$before $1"
+	#elif test "$1" = "-l" -o "$1" = "-h" -o "$1" = "-s" -o "$1" = "-t"; then
+	#    argument=yes
+	#    before="$before $1"
+	else
+	    flag=""
+	fi
+    fi
+    if test -z "$flag"; then
+	after="$after $1"
+    fi
+    shift
+done
+
+#echo "${YAP} $before -- $after"
+
+QUOTE="'"
+
+echo "consult(${QUOTE}${PLFILE}${QUOTE}). start(0,20). halt."
+echo "consult(${QUOTE}${PLFILE}${QUOTE}). start(0,20). halt." | ${YAP} $before -- $after
+
diff --git a/library/mpi/examples/demo1_mpe.pl b/library/mpi/examples/demo1_mpe.pl
new file mode 100644
index 000000000..b8cffd4af
--- /dev/null
+++ b/library/mpi/examples/demo1_mpe.pl
@@ -0,0 +1,140 @@
+%% demo1_mpe.pl
+%%
+%% This file was originally created on 3/7/2003
+%% by Stasinos Konstantopoulos <konstant@let.rug.nl>
+%% as part of the YAP Prolog distribution.
+%%
+%% This file is in the Public Domain.
+
+%% demo1_mpe.pl is the same as demo1.pl, except
+%% that MPE is used to log MPI activity.
+
+% make the `floor' operator return integer values
+:- set_prolog_flag(language, iso).
+
+
+%%
+%% This the calculation that needs to be performed, in this case
+%% the sum of [From..To]
+%%
+
+calc(From, From, Acc, Res) :- !,
+	Res is Acc + From.
+calc(From, To, Acc, Res) :- !,
+	Acc1 is Acc + To,
+	To1 is To - 1,
+	calc(From, To1, Acc1, Res).
+
+%%
+%% The master node sends teh task definition to
+%% the workers and then collects the results
+%%
+
+do(0, NumProc):-
+	!,
+
+	% processing state
+	mpe_create_event(Ev1),
+	mpe_create_event(Ev2),
+	mpe_create_state(Ev1,Ev2,processing,green),
+
+	% bcasting state
+	mpe_create_event(Ev3),
+	mpe_create_event(Ev4),
+	mpe_create_state(Ev3,Ev4,bcasting,red),
+
+	% sending/recving state
+	mpe_create_event(Ev5),
+	mpe_create_event(Ev6),
+	mpe_create_state(Ev5,Ev6,'sending/receiving',brown),
+
+	% pretend that the other nodes do not have
+	% access to the task definition.
+	% retrieve it and broadcast it.
+	get_value(from, From),
+	get_value(to, To),
+	mpe_log(Ev3,0,event3),
+	mpi_bcast(msg(From,To), 0),
+	mpe_log(Ev4,0,event4),
+
+	% loop to collect and accumulate partial results
+	% from the workers
+	set_value(n, NumProc),
+	set_value(acc, 0),
+	repeat,
+	  mpe_log(Ev5,0,event5),
+	  mpi_receive(T, Source, Tag),
+	  mpe_log(Ev6,0,event6),
+	  format( '0: Proc ~q said: ~q (Tag: ~q)~n', [Source,T,Tag] ),
+	  % accumulate results
+	  get_value(acc, Acc),
+	  NewAcc is Acc + T,
+	  set_value(acc, NewAcc),
+	  % processors still left
+	  get_value(n, Counter),
+	  NewCounter is Counter - 1,
+	  set_value(n, NewCounter),
+	  NewCounter == 1,
+	!,
+	format('0: Sum(~q..~q) = ~q.~n', [From,To,NewAcc]).
+
+
+%%
+%% The workers hear from the master what needs to
+%% be done, do the work and then send the results back.
+%%
+
+do(Rank, NumProc):-
+	!,
+
+	% processing state
+	mpe_create_event(Ev1),
+	mpe_create_event(Ev2),
+	mpe_create_state(Ev1,Ev2,processing,green),
+
+	% bcasting state
+	mpe_create_event(Ev3),
+	mpe_create_event(Ev4),
+	mpe_create_state(Ev3,Ev4,bcasting,red),
+
+	% sending/recving state
+	mpe_create_event(Ev5),
+	mpe_create_event(Ev6),
+	mpe_create_state(Ev5,Ev6,'sending/receiving',brown),
+
+	% catch the task broadcast
+	mpe_log(Ev3,0,event3),
+	mpi_bcast(Msg, 0),
+	mpe_log(Ev4,0,event4),
+	Msg = msg(From,To),
+	format( '~q: All are calculating ~q..~q.~n', [Rank,From,To] ),
+	MyFrom is floor(To * (Rank - 1) / (NumProc - 1)) + From,
+	MyTo is floor(To * Rank / (NumProc - 1)) + From - 1,
+	format( '~q: I am calculating ~q..~q.~n', [Rank,MyFrom,MyTo] ),
+	% do the job
+	mpe_log(Ev1,0,event1),
+	calc( MyFrom, MyTo, 0, Result ),
+	mpe_log(Ev2,0,event2),
+	format( '~q: sending ~q to 0. (Tag: 1)~n', [Rank,Result] ),
+	% send back the results
+	mpe_log(Ev5,0,event5),
+	mpi_send(Result, 0, 1),
+	mpe_log(Ev6,0,event6).
+
+
+%%
+%% This is the entry point
+%%
+
+start(From, To):-
+	% store the task definition
+	set_value(from, From),
+	set_value(to, To),
+
+	mpi_open(Rank, NumProc, ProcName),
+	format('Rank: ~q NumProc: ~q, ProcName: ~q~n', [Rank,NumProc,ProcName]),
+	mpe_open,
+	do(Rank, NumProc),
+	format( 'Rank ~q finished!~n', [Rank] ),
+	mpe_close( demo1_mpe ),
+	mpi_close.
diff --git a/library/mpi/examples/mk_gmiconf.awk b/library/mpi/examples/mk_gmiconf.awk
new file mode 100644
index 000000000..e5d8d7779
--- /dev/null
+++ b/library/mpi/examples/mk_gmiconf.awk
@@ -0,0 +1,19 @@
+#!/bin/awk
+# mk_gmiconf.awk
+#
+# This file was originally created on 3/7/2003
+# by Stasinos Konstantopoulos konstant@let.rug.nl
+# as part of the YAP Prolog distribution.
+#
+# This file is in the Public Domain.
+
+
+! /^$/ { 
+  if( $0 in BUFF ) {
+    BUFF[$0] = BUFF[$0] + 2;
+  }
+  else {
+    BUFF[$0] = 2;
+  }
+  print $0 " " BUFF[$0];
+}