X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=blobdiff_plain;f=octave_packages%2Fopenmpi_ext-1.0.2%2Fmc_example.m;fp=octave_packages%2Fopenmpi_ext-1.0.2%2Fmc_example.m;h=ea4967ff72dfb7b7a28136248ed4a990713d98c9;hp=0000000000000000000000000000000000000000;hb=c880e8788dfc484bf23ce13fa2787f2c6bca4863;hpb=1705066eceaaea976f010f669ce8e972f3734b05

diff --git a/octave_packages/openmpi_ext-1.0.2/mc_example.m b/octave_packages/openmpi_ext-1.0.2/mc_example.m
new file mode 100644
index 0000000..ea4967f
--- /dev/null
+++ b/octave_packages/openmpi_ext-1.0.2/mc_example.m
@@ -0,0 +1,53 @@
+## Copyright (C) 2009  Michael Cree <michael.creel@uab.es>
+## under the terms of the GNU General Public License.
+##
+## This program is free software; you can redistribute it and/or modify
+## it under the terms of the GNU General Public License as published by
+## the Free Software Foundation; either version 2 of the License, or
+## (at your option) any later version.
+##
+## This program is distributed in the hope that it will be useful,
+## but WITHOUT ANY WARRANTY; without even the implied warranty of
+## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+## GNU General Public License for more details.
+##
+## You should have received a copy of the GNU General Public License
+## along with this program; If not, see <http://www.gnu.org/licenses/>.
+
+
+# mc_example: shows how Monte Carlo can be done using mpi, Does Monte
+# Carlo on the OLS estimator. Uses montecarlo.m
+#
+# USAGE: from the command prompt, not the octave prompt, execute
+# orterun -np 3 octave --eval mc_example
+1;
+function betahat = olswrapper(args)
+	n = args{1};
+	theta = args{2};
+	x = [ones(n,1) randn(n,1)];
+	y = x*theta + randn(n,1);
+	betahat = ols(y,x);
+  	betahat = betahat';
+endfunction
+
+
+n = 30;
+theta = [1;1];
+
+reps = 10000;
+f = "olswrapper";
+args = {n, theta};
+outfile = "mc_output";
+n_pooled = 100;
+verbose = false;
+tic;
+ montecarlo(f, args, reps, outfile, n_pooled, false, verbose);
+t=toc
+
+# info = MPI_Init 
+# tic;
+# montecarlo(f, args, reps, outfile, n_pooled, verbose);
+# disp("I am here!");
+# t=toc
+#  info = MPI_Finalize() 
+