AIJ and BAIJ convertion

Satish Balay balay at
Tue May 5 11:09:02 CDT 2009

On Tue, 5 May 2009, Hong Zhang wrote:

> On Tue, 5 May 2009, SUN Chun wrote:
> > Thanks Hong and Satish,
> > 
> > Unfortunately I did use MatCreateMPIAIJ. It's difficult to change that part
> > of my code.
> > Plus I have matrices dumped out with AIJ format and I want to read it as
> > BAIJ. It seems that I have no option other than MatCreateMPIBAIJ then
> > MatSetBlockSize then add entries one by one...?

> You can call
> ierr = MatLoad(fd,MATBAIJ,&newbaijmat);CHKERRQ(ierr);
> and run your code with '-matload_block_size 6'.
> In this way, a new baij matrix is created with bs=6.

Also the call to MatCreateMPIAIJ() can be substituted with calls to
MatCreate(), MatSetType(MATMPIAIJ),MatSetSizes() etc..] without
changing the rest of the code..


More information about the petsc-users mailing list