On Wed, Sep 3, 2014 at 9:46 AM, Klaij, Christiaan <[email protected]> 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 <[email protected]> > *Sent:* Wednesday, September 03, 2014 4:44 PM > *To:* Klaij, Christiaan > *Cc:* [email protected] > *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran > > On Wed, Sep 3, 2014 at 9:43 AM, Klaij, Christiaan <[email protected]> > 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 <[email protected]> >> *Sent:* Wednesday, September 03, 2014 4:38 PM >> *To:* Klaij, Christiaan >> *Cc:* [email protected] >> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran >> >> On Wed, Sep 3, 2014 at 9:23 AM, Klaij, Christiaan <[email protected]> >> 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 <[email protected]> >>> *Sent:* Wednesday, September 03, 2014 2:12 PM >>> *To:* Klaij, Christiaan >>> *Cc:* [email protected] >>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran >>> >>> On Wed, Sep 3, 2014 at 7:00 AM, Klaij, Christiaan <[email protected]> >>> 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 <[email protected]> >>>> *Sent:* Tuesday, September 02, 2014 5:42 PM >>>> *To:* Klaij, Christiaan >>>> *Cc:* [email protected] >>>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran >>>> >>>> On Tue, Sep 2, 2014 at 2:08 AM, Klaij, Christiaan <[email protected]> >>>> 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 >>>>> >>>>> dr. ir. Christiaan Klaij >>>>> >>>>> CFD Researcher >>>>> Research & Development >>>>> >>>>> >>>>> >>>>> *MARIN* >>>>> >>>>> >>>>> 2, Haagsteeg E [email protected] P.O. Box 28 T +31 317 49 39 11 6700 >>>>> AA Wageningen F +31 317 49 32 45 T +31 317 49 33 44 The Netherlands I >>>>> 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:* [email protected] >>>>> *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 <[email protected]> >>>>> *Sent:* Friday, August 29, 2014 4:14 PM >>>>> *To:* Klaij, Christiaan >>>>> *Cc:* [email protected] >>>>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran >>>>> >>>>> On Fri, Aug 29, 2014 at 8:55 AM, Klaij, Christiaan <[email protected] >>>>> > 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:[email protected] >>>>>> T +31 317 49 33 44 >>>>>> >>>>>> >>>>>> MARIN >>>>>> 2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands >>>>>> T +31 317 49 39 11, F +31 317 49 32 45, I 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
