[petsc-users] PCFieldSplitSetSchurPre in fortran

Wed Sep 3 07:12:00 CDT 2014

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