[petsc-users] PETSC_VIEWER_DRAW_(Comm) Fortran interface missing?

[Tim Steinhoff]

Yeah, thanks Satish.
PETSC_VIEWER_DRAW_() works perfectly and I could not reproduce my
above mentioned issue with PETSC_VIEWER_STDOUT_(), so probably I
messed something up.
Thanks again for your great support,
Volker

2018-03-06 18:14 GMT+01:00 Satish Balay <balay at mcs.anl.gov>:
> Ok - I've added ftn interface for PETSC_VIEWER_DRAW_() [similar to
> PETSC_VIEWER_STDOUT_()] in branch balay/ftn-VIEWER_DRAW/maint.
>
> You can give it a try and see if it works.
>
> Satish
>
> On Tue, 6 Mar 2018, Satish Balay wrote:
>
>> Hm - I do see the fortran interface code for PETSC_VIEWER_STDOUT_()
>>
>> https://bitbucket.org/petsc/petsc/src/049452fa31a1384b9850ba16ebcec79c99e4198c/src/sys/classes/viewer/impls/ascii/ftn-custom/zvcreatef.c?at=master&fileviewer=file-view-default
>>
>> and it does appear to work for me.
>>
>> balay at asterix /home/balay/petsc/src/mat/examples/tests (master *=)
>> $ git diff
>> diff --git a/src/mat/examples/tests/ex16f90.F90 b/src/mat/examples/tests/ex16f90.F90
>> index 6bb83b4644..fee44aea95 100644
>> --- a/src/mat/examples/tests/ex16f90.F90
>> +++ b/src/mat/examples/tests/ex16f90.F90
>> @@ -15,6 +15,7 @@
>>        PetscInt one
>>        PetscScalar  v(1)
>>        PetscScalar, pointer :: array(:,:)
>> +      PetscViewer vw
>>
>>
>>        call PetscInitialize(PETSC_NULL_CHARACTER,ierr)
>> @@ -49,7 +50,8 @@
>>  !
>>  !       Print the matrix to the screen
>>  !
>> -      call MatView(A,PETSC_VIEWER_STDOUT_WORLD,ierr);CHKERRA(ierr)
>> +      vw = PETSC_VIEWER_STDOUT_(PETSC_COMM_WORLD)
>> +      call MatView(A,vw,ierr);CHKERRA(ierr)
>>
>>
>>  !
>> balay at asterix /home/balay/petsc/src/mat/examples/tests (master *=)
>> $ mpiexec -n 2 ./ex16f90
>> Mat Object: 2 MPI processes
>>   type: mpidense
>> 9.0000000000000000e+00 4.5000000000000000e+00
>> 4.5000000000000000e+00 3.0000000000000000e+00
>> 3.0000000000000000e+00 2.2500000000000000e+00
>> balay at asterix /home/balay/petsc/src/mat/examples/tests (master *=)
>> $
>>
>>
>> Satish
>>
>>
>> On Tue, 6 Mar 2018, Tim Steinhoff wrote:
>>
>> > Hi Barry,
>> >
>> > thanks for your fast response. Yes, I'm fine with that.
>> > Indeed, my code with the routine PETSC_VIEWER_STDOUT_(comm) was
>> > compilable, but turned out it did not work properly at runtime
>> > (hangs).
>> > When using PETSC_VIEWER_DRAW_(comm), the compiler immediately gives an
>> > error and I was confused.
>> >
>> > Thanks,
>> >
>> > Volker
>> >
>> > By the way, the output parameter description of
>> > PetscViewerASCIIGetStdout is missing in the documentation:
>> > http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Viewer/PetscViewerASCIIGetStdout.html
>> >
>> >
>> > 2018-03-06 16:26 GMT+01:00 Smith, Barry F. <bsmith at mcs.anl.gov>:
>> > >
>> > >
>> > >    Tim,
>> > >
>> > >      The PETSC_VIEWER_STDOUT_(comm)  construct doesn't work in Fortran. Only the
>> > > PETSC_VIEWER_STDOUT_WORLD or PETSC_VIEWER_STDOUT_SELF work. This is also true for
>> > > draw see src/ksp/ksp/examples/tutorials/ex100f.F90
>> > >
>> > >    Note that you can use PetscViewerDrawOpen() to get a draw to open on any communicator you like.
>> > >
>> > >
>> > >    Does this answer your question or did I misunderstand?
>> > >
>> > >    Barry
>> > >
>> > >
>> > >
>> > >
>> > >> On Mar 6, 2018, at 9:02 AM, Tim Steinhoff <kandanovian at gmail.com> wrote:
>> > >>
>> > >> Hi all,
>> > >>
>> > >> it seems like the routine PETSC_VIEWER_DRAW_(comm)  is not available
>> > >> in Fortran, while PETSC_VIEWER_STDOUT_(comm) is.
>> > >> Is there anything I can do about it? (I am using the current maint
>> > >> branch of PETSc)
>> > >>
>> > >> Thank you very much and kind regards,
>> > >>
>> > >> Volker
>> > >
>> >
>>
>>
>