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>