[petsc-users] Question about KSP
Elias Karabelas
elias.karabelas at medunigraz.at
Thu May 19 07:43:43 CDT 2016
Well I don't bother to much about which PC is use. I just wanted to look
what's possible. One can of course come up with some kind of augmented
Lagrange formulation, but then you end up with a saddle point problem.
On 19.05.2016 14:41, Matthew Knepley wrote:
> On Thu, May 19, 2016 at 7:39 AM, Elias Karabelas
> <elias.karabelas at medunigraz.at <mailto:elias.karabelas at medunigraz.at>>
> wrote:
>
> Ok maybe I should go a little bit more into detail. I have some
> stokes problem in a bifurcating Y-tube. Now I would like to
> enforce some prescribed flux splits at the outlets. After some
> literature research i found a Nitsche method to do that. And if
> you consider the additional billinear forms you end up with some
> nonlocal coupling terms at the outlets (which are those vectors I
> talked earlier)
>
>
> Keep the Cc.
>
> I can see that you might have that formulation, but how you would
> incorporate that into AMG, I have no idea.
>
> Matt
>
>
> On 19.05.2016 14:31, Matthew Knepley wrote:
>> On Thu, May 19, 2016 at 7:29 AM, Elias Karabelas
>> <elias.karabelas at medunigraz.at
>> <mailto:elias.karabelas at medunigraz.at>> wrote:
>>
>> maybe something with PCSHELL?
>>
>>
>> No, I mean that I do not understand what algorithm you might use
>> (apart from the implementation).
>>
>> Matt
>>
>>
>> On 19.05.2016 14:28, Matthew Knepley wrote:
>>> On Thu, May 19, 2016 at 2:48 AM, Elias Karabelas
>>> <elias.karabelas at medunigraz.at
>>> <mailto:elias.karabelas at medunigraz.at>> wrote:
>>>
>>> Dear all,
>>>
>>> I have a question about preconditioned solvers. So, I
>>> have a Sparsematrix, say A, and now for some reason I
>>> would like to add some rank-one term u v^T to that matrix.
>>> As far as I know about Petsc, I can define the action of
>>> this matrix with MatShell. But is it possible to adapt
>>> the preconditioner (like a AMG) of my KSP to handle that
>>> kind of MatShell?
>>>
>>>
>>> I do not know how you would do that. Its possible other
>>> people have ideas.
>>>
>>> Thanks,
>>>
>>> Matt
>>>
>>> Kind Regards
>>> Elias
>>>
>>> --
>>> Dr Elias Karabelas
>>>
>>> Medical University of Graz
>>> Institute of Biophysics
>>> Harrachgasse 21/IV
>>> 8010 Graz, Austria
>>>
>>> Phone: +43 316 380 7759 <tel:%2B43%20316%20380%207759>
>>> Email: elias.karabelas at medunigraz.at
>>> <mailto:elias.karabelas at medunigraz.at>
>>> Web :
>>> http://forschung.medunigraz.at/fodok/staff?name=EliasKarabelas
>>>
>>>
>>>
>>>
>>> --
>>> 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
>>
>> --
>> Dr Elias Karabelas
>>
>> Medical University of Graz
>> Institute of Biophysics
>> Harrachgasse 21/IV
>> 8010 Graz, Austria
>>
>> Phone:+43 316 380 7759 <tel:%2B43%20316%20380%207759>
>> Email:elias.karabelas at medunigraz.at
>> <mailto:elias.karabelas at medunigraz.at>
>> Web :http://forschung.medunigraz.at/fodok/staff?name=EliasKarabelas
>>
>>
>>
>>
>> --
>> 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
>
> --
> Dr Elias Karabelas
>
> Medical University of Graz
> Institute of Biophysics
> Harrachgasse 21/IV
> 8010 Graz, Austria
>
> Phone:+43 316 380 7759 <tel:%2B43%20316%20380%207759>
> Email:elias.karabelas at medunigraz.at <mailto:elias.karabelas at medunigraz.at>
> Web :http://forschung.medunigraz.at/fodok/staff?name=EliasKarabelas
>
>
>
>
> --
> 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
--
Dr Elias Karabelas
Medical University of Graz
Institute of Biophysics
Harrachgasse 21/IV
8010 Graz, Austria
Phone: +43 316 380 7759
Email: elias.karabelas at medunigraz.at
Web : http://forschung.medunigraz.at/fodok/staff?name=EliasKarabelas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160519/dd4938c6/attachment-0001.html>
More information about the petsc-users
mailing list