On Sun, Jul 22, 2012 at 9:22 PM, Juha Jäykkä <span dir="ltr"><<a href="mailto:juhaj@iki.fi" target="_blank">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">
Hi list!<br>
<br>
Petsc3.2-p7, make test fails and at a closer look,<br>
src/snes/examples/tutorials/ex19 gives:<br></blockquote><div><br></div><div>Is this complex? There was a problem with the complex dot product for some</div><div>BLAS. If this is the problem, you can reconfigure with --download-f-blas-lapack.</div>
<div>Also, upgrading should fix it.</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">
orterun -n 1 ./ex19 -dmmg_nlevels 4 -snes_monitor_short -<br>
on_error_attach_debugger<br>
lid velocity = 0.0016, prandtl # = 1, grashof # = 1<br>
[0]PETSC ERROR:<br>
------------------------------------------------------------------------<br>
[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably<br>
memory access out of range<br>
[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br>
[0]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/petsc-" target="_blank">http://www.mcs.anl.gov/petsc/petsc-</a><br>
as/documentation/faq.html#valgrind[0]PETSC ERROR: or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a><br>
on GNU/linux and Apple Mac OS X to find memory corruption errors<br>
[0]PETSC ERROR: likely location of problem given in stack below<br>
[0]PETSC ERROR: --------------------- Stack Frames<br>
------------------------------------<br>
[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br>
[0]PETSC ERROR: INSTEAD the line number of the start of the function<br>
[0]PETSC ERROR: is given.<br>
[0]PETSC ERROR: [0] VecNorm_Seq line 236 src/vec/vec/impls/seq/bvec2.c<br>
[0]PETSC ERROR: [0] VecNormBegin line 479 src/vec/vec/utils/comb.c<br>
[0]PETSC ERROR: [0] SNESSolve_LS line 142 src/snes/impls/ls/ls.c<br>
[0]PETSC ERROR: [0] SNESSolve line 2647 src/snes/interface/snes.c<br>
[0]PETSC ERROR: [0] DMMGSolveSNES line 538 src/snes/utils/damgsnes.c<br>
[0]PETSC ERROR: [0] DMMGSolve line 303 src/snes/utils/damg.c<br>
[0]PETSC ERROR: User provided function() line 0 in unknown directory unknown<br>
file<br>
--------------------------------------------------------------------------<br>
<br>
Any idea where this could come from? Worst of all, if I run this in debugger,<br>
it never sigsegvs, likewise under valgrind. Sometimes (about 1/3) it also<br>
succeeds without a debugger, but even when not segfaulting, it gives wrong<br>
results (so I assume not segfaulting is just good (or bad?) luck):<br>
<br>
2,5c2,3<br>
< 0 SNES Function norm 0.0406612<br>
< 1 SNES Function norm 3.33175e-06<br>
< 2 SNES Function norm 1.092e-11<br>
< Number of Newton iterations = 2<br>
---<br>
> 0 SNES Function norm < 1.e-11<br>
> Number of Newton iterations = 0<br>
7,10c5,6<br>
< 0 SNES Function norm 0.0406612<br>
< 1 SNES Function norm 3.33175e-06<br>
< 2 SNES Function norm 1.092e-11<br>
< Number of Newton iterations = 2<br>
---<br>
> 0 SNES Function norm < 1.e-11<br>
> Number of Newton iterations = 0<br>
<br>
Any help appreciated, thanks.<br>
<span class="HOEnZb"><font color="#888888">Juha<br>
<br>
--<br>
-----------------------------------------------<br>
| Juha Jäykkä, <a href="mailto:juhaj@iki.fi">juhaj@iki.fi</a> |<br>
| <a href="http://koti.kapsi.fi/~juhaj/" target="_blank">http://koti.kapsi.fi/~juhaj/</a> |<br>
-----------------------------------------------<br>
</font></span></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>