<div dir="ltr"><div class="gmail_quote">On Tue, Mar 29, 2011 at 08:07, Kontsantinos Kontzialis <span dir="ltr">&lt;<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":zj"> When should I call MatMFFDGeth in order to get the differencing parameter?<br>
I call it after SNESSolve but I get zero for h.</div></blockquote></div><br><div>I suspect you are calling this function on a matrix of the wrong type (and, because this function does not check types properly, returning nonsense bytes that happen to represent a very small floating point value). Try calling</div>
<div><br></div><div><div>    ierr = SNESGetJacobian(snes,&amp;mat,0,0,0);CHKERRQ(ierr);</div><div>    ierr = MatMFFDGetH(mat,&amp;h);CHKERRQ(ierr);</div></div><div><br></div></div>