[petsc-dev] how come PETSc doesn't contain the SPIKE algorithm?

[Hong Zhang]

Hey, SPIKE_algorithm is a block version of the algorithm PPT and PPD
we published almost 20 years ago
http://www.cs.iit.edu/~scs/psfiles/efficient92.pdf
I did this work at graduate school; the paper took 3+ years to appear ...
Surprisingly, the SPIKE_algorithm uses the same notations as we did in
the paper, V's, W's ...

4-5 years ago, a Columbia physic graduate student implemented PPT using PETSc.
I helped him improved it (see attached report). We were at the point
to move to the block version for its optimization.
Then I was drafted to other projects, and the student turned to physic
experiments. The work was
interrupted. I still keep the codes of this work.
I would love to resume it and get it included in petsc library.

Hong


On Wed, Dec 7, 2011 at 8:17 PM, Dave Nystrom
<Dave.Nystrom at tachyonlogic.com> wrote:
> Barry Smith writes:
>  > http://en.wikipedia.org/wiki/SPIKE_algorithm
>
> Looks great to me.  I have a bunch of banded systems in my code to solve.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FastPoissonSolver.pdf
Type: application/pdf
Size: 459278 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20111207/b2b65ab2/attachment.pdf>