PETSc, hangs, on, PetscInitialize

[Sanjay Govindjee]

As the author of this bit of code, I can suggest the following issue.  
FEAP  relies on an mpi implementation
that send standard fortran writes, things like write(*,*),
to the screen when issued, but this is in no way gaurenteed by the mpi 
standard.  FEAP is probably sitting at a prompt waiting for your input 
but you simply can not see it.  I have experienced this
problem on several new installations.  There are two solutions: (1) Run 
in a strictly batch mode, making sure
that your feapname file is setup correctly, and that you have a means to 
read inputs from a file for the prompts that will appear (in this case a 
file with the letter y in it to accept the file names read from 
feapname).  (2) re-write all the old non-parallel IO in the code 
(something I am thinking about but which is a lot of work to set it up 
so that it works\
in serial mode without petsc and in parallel mode with it). 

The fact that it looks like it is hanging in the initialize routine is 
probably a red-herring. 

-sg

Barry Smith wrote:
>   I would run in the debugger and then stop it after it has
> hanged a long time (control c in gdb or dbx) then do where to
> determine what it is doing. Send the output from where to 
> petsc-maint at mcs.anl.gov if you need to.
>
>
>    Barry
>
> On Fri, 15 Jun 2007, Marcus Heck wrote:
>
>   
>> Hi,
>>
>> I am runing PETSc on the FE-solver FEAP. If I run a normal FE analysis
>> everithing is fine!
>>
>> Then I want do run a coupled FSI computation. So I initialize PETSc and
>> MpCCI (my coupling library) or the other way round (MpCCI, PETSc). But
>> PETSc enters the PetscInitialize and never comes out. Does anyone has an
>> idea.
>>
>> Thanks a lot
>> Marcus
>>
>> Here some infos about PETSc and my code:
>> **************
>> PETSc config:
>> ./config/configure.py --prefix=/amd/software/linux/petsc-2.3.2-p10 --with-cc=gcc --with-fc=g77 --LD_SHARED= g++ --LDFLAGS="-O1" --with-mpi-dir=/usr/lib/mpich --with-mpi-shared=1
>>
>>
>> Compile and linking command in both cases are identical, except the
>> MpCCI library includes.
>>
>> If there is any idea let me know please
>> Thanks a lot in advance
>> Marcus
>>
>> **************
>> CODE inside FEAP75.f
>>       write(*,*) 'FEAP starting MpCCI setup and interface definition.'
>>       call CCI_Init( iError )
>>       write(*,*) 'FEAP starting MpCCI done'
>>       write(*,*)'pcontr(): initializing PETSC'
>>       call PetscInitialize(PETSC_NULL_CHARACTER,petsc_ierr)
>>       write(*,*)'pcontr(): initializing PETSC done'
>>       call pstart()
>>
>> **************
>> CODE OUTPUT:
>>  FEAP starting MpCCI setup and interface definition.
>> ===========================================================================
>> MpCCI 3.0.3 (for MPICH 1.2.5) - compilation date Dec 22 2004 20:08:09
>> FhG-SCAI, Fraunhofer Institute for Algorithms and Scientific Computing
>> Schloss Birlinghoven, D-53754 Sankt Augustin, Germany
>> Coupled computation started at 15-6-2007 11:16:16.
>> Running on host 'ares'.
>> ===========================================================================
>>
>>                    ************************************
>>                    *    This MpCCI 3.0 Professional   *
>>                    *         is licensed for          *
>>                    * TU Darmstadt, FB Maschinenbau, M *
>>                    *    until date 28-feb-2008.       *
>>                    ************************************
>>
>>  FEAP starting MpCCI done
>>  pcontr(): initializing PETSC
>>
>> **************
>> strace on the pid tells:
>> 3p$ strace -p17164
>> Process 17164 attached - interrupt to quit
>> select(5, [3 4], NULL, NULL, {6, 266000}) = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0})    = 0 (Timeout)
>> select(5, [3 4], NULL, NULL, {9, 0}
>> -- 8< --- snip -- 8< ---
>>
>> **************
>> Totalview telals following:
>>
>> main (f77) -> perscinitialize_  (C) -> PETScParseFortranArgs_Private (C)
>> -> PMPI_Bcast -> inta_Bcast -> ...
>>
>> Please see screen shot.
>>
>>
>>     
>
>