On Sat, Dec 17, 2011 at 4:55 PM, Juha Jäykkä <span dir="ltr"><<a href="mailto:juhaj@iki.fi">juhaj@iki.fi</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">
> This is memory overwriting. Use valgrind to find it.<br>
<br>
Valgrind gives it a clean bill of health:<br>
<br>
==30174== ERROR SUMMARY: 0 errors from 0 contexts (suppressed: 153 from 1)<br>
==30174== malloc/free: in use at exit: 0 bytes in 0 blocks.<br>
==30174== malloc/free: 0 allocs, 0 frees, 0 bytes allocated.<br>
==30174== For counts of detected errors, rerun with: -v<br>
==30174== All heap blocks were freed -- no leaks are possible.<br>
<br>
What am I actually even looking for? Myself overwriting stuff or a PETSc bug?<br>
I do not quite understand how PETSc even could manage to overwrite unrelated<br>
vectors. Or do you suspect I overwrite the vector structure at somepoint and<br>
that confuses PETSc later on?<br></blockquote><div><br></div><div>I guarantee you this is not a bug. This stuff is used by thousands of people every day.</div><div>Its very simple and old code. There is a misunderstanding somewhere in your code</div>
<div>about how this mechanism works. I suggest stripping down the code until there are</div><div>just two vectors and you alter one. If you send me that, I will go through it to figure</div><div>out what is happening.</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">
Cheers,<br>
<div class="HOEnZb"><div class="h5">Juha<br>
<br>
--<br>
-----------------------------------------------<br>
| Juha Jäykkä, <a href="mailto:juhaj@iki.fi">juhaj@iki.fi</a> |<br>
| <a href="http://www.maths.leeds.ac.uk/~juhaj" target="_blank">http://www.maths.leeds.ac.uk/~juhaj</a> |<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>