On Thu, Jan 30, 2020 at 6:04 PM Jeremy Theler <[email protected]> wrote:
>
> On Thu, 2020-01-30 at 21:10 +0000, Smith, Barry F. wrote:
>
> > MatGetSubMatrix() and then do the product on the sub matrix then
> > VecSum
> >
>
> Ok, I have it working in a single-processor and throws the expected
> value. Yet I have a segfault in parallel when I ask for the IS
> corresponding to the rows. When I call this instruction in parallel I
> get a segfault (I can post the full debug output if needed):
>
> ISCreateStride(PETSC_COMM_WORLD, size_local, first_row, 1, &set_cols);
>
> I tried also
>
> ISCreateStride(PETSC_COMM_WORLD, size_global, 0, 1, &set_cols));
>
> but it also fails with the same segfault.
>
> What am I getting wrong?
>
Show the entire error, including stack.
Thanks,
Matt
> --
> jeremy
>
>
>
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener
https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/>
