[petsc-users] PCFieldSplitSetSchurPre in fortran

Matthew Knepley knepley at gmail.com
Wed Sep 3 09:44:40 CDT 2014 

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140903/64a1fb37/attachment.html>
-------------- 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/64a1fb37/attachment.jpe>
-------------- 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/64a1fb37/attachment-0001.jpe>

More information about the petsc-users mailing list