<div dir="ltr">Hi all,<div><br></div><div>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:</div>
<div><br></div><div><div> "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)."</div></div><div><br></div><div>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.</div>
<div><br></div><div>  if (bs != 1) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_SUP,"block size %D > 1 is not supported yet",bs);<br></div><div><br></div><div>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.</div>
<div><br></div><div>Thanks,</div><div>Travis</div></div>