[petsc-users] How to use the matrix format of block compressed row (MATBAIJ)

[Li, Zhisong (lizs)]

 Hi, Petsc Team,

The PETSC user manual says that the matrix format of block compressed row (BAIJ) will be more efficient than the default AIJ format for problems with multiple DOF. But the description of using BAIJ matrix is very limited in the manual. When trying to replace AIJ with BAIJ in the code of iterative solver, I often meet the errors " No support for this operation for this object type!", " Matrix format mpibaij does not have a built-in PETSc direct solver! ". I think something must be missing here when we change the format.

And how to prelocate memory for BAIJ matrices when we have multiple DOF? Where can we find detailed statements or example codes of using BAIJ with preconditioners in iterative computation?

BTW, When we use DA to organize the data for parallel computing, is it more natural to use DAGetMatrix() to create the matrix than using MatCreateMPIAIJ() ?


Thank you.


Zhisong Li

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20100628/05b6cac8/attachment.htm>