<div dir="ltr"><div dir="ltr">On Thu, Jan 30, 2020 at 6:04 PM Jeremy Theler <<a href="mailto:jeremy@seamplex.com">jeremy@seamplex.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><br>
On Thu, 2020-01-30 at 21:10 +0000, Smith, Barry F. wrote:<br>
<br>
>   MatGetSubMatrix() and then do the product on the sub matrix then<br>
> VecSum<br>
> <br>
<br>
Ok, I have it working in a single-processor and throws the expected<br>
value. Yet I have a segfault in parallel when I ask for the IS<br>
corresponding to the rows. When I call this instruction in parallel I<br>
get a segfault (I can post the full debug output if needed):<br>
<br>
ISCreateStride(PETSC_COMM_WORLD, size_local, first_row, 1, &set_cols);<br>
<br>
I tried also<br>
<br>
ISCreateStride(PETSC_COMM_WORLD, size_global, 0, 1, &set_cols));<br>
<br>
but it also fails with the same segfault.<br>
<br>
What am I getting wrong?<br></blockquote><div><br></div><div>Show the entire error, including stack.</div><div><br></div><div>  Thanks,</div><div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
--<br>
jeremy<br>
<br>
<br>
<br>
<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>