Menu
▾
▴
[Octave-cvsupdate] SF.net SVN: octave:[6476] trunk/octave-forge/extra
|
From: <rik...@us...> - 2009-11-16 15:40:02
|
Revision: 6476
http://octave.svn.sourceforge.net/octave/?rev=6476&view=rev
Author: rikcorradini
Date: 2009-11-16 15:39:52 +0000 (Mon, 16 Nov 2009)
Log Message:
-----------
Added Paths:
-----------
trunk/octave-forge/extra/openmpi_ext/
trunk/octave-forge/extra/openmpi_ext/COPYING
trunk/octave-forge/extra/openmpi_ext/DESCRIPTION
trunk/octave-forge/extra/openmpi_ext/INDEX
trunk/octave-forge/extra/openmpi_ext/doc/
trunk/octave-forge/extra/openmpi_ext/doc/README
trunk/octave-forge/extra/openmpi_ext/inst/
trunk/octave-forge/extra/openmpi_ext/inst/hello2dimmat.m
trunk/octave-forge/extra/openmpi_ext/inst/hellosparsemat.m
trunk/octave-forge/extra/openmpi_ext/inst/helloworld.m
trunk/octave-forge/extra/openmpi_ext/src/
trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_rank.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_size.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalize.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalized.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Init.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Initialized.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Recv.cc
trunk/octave-forge/extra/openmpi_ext/src/MPI_Send.cc
trunk/octave-forge/extra/openmpi_ext/src/Makefile
Added: trunk/octave-forge/extra/openmpi_ext/COPYING
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/COPYING (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/COPYING 2009-11-16 15:39:52 UTC (rev 6476)
@@ -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/extra/openmpi_ext/DESCRIPTION
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/DESCRIPTION (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/DESCRIPTION 2009-11-16 15:39:52 UTC (rev 6476)
@@ -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), openmpi1.3.3
+Autoload: ?
+License: GPL version 2 or later
+Url: http://octave.sf.net
Added: trunk/octave-forge/extra/openmpi_ext/INDEX
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/INDEX (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/INDEX 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1 @@
+TO be filled
Added: trunk/octave-forge/extra/openmpi_ext/doc/README
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/doc/README (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/doc/README 2009-11-16 15:39:52 UTC (rev 6476)
@@ -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/extra/openmpi_ext/inst/hello2dimmat.m
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/inst/hello2dimmat.m (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/inst/hello2dimmat.m 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,33 @@
+if not(MPI_Initialized)
+ info = MPI_Init();
+end
+
+ [info my_rank] = MPI_Comm_rank();
+ [info p ] = MPI_Comm_size();
+ 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);
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [info messager] = MPI_Recv(source,mytag);
+# You could also save each result and make comparisons if you don't trust MPI
+ disp("Rank 0 is the master receiving ... :");
+ if (info == true)
+ disp('OK!');
+ endif
+ endfor
+ end
+
+
+ if not(MPI_Finalized)
+ info = MPI_Finalize();
+ end
Added: trunk/octave-forge/extra/openmpi_ext/inst/hellosparsemat.m
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/inst/hellosparsemat.m (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/inst/hellosparsemat.m 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,46 @@
+if not(MPI_Initialized)
+ info = MPI_Init();
+end
+
+ [info my_rank] = MPI_Comm_rank();
+ [info p ] = MPI_Comm_size();
+# 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);
+ 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..");
+ [info messager] = MPI_Recv(source,mytag);
+ disp("Rank 0 is the master receiving ... :");
+ if (messager/message)
+ disp('OK!');
+ endif
+ messager
+ endfor
+ end
Added: trunk/octave-forge/extra/openmpi_ext/inst/helloworld.m
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/inst/helloworld.m (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/inst/helloworld.m 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,56 @@
+## 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
+#h = genpath("/home/user/octave/openmpi_ext/");
+#addpath(h);
+#clear h;
+# 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
+
+
+if not(MPI_Initialized)
+ info = MPI_Init();
+end
+
+ [info my_rank] = MPI_Comm_rank();
+ [info p ] = MPI_Comm_size();
+ # Could be any number
+ mytag = 48;
+
+
+ message="";
+ if (my_rank != 0)
+ message = sprintf('Greetings from process: %d!',my_rank);
+ dest = 0;
+# rankvect is the vector containing the list of rank destination process
+ rankvect = 0;
+ [info] = MPI_Send(message,rankvect,mytag);
+ info
+ else
+ for source = 1:p-1
+ disp("We are at rank 0 that is master etc..");
+ [info message] = MPI_Recv(source,mytag);
+ printf('%s\n', message);
+ endfor
+ end
+
+
+ if not(MPI_Finalized)
+ info = MPI_Finalize();
+ end
\ No newline at end of file
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_rank.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_rank.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_rank.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,59 @@
+// 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_Comm_rank
+/*
+ * ----------------------------------------------------
+ * Determines the rank of the calling process in the communicator
+ * [info rank] = MPI_Comm_rank (comm)
+ * ----------------------------------------------------
+ */
+
+#include <octave/oct.h>
+#include "mpi.h"
+
+DEFUN_DLD(NAME, args, nargout,
+"MPI_Comm_rank Determines rank of calling process in communicator\n\
+\n\
+ [info rank] = MPI_Comm_rank (comm)\n\
+\n\
+ comm (int) communicator handle. MPI_COMM_NULL not valid\n\
+ rank (int) rank of the calling process in group of comm\n\
+\n\
+ 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\
+\n\
+ SEE ALSO: MPI_Comm_size\n\
+ comms\n\
+\n\
+")
+
+{
+ octave_value_list results;
+ int my_rank;
+ int info = MPI_Comm_rank(MPI_COMM_WORLD,&my_rank);
+ results(0) = info;
+ results(1) = my_rank;
+
+ /* [info rank] = MPI_Comm_rank (comm) */
+
+ return results;
+}
\ No newline at end of file
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_size.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_size.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Comm_size.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,59 @@
+// 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_Comm_size
+/*
+ * ----------------------------------------------------
+ * Determines the size of the calling process in the communicator
+ * [info rank] = MPI_Comm_size (comm)
+ * ----------------------------------------------------
+ */
+
+#include <octave/oct.h>
+#include "mpi.h"
+
+DEFUN_DLD(NAME, args, nargout,
+"MPI_Comm_size Determines rank of calling process in communicator\n\
+\n\
+ [info rank] = MPI_Comm_size (comm)\n\
+\n\
+ comm (int) communicator handle. MPI_COMM_NULL not valid\n\
+ size (int) size of the calling process in group of comm\n\
+\n\
+ 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\
+\n\
+ SEE ALSO: MPI_Comm_rank\n\
+ comms\n\
+\n\
+")
+
+{
+ octave_value_list results;
+ int my_size;
+ int info = MPI_Comm_size(MPI_COMM_WORLD,&my_size);
+ results(0) = info;
+ results(1) = my_size;
+
+ /* [info rank] = MPI_Comm_size (comm) */
+
+ return results;
+}
+
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalize.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalize.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalize.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,49 @@
+// 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,
+"MPI_Finalize Terminates MPI execution environment\n\
+\n\
+ info = MPI_Finalize\n\
+\n\
+ info(int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 17 MPI_ERR_INTERN This is fatal. Please send bug report to OPENMPI\n\
+ 16 MPI_ERR_OTHER A different thread attempts to finalize MPI than\n\
+ the thread that initialized MPI\n\
+\n\
+ SEE ALSO: MPI_Init, MPI_Initialized, MPI_Finalized\n\
+ misc\n\
+\n\
+")
+{
+
+ int info = MPI_Finalize();
+
+ return octave_value(info);
+}
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalized.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalized.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Finalized.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,57 @@
+// 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...>
+
+// 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/>.
+
+#define NAME MPI_Finalized
+/*
+ * ----------------------------------------------------
+ * Indicates whether MPI_Finalize has completed
+ * [info flag] = MPI_Finalized
+ * ----------------------------------------------------
+ */
+#include "mpi.h"
+#include <octave/oct.h>
+
+DEFUN_DLD(NAME, args, nargout,
+"MPI_Finalized Indicates whether MPI_Finalize has completed\n\
+\n\
+ [info flag] = MPI_Finalized\n\
+\n\
+ flag(int) 0 false\n\
+ 1 true\n\
+\n\
+ info(int) return code\n\
+ 0 MPI_SUCCESS This function always returns MPI_SUCCESS\n\
+\n\
+ SEE ALSO: MPI_Init, MPI_Finalize, MPI_Initialized\n\
+ misc\n\
+\n\
+")
+{
+ octave_value_list results;
+ int flag;
+
+
+
+ int info = MPI_Finalized(&flag);
+ results(0) = info;
+ results(1) = flag;
+ return results;
+
+ /* [info flag] = MPI_Finalized */
+}
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Init.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Init.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Init.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,71 @@
+// 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/>.
+
+
+
+
+
+
+/*
+ * ----------------------------------------------------
+ * Initialize the MPI execution environment
+ * info = MPI_Init [ ( 'arg' [, 'arg']... ) ]
+ * ----------------------------------------------------
+ */
+
+#define NAME MPI_Init
+
+#include "mpi.h" // mpi.h, oct.h
+#include <octave/oct.h>
+
+DEFUN_DLD(NAME, args, nargout,
+"MPI_Init Initialize the MPI execution environment\n\
+\n\
+ info = MPI_Init [ ( 'arg' [, 'arg']... ) ]\n\
+\n\
+ arg (str) arguments to pass\n\
+ info(int) return code\n\
+ 0 MPI_SUCCESS No error\n\
+ 16 MPI_ERR_OTHER Attempt was made to call MPI_INIT a second time\n\
+ MPI_INIT may only be called once in a program\n\
+ \n\
+\n\
+ SEE ALSO: MPI_Finalize, MPI_Initialized, MPI_Finalized\n\
+ misc\n\
+\n\
+")
+{
+ int nargin = args.length();
+ for (int i=0; i<nargin; i++){
+ if( ! args(i).is_string() ) {
+ error("MPI_Init: args must be strings");
+ return octave_value (MPI_ERR_ARG); // error returns nothing
+ }
+ }
+
+ string_vector argvec = args.make_argv("MPI_Init");
+ char **argve= argvec.c_str_vec();
+ char **argv =&argve[1];
+
+// printf("args: "); for (int i=0; i<nargin; i++) printf("%s ",argv[i]);
+// printf("\n");
+
+ int info = MPI_Init(&nargin, &argv);
+ free(argve);
+ return octave_value (info);
+}
+
\ No newline at end of file
Added: trunk/octave-forge/extra/openmpi_ext/src/MPI_Recv.cc
===================================================================
--- trunk/octave-forge/extra/openmpi_ext/src/MPI_Recv.cc (rev 0)
+++ trunk/octave-forge/extra/openmpi_ext/src/MPI_Recv.cc 2009-11-16 15:39:52 UTC (rev 6476)
@@ -0,0 +1,1432 @@
+// 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/>.
+/*
+ * Receives most Octave datatypes into contiguous memory
+ * using derived datatypes
+ * info = MPI_Send(var,rank)
+ */
+#include "mpi.h"
+#include <octave/oct.h>
+#include <ov-cell.h> // avoid errmsg "cell -- incomplete datatype"
+#include <oct-map.h> // avoid errmsg "Oct.map -- invalid use undef type"
+
+// tested on Octave 3.2.3
+// Octave 3.2.X
+
+enum ov_t_id
+{
+
+ov_unknown=0, // t_id=0
+ov_cell, // t_id=1
+ov_scalar, // t_id=2
+ov_complex_scalar, // t_id=3
+ov_matrix, // t_id=4
+ov_diagonal_matrix, // t_id=5
+ov_complex_matrix, // t_id=6
+ov_complex_diagonal_matrix, // t_id=7
+ov_range, // t_id=8
+ov_bool, // t_id=9
+ov_bool_matrix, // t_id=10
+ov_char_matrix, // t_id=11
+ov_string, // t_id=12
+ov_sq_string, // t_id=13
+ov_int8_scalar, // t_id=14
+ov_int16_scalar, // t_id=15
+ov_int32_scalar, // t_id=16
+ov_int64_scalar, // t_id=17
+ov_uint8_scalar, // t_id=18
+ov_uint16_scalar, // t_id=19
+ov_uint32_scalar, // t_id=20
+ov_uint64_scalar, // t_id=21
+ov_int8_matrix, // t_id=22
+ov_int16_matrix, // t_id=23
+ov_int32_matrix, // t_id=24
+ov_int64_matrix, // t_id=25
+ov_uint8_matrix, // t_id=26
+ov_uint16_matrix, // t_id=27
+ov_uint32_matrix, // t_id=28
+ov_uint64_matrix, // t_id=29
+ov_sparse_bool_matrix, // t_id=30
+
+ov_sparse_matrix, // t_id=31
+ov_sparse_complex_matrix,
+ov_struct,
+ov_class,
+ov_list,
+ov_cs_list,
+ov_magic_colon,
+ov_built_in_function,
+ov_user_defined_function,
+ov_dynamically_linked_function,
+ov_function_handle,
+ov_inline_function,
+ov_float_scalar,
+ov_float_complex_scalar,
+ov_float_matrix,
+ov_float_diagonal_matrix,
+ov_float_complex_matrix,
+ov_float_complex_diagonal_matrix,
+ov_permutation_matrix,
+ov_null_matrix,
+ov_null_string,
+ov_null_sq_string,
+};
+
+
+
+/*----------------------------------*/ /* forward declaration */
+
+
+
+int recv_class( octave_value &ov, int source, int mytag); /* along the datatype */
+/*----------------------------------*/ /* to send any octave_value */
+
+
+int recv_string( octave_value &ov, int source, int mytag){ /* directly MPI_Send it, */
+/*-----------------------------*/ /* it's just a string value */
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 2);
+tanktag[0]=mytag;
+tanktag[1]=mytag+1;
+tanktag[2]=mytag+2;
+
+ int info,nitem;
+ MPI_Status stat;
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ OCTAVE_LOCAL_BUFFER(char,mess,nitem+1);
+ if (info !=MPI_SUCCESS) return info;
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_CHAR, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+
+ info = MPI_Recv(mess, 1,fortvec, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("Flag for string received %i \n",info);
+ std::string cpp_string;
+ cpp_string = mess;
+ ov = cpp_string;
+ if (info !=MPI_SUCCESS) return info;
+ return(MPI_SUCCESS);
+}
+
+int recv_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 5);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+
+ dv(i) = dimV[i] ;
+// printf("I am printing dimvector %i \n",dimV[i]);
+ }
+ NDArray myNDA(dv);
+// double *p = myNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(double,LBNDA,nitem);
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_DOUBLE, &fortvec);
+ MPI_Type_commit(&fortvec);
+// printf("I am printing dimvector %i \n",tanktag[4]);
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ myNDA(i)=LBNDA[i];
+ }
+ ov = myNDA;
+// printf("info for receiving scalar matrix is = %i \n",info);
+
+}
+
+int recv_float_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 5);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+// printf("I am printing dimvector %i \n",dimV[i]);
+ dv(i) = dimV[i] ;
+ }
+ FloatNDArray myNDA(dv);
+ float *p = myNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(float,LBNDA,nitem);
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_FLOAT, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ myNDA(i)=LBNDA[i];
+ }
+ ov = myNDA;
+// printf("info for receiving scalar matrix is = %i \n",info);
+
+}
+
+
+
+int recv_complex_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 6);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+tanktag[5] = mytag+5;
+
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+// printf("I am printing dimvector %i \n",i);
+ dv(i) = dimV[i] ;
+ }
+ NDArray rmyNDA(dv); // for the real part
+ NDArray imyNDA(dv); // for the img part
+
+ ComplexNDArray oa (dv); /* Create->Delete on ret?!? */
+
+ double *p = rmyNDA.fortran_vec();
+ double *ip = imyNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(double,LBNDA,nitem);
+ OCTAVE_LOCAL_BUFFER(double,CLBNDA,nitem);
+
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_DOUBLE, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+ info = MPI_Recv((CLBNDA), 1,fortvec, source, tanktag[5] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ rmyNDA(i)=LBNDA[i];
+ imyNDA(i)=CLBNDA[i];
+ oa(i) = Complex (rmyNDA (i), imyNDA (i));
+ }
+ ov = oa;
+
+return(MPI_SUCCESS);
+}
+
+
+int recv_float_complex_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 6);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+tanktag[5] = mytag+5;
+
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+// printf("I am printing dimvector %i \n",i);
+ dv(i) = dimV[i] ;
+ }
+ FloatNDArray rmyNDA(dv); // for the real part
+ FloatNDArray imyNDA(dv); // for the img part
+
+ FloatComplexNDArray oa (dv); /* Create->Delete on ret?!? */
+
+ float *p = rmyNDA.fortran_vec();
+ float *ip = imyNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(float,LBNDA,nitem);
+ OCTAVE_LOCAL_BUFFER(float,CLBNDA,nitem);
+
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_DOUBLE, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+ info = MPI_Recv((CLBNDA), 1,fortvec, source, tanktag[5] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ rmyNDA(i)=LBNDA[i];
+ imyNDA(i)=CLBNDA[i];
+ std::complex<float> c = real(rmyNDA(i))+imag(imyNDA(i));
+ oa(i) = c;
+ }
+ ov = oa;
+
+return(MPI_SUCCESS);
+}
+
+
+
+
+
+
+int recv_sp_mat(octave_value &ov,int source, int mytag){
+int info;
+/* int nr = m.rows();
+ int nc = m.cols ();
+ int nnzm = m.nzmax ();*/
+// 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[5] ----> vector of number of non zero elements
+OCTAVE_LOCAL_BUFFER(int, tanktag,6);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+tanktag[5] = mytag+5;
+
+MPI_Status stat;
+
+OCTAVE_LOCAL_BUFFER(int,s,3);
+
+// Create a contiguous derived datatype
+MPI_Datatype sintsparse;
+MPI_Type_contiguous(3,MPI_INT, &sintsparse);
+MPI_Type_commit(&sintsparse);
+
+// receive the sintsparse vector named s
+info = MPI_Recv(s, 1, sintsparse, source, tanktag[1], MPI_COMM_WORLD, &stat);
+if (info !=MPI_SUCCESS) return info;
+
+// int appo1, appo2, appo3;
+// appo1 = s[0];
+// appo2 = s[1];
+// appo3 = s[2];
+//
+// printf("appo 1 is %i \n",appo1);
+// printf("appo 2 is %i \n",appo2);
+// printf("appo 3 is %i \n",appo3);
+
+
+
+SparseMatrix m(s[0],s[1],s[2]);
+
+// Create a contiguous derived datatype for row index
+OCTAVE_LOCAL_BUFFER(int,sridx,s[2]);
+OCTAVE_LOCAL_BUFFER(double,sdata,s[2]);
+// Create a contiguous derived datatype for column index
+OCTAVE_LOCAL_BUFFER(int,scidx,s[1]+1);
+MPI_Datatype rowindex;
+MPI_Type_contiguous(s[2],MPI_INT, &rowindex);
+MPI_Type_commit(&rowindex);
+MPI_Datatype datavect;
+MPI_Type_contiguous(s[2],MPI_DOUBLE, &datavect);
+MPI_Type_commit(&datavect);
+MPI_Datatype columnindex;
+MPI_Type_contiguous(s[1]+1,MPI_INT, &columnindex);
+MPI_Type_commit(&columnindex);
+
+// Now receive the two vectors
+
+// receive the vector with row indexes
+ info = MPI_Recv(sridx,1,rowindex,source,tanktag[2],MPI_COMM_WORLD,&stat);
+// printf("Hope everything is fine here with ridx %i =\n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+// receive the vector with column indexes
+ info = MPI_Recv(scidx,1,columnindex,source,tanktag[3],MPI_COMM_WORLD, &stat);
+// printf("Hope everything is fine here with scidx %i =\n",info);
+ if (info !=MPI_SUCCESS) return info;
+
+ info = MPI_Recv(sdata,1,datavect,source,tanktag[4],MPI_COMM_WORLD,&stat);
+// printf("Hope everything is fine here with data vect %i =\n",info);
+
+// int appoc;
+NDArray buf (dim_vector (m.capacity(), 1));
+for (octave_idx_type i = 0; i < s[2]; i++)
+{
+ m.ridx(i) = sridx[i];
+ m.data(i) = sdata[i];
+ buf(i) = sdata[i];
+// appoc = m.ridx(i);
+// printf("apporow is %i \n",appoc);
+// printf("appodata is %d \n",sdata[i]);
+}
+for (octave_idx_type i = 0; i < s[1]+1; i++)
+{
+ m.cidx(i) = scidx[i];
+// appoc=m.cidx(i);
+// printf("appocol is %i \n",appoc);
+}
+
+
+ov = m;
+return(info);
+
+
+}
+
+
+//here ! int matrixes
+int recv_int8_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 5);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[4] = mytag+4;
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+// printf("I am printing dimvector %i \n",i);
+ dv(i) = dimV[i] ;
+ }
+ int8NDArray myNDA(dv);
+// double *p = myNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(octave_int8,LBNDA,nitem);
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_BYTE, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ myNDA(i)=LBNDA[i];
+ }
+ ov = myNDA;
+// printf("info for receiving scalar matrix is = %i \n",info);
+
+}
+
+int recv_int16_matrix(octave_value &ov,int source, int mytag){
+
+OCTAVE_LOCAL_BUFFER(int, tanktag, 5);
+tanktag[0] = mytag;
+tanktag[1] = mytag+1;
+tanktag[2] = mytag+2;
+tanktag[3] = mytag+3;
+tanktag[3] = mytag+4;
+ int info;
+ int nitem,nd;
+ MPI_Status stat;
+ dim_vector dv;
+
+// nitem is the total number of elements
+ info = MPI_Recv((&nitem), 1,MPI_INT, source, tanktag[1] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of elements %i \n",nitem);
+ if (info !=MPI_SUCCESS) return info;
+// ndims is number of dimensions
+ info = MPI_Recv((&nd), 1,MPI_INT, source, tanktag[2] , MPI_COMM_WORLD,&stat);
+// printf("I have received number of dimensions %i \n",nd);
+ if (info !=MPI_SUCCESS) return info;
+// Now create contiguos datatype for dim vector
+ dv.resize(nd);
+ OCTAVE_LOCAL_BUFFER(int,dimV,nd);
+ MPI_Datatype dimvec;
+ MPI_Type_contiguous(nd,MPI_INT, &dimvec);
+ MPI_Type_commit(&dimvec);
+
+ info = MPI_Recv((dimV), 1,dimvec, source, tanktag[3] , MPI_COMM_WORLD,&stat);
+ if (info !=MPI_SUCCESS) return info;
+
+// Now reverse the content of int vector into dim vector
+ for (octave_idx_type i=0; i<nd; i++)
+ {
+// printf("I am printing dimvector %i \n",i);
+ dv(i) = dimV[i] ;
+ }
+ int16NDArray myNDA(dv);
+// double *p = myNDA.fortran_vec();
+ OCTAVE_LOCAL_BUFFER(octave_int16,LBNDA,nitem);
+// printf("BUFFER CREATED \n");
+ // Now create the contiguous derived datatype
+ MPI_Datatype fortvec;
+ MPI_Type_contiguous(nitem,MPI_SHORT, &fortvec);
+ MPI_Type_commit(&fortvec);
+
+ info = MPI_Recv((LBNDA), 1,fortvec, source, tanktag[4] , MPI_COMM_WORLD,&stat);
+// printf("info for receiving data is = %i \n",info);
+ if (info !=MPI_SUCCESS) return info;
+ for (octave_idx_type i=0; i<nitem; i++)
+ {
+// *LBNDA = *p;
+// LBNDA++;
+// p++;
+ myNDA(i)=LBNDA[i];
+ }
+ ov = myNDA;
+// printf("info for receiving scalar matrix is = %i \n",info);
+ return(info);
+}
+
+
+int recv_int32_matrix(octave_value &ov,int source, int mytag){
+
+OCTA...
[truncated message content] |
