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)
