[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Schur complement system
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Schur complement system
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Fri, 7 Mar 2008 09:25:02 -0600
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=2CjcE57VO5rMy84LmIdRcTSgEvZAe7/V7HAyd8zVuOk=; b=IuneSgT+Apa21YxjqsmFHc9cfhhNYEq1X66Ba13s+C+zz0qKT2ncIQO6tGcg5m02iASX//Yrfc6GUrwIheZ1OcODmUi+9OCaz+u971yNeFJwnoR4fTZRHdUQnEwP2DvX527B/0dpkeIcnGC4EU3c+b50BSwLLTZ58qJkRAPqS1U=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=OmdoxkfEHb+9quc2YD/1D640oJnvP0GYrkqFhgBpeKu7ywcTvzj43V7hGpR82me8squYpvuXIqYNbApGMVGrL8+l6OuC7/3AuG/f4RjKRVl1oJEURP+/biq0Y2t9AgLBAftKw1nTq+g/nml2iFEBap+GHtRM5Jy/zvTA6McYQaU=
- In-reply-to: <Pine.LNX.4.64.0803071249030.25330@xeon.ee.ethz.ch>
- References: <Pine.LNX.4.64.0803071249030.25330@xeon.ee.ethz.ch>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
Yes, you can certainly solve this system. Normally, you would just create a
simple MATSHELL to apply the Schur complement, which does two MatMult()s
and one KSPSolve(). We are integrating specific support for this, but it will
not be released until the summer I think.
Thanks,
Matt
On Fri, Mar 7, 2008 at 5:50 AM, Kathrin Burckhardt
<tribur@xxxxxxxxxxxxxxxxx> wrote:
> Dear nice people,
>
> Did you ever try to solve a Schur complement system using PETSc?
>
>
--
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