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

Barry Smith bsmith at mcs.anl.gov
Wed Dec 7 21:23:47 CST 2011


  Time to update the Wikipedia page for spike

On Dec 7, 2011, at 8:39 PM, Hong Zhang wrote:

> 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.
> <FastPoissonSolver.pdf>




More information about the petsc-dev mailing list