[petsc-dev] Questions about MatCreateSeqBAIJWithArrays

[Travis Austin]

Hi all,

I want to use MatCreateSeqBAIJWithArrays to set a sequential block matrix
without allocating new arrays.  I believe I understand how PETSc stores the
matrix in the block format so I can transform my current data appropriately
but I am confused by the following comments in the code (baij.c:line 3490
in most recent git version).  I've also pasted this comment below:

 "When block size is greater than 1 the matrix values must be stored using
the BAIJ storage format (see the BAIJ code to determine this)."

I then see this in the code (baij.c: line 3507) and have dealt with this
when trying to use a block size larger than 1.

  if (bs != 1) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_SUP,"block size %D > 1 is
not supported yet",bs);

I'm wondering if this is just an oversight and it really works.  Can
someone provide me an update on MatCreateSeqBAIJWithArrays?  If it doesn't
work, I'll probably go ahead and try and make it work.

Thanks,
Travis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20140404/844897c0/attachment.html>