<div class="gmail_quote">On Wed, Aug 24, 2011 at 09:33, Wienand Drenth <span dir="ltr">&lt;<a href="mailto:w.drenth@gmail.com">w.drenth@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":2d">The program compiles without troubles, but running it gives me an error. It is, I think, related to the MPI_SUM operation:<br>
[walrus:463] *** An error occurred in MPI_Scan: the reduction operation MPI_SUM is not defined on the MPI_BYTE datatype<br>[walrus:463] *** on communicator MPI_COMM_WORLD<br>[walrus:463] *** MPI_ERR_OP: invalid reduce operation<br>

[walrus:463] *** MPI_ERRORS_ARE_FATAL (your MPI job will now abort)<br><br>In addition, I needed to declare MPIU_INT explicitly as integer; this was not required for the C program.</div></blockquote></div><br><div>I believe MPIU_INT (which is normally made to match PetscInt) is not defined in Fortran. You can use MPI_INTEGER4 or MPI_INTEGER8 as needed.</div>
<div><br></div><div>Perhaps PETSc should define MPIU_INT, MPIU_SCALAR, etc in Fortran too?</div>