[petsc-users] PCFieldSplitSetSchurPre in fortran

[Matthew Knepley]

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