[petsc-dev] print statements in fortran interface functions

Jacob Faibussowitsch jacob.fai at gmail.com
Mon Mar 16 20:50:02 CDT 2020


Hello Mark,

If you are using PetscInfo to print, PetscInfo requires that you pass “-info” as a command line flag in order to print anything. Sounds simple enough to miss, did you run your code with it? i.e. ./yourCode -info? 

For what its worth, PetscPrintf on PETSC_COMM_WORLD will always print but only from proc 0. PetscSynchronizedPrintf prints from all in order. Might be more useful for debugging than PetscInfo!

Best regards,

Jacob Faibussowitsch
(Jacob Fai - booss - oh - vitch)
Cell: (312) 694-3391

> On Mar 16, 2020, at 7:54 PM, Mark Adams <mfadams at lbl.gov> wrote:
> 
> I am trying to debug an application code that works with v3.7 but fails with master. The code works for "normal" solvers but for a solver that uses FieldSplit it fails. It looks like vectors are not getting created from MatCreateVecs with a matrix that is a MatNest (I can't run the code). 
> 
> I have put print statements in the interface files to try to see what is going on but the print statements seem to disappear somehow. (I added the same print statements in my installation of v3.7 to see what is different) I was using PetscPrint and then found and example of a Fortran PetscInfo and switched to that. I stepped through the code with DDT and see the print statements but nothing goes to the output. I did this DDT stuff with ex73f90t, which has a MatCreateNest.
> 
> I am completely baffled. Does anyone have any ideas?
> 
> Thanks,
> Mark

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20200316/88499e84/attachment.html>


More information about the petsc-dev mailing list