<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:SimSun;
        panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
        {font-family:SimSun;
        panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:SimSun;
        panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri","sans-serif";
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.25in 1.0in 1.25in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal>Dear MPICH2 users and developers,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am recently learning MPI and reading the source code of a MPI implemented open software. One of the functions which does something strange is MPI_Allreduce.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am working on Windows 7 pro 64bit machine with MinGW (gcc version 4.6.2, 32bit) and using MPICH2 1.4.1-p1 32bit version downloaded from the MPICH2 site. The code can be compiled without any problem, but however it failed when running (maybe invalid memory access?). There must be some problem with the windows version MPI_Allreduce since it works fine if I remove that line. And it also works if I make the matrix smaller. I tried it on a Ubuntu machine with same version MPI as well. No problem in Linux.<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>To make the question clear, I added MPI_Allreduce into a hello world code. The code is written in F90. I haven't test the c version of it but I think they should be very similar (differed by the function name and the error parameter).<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here is the command that I used to compile:<o:p></o:p></p><p class=MsoNormal>gfortran hello1.f90 -g -o hello.exe -IC:\MPICH2_x86\include -LC:\MPICH2_x86\lib -lfmpich2g<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here is the source code:<o:p></o:p></p><p class=MsoNormal>!------------------code begin----------------------<o:p></o:p></p><p class=MsoNormal> program main<o:p></o:p></p><p class=MsoNormal>  include 'mpif.h'<o:p></o:p></p><p class=MsoNormal>  character * (MPI_MAX_PROCESSOR_NAME) processor_name<o:p></o:p></p><p class=MsoNormal>  integer myid, numprocs, namelen, rc,ierr<o:p></o:p></p><p class=MsoNormal>  integer, allocatable :: mat1(:, :, :)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>  call MPI_INIT( ierr )<o:p></o:p></p><p class=MsoNormal>  call MPI_COMM_RANK( MPI_COMM_WORLD, myid, ierr )<o:p></o:p></p><p class=MsoNormal>  call MPI_COMM_SIZE( MPI_COMM_WORLD, numprocs, ierr )<o:p></o:p></p><p class=MsoNormal>  call MPI_GET_PROCESSOR_NAME(processor_name, namelen, ierr)<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>  allocate(mat1(-36:36, -36:36, -36:36))<o:p></o:p></p><p class=MsoNormal>  mat1(:,:,:) = 0<o:p></o:p></p><p class=MsoNormal>  call MPI_Bcast( mat1(-36, -36, -36), 389017, MPI_INT, 0, MPI_COMM_WORLD, ierr )<o:p></o:p></p><p class=MsoNormal>  call MPI_Allreduce(MPI_IN_PLACE, mat1(-36, -36, -36), 389017, MPI_INTEGER, MPI_BOR, MPI_COMM_WORLD, ierr)<o:p></o:p></p><p class=MsoNormal>  print *,"MPI_Allreduce done!!!"<o:p></o:p></p><p class=MsoNormal>  print *,"Hello World! Process ", myid, " of ", numprocs, " on ", processor_name<o:p></o:p></p><p class=MsoNormal>  call MPI_FINALIZE(rc)<o:p></o:p></p><p class=MsoNormal>  end<o:p></o:p></p><p class=MsoNormal>!------------------code end----------------------<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>When I use gdb (comes with MinGW) to check (gdb hello.exe then backtrace). I got something meaningless (or seems to be for myself):<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Program received signal SIGSEGV, Segmentation fault.<o:p></o:p></p><p class=MsoNormal>[Switching to Thread 16316.0x4fd0]<o:p></o:p></p><p class=MsoNormal>0x01c03100 in mpich2nemesis!PMPI_Wtime ()<o:p></o:p></p><p class=MsoNormal>   from C:\Windows\system32\mpich2nemesis.dll<o:p></o:p></p><p class=MsoNormal>(gdb) backtrace<o:p></o:p></p><p class=MsoNormal>#0  0x01c03100 in mpich2nemesis!PMPI_Wtime ()<o:p></o:p></p><p class=MsoNormal>   from C:\Windows\system32\mpich2nemesis.dll<o:p></o:p></p><p class=MsoNormal>#1  0x0017be00 in ?? ()<o:p></o:p></p><p class=MsoNormal>#2  0x00000000 in ?? ()<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Does this actually mean there are something wrong with the windows version MPI library?<o:p></o:p></p><p class=MsoNormal>What will be the solution to make it work?<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Thanks.<o:p></o:p></p></div></body></html>