[petsc-dev] ASM request

Jed Brown jed at 59A2.org
Tue Mar 16 13:26:27 CDT 2010

On Tue, 16 Mar 2010 12:51:53 -0500, Barry Smith <bsmith at mcs.anl.gov> wrote:
> If you use ILU(0) and the inplace option then it only needs one copy  
> of the blocks instead of the two.

With block Jacobi, it would be possible to avoid copying out the
sub-blocks by assembling into a matrix with the blocks already defined
(I think that Dmitry's MatFwk will allow this).  Similarly for
Neumann-Neumann with many blocks per process.

It would also be possible to assemble matrices with overlapping
subdomains, and define MatMult suitably (dropping redundant dofs).  I
don't know whether Dmitry has plans to support this.


More information about the petsc-dev mailing list