[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: ***LOOPING MAIL*** PetscInitialize problem?

  Are you sure that mpif.h is actually being included
by include/finclude/petsc.h?

   Barry


On Tue, 14 Mar 2006, ziemer@lncc.br wrote:


This is part of the source code where the problem happens:

Thank you for your patience.

*****************************************************************
PROGRAM Premain
C     ------------------------------------------------------------------
!     Use ModuloLectura
     IMPLICIT REAL*8 (A-H,O-Z)
!
!     Parallel begin
!     petsc include files
#include "include/finclude/petsc.h"
!     Parallel end
!
     Logical Sym, NonLinearProblem,Iterative
     INTEGER ERR_ALLOC
     Parameter (Lch=16)
     Character Str*120
     double precision noMPI_Wtime
C
!     INTERFACE ReadMat
!     MODULE PROCEDURE ReadMat
!     END INTERFACE
!     INTERFACE ReadCoupling
!     MODULE PROCEDURE ReadCoupling
!     END INTERFACE
!     INTERFACE Dimen
!     MODULE PROCEDURE Dimen
!     END INTERFACE
C
     Allocatable X         (:)
     Allocatable Ie__Param (:)
     Allocatable Ie_JParam (:)
     Allocatable Param     (:)

!     Parallel begin
     call PetscInitialize (PETSC_NULL_CHARACTER, iError)
     if (iError.ne.0) then
        Write (6,*) "Can't initialize PETSC"
        Stop
     End if
!     Ask total number of proccesses
!
!     PROBLEM BEGINS IN THE NEXT FUNCTION - MPI_COMM_WORLD IS NOT SET!!
     call MPI_Comm_size (MPI_COMM_WORLD, iSize, iError)
!     ID of this proccess
     call MPI_Comm_rank (MPI_COMM_WORLD, iRank, iError)
!
**********************************************************


Quoting Satish Balay <balay@mcs.anl.gov>:


You should be able to use either MPI_COMM_WORLD or PETSC_COMM_WORLD or
any other communicator. There must be something else wrong in your
code. Can't comment without looking at the code.

Satish

On Mon, 13 Mar 2006, ziemer@lncc.br wrote:


Hi,

The Petsc installation examples run ok.

The strange thing is if I change the parameter MPI_COMM_WORLD of the

function

MPI_Comm_size() to PTESC_COMM_WORLD this function works, but at same time

this

same strategy does not work with the MPI_Bcast() function.

Thank you very much for your nice help.

Paulo Ziemer.


Quoting Satish Balay <balay@mcs.anl.gov>:


According to the error message - the communicator is invalid.  You
might want to check that..

You should be able to call MPI_Comm_size() after PetscInitialize()
Check src/sys/examples/tutorials/ex1.c for an example of this usage.

BTW: Do PETSc examples run fine with your install?

Satish


On Mon, 13 Mar 2006, ziemer@lncc.br wrote:


Hello,

I would be very glad if someone could give me some direction:

After the MPI_Comm_size function is called, I got this error:

aborting job:
Fatal error in MPI_Comm_size: Invalid communicator, error stack:
MPI_Comm_size(110): MPI_Comm_size(comm=0x1, size=0xbfffdf78) failed
MPI_Comm_size(69): Invalid communicator

It seems that the MPI environment is not being correctly created by the

function

PetscInitialize.


Does anyone have any idea?


Best Regards.

----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.








----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.








----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.