[petsc-users] Product of matrix row times a vector

[Matthew Knepley]

On Thu, Jan 30, 2020 at 6:04 PM Jeremy Theler <jeremy at seamplex.com> 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/>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200130/f872d41c/attachment.html>