On Mon, Sep 7, 2009 at 2:01 AM, Klaij, Christiaan <span dir="ltr">&lt;<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a>&gt;</span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Matt,<br>
<br>
Thanks for your answer!<br>
<br>
dummy is defined as a PETSc integer inside MyKSPConverged, so why can&#39;t I set it to 0? I know it&#39;s not useful because it&#39;s just a dummy and nothing happens with it, but setting it shouldn&#39;t crash the code, right? Also, if setting it changes the basic definition of PETSC_NULL_OBJECT shouldn&#39;t you get a SEGV as well, or at least some kind of trouble?<br>
</blockquote><div><br>dummy is something <b>you passed in</b> earlier in the code, and you passed in something that should not change. I<br>do not know why you are getting an SEGV (I do not get one), but I do know you should not change this for the<br>
code to function correctly.<br><br>  Matt<br> </div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Chris<br>
<br>
<br>
<br>
Message: 5<br>
Date: Fri, 4 Sep 2009 10:49:25 -0500<br>
From: Matthew Knepley &lt;<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>&gt;<br>
Subject: Re: dummy in MyKSPConverged<br>
To: PETSc users list &lt;<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>&gt;<br>
Message-ID:<br>
        &lt;<a href="mailto:a9f269830909040849t7a8dd52dwc7fa33c6a0d313f7@mail.gmail.com">a9f269830909040849t7a8dd52dwc7fa33c6a0d313f7@mail.gmail.com</a>&gt;<br>
Content-Type: text/plain; charset=&quot;iso-8859-1&quot;<br>
<div><div></div><div class="h5"><br>
On Fri, Sep 4, 2009 at 10:35 AM, Klaij, Christiaan &lt;<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a>&gt; wrote:<br>
<br>
&gt;   Hello,<br>
&gt;<br>
&gt; With the new version 3.0.0-p8, if I assign a value to the dummy in<br>
&gt; MyKSPConverged at the end of src/ksp/ksp/examples/tutorials/ex2f.F, I get a<br>
&gt; segmentation fault (this didn&#39;t happen with version 2.3.3-p13). Any idea<br>
&gt; why?<br>
&gt;<br>
<br>
1) I do not get a SEGV<br>
<br>
2) You should not be setting this at all. You pass in PETSC_NULL_OBJECT,<br>
which means you are changing the definition of that<br>
    basic thing in PETSc. This can produce undefined results.<br>
<br>
  Matt<br>
<br>
<br>
&gt; $ tail ex2f.F<br>
&gt;       else<br>
&gt;         flag = 0<br>
&gt;       endif<br>
&gt;       ierr = 0<br>
&gt;       dummy = 0<br>
&gt;<br>
&gt;       end<br>
&gt; $ ./ex2f -my_ksp_convergence<br>
&gt; [0]PETSC ERROR:<br>
&gt; ------------------------------------------------------------------------<br>
&gt; [0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation,<br>
&gt; probably memory access out of range<br>
&gt; [0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br>
&gt; [0]PETSC ERROR: or see<br>
&gt; <a href="http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal%5B0%5DPETSC" target="_blank">http://www.mcs.anl.gov/petsc/petsc-as/documentation/troubleshooting.html#Signal[0]PETSC</a><br>
&gt; ERROR: or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a> on linux or man libgmalloc on Apple to<br>
&gt; find memory corruption errors<br>
&gt; [0]PETSC ERROR: likely location of problem given in stack below<br>
&gt; [0]PETSC ERROR: ---------------------  Stack Frames<br>
&gt; ------------------------------------<br>
&gt; [0]PETSC ERROR: Note: The EXACT line numbers in the stack are not<br>
&gt; available,<br>
&gt; [0]PETSC ERROR:       INSTEAD the line number of the start of the function<br>
&gt; [0]PETSC ERROR:       is given.<br>
&gt; [0]PETSC ERROR: [0] GMREScycle line 133 src/ksp/ksp/impls/gmres/gmres.c<br>
&gt; [0]PETSC ERROR: [0] KSPSolve_GMRES line 228 src/ksp/ksp/impls/gmres/gmres.c<br>
&gt; [0]PETSC ERROR: [0] KSPSolve line 308 src/ksp/ksp/interface/itfunc.c<br>
&gt; [0]PETSC ERROR: --------------------- Error Message<br>
&gt; ------------------------------------<br>
&gt; [0]PETSC ERROR: Signal received!<br>
&gt; [0]PETSC ERROR:<br>
&gt; ------------------------------------------------------------------------<br>
&gt; [0]PETSC ERROR: Petsc Release Version 3.0.0, Patch 8, Fri Aug 21 14:02:12<br>
&gt; CDT 2009<br>
&gt; [0]PETSC ERROR: See docs/changes/index.html for recent updates.<br>
&gt; [0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.<br>
&gt; [0]PETSC ERROR: See docs/index.html for manual pages.<br>
&gt; [0]PETSC ERROR:<br>
&gt; ------------------------------------------------------------------------<br>
<br>
dr. ir. Christiaan Klaij<br>
CFD Researcher<br>
Research &amp; Development<br>
</div></div>mailto:<a href="mailto:C.Klaij@marin.nl">C.Klaij@marin.nl</a><br>
<div class="im">T +31 317 49 33 44<br>
<br>
</div>MARIN<br>
2, Haagsteeg, P.O. Box 28, 6700 AA Wageningen, The Netherlands<br>
T +31 317 49 39 11, F +31 317 49 32 45, I <a href="http://www.marin.nl/" target="_blank">http://www.marin.nl/</a><br>
<a href="http://www.marin.nl/web/News/News-items/RD-seminar-on-September-22-few-places-left.htm" target="_blank">http://www.marin.nl/web/News/News-items/RD-seminar-on-September-22-few-places-left.htm</a> R&amp;D seminar on September 22: few places left<br>

<div><div></div><div class="h5"><br>
This e-mail may be confidential, privileged and/or protected by copyright. If you are not the intended recipient, you should return it to the sender immediately and delete your copy from your system.<br>
<br>
</div></div></blockquote></div><br><br clear="all"><br>-- <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>