<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Apr 4, 2014 at 12:38 PM, Travis Austin <span dir="ltr"><<a href="mailto:austint73@gmail.com" target="_blank">austint73@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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></blockquote><div><br></div><div>It looks like whoever coded it only tested this with bs = 1. Feel free to take that out and run. Even better would be</div><div>to write a little test as a PETSc example. That would make us pus hit into master faster.</div>
<div><br></div><div> Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div>Thanks,</div><div>Travis</div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>