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
