Menu
▾
▴
[Octave-cvsupdate] SF.net SVN: octave:[7052] trunk/octave-forge/main
|
From: <rik...@us...> - 2010-03-12 13:07:08
|
Revision: 7052
http://octave.svn.sourceforge.net/octave/?rev=7052&view=rev
Author: rikcorradini
Date: 2010-03-12 13:06:46 +0000 (Fri, 12 Mar 2010)
Log Message:
-----------
First working version for main folder
Added Paths:
-----------
trunk/octave-forge/main/openmpi_ext/
trunk/octave-forge/main/openmpi_ext/COPYING
trunk/octave-forge/main/openmpi_ext/DESCRIPTION
trunk/octave-forge/main/openmpi_ext/INDEX
trunk/octave-forge/main/openmpi_ext/doc/
trunk/octave-forge/main/openmpi_ext/doc/README
trunk/octave-forge/main/openmpi_ext/inst/
trunk/octave-forge/main/openmpi_ext/inst/Pi.m
trunk/octave-forge/main/openmpi_ext/inst/allnodes
trunk/octave-forge/main/openmpi_ext/inst/hello2dimmat.m
trunk/octave-forge/main/openmpi_ext/inst/hellocell.m
trunk/octave-forge/main/openmpi_ext/inst/hellosparsemat.m
trunk/octave-forge/main/openmpi_ext/inst/hellostruct.m
trunk/octave-forge/main/openmpi_ext/inst/helloworld.m
trunk/octave-forge/main/openmpi_ext/inst/mc_example.m
trunk/octave-forge/main/openmpi_ext/inst/montecarlo.m
trunk/octave-forge/main/openmpi_ext/src/
trunk/octave-forge/main/openmpi_ext/src/MPI_Barrier.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Load.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Test.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_rank.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_size.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Finalize.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Finalized.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Init.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Initialized.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Iprobe.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Op_Load.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Op_Test.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Probe.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Recv.cc
trunk/octave-forge/main/openmpi_ext/src/MPI_Send.cc
trunk/octave-forge/main/openmpi_ext/src/Makefile
trunk/octave-forge/main/openmpi_ext/src/Makefile.OPENMPI1.3.3
trunk/octave-forge/main/openmpi_ext/src/simple.h
trunk/octave-forge/main/openmpi_ext/src/simpleop.h
Added: trunk/octave-forge/main/openmpi_ext/COPYING
===================================================================
--- trunk/octave-forge/main/openmpi_ext/COPYING (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/COPYING 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,337 @@
+ GNU GENERAL PUBLIC LICENSE
+ Version 2, June 1991
+
+ Copyright (C) 1989, 1991 Free Software Foundation, Inc. <http://fsf.org/>
+ Everyone is permitted to copy and distribute verbatim copies
+ of this license document, but changing it is not allowed.
+
+ Preamble
+
+ The licenses for most software are designed to take away your
+freedom to share and change it. By contrast, the GNU General Public
+License is intended to guarantee your freedom to share and change free
+software--to make sure the software is free for all its users. This
+General Public License applies to most of the Free Software
+Foundation's software and to any other program whose authors commit to
+using it. (Some other Free Software Foundation software is covered by
+the GNU Library General Public License instead.) You can apply it to
+your programs, too.
+
+ When we speak of free software, we are referring to freedom, not
+price. Our General Public Licenses are designed to make sure that you
+have the freedom to distribute copies of free software (and charge for
+this service if you wish), that you receive source code or can get it
+if you want it, that you can change the software or use pieces of it
+in new free programs; and that you know you can do these things.
+
+ To protect your rights, we need to make restrictions that forbid
+anyone to deny you these rights or to ask you to surrender the rights.
+These restrictions translate to certain responsibilities for you if you
+distribute copies of the software, or if you modify it.
+
+ For example, if you distribute copies of such a program, whether
+gratis or for a fee, you must give the recipients all the rights that
+you have. You must make sure that they, too, receive or can get the
+source code. And you must show them these terms so they know their
+rights.
+
+ We protect your rights with two steps: (1) copyright the software, and
+(2) offer you this license which gives you legal permission to copy,
+distribute and/or modify the software.
+
+ Also, for each author's protection and ours, we want to make certain
+that everyone understands that there is no warranty for this free
+software. If the software is modified by someone else and passed on, we
+want its recipients to know that what they have is not the original, so
+that any problems introduced by others will not reflect on the original
+authors' reputations.
+
+ Finally, any free program is threatened constantly by software
+patents. We wish to avoid the danger that redistributors of a free
+program will individually obtain patent licenses, in effect making the
+program proprietary. To prevent this, we have made it clear that any
+patent must be licensed for everyone's free use or not licensed at all.
+
+ The precise terms and conditions for copying, distribution and
+modification follow.
+
+ GNU GENERAL PUBLIC LICENSE
+ TERMS AND CONDITIONS FOR COPYING, DISTRIBUTION AND MODIFICATION
+
+ 0. This License applies to any program or other work which contains
+a notice placed by the copyright holder saying it may be distributed
+under the terms of this General Public License. The "Program", below,
+refers to any such program or work, and a "work based on the Program"
+means either the Program or any derivative work under copyright law:
+that is to say, a work containing the Program or a portion of it,
+either verbatim or with modifications and/or translated into another
+language. (Hereinafter, translation is included without limitation in
+the term "modification".) Each licensee is addressed as "you".
+
+Activities other than copying, distribution and modification are not
+covered by this License; they are outside its scope. The act of
+running the Program is not restricted, and the output from the Program
+is covered only if its contents constitute a work based on the
+Program (independent of having been made by running the Program).
+Whether that is true depends on what the Program does.
+
+ 1. You may copy and distribute verbatim copies of the Program's
+source code as you receive it, in any medium, provided that you
+conspicuously and appropriately publish on each copy an appropriate
+copyright notice and disclaimer of warranty; keep intact all the
+notices that refer to this License and to the absence of any warranty;
+and give any other recipients of the Program a copy of this License
+along with the Program.
+
+You may charge a fee for the physical act of transferring a copy, and
+you may at your option offer warranty protection in exchange for a fee.
+
+ 2. You may modify your copy or copies of the Program or any portion
+of it, thus forming a work based on the Program, and copy and
+distribute such modifications or work under the terms of Section 1
+above, provided that you also meet all of these conditions:
+
+ a) You must cause the modified files to carry prominent notices
+ stating that you changed the files and the date of any change.
+
+ b) You must cause any work that you distribute or publish, that in
+ whole or in part contains or is derived from the Program or any
+ part thereof, to be licensed as a whole at no charge to all third
+ parties under the terms of this License.
+
+ c) If the modified program normally reads commands interactively
+ when run, you must cause it, when started running for such
+ interactive use in the most ordinary way, to print or display an
+ announcement including an appropriate copyright notice and a
+ notice that there is no warranty (or else, saying that you provide
+ a warranty) and that users may redistribute the program under
+ these conditions, and telling the user how to view a copy of this
+ License. (Exception: if the Program itself is interactive but
+ does not normally print such an announcement, your work based on
+ the Program is not required to print an announcement.)
+
+These requirements apply to the modified work as a whole. If
+identifiable sections of that work are not derived from the Program,
+and can be reasonably considered independent and separate works in
+themselves, then this License, and its terms, do not apply to those
+sections when you distribute them as separate works. But when you
+distribute the same sections as part of a whole which is a work based
+on the Program, the distribution of the whole must be on the terms of
+this License, whose permissions for other licensees extend to the
+entire whole, and thus to each and every part regardless of who wrote it.
+
+Thus, it is not the intent of this section to claim rights or contest
+your rights to work written entirely by you; rather, the intent is to
+exercise the right to control the distribution of derivative or
+collective works based on the Program.
+
+In addition, mere aggregation of another work not based on the Program
+with the Program (or with a work based on the Program) on a volume of
+a storage or distribution medium does not bring the other work under
+the scope of this License.
+
+ 3. You may copy and distribute the Program (or a work based on it,
+under Section 2) in object code or executable form under the terms of
+Sections 1 and 2 above provided that you also do one of the following:
+
+ a) Accompany it with the complete corresponding machine-readable
+ source code, which must be distributed under the terms of Sections
+ 1 and 2 above on a medium customarily used for software interchange; or,
+
+ b) Accompany it with a written offer, valid for at least three
+ years, to give any third party, for a charge no more than your
+ cost of physically performing source distribution, a complete
+ machine-readable copy of the corresponding source code, to be
+ distributed under the terms of Sections 1 and 2 above on a medium
+ customarily used for software interchange; or,
+
+ c) Accompany it with the information you received as to the offer
+ to distribute corresponding source code. (This alternative is
+ allowed only for noncommercial distribution and only if you
+ received the program in object code or executable form with such
+ an offer, in accord with Subsection b above.)
+
+The source code for a work means the preferred form of the work for
+making modifications to it. For an executable work, complete source
+code means all the source code for all modules it contains, plus any
+associated interface definition files, plus the scripts used to
+control compilation and installation of the executable. However, as a
+special exception, the source code distributed need not include
+anything that is normally distributed (in either source or binary
+form) with the major components (compiler, kernel, and so on) of the
+operating system on which the executable runs, unless that component
+itself accompanies the executable.
+
+If distribution of executable or object code is made by offering
+access to copy from a designated place, then offering equivalent
+access to copy the source code from the same place counts as
+distribution of the source code, even though third parties are not
+compelled to copy the source along with the object code.
+
+ 4. You may not copy, modify, sublicense, or distribute the Program
+except as expressly provided under this License. Any attempt
+otherwise to copy, modify, sublicense or distribute the Program is
+void, and will automatically terminate your rights under this License.
+However, parties who have received copies, or rights, from you under
+this License will not have their licenses terminated so long as such
+parties remain in full compliance.
+
+ 5. You are not required to accept this License, since you have not
+signed it. However, nothing else grants you permission to modify or
+distribute the Program or its derivative works. These actions are
+prohibited by law if you do not accept this License. Therefore, by
+modifying or distributing the Program (or any work based on the
+Program), you indicate your acceptance of this License to do so, and
+all its terms and conditions for copying, distributing or modifying
+the Program or works based on it.
+
+ 6. Each time you redistribute the Program (or any work based on the
+Program), the recipient automatically receives a license from the
+original licensor to copy, distribute or modify the Program subject to
+these terms and conditions. You may not impose any further
+restrictions on the recipients' exercise of the rights granted herein.
+You are not responsible for enforcing compliance by third parties to
+this License.
+
+ 7. If, as a consequence of a court judgment or allegation of patent
+infringement or for any other reason (not limited to patent issues),
+conditions are imposed on you (whether by court order, agreement or
+otherwise) that contradict the conditions of this License, they do not
+excuse you from the conditions of this License. If you cannot
+distribute so as to satisfy simultaneously your obligations under this
+License and any other pertinent obligations, then as a consequence you
+may not distribute the Program at all. For example, if a patent
+license would not permit royalty-free redistribution of the Program by
+all those who receive copies directly or indirectly through you, then
+the only way you could satisfy both it and this License would be to
+refrain entirely from distribution of the Program.
+
+If any portion of this section is held invalid or unenforceable under
+any particular circumstance, the balance of the section is intended to
+apply and the section as a whole is intended to apply in other
+circumstances.
+
+It is not the purpose of this section to induce you to infringe any
+patents or other property right claims or to contest validity of any
+such claims; this section has the sole purpose of protecting the
+integrity of the free software distribution system, which is
+implemented by public license practices. Many people have made
+generous contributions to the wide range of software distributed
+through that system in reliance on consistent application of that
+system; it is up to the author/donor to decide if he or she is willing
+to distribute software through any other system and a licensee cannot
+impose that choice.
+
+This section is intended to make thoroughly clear what is believed to
+be a consequence of the rest of this License.
+
+ 8. If the distribution and/or use of the Program is restricted in
+certain countries either by patents or by copyrighted interfaces, the
+original copyright holder who places the Program under this License
+may add an explicit geographical distribution limitation excluding
+those countries, so that distribution is permitted only in or among
+countries not thus excluded. In such case, this License incorporates
+the limitation as if written in the body of this License.
+
+ 9. The Free Software Foundation may publish revised and/or new versions
+of the General Public License from time to time. Such new versions will
+be similar in spirit to the present version, but may differ in detail to
+address new problems or concerns.
+
+Each version is given a distinguishing version number. If the Program
+specifies a version number of this License which applies to it and "any
+later version", you have the option of following the terms and conditions
+either of that version or of any later version published by the Free
+Software Foundation. If the Program does not specify a version number of
+this License, you may choose any version ever published by the Free Software
+Foundation.
+
+ 10. If you wish to incorporate parts of the Program into other free
+programs whose distribution conditions are different, write to the author
+to ask for permission. For software which is copyrighted by the Free
+Software Foundation, write to the Free Software Foundation; we sometimes
+make exceptions for this. Our decision will be guided by the two goals
+of preserving the free status of all derivatives of our free software and
+of promoting the sharing and reuse of software generally.
+
+ NO WARRANTY
+
+ 11. BECAUSE THE PROGRAM IS LICENSED FREE OF CHARGE, THERE IS NO WARRANTY
+FOR THE PROGRAM, TO THE EXTENT PERMITTED BY APPLICABLE LAW. EXCEPT WHEN
+OTHERWISE STATED IN WRITING THE COPYRIGHT HOLDERS AND/OR OTHER PARTIES
+PROVIDE THE PROGRAM "AS IS" WITHOUT WARRANTY OF ANY KIND, EITHER EXPRESSED
+OR IMPLIED, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF
+MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE ENTIRE RISK AS
+TO THE QUALITY AND PERFORMANCE OF THE PROGRAM IS WITH YOU. SHOULD THE
+PROGRAM PROVE DEFECTIVE, YOU ASSUME THE COST OF ALL NECESSARY SERVICING,
+REPAIR OR CORRECTION.
+
+ 12. IN NO EVENT UNLESS REQUIRED BY APPLICABLE LAW OR AGREED TO IN WRITING
+WILL ANY COPYRIGHT HOLDER, OR ANY OTHER PARTY WHO MAY MODIFY AND/OR
+REDISTRIBUTE THE PROGRAM AS PERMITTED ABOVE, BE LIABLE TO YOU FOR DAMAGES,
+INCLUDING ANY GENERAL, SPECIAL, INCIDENTAL OR CONSEQUENTIAL DAMAGES ARISING
+OUT OF THE USE OR INABILITY TO USE THE PROGRAM (INCLUDING BUT NOT LIMITED
+TO LOSS OF DATA OR DATA BEING RENDERED INACCURATE OR LOSSES SUSTAINED BY
+YOU OR THIRD PARTIES OR A FAILURE OF THE PROGRAM TO OPERATE WITH ANY OTHER
+PROGRAMS), EVEN IF SUCH HOLDER OR OTHER PARTY HAS BEEN ADVISED OF THE
+POSSIBILITY OF SUCH DAMAGES.
+
+ END OF TERMS AND CONDITIONS
+
+ How to Apply These Terms to Your New Programs
+
+ If you develop a new program, and you want it to be of the greatest
+possible use to the public, the best way to achieve this is to make it
+free software which everyone can redistribute and change under these terms.
+
+ To do so, attach the following notices to the program. It is safest
+to attach them to the start of each source file to most effectively
+convey the exclusion of warranty; and each file should have at least
+the "copyright" line and a pointer to where the full notice is found.
+
+ <one line to give the program's name and a brief idea of what it does.>
+ Copyright (C) <year> <name of author>
+
+ 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/>.
+
+Also add information on how to contact you by electronic and paper mail.
+
+If the program is interactive, make it output a short notice like this
+when it starts in an interactive mode:
+
+ Gnomovision version 69, Copyright (C) year name of author
+ Gnomovision comes with ABSOLUTELY NO WARRANTY; for details type `show w'.
+ This is free software, and you are welcome to redistribute it
+ under certain conditions; type `show c' for details.
+
+The hypothetical commands `show w' and `show c' should show the appropriate
+parts of the General Public License. Of course, the commands you use may
+be called something other than `show w' and `show c'; they could even be
+mouse-clicks or menu items--whatever suits your program.
+
+You should also get your employer (if you work as a programmer) or your
+school, if any, to sign a "copyright disclaimer" for the program, if
+necessary. Here is a sample; alter the names:
+
+ Yoyodyne, Inc., hereby disclaims all copyright interest in the program
+ `Gnomovision' (which makes passes at compilers) written by James Hacker.
+
+ <signature of Ty Coon>, 1 April 1989
+ Ty Coon, President of Vice
+
+This General Public License does not permit incorporating your program into
+proprietary programs. If your program is a subroutine library, you may
+consider it more useful to permit linking proprietary applications with the
+library. If this is what you want to do, use the GNU Library General
+Public License instead of this License.
Added: trunk/octave-forge/main/openmpi_ext/DESCRIPTION
===================================================================
--- trunk/octave-forge/main/openmpi_ext/DESCRIPTION (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/DESCRIPTION 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,11 @@
+Name: openmpi_ext
+Version: 1.0.0
+Date: 2009-11-16
+Author: Riccardo Corradini <ric...@ya...>
+Maintainer: Riccardo Corradini <ric...@ya...>
+Title: openmpi_ext
+Description: MPI functions for parallel computing using simple MPI Derived Datatypes.
+Depends: octave (>= 3.2.3)
+Autoload: yes
+License: GPL version 2 or later
+Url: http://octave.sf.net
Added: trunk/octave-forge/main/openmpi_ext/INDEX
===================================================================
--- trunk/octave-forge/main/openmpi_ext/INDEX (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/INDEX 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1 @@
+TO be filled
Added: trunk/octave-forge/main/openmpi_ext/doc/README
===================================================================
--- trunk/octave-forge/main/openmpi_ext/doc/README (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/doc/README 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,44 @@
+The code is general-purpose, but I would like to use it for econometrics.
+So the first step will be to install the following tarball from
+http://www.open-mpi.org/software/ompi/v1.3/downloads/openmpi-1.3.3.tar.bz2
+possibly in a multi-core computer to run my simple examples
+and configure it for instance in the following way (/home/user is your $HOME)
+./configure --enable-mpirun-prefix-by-default --enable-heterogeneous --prefix=/home/user/openmpi-1.3.3/ --enable-static
+
+and modify .bashrc in your home
+ OMPIBIN=`$ompi_info -path bindir -parsable | cut -d: -f3`
+ OMPILIB=`$ompi_info -path libdir -parsable | cut -d: -f3`
+ OMPISCD=`$ompi_info -path sysconfdir -parsable | cut -d: -f3`
+
+
+
+export PATH=$OMPIBIN:$PATH
+
+export LD_LIBRARY_PATH=:$OMPILIB:$LD_LIBRARY_PATH
+
+unset ompi_info OMPIBIN OMPILIB OMPISCD
+
+If you want to install it on a simple toy network, just assign a static ip address on every linux computer and set up
+an ssh connection with no password (see for instance http://linuxproblem.org/art_9.html ) and then install openmpi and octave always with the same versions and with the same info on .bashrc for the same user.
+
+After this type in a terminal mpiCC --showme
+In my case I will have something like
+
+g++ -I/home/user/openmpi-1.3.3/include -pthread -L/home/user/openmpi-1.3.3/lib -lmpi_cxx -lmpi -lopen-rte -lopen-pal -ldl -Wl,--export-dynamic -lnsl -lutil -lm -ldl
+
+This will be useful for mkoctfile
+for instance for MPI_Init.cc we shall have
+mkoctfile -I/home/user/openmpi-1.3.3/include -lpthread -L/home/user/openmpi-1.3.3/lib -lmpi_cxx -lmpi -lopen-rte -lopen-pal -ldl -lnsl -lutil -lm -ldl MPI_Init.cc
+
+
+
+The m files just contain some very simple examples
+More complex examples will be provided in the next future.
+See also
+http://static.msi.umn.edu/tutorial/scicomp/general/MPI/content6.html
+to understand the logic of MPI Derived Datatypes and how could they be easily handled by openmpi_ext package.
+Bests regards
+Riccardo Corradini
+
+
+
Added: trunk/octave-forge/main/openmpi_ext/inst/Pi.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/Pi.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/Pi.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,113 @@
+## Copyright (C) 2004-2007 Javier Fernández Baldomero, Mancia Anguita López
+## This code has been adjusted for octave3.2.3 and octave 3.3.50+ in
+## 2009 by Riccardo Corradini <ric...@ya...>
+##
+## 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/>.
+
+
+# Please add the oct files openmpi_ext folder
+# For instance addpath("../src");
+# mpirun -np 5 octave -q --eval "Pi(2E7,'s')"
+
+function Pi(N,mod)
+addpath("../src");
+# Pi: Classic PI computation by numeric integration of arctan'(x) in [0..1]
+#
+# Pi [ ( N [ ,mod ] ) ]
+#
+# N [1E7] #subdivisions of the [0, 1] interval
+# mod ['s'] communication modality: (s)end (r)educe
+#
+# printed results struct contains
+# pi estimated pi value
+# err error
+# time from argument xmit to pi computed
+#
+
+
+##########
+# ArgChk #
+##########
+if nargin<1, N=1E7; end
+if nargin<2, mod='s'; end
+if nargin>2, usage("Pi(N,mod)"); end # let all ranks complain
+flag=0; # code much simpler
+flag=flag || ~isscalar(N) || ~isnumeric(N);
+flag=flag | fix(N)~=N | N<1;
+ mod=lower(mod); mods='sr';
+flag=flag | isempty(findstr(mod, mods)); # let them all error out
+if flag, usage("Pi( <int> N>0, <char> mod=='s|r' )"); end
+
+##################
+# Results struct #
+##################
+results.pi =0;
+results.err =0;
+results.time =0;
+
+
+############
+# PARALLEL # initialization, include MPI_Init time in measurement
+############
+ T=clock; #
+############
+ MPI_ANY_SOURCE = -1;
+ MPI_Init();
+ MPI_COMM_WORLD = MPI_Comm_Load("NEWORLD");
+ rnk = MPI_Comm_rank (MPI_COMM_WORLD); # let it abort if it fails
+ siz = MPI_Comm_size (MPI_COMM_WORLD);
+
+ SLV = logical(rnk); # handy shortcuts, master is rank 0
+ MST = ~ SLV; # slaves are all other
+
+############
+# PARALLEL # computation (depends on rank/size)
+############ # vectorized code, equivalent to
+ width=1/N; lsum=0; # for i=rnk:siz:N-1
+ i=rnk:siz:N-1; # x=(i+0.5)*width;
+ x=(i+0.5)*width; # lsum=lsum+4/(1+x^2);
+ lsum=sum(4./(1+x.^2)); # end
+
+############
+# PARALLEL # reduction and finish
+############
+switch mod
+ case 's', TAG=7; # Any tag would do
+ if SLV # All slaves send result back
+ MPI_Send(lsum, 0,TAG,MPI_COMM_WORLD);
+ else # Here at master
+ Sum =lsum; # save local result
+ for slv=1:siz-1 # collect in any order
+ lsum = MPI_Recv(MPI_ANY_SOURCE,TAG,MPI_COMM_WORLD);
+ Sum+=lsum; # and accumulate
+ end # order: slv or MPI_ANY_SOURCE
+ end
+ case 'r',
+ disp("not yet implemented");
+# Sum=0;
+# reduction master = rank 0 @ WORLD
+# MPI_Reduce(lsum,Sum, MPI_SUM, 0,MPI_COMM_WORLD);
+end
+
+MPI_Finalize();
+
+if MST
+ Sum = Sum/N ; # better at end: don't loose resolution
+################################# # stopwatch measurement
+results.time = etime(clock,T); # # but only at master after PI computed
+################################# # all them started T=clock;
+results.err = Sum-pi;
+results.pi = Sum # ;
+
+end
\ No newline at end of file
Added: trunk/octave-forge/main/openmpi_ext/inst/allnodes
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/allnodes (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/allnodes 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,7 @@
+#!/bin/bash
+# this script will run an Octave script using MPI with a given hostfile,
+# using a given number of ranks. Edit it to set the hostfile and number of nodes
+# then use it as follows: allnodes <your_script_name.m>
+HOSTFILE="/home/user/tmp/bhosts"
+NUMBER_OF_NODES="33"
+mpirun --hostfile $HOSTFILE -np $NUMBER_OF_NODES octave -q --eval $1
Added: trunk/octave-forge/main/openmpi_ext/inst/hello2dimmat.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/hello2dimmat.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/hello2dimmat.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,36 @@
+# Please add the oct files openmpi_ext folder
+# For instance
+ addpath("../src");
+ MPI_SUCCESS =0;
+ MPI_Init();
+
+ # the string NEWORLD is just a label could be whatever you want
+ CW = MPI_Comm_Load("NEWORLD");
+ my_rank = MPI_Comm_rank(CW);
+ p = MPI_Comm_size(CW);
+ mytag = 48;
+
+
+
+ if (my_rank != 0)
+# Generate a random matrix
+ message=rand(90,90);
+# load message
+# rankvect is the vector containing the list of rank destination process
+ rankvect = 0;
+ [info] = MPI_Send(message,rankvect,mytag,CW);
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [messager, info] = MPI_Recv(source,mytag,CW);
+
+# You could also save each result and make comparisons if you don't trust MPI
+ disp("Rank 0 is the master receiving ... :");
+ if (info == MPI_SUCCESS)
+ disp('OK!');
+ endif
+ endfor
+ end
+
+
+ MPI_Finalize();
Added: trunk/octave-forge/main/openmpi_ext/inst/hellocell.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/hellocell.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/hellocell.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,55 @@
+## Copyright (C) 2009 Riccardo Corradini <ric...@ya...>
+## 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/>.
+
+
+# Please add the oct files openmpi_ext folder
+# For instance
+addpath("../src");
+# if you have 4 cores or a network of 4 computers with a ssh connection with no password and same openmpi 1.3.3 installation
+# type at the terminal mpirun -np 4 octave --eval hellocell
+
+
+ MPI_Init();
+ # the string NEWORLD is just a label could be whatever you want
+ CW = MPI_Comm_Load("NEWORLD");
+
+
+
+ my_rank = MPI_Comm_rank(CW);
+ p = MPI_Comm_size(CW);
+ # TAG is very important to identify the message
+ TAG = 1;
+
+
+ message="";
+ if (my_rank != 0)
+ message = {magic(3) 17 'fred'; ...
+ 'AliceBettyCarolDianeEllen' 'yp' 42; ...
+ {1} 2 3};
+ rankvect = 0;
+ [info] = MPI_Send(message,rankvect,TAG,CW);
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [messager, info] = MPI_Recv(source,TAG,CW);
+ info
+ messager
+ endfor
+ end
+
+
+ MPI_Finalize();
+
Added: trunk/octave-forge/main/openmpi_ext/inst/hellosparsemat.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/hellosparsemat.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/hellosparsemat.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,50 @@
+# Please add the oct files openmpi_ext folder
+# For instance
+ addpath("../src");
+ MPI_Init();
+ # the string NEWORLD is just a label could be whatever you want
+ CW = MPI_Comm_Load("NEWORLD");
+ my_rank = MPI_Comm_rank(CW);
+ p = MPI_Comm_size(CW);
+# tag[0] ----> type of octave_value
+# tag[1] ----> array of three elements 1) num of rows 2) number of columns 3) number of non zero elements
+# tag[2] ----> vector of rowindex
+# tag[3] ----> vector of columnindex
+# tag[4] ----> vector of non zero elements
+# These tags will be generated after mytag by the MPI_Send and MPI_Recv (see source code)
+
+ mytag = 48;
+
+
+
+
+# This is just to fill the sparse matrix
+ M=5;
+ N=5;
+ D=0.9;
+ message = sprand (M, N, D);
+# load message
+
+
+
+ if (my_rank != 0)
+ dest = 0;
+# rankvect is the vector containing the list of rank destination process
+ rankvect(1,1) = 0;
+ [info] = MPI_Send(message,rankvect,mytag,CW);
+ disp("This is flag for sending the message --")
+ info
+ else
+ for source = 1:p-1
+ messager='';
+ disp("We are at rank 0 that is master etc..");
+ [messager, info] = MPI_Recv(source,mytag,CW);
+ disp("Rank 0 is the master receiving ... :");
+ if (messager/message)
+ disp('OK!');
+ endif
+ messager
+ endfor
+ end
+
+ MPI_Finalize();
Added: trunk/octave-forge/main/openmpi_ext/inst/hellostruct.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/hellostruct.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/hellostruct.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,51 @@
+## Copyright (C) 2009 Riccardo Corradini <ric...@ya...>
+## 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/>.
+
+
+# Please add the oct files openmpi_ext folder
+# For instance
+addpath("../src");
+# if you have 4 cores or a network of 4 computers with a ssh connection with no password and same openmpi 1.3.3 installation
+# type at the terminal mpirun -np 4 octave --eval hellostruct
+
+
+ MPI_Init();
+ # the string NEWORLD is just a label could be whatever you want
+ CW = MPI_Comm_Load("NEWORLD");
+
+
+
+ my_rank = MPI_Comm_rank(CW);
+ p = MPI_Comm_size(CW);
+ # TAG is very important to identify the message
+ TAG = 1;
+
+
+ message="";
+ if (my_rank != 0)
+ message = struct('f1', {1 3; 2 4}, 'f2', 25);
+ # Could be a vector containing the list of ranks identifiers;
+ rankvect = 0;
+ [info] = MPI_Send(message,rankvect,TAG,CW);
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [message, info] = MPI_Recv(source,TAG,CW);
+ message
+ endfor
+ end
+ MPI_Finalize();
+
Added: trunk/octave-forge/main/openmpi_ext/inst/helloworld.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/helloworld.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/helloworld.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,52 @@
+## Copyright (C) 2009 Riccardo Corradini <ric...@ya...>
+## 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/>.
+
+
+# Please add the oct files openmpi_ext folder
+# For instance
+addpath("../src");
+# if you have 4 cores or a network of 4 computers with a ssh connection with no password and same openmpi 1.3.3 installation
+# type at the terminal mpirun -np 4 octave --eval helloworld
+
+
+ MPI_Init();
+ # the string NEWORLD is just a label could be whatever you want
+ CW = MPI_Comm_Load("NEWORLD");
+
+
+
+ my_rank = MPI_Comm_rank(CW);
+ p = MPI_Comm_size(CW);
+ # Could be any number
+ TAG=1;
+
+
+ message="";
+ if (my_rank != 0)
+ message = sprintf('Greetings from process: %d!',my_rank);
+ # rankvect is the vector containing the list of rank destination process
+ rankvect = 0;
+ [info] = MPI_Send(message,rankvect,TAG,CW);
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [message, info] = MPI_Recv(source,TAG,CW);
+ printf('%s\n', message);
+ endfor
+ end
+
+ MPI_Finalize();
+
Added: trunk/octave-forge/main/openmpi_ext/inst/mc_example.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/mc_example.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/mc_example.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,50 @@
+## Copyright (C) 2009 Michael Cree <mic...@ua...>
+## 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;
+addpath("../src");
+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 = 1000;
+f = "olswrapper";
+args = {n, theta};
+outfile = "mc_output";
+n_pooled = 10;
+verbose = true;
+
+# montecarlo(f, args, reps, outfile, n_pooled, false, verbose);
+
+if not(MPI_Initialized) MPI_Init; endif
+montecarlo(f, args, reps, outfile, n_pooled, verbose);
+if not(MPI_Finalized) MPI_Finalize; endif
Added: trunk/octave-forge/main/openmpi_ext/inst/montecarlo.m
===================================================================
--- trunk/octave-forge/main/openmpi_ext/inst/montecarlo.m (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/inst/montecarlo.m 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,135 @@
+# Copyright (C) 2006, 2009 Michael Creel <mic...@ua...>
+# 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, write to the Free Software
+# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
+
+# montecarlo.m: generates a specified number of replications of a function's
+# output and writes them to a user-specified output file.
+#
+# USAGE: montecarlo(f,f_args,reps,outfile,n_pooled,n_returns,usempi, verbose)
+#
+# IMPORTANT: f should return a row vector of output from feval(f,f_args)
+#
+# For normal evaluation on one core, only the first 4 arguments are required.
+# * Arg 1: (required) the function that generates a row vector of output
+# * Arg 2: (required) the arguments of the function, in a cell
+# * Arg 3: (required) the number of replications to generate
+# * Arg 4: (required) the output file name
+# * Arg 5 (optional) number of replications to be pooled together between writes
+# * Arg 6 (optional) verbose: 1 for on, 0 for off
+#
+# If using MPI, you should run using ranks equal to number of cores plus 1,
+# and should make sure that the core running the frontend is also the one that
+# has the second rank. That way the core the frontend is on will also do work.
+
+function n_received = montecarlo(f,f_args,reps,outfile,n_pooled,verbose)
+
+ t0 = clock(); # initialize timing
+
+ # defaults for optional arguments
+ if (nargin < 6) verbose = false; endif
+ if (nargin < 5) n_pooled = 1; endif;
+
+ if MPI_Initialized # check if doing this parallel or serial
+ use_mpi = true;
+ CW = MPI_Comm_Load("NEWORLD");
+ is_node = MPI_Comm_rank(CW);
+ nodes = MPI_Comm_size(CW);
+ mytag = 48;
+ else
+ use_mpi = false;
+ is_node = 0;
+ endif
+
+ if is_node # compute nodes
+ more_please = 1;
+ while more_please
+ for i = 1:n_pooled
+ contrib = feval(f, f_args);
+ contribs(i,:) = contrib;
+ endfor
+ MPI_Send(contribs, 0, mytag, CW);
+ # check if we're done
+ if (MPI_Iprobe(0, is_node, CW)) # check for ping from rank 0
+ junk = MPI_Recv(0, is_node, CW);
+ break;
+ endif
+ endwhile
+ else # frontend
+ received = 0;
+ done = false;
+ while received < reps
+ if use_mpi
+ # retrieve results from compute nodes
+ for i = 1:nodes-1
+ # compute nodes have results yet?
+ ready = false;
+ ready = MPI_Iprobe(i, mytag, CW); # check if message pending
+ if ready
+ # get it if it's there
+ contribs = MPI_Recv(i, mytag, CW);
+ need = reps - received;
+ received = received + n_pooled;
+ # truncate?
+ if n_pooled >= need
+ contribs = contribs(1:need,:);
+ done = true;
+ endif
+ # write to output file
+ FN = fopen (outfile, "a");
+ if (FN < 0) error ("montecarlo: couldn't open output file %s", outfile); endif
+ t = etime(clock(), t0);
+ for j = 1:rows(contribs)
+ fprintf(FN, "%f ", i, t, contribs(j,:));
+ fprintf(FN, "\n");
+ endfor
+ fclose(FN);
+ if verbose printf("\nContribution received from node%d. Received so far: %d\n", i, received); endif
+ endif
+ # tell compute nodes to stop loop
+ if done
+ for i = 1:(nodes-1)
+ ready = MPI_Iprobe(i, mytag, CW); # get last messages
+ if ready contribs = MPI_Recv(i, mytag, CW); endif
+ MPI_Send(" ",i,i,CW); # send out message to stop
+ endfor
+ break;
+ endif
+ endfor
+ else
+ for i = 1:n_pooled
+ contrib = feval(f, f_args);
+ contribs(i,:) = contrib;
+ endfor
+ need = reps - received;
+ received = received + n_pooled;
+ # truncate?
+ if n_pooled >= need
+ contribs = contribs(1:need,:);
+ endif
+ # write to output file
+ FN = fopen (outfile, "a");
+ if (FN < 0) error ("montecarlo: couldn't open output file %s", outfile); endif
+ t = etime(clock(), t0);
+ for j = 1:rows(contribs)
+ fprintf(FN, "%f ", 0, t, contribs(j,:));
+ fprintf(FN, "\n");
+ endfor
+ fclose(FN);
+ if verbose printf("\nContribution received from node 0. Received so far: %d\n", received); endif
+ endif
+ endwhile
+ endif
+endfunction
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Barrier.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Barrier.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Barrier.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,85 @@
+// Copyright (C) 2004-2007 Javier Fernández Baldomero, Mancia Anguita López
+// This code has been adjusted for octave3.2.3 and more in
+// 2009 by Riccardo Corradini <ric...@ya...>
+// Copyright (C) 2009 VZLU Prague
+
+//
+// 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/>.
+
+#define NAME MPI_Barrier
+/*
+ * ----------------------------------------------------
+ * Blocks until all processes in the communicator have reached this routine
+ * [info ] = MPI_Barrier (comm)
+ * ----------------------------------------------------
+ */
+
+
+#include "simple.h"
+DEFUN_DLD(NAME, args, ,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} @var{INFO} = MPI_Barrier (@var{COMM})\n\
+Blocks until all processes in the communicator have reached this routine.\n\
+If @var{COMM} octave comunicator object loaded with MPI_Comm_Load is omitted \n\
+returns an error. \n\
+ @example\n\
+ @group\n\
+ @var{INFO} (int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 5 MPI_ERR_COMM Invalid communicator (NULL?)\n\
+ 13 MPI_ERR_ARG Invalid argument (typically a NULL pointer?)\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+{
+
+ octave_value results;
+ int nargin = args.length ();
+ if (nargin != 1)
+ {
+ error ("expecting 1 input argument");
+ return results;
+ }
+
+ if (!simple_type_loaded)
+ {
+ simple::register_type ();
+ simple_type_loaded = true;
+ mlock ();
+ }
+
+ if((args.length() != 1 )
+ || args(0).type_id()!=simple::static_type_id()){
+
+ error("Please enter octave comunicator object!");
+ return octave_value(-1);
+ }
+
+ const octave_base_value& rep = args(0).get_rep();
+ const simple& B = ((const simple &)rep);
+ MPI_Comm comm = ((const simple&) B).comunicator_value ();
+ if (! error_state)
+ {
+ int my_size;
+ int info = MPI_Barrier (comm);
+
+ results = info;
+ }
+ else
+ print_usage ();
+ comm= NULL;
+ /* [info] = MPI_Barrier (comm) */
+
+ return results;
+}
+
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Load.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Load.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Load.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,56 @@
+// Copyright (C) 2009 Riccardo Corradini <ric...@ya...>
+// under the terms of the GNU General Public License.
+// Copyright (C) 2009 VZLU Prague
+//
+// 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/>.
+#include "simple.h"
+
+DEFUN_DLD(MPI_Comm_Load, args, ,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} @var{COMM} = MPI_Comm_Load (@var{DESCRIPTION})\n\
+Return @var{COMM} the MPI_Communicator object whose description is @var{DESCRIPTION}, as a string.\n\
+The default value will be MPI_COMM_WORLD. \n\
+If @var{DESCRIPTION} is omitted, return anyway an MPI_COMM_WORLD comunicator object \n\
+with no decription.\n\
+For\n\
+example,\n\
+\n\
+@example\n\
+@group\n\
+MPI_Init();\n\
+X = MPI_Comm_Load(\"description\"); \n\
+whos X\n\
+MPI_Finalize();\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+{
+ if (!simple_type_loaded)
+ {
+ simple::register_type ();
+ simple_type_loaded = true;
+ mlock ();
+ }
+
+ octave_value retval;
+ if (args.length () != 1 || !args (0).is_string ())
+ {
+ error ("MPI_Comm_Load: first argument must be a string");
+ return retval;
+ }
+
+ const std::string name = args (0).string_value ();
+ retval = new simple (name,MPI_COMM_WORLD);
+
+ return retval;
+}
\ No newline at end of file
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Test.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Test.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_Test.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,40 @@
+#include "simple.h"
+DEFUN_DLD(MPI_Comm_Test, args, ,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} @var{DESCRIPTION} = MPI_Comm_Test (@var{COMM})\n\
+Return @var{DESCRIPTION} string description of the MPI_Communicator @var{COMM}\n\
+For\n\
+example,\n\
+\n\
+@example\n\
+@group\n\
+MPI_Init();\n\
+X = MPI_Comm_Load(\"description\"); \n\
+whos X\n\
+MPI_Comm_Test(X) \n\
+@result{} \"description\"\n\
+MPI_Finalize();\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+{
+ if (!simple_type_loaded)
+ {
+ simple::register_type ();
+ simple_type_loaded = true;
+ mlock ();
+ }
+
+ octave_value retval;
+ if(args.length() != 1
+ || args(0).type_id()!=simple::static_type_id()){
+
+ error("usage: MPI_Comm_Test(octave_comunicator_object)");
+ return octave_value(-1);
+ }
+ const octave_base_value& rep = args(0).get_rep();
+ const simple& b = ((const simple &)rep);
+ octave_stdout << "MPI_Comm_Test has " << b.name_value() << " output arguments.\n";
+ MPI_Comm res = b.comunicator_value();
+ retval = b.name_value();
+ return retval;
+}
\ No newline at end of file
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_rank.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_rank.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_rank.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,88 @@
+// Copyright (C) 2004-2007 Javier Fernández Baldomero, Mancia Anguita López
+// This code has been adjusted for octave3.2.3 and more in
+// 2009 by Riccardo Corradini <ric...@ya...>
+// Copyright (C) 2009 VZLU Prague
+
+//
+// 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/>.
+
+
+#define NAME MPI_Comm_rank
+/*
+ * ----------------------------------------------------
+ * Determines the rank of the calling process in the communicator
+ * [info rank] = MPI_Comm_rank (comm)
+ * ----------------------------------------------------
+ */
+
+#include "simple.h"
+DEFUN_DLD(NAME, args,nargout ,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} [@var{RANK} @var{INFO}] = MPI_Comm_rank (@var{COMM})\n\
+Determines rank of calling process in communicator.\n\
+If @var{COMM} octave comunicator object loaded with MPI_Comm_Load is omitted \n\
+returns an error. \n\
+ @example\n\
+ @group\n\
+ @var{RANK} rank of the calling process in group of communicator\n\
+ @var{INFO} (int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 5 MPI_ERR_COMM Invalid communicator (NULL?)\n\
+ 13 MPI_ERR_ARG Invalid argument (typically a NULL pointer?)\n\
+SEE ALSO: MPI_Comm_size\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+
+{
+ octave_value_list results;
+ int nargin = args.length ();
+ if (nargin != 1)
+ {
+ error ("expecting 1 input argument");
+ return results;
+ }
+
+ if (!simple_type_loaded)
+ {
+ simple::register_type ();
+ simple_type_loaded = true;
+ mlock ();
+ }
+
+ if((args.length() != 1 )
+ || args(0).type_id()!=simple::static_type_id()){
+
+ error("Please enter octave comunicator object!");
+ return octave_value(-1);
+ }
+
+ const octave_base_value& rep = args(0).get_rep();
+ const simple& B = ((const simple &)rep);
+ MPI_Comm comm = ((const simple&) B).comunicator_value ();
+ if (! error_state)
+ {
+ int my_rank;
+ int info = MPI_Comm_rank (comm, &my_rank);
+ if (nargout > 1)
+ results(1) = info;
+ results(0) = my_rank;
+ }
+ else
+ print_usage ();
+ comm= NULL;
+ /* [rank info] = MPI_Comm_rank (comm) */
+
+ return results;
+
+}
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_size.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_size.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Comm_size.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,89 @@
+// Copyright (C) 2004-2007 Javier Fernández Baldomero, Mancia Anguita López
+// This code has been adjusted for octave3.2.3 and more in
+// 2009 by Riccardo Corradini <ric...@ya...>
+// Copyright (C) 2009 VZLU Prague
+
+//
+// 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/>.
+
+#define NAME MPI_Comm_size
+/*
+ * ----------------------------------------------------
+ * Determines the size of the calling process in the communicator
+ * [info rank] = MPI_Comm_size (comm)
+ * ----------------------------------------------------
+ */
+
+#include "simple.h"
+DEFUN_DLD(NAME, args,nargout ,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} [@var{RANK} @var{INFO}] = MPI_Comm_size (@var{COMM})\n\
+Determines size of calling process in communicator.\n\
+If @var{COMM} octave comunicator object loaded with MPI_Comm_Load is omitted \n\
+returns an error. \n\
+ @example\n\
+ @group\n\
+ @var{exprank} rank of the calling process in group of communicator\n\
+ @var{exprinfo} (int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 5 MPI_ERR_COMM Invalid communicator (NULL?)\n\
+ 13 MPI_ERR_ARG Invalid argument (typically a NULL pointer?)\n\
+SEE ALSO: MPI_Comm_rank\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+
+{
+ octave_value_list results;
+ int nargin = args.length ();
+ if (nargin != 1)
+ {
+ error ("expecting 1 input argument");
+ return results;
+ }
+
+ if (!simple_type_loaded)
+ {
+ simple::register_type ();
+ simple_type_loaded = true;
+ mlock ();
+ }
+
+ if((args.length() != 1 )
+ || args(0).type_id()!=simple::static_type_id()){
+
+ error("Please enter octave comunicator object!");
+ return octave_value(-1);
+ }
+
+ const octave_base_value& rep = args(0).get_rep();
+ const simple& B = ((const simple &)rep);
+ MPI_Comm comm = ((const simple&) B).comunicator_value ();
+
+
+ if (! error_state)
+ {
+ int my_size;
+ int info = MPI_Comm_size (comm, &my_size);
+ if (nargout > 1)
+ results(1) = info;
+ results(0) = my_size;
+ }
+ else
+ print_usage ();
+ comm= NULL;
+ /* [size info] = MPI_Comm_size (comm) */
+
+ return results;
+}
+
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Finalize.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Finalize.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Finalize.cc 2010-03-12 13:06:46 UTC (rev 7052)
@@ -0,0 +1,48 @@
+// Copyright (C) 2004-2007 Javier Fernández Baldomero, Mancia Anguita López
+// This code has been adjusted for octave3.2.3 and more in
+// 2009 by Riccardo Corradini <ric...@ya...>
+
+//
+// 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/>.
+
+#define NAME MPI_Finalize
+/*
+ * ----------------------------------------------------
+ * Terminates MPI execution environment
+ * info = MPI_Finalize
+ * ----------------------------------------------------
+ */
+#include "mpi.h"
+#include <octave/oct.h>
+
+DEFUN_DLD(NAME, args, nargout,"-*- texinfo -*-\n\
+@deftypefn {Loadable Function} {} @var{INFO} = MPI_Finalize()\n\
+ Terminates MPI execution environment\n\
+\n\
+ @example\n\
+ @group\n\
+ @var{INFO} (int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 5 MPI_ERR_COMM Invalid communicator (NULL?)\n\
+ 13 MPI_ERR_ARG Invalid argument (typically a NULL pointer?)\n\
+SEE ALSO: MPI_Init\n\
+@end group\n\
+@end example\n\
+@end deftypefn")
+{
+
+ int info = MPI_Finalize();
+
+ return octave_value(info);
+}
Added: trunk/octave-forge/main/openmpi_ext/src/MPI_Finalized.cc
===================================================================
--- trunk/octave-forge/main/openmpi_ext/src/MPI_Finalized.cc (rev 0)
+++ trunk/octave-forge/main/openmpi_ext/src/MPI_Finalized.cc 2010-03-12 13:06:46...
[truncated message content] |
