[aspect-devel] ASPECT with PETSc

[Matthew Knepley]

On Sun, Mar 11, 2012 at 8:20 AM, Timo Heister <heister at math.tamu.edu> wrote:

> > [Timo: Question -- do we want this anyway? If so, please apply to the
> > repository since I will be offline for the next couple of days.]
>
> I think so, yes. This was on my "10% time to be done on weekends" list.
>
> > The second part you talked about, namely not building the matrix as
> > individual blocks but instead as one matrix from which you let PETSc pick
> > and choose is likely going to be much more difficult. I'll try to explain
> > how this currently works when we meet in 20 minutes down at the AGU
> > building. For now I've gotta run :-)
>
> Matt, you want to use PETSc solvers for the whole linear system? It
> shouldn't be too hard to make this change. To explain: right now a
> PETScWrappers::BlockSparseMatrix is a two-dimensional array of PETSc
> matrices (same with vectors). Of course, we are using this block
> structure in many places.


Yes, do no blocking. Just give the whole matrix, do not hardcode the
blocking
or the matrix type, or the solver type or the preconditioner. That is
design from
the 1980s.

   Matt


>
> --
> Timo Heister
> http://www.math.tamu.edu/~heister/
>



-- 
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://geodynamics.org/pipermail/aspect-devel/attachments/20120311/7166a7f0/attachment.htm