I did a make allclean, make and install. That solved the two problems. So I'm guessing these values are taken from the petscksp module or its predecessors instead of being taken directly from the finclude/petscpc.h file. Case closed.
Chris ________________________________ From: Matthew Knepley <knep...@gmail.com> Sent: Wednesday, September 03, 2014 5:02 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Wed, Sep 3, 2014 at 9:46 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: print *, PC_FIELDSPLIT_SCHUR_PRE_USER gives: 2. (changing to *_SELFP gives: This name does not have a type, and must have an explicit type. [PC_FIELDSPLIT_SCHUR_PRE_SELFP]) I have next:/PETSc3/petsc/petsc-pylith$ find include/finclude/ -type f | xargs grep USER find include/finclude/ -type f | xargs grep USER include/finclude//petscerrordef.h:#define PETSC_ERR_USER 83 include/finclude//petscpc.h: PetscEnum PC_FIELDSPLIT_SCHUR_PRE_USER include/finclude//petscpc.h: parameter (PC_FIELDSPLIT_SCHUR_PRE_USER=3) include/finclude//petsctao.h: PetscEnum TAO_CONVERGED_USER include/finclude//petsctao.h: PetscEnum TAO_DIVERGED_USER include/finclude//petsctao.h: parameter ( TAO_CONVERGED_USER = 8) include/finclude//petsctao.h: parameter ( TAO_DIVERGED_USER = -8) Where could it possibly be getting this value from? I think from some compiled module which you need to rebuild. Matt MARIN news: MARIN Report 112: Hydro-structural<http://www.marin.nl/web/News/News-items/MARIN-Report-112-Hydrostructural.htm> This e-mail may be confidential, privileged and/or protected by copyright. If you are not the intended recipient, you should return it to the sender immediately and delete your copy from your system. ________________________________ From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Sent: Wednesday, September 03, 2014 4:44 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Wed, Sep 3, 2014 at 9:43 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: I'm sorry, how do I do that? print it Matt Chris MARIN news: Applied Hydrodynamics of Floating Offshore Structures course, Oct 8 - 10, Houston<http://www.marin.nl/web/News/News-items/Applied-Hydrodynamics-of-Floating-Offshore-Structures-course-Oct-8-10-Houston.htm> This e-mail may be confidential, privileged and/or protected by copyright. If you are not the intended recipient, you should return it to the sender immediately and delete your copy from your system. ________________________________ From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Sent: Wednesday, September 03, 2014 4:38 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Wed, Sep 3, 2014 at 9:23 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: Matt, Thanks, after applying the fix to my petsc-3.5.1 install, the small Fortran program works as expected. Now, I would like to change the fortran strategy to the option "3) Using Fortran modules". So, in the small fortran program I replace these seven lines #include <finclude/petscsys.h> #include <finclude/petscis.h> #include <finclude/petscvec.h> #include <finclude/petscmat.h> #include <finclude/petscpc.h> #include <finclude/petscksp.h> #include <finclude/petscviewer.h> by the two following lines use petscksp #include <finclude/petsckspdef.> This still compiles but I get the two old problem back... I have no idea why. Can you look at the values of those enumerations? Matt Chris ________________________________ From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Sent: Wednesday, September 03, 2014 2:12 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Wed, Sep 3, 2014 at 7:00 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: Matt, Thanks for the fix. If I understand correctly, in an existing install of petsc-3.5.1, I would only need to replace the file "finclude/petscpc.h" by the new file for the fix to work? (instead of downloading dev, configuring, installing on various machines). Yes Matt Chris MARIN news: Bas Buchner speaker at Lowpex conference at SMM Hamburg<http://www.marin.nl/web/News/News-items/Bas-Buchner-speaker-at-Lowpex-conference-at-SMM-Hamburg.htm> This e-mail may be confidential, privileged and/or protected by copyright. If you are not the intended recipient, you should return it to the sender immediately and delete your copy from your system. ________________________________ From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Sent: Tuesday, September 02, 2014 5:42 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Tue, Sep 2, 2014 at 2:08 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: Matt, Attached is a small Fortran code that replicates the second problem. This was a Fortran define problem. I fixed it on next https://bitbucket.org/petsc/petsc/branch/knepley/fix-pc-fieldsplit-fortran and it will be in maint and master tomorrow. Thanks, Matt Chris [cid:image76463f.JPG@e74baedb.4a899186][cid:imageae8e3c.JPG@1f0a6dd3.4e92756a] dr. ir. Christiaan Klaij CFD Researcher Research & Development MARIN 2, Haagsteeg E c.kl...@marin.nl<mailto:c.kl...@marin.nl> P.O. Box 28 T +31 317 49 39 11<tel:%2B31%20317%2049%2039%2011> 6700 AA Wageningen F +31 317 49 32 45<tel:%2B31%20317%2049%2032%2045> T +31 317 49 33 44<tel:%2B31%20317%2049%2033%2044> The Netherlands I www.marin.nl<http://www.marin.nl> MARIN news: MARIN at SMM, Hamburg, September 9-12<http://www.marin.nl/web/News/News-items/MARIN-at-SMM-Hamburg-September-912.htm> This e-mail may be confidential, privileged and/or protected by copyright. If you are not the intended recipient, you should return it to the sender immediately and delete your copy from your system. ________________________________ From: Klaij, Christiaan Sent: Friday, August 29, 2014 4:42 PM To: Matthew Knepley Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: RE: [petsc-users] PCFieldSplitSetSchurPre in fortran Matt, The small test code (ex70) is in C and it works fine, the problem happens in a big Fortran code. I will try to replicate the problem in a small Fortran code, but that will take some time. Chris ________________________________ From: Matthew Knepley <knep...@gmail.com<mailto:knep...@gmail.com>> Sent: Friday, August 29, 2014 4:14 PM To: Klaij, Christiaan Cc: petsc-users@mcs.anl.gov<mailto:petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PCFieldSplitSetSchurPre in fortran On Fri, Aug 29, 2014 at 8:55 AM, Klaij, Christiaan <c.kl...@marin.nl<mailto:c.kl...@marin.nl>> wrote: I'm trying PCFieldSplitSetSchurPre with PC_FIELDSPLIT_SCHUR_PRE_SELFP in petsc-3.5.1 using fortran. The first problem is that PC_FIELDSPLIT_SCHUR_PRE_SELFP seems to be missing in fortran, I get the compile error: This name does not have a type, and must have an explicit type. [PC_FIELDSPLIT_SCHUR_PRE_SELFP] while compilation works fine with _A11, _SELF and _USER. Mark Adams has just fixed this. The second problem is that the call doesn't seem to have any effect. For example, I have CALL PCFieldSplitSetSchurPre(pc,PC_FIELDSPLIT_SCHUR_PRE_USER,aa,ierr) CALL PCFieldSplitSetSchurFactType(pc,PC_FIELDSPLIT_SCHUR_FACT_LOWER,ierr) This compiles and runs, but ksp_view tells me PC Object:(sys_) 3 MPI processes type: fieldsplit FieldSplit with Schur preconditioner, factorization LOWER Preconditioner for the Schur complement formed from A11 So changing the factorization from the default FULL to LOWER did work, but changing the preconditioner from A11 to USER didn't. I've also tried to run directly from the command line using -sys_pc_fieldsplit_schur_precondition user -sys_ksp_view This works in the sense that I don't get the "WARNING! There are options you set that were not used!" message, but still ksp_view reports A11 instead of user provided matrix. Can you send a small test code, since I use this everyday here and it works. Thanks, Matt Chris dr. ir. Christiaan Klaij CFD Researcher Research & Development E mailto:c.kl...@marin.nl<mailto:c.kl...@marin.nl> T +31 317 49 33 44<tel:%2B31%20317%2049%2033%2044> MARIN 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands T +31 317 49 39 11<tel:%2B31%20317%2049%2039%2011>, F +31 317 49 32 45<tel:%2B31%20317%2049%2032%2045>, I www.marin.nl<http://www.marin.nl> -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener -- What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead. -- Norbert Wiener