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
