[Dune] dune parallel matrix format, example

[Markus Blatt]

On Mon, Jul 18, 2016 at 10:19:56AM +0200, Markus Blatt wrote:
> 
> Now you have your complete parallel index set with a global index and
> a marker for each local index. Let l(i) be the function that
> determines the local index of the global (PETSc) index i, them your
> local matrix can be set up, if you insert an entry l(i),l(j) for each
> PETSc sparse matrix entry (i,j). For all column indices g not in {n_i,
> ..., n_{i+1}-1}, you insert a row l(g) with a 1 on the diagonal and
> zeros elsewhere. Please make sure that for a symmetric sparse matrix
> the sparsity pattern should be symetric, too.
> 

Small addition: Now you can setup a parallel linear Operator by
passing the local sparse matrix and the communication object to
OverlappingSchwarzOperator. Or even write your own one.

Markus

-- 
Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de
Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany
Tel.: +49 (0) 160 97590858
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 819 bytes
Desc: Digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20160718/7a69c46b/attachment.sig>