[OMPI users] Library interposing works with 1.6.5 but not with 2.0.1 (Fortran)?

2017-01-16 Thread Alef Farah
Hi, I contribute to a tracing library which uses PMPI. It's loaded with LD_PRELOAD so it interposes libmpi, intercepting MPI_ calls. Since we upgraded from OpenMPI 1.6.5 to OpenMPI 2.0.1 it seems to have stopped intercepting the calls from Fortran application, although it continues to work with C

Re: [OMPI users] Rounding errors and MPI

2017-01-16 Thread Oscar Mojica
Thanks guys for your answers. Actually, the optimization was not disabled, and that was the problem, compiling it with -o0 solves it. Sorry. Oscar Mojica Geologist Ph.D. in Geophysics SENAI CIMATEC Supercomputing Center Lattes: http://lattes.cnpq.br/0796232840554652

Re: [OMPI users] Rounding errors and MPI

2017-01-16 Thread Yann Jobic
Hi, Is there an overlapping section in the MPI part ? Otherwise, please check : - declaration type of all the variables (consistency) - correct initialization of the array "wave" (to zero) - maybe use temporary variables like real size1,size2,factor size1 = dx+dy size2 = dhx+dhy factor = dt*size

Re: [OMPI users] Rounding errors and MPI

2017-01-16 Thread Tim Prince via users
You might try inserting parentheses so as to specify your preferred order of evaluation. If using ifort, you would need -assume protect-parens . Sent via the ASUS PadFone X mini, an AT&T 4G LTE smartphone Original Message From:Oscar Mojica Sent:Mon, 16 Jan 2017 08:28:05 -0500

[OMPI users] Rounding errors and MPI

2017-01-16 Thread Oscar Mojica
Hello everybody I'm having a problem with a parallel program written in fortran. I have a 3D array which is divided in two in the third dimension so thats two processes perform some operations with a part of the cube, using a subroutine. Each process also has the complete cube. Before each proc