MatFeti

Matthew Knepley knepley at gmail.com
Tue Sep 22 13:22:50 CDT 2009


On Tue, Sep 22, 2009 at 1:19 PM, Dima Karpeyev <karpeev at gmail.com> wrote:

> It is a preconditioner, rather than a matrix class, isn't it?
>

You have to go back to the Mat level, because there is no way to
"unassemble" a
Mat. This is the same old bullshit we deal with every time we try to move
beyond
a single linear algebra problem. Somehow, discretization needs to come in.

  Matt


> Dmitry.
>
>
> On Tue, Sep 22, 2009 at 1:17 PM, Matthew Knepley <knepley at gmail.com>wrote:
>
>> That is Oliver's old code. It is really complicated. I believe even Oliver
>> thinks that any new implementation should start from scratch, using his
>> current PETSc code as a reference.
>>
>>   Matt
>>
>>
>> On Tue, Sep 22, 2009 at 1:15 PM, Dmitry Karpeev <karpeev at mcs.anl.gov>wrote:
>>
>>> I've found something that looks like somewhat outdated code (e.g., it
>>> uses SLESxxx) for FETI-DP:
>>> src/mat/impls/feti
>>>
>>> It appears that this matrix class actually implements a FETI-DP based
>>> preconditioner.
>>> Is that right?
>>> Perhaps some (or all) of it can be rolled into a more general
>>> PCFieldSplit-type framework?
>>> Any ideas?
>>>
>>> Thanks.
>>> Dmitry.
>>>
>>
>>
>>
>> --
>> 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-dev/attachments/20090922/21c5561c/attachment.html>


More information about the petsc-dev mailing list