<div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 31, 2015 at 4:13 PM, Matthew Knepley <span dir="ltr"><<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><span class="">On Tue, Mar 31, 2015 at 3:05 PM, Mark Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Humm, that is what I see (appended).  I put a print in PETSC and see this but I still have PETSC_COMM_SELF==0 in my fortran code.  <br></div></div></blockquote><div><br></div></span><div>Is it possible that you do not have the right #include? Can you put exactly the same includes as ex61f</div></div></div></div></blockquote><div><br></div><div>Yea, I tried that. and changed it to .h90<br><br></div><div>BTW, I spent most of Friday finding that if I did not include the .h90 file ISGetIndicesF90 would segv at the call site.<br><br><br></div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_extra"><div class="gmail_quote"><div><br></div><div>  Matt</div><div><div class="h5"><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div></div>Strange, the compile lines look close enough and the example ex61f sees a good PETSC_COMM_SELF ...<br><div><br></div><div>Mark<br></div><div><div><div><br>(gdb) where<br>#0  PetscInitializeFortran () at /Users/markadams/Codes/petsc/src/sys/objects/ftn-custom/zstartf.c:48<br>#1  0x000000010027d985 in petscinitialize_ (filename=0x101eaab40 "", ierr=0x7fff5fbff18c, len=80) at /Users/markadams/Codes/petsc/src/sys/objects/ftn-custom/zstart.c:398<br>#2  0x00000001000015b2 in tpetsc () at ex61f.F90:111<br>#3  0x0000000100002a3d in main (argc=1, argv=0x7fff5fbff968 "/Users/markadams/Codes/petsc/src/ksp/ksp/examples/tutorials/ex61f") at ex61f.F90:60<br>Current language:  auto; currently minimal<br>(gdb) n<br>50        if (PETSC_COMM_WORLD) c1 =  MPI_Comm_c2f(PETSC_COMM_WORLD);<br>(gdb) p c1<br>$1 = 0<br>(gdb) n <br>51        c2 =  MPI_Comm_c2f(PETSC_COMM_SELF);<br>(gdb) p c1<br>$2 = 1140850688<br>(gdb) <br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Mar 31, 2015 at 1:32 PM, Satish Balay <span dir="ltr"><<a href="mailto:balay@mcs.anl.gov" target="_blank">balay@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">1. Check if this is reproduceable with a petsc example for your build.<br>
<br>
2. if you run in the debugger - and place a breakpoint in<br>
PetscInitializeFortran() - you should see the following trace.<br>
<br>
(gdb) where<br>
#0  PetscInitializeFortran () at /home/balay/petsc/src/sys/objects/ftn-custom/zstartf.c:48<br>
#1  0x00007ffff66c396d in petscinitialize_ (filename=0x603640 <petscfortran1_> "",<br>
    ierr=0x7fffffffe0bc, len=80) at /home/balay/petsc/src/sys/objects/ftn-custom/zstart.c:398<br>
#2  0x0000000000401a42 in MAIN__ () at ex2f.F:102<br>
(gdb) n<br>
50      if (PETSC_COMM_WORLD) c1 =  MPI_Comm_c2f(PETSC_COMM_WORLD);<br>
(gdb) p c1<br>
$1 = 0<br>
(gdb) n<br>
51      c2 =  MPI_Comm_c2f(PETSC_COMM_SELF);<br>
(gdb) p c1<br>
$2 = 1140850688<br>
<span><font color="#888888"><br>
<br>
Satish<br>
</font></span><div><div><br>
On Tue, 31 Mar 2015, Mark Adams wrote:<br>
<br>
> I'm not sure where we are.  My code works on Edison but fails on OSX.<br>
> Should I call PetscInitializeFortran?  I do not have a C main.<br>
><br>
><br>
> On Tue, Mar 31, 2015 at 12:14 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
> ><br>
> > PetscErrorCode PetscInitializeFortran(void)<br>
> > {<br>
> >   MPI_Fint c1=0,c2=0;<br>
> ><br>
> >   if (PETSC_COMM_WORLD) c1 =  MPI_Comm_c2f(PETSC_COMM_WORLD);<br>
> >   c2 =  MPI_Comm_c2f(PETSC_COMM_SELF);<br>
> >   petscsetcommonblock_(&c1,&c2);<br>
> ><br>
> > Note that this routine is called automatically by the Fortran<br>
> > petscinitialize function.<br>
> ><br>
> >   Barry<br>
> ><br>
> > > On Mar 31, 2015, at 10:19 AM, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>> wrote:<br>
> > ><br>
> > > On Tue, Mar 31, 2015 at 10:13 AM, Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br>
> > > No, I don't have a C main program (like the doc says).<br>
> > ><br>
> > > I can see COMM_WORLD set in the Fortran initialize wrapper, but where is<br>
> > COMM_SELF set?<br>
> > ><br>
> > >    MAtt<br>
> > ><br>
> > > On Tue, Mar 31, 2015 at 11:05 AM, Matthew Knepley <<a href="mailto:knepley@gmail.com" target="_blank">knepley@gmail.com</a>><br>
> > wrote:<br>
> > > On Tue, Mar 31, 2015 at 5:44 AM, Mark Adams <<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>> wrote:<br>
> > > I'm finding that my fortran code on OSX has PETSC_COMM_SELF == 0.<br>
> > ><br>
> > > I can not find where PETSC_COMM_SELF is created.<br>
> > ><br>
> > > Its set in PetscInitializeFortran(). Did you call that?<br>
> > ><br>
> > >   Matt<br>
> > ><br>
> > ><br>
> > > Mark<br>
> > ><br>
> > ><br>
> > ><br>
> > > --<br>
> > > What most experimenters take for granted before they begin their<br>
> > experiments is infinitely more interesting than any results to which their<br>
> > experiments lead.<br>
> > > -- Norbert Wiener<br>
> > ><br>
> > ><br>
> > ><br>
> > ><br>
> > > --<br>
> > > What most experimenters take for granted before they begin their<br>
> > experiments is infinitely more interesting than any results to which their<br>
> > experiments lead.<br>
> > > -- Norbert Wiener<br>
> ><br>
> ><br>
><br>
<br>
</div></div></blockquote></div><br></div>
</blockquote></div></div></div><div><div class="h5"><br><br clear="all"><div><br></div>-- <br><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div></div></div>
</blockquote></div><br></div></div>