Banded Tridiagonal Matrices in PETSc

[Tim Stitt]

Thanks Barry, that is all I need to know.

I can stick with my current implementation then...great stuff.

Barry Smith wrote:
>
>    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)
>>
>


-- 
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)