I don’t know if this is related, but the following example triggers a
[0]PETSC ERROR: Argument out of range
[0]PETSC ERROR: New nonzero at (0,1) caused a malloc
While everything works OK if I switch from MPISBAIJ to MPIBAIJ.
Thanks,
Pierre
On 2 Aug 2019, at 1:49 AM, Smith, Barry F. <[email protected]> wrote:
Yes it is a bug, working on it now.
On Aug 1, 2019, at 9:13 AM, Pierre Jolivet via petsc-dev <[email protected]> wrote:
Hello,
The attached example is a little confusing for me.
How come I don’t get the same matrix out-of-the-box?
For me, the “correct” matrix is the SeqSBAIJ, how can I get MatMPISBAIJSetPreallocationCSR to assemble the same matrix?
Do I have to resort to first assembling SeqSBAIJ matrices and then calling MatCreateMPIMatConcatenateSeqMat if I don’t want to play around with my input arrays?
Thanks,
Pierre
<foo.c>
ex232.c