[petsc-users] Column #j is wrong in parallel from message "Inserting a new nonzero (i, j) into matrix"

[Matthew Knepley]

On Wed, Mar 25, 2015 at 1:03 PM, Eric Chamberland <
Eric.Chamberland at giref.ulaval.ca> wrote:

> Hi,
>
> while looking for where in the world do I insert the (135,9) entry in my
> matrix, I have discovered that the column # shown is wrong in parallel!
>

We have talked about this before. It is certainly not easy. The error is
given by the embedded serial matrix, whereas you want
the parallel column. It would be good if we could restructure to give that
back.

  Thanks,

      Matt


> I am using PETsc 3.5.3.
>
> The full error message is:
>
> [0]PETSC ERROR: MatSetValues_MPIAIJ() line 564 in
> /home/mefpp_ericc/petsc-3.5.3/src/mat/impls/aij/mpi/mpiaij.c Inserting a
> new nonzero (135, 9) into matrix
>
> This line code is a call to a #defined macro:
>
> MatSetValues_SeqAIJ_B_Private(row,col,value,addv);
>
> where the "col" parameter is not equal to "in[j]"!!!
>
> in gdb, printing "in[j]" gave me:
>
> print in[j]
> $6 = 537
>
> while "col" is:
>
> print col
> $7 = 9
>
> So, I expected to have a message telling me that (135,537) and not (135,9)
> is a new entry matrix!!!
>
> Would it be a big work to fix this so that the col # displayed is correct?
>
> Thanks!
>
> Eric
>



-- 
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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150325/ccb10655/attachment.html>