Hei, thanks for the idea, that solved my problem! I noticed that I was writing out of bounds in an earlier loop, which then resulted in the behavior observed below.
Regards, Roland Am 17.01.21 um 19:40 schrieb Stefano Zampini: > Valgrind https://valgrind.org/ <https://valgrind.org/> > > Il Dom 17 Gen 2021, 21:28 Roland Richter <[email protected] > <mailto:[email protected]>> ha scritto: > > Dear all, > > I am currently encountering issues in my program when calling > MatDiagonalScale(matrix, NULL, vector) in an MPI-context for at > least two threads. The backtrace from gdb (and the error output) > tells me that the problem lies apparently in > PetscSFBcastAndOpBegin(), with the exact error listed as > /petsc/src/vec/is/sf/interface/sf.c Wrong type of object: > Parameter # 1 > / > > When printing the involved matrix and vector, everything looks as > expected (the matrix is split over the processes equally with > continuous rows, the vector is distributed over both processes). > > Is the issue that my vector is distributed over several processes? > Or is there another issue I am not seeing here? > > I tried to narrow the problem down by writing a MWE, but was not > able to reproduce it in a smaller scale. Therefore, are there > other approaches I could try, too? > > Thanks! > > Regards, > > Roland Richter >
