Banded Tridiagonal Matrices in PETSc

[Barry Smith]

    There is a format Bdiag that stores by "banded diagonal". You will  
find that
this performs slower than then the AIJ format.

   If your matrix has constant values along the "diagaonals" then you  
will benefit from
using a MatShell and writing custom code. if the values along the  
diagonals are
not constant you will not do any better than AIJ anyways.

    Barry

On Nov 22, 2007, at 11:16 AM, Tim Stitt wrote:

> Hi,
>
> I was just wondering if PETSc has any special provision for banded  
> tridiagonal complex matrices when used in conjunction with  
> KSPSolve(). Are there any special PETSc matrix types or  
> factorisation/solver methods that benefit more from this matrix form?
>
> Currently I am just using standard AIJ representation in my serial/ 
> parallel codes.
>
> I would be grateful for any thoughts.
>
> Thanks,
>
> Tim.
>
> -- 
> Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie>
> HPC Application Consultant - ICHEC (www.ichec.ie)
>
> Dublin Institute for Advanced Studies
> 5 Merrion Square - Dublin 2 - Ireland
>
> +353-1-6621333 (tel) / +353-1-6621477 (fax)
>