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
