[aspect-devel] ASPECT with PETSc

[Wolfgang Bangerth]

On 03/11/2012 08:20 AM, Timo Heister 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.

OK, committed.


> 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, Matt and I talked about how this works and what would need to 
change to have access to all matrix entries. We agreed that a first step 
would be to make the current scheme work with PETSc.

Cheers
  W.

------------------------------------------------------------------------
Wolfgang Bangerth               email:            bangerth at math.tamu.edu
                                 www: http://www.math.tamu.edu/~bangerth/