[petsc-users] PCFieldSplitSetSchurPre in fortran
Matthew Knepley
knepley at gmail.com
Wed Sep 3 10:02:12 CDT 2014
On Wed, Sep 3, 2014 at 9:46 AM, Klaij, Christiaan <C.Klaij at 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 <knepley at gmail.com>
> *Sent:* Wednesday, September 03, 2014 4:44 PM
> *To:* Klaij, Christiaan
> *Cc:* petsc-users at mcs.anl.gov
> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
>
> On Wed, Sep 3, 2014 at 9:43 AM, Klaij, Christiaan <C.Klaij at 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 <knepley at gmail.com>
>> *Sent:* Wednesday, September 03, 2014 4:38 PM
>> *To:* Klaij, Christiaan
>> *Cc:* petsc-users at mcs.anl.gov
>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
>>
>> On Wed, Sep 3, 2014 at 9:23 AM, Klaij, Christiaan <C.Klaij at 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 <knepley at gmail.com>
>>> *Sent:* Wednesday, September 03, 2014 2:12 PM
>>> *To:* Klaij, Christiaan
>>> *Cc:* petsc-users at mcs.anl.gov
>>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
>>>
>>> On Wed, Sep 3, 2014 at 7:00 AM, Klaij, Christiaan <C.Klaij at 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 <knepley at gmail.com>
>>>> *Sent:* Tuesday, September 02, 2014 5:42 PM
>>>> *To:* Klaij, Christiaan
>>>> *Cc:* petsc-users at mcs.anl.gov
>>>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
>>>>
>>>> On Tue, Sep 2, 2014 at 2:08 AM, Klaij, Christiaan <C.Klaij at 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
>>>>>
>>>>> dr. ir. Christiaan Klaij
>>>>>
>>>>> CFD Researcher
>>>>> Research & Development
>>>>>
>>>>>
>>>>>
>>>>> *MARIN*
>>>>>
>>>>>
>>>>> 2, Haagsteeg E C.Klaij at marin.nl 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:* petsc-users at 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 <knepley at gmail.com>
>>>>> *Sent:* Friday, August 29, 2014 4:14 PM
>>>>> *To:* Klaij, Christiaan
>>>>> *Cc:* petsc-users at mcs.anl.gov
>>>>> *Subject:* Re: [petsc-users] PCFieldSplitSetSchurPre in fortran
>>>>>
>>>>> On Fri, Aug 29, 2014 at 8:55 AM, Klaij, Christiaan <C.Klaij at 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.Klaij at marin.nl
>>>>>> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140903/2d6c7eb0/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image76463f.JPG
Type: image/jpeg
Size: 1069 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140903/2d6c7eb0/attachment-0002.jpe>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: imageae8e3c.JPG
Type: image/jpeg
Size: 1622 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140903/2d6c7eb0/attachment-0003.jpe>
More information about the petsc-users
mailing list