On Thu, Dec 8, 2011 at 7:51 AM, Thomas Witkowski <span dir="ltr"><<a href="mailto:Thomas.Witkowski@tu-dresden.de">Thomas.Witkowski@tu-dresden.de</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Maybe some more note on my problem: I've thounds of matrices in my code, where doing this works fine! The only difference here is, that row 1129 of the matrix is owned by rank 0, but all of its entries are created on rank 1. For all of the other matrices in my codes, at least the diagonal element is created on the rank that owns the matrix row.<br>
</blockquote><div><br></div><div>Did you remember to call MatAssemblyBegin/End()?</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Thomas<br>
<br>
Zitat von Thomas Witkowski <<a href="mailto:Thomas.Witkowski@tu-dresden.de" target="_blank">Thomas.Witkowski@tu-dresden.<u></u>de</a>>:<div class="HOEnZb"><div class="h5"><br>
<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I've some strange problem: I create a mpiaij matrix, fill it with some<br>
values and want to make use of it in kspsolve. But I get the error<br>
<br>
[15]PETSC ERROR: Object is in wrong state!<br>
[15]PETSC ERROR: Matrix is missing diagonal entry 1129!<br>
<br>
I've checked the following things:<br>
- There is a MatSetValues call that contains the entry 1129/1129 with<br>
some value<br>
- MatView shows that there exists the value 1129/1129<br>
- MatAssemblyBegin/End is called correctly<br>
- valgrind says: 0 errors<br>
- ddt also cannot find anything unusual<br>
<br>
I use petsc 3.2-p5, its compiled in debug mode. What could I check to<br>
solve this problem? Any hints?<br>
<br>
Thomas<br>
</blockquote>
<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>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<br>