[petsc-dev] petsc-dev fails test with debugging=0

[Jed Brown]

Karl Rupp <rupp at mcs.anl.gov> writes:
> It depends on the amount of stability one wants to put it. If Gauss 
> without pivoting is sufficient, Block-Jacobi shouldn't be too hard for 
> the block size of interest (9-by-9). 

This is called pbjacobi.  I think that is a better choice to put on the
GPU.

> However, I suppose we do want robust solves, so it's better to use
> some explicit solution formulas.  Cramer's rule is okay up to 4-by-4
> blocks or maybe 5-by-5, but likely to be a mess for 9-by-9.

PBJacobi inverts all the diagonals.  I suggest just doing that on the
CPU and then transfer the inverses to the GPU.  You'll need a custom
kernel on the GPU.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140119/c9600e5b/attachment.sig>