<br><br><div class="gmail_quote">On Wed, Nov 23, 2011 at 4:33 PM, Jed Brown <span dir="ltr"><<a href="mailto:jedbrown@mcs.anl.gov">jedbrown@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div class="gmail_quote"><div class="im">On Wed, Nov 23, 2011 at 01:06, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com" target="_blank">behzad.baghapour@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div> for( c=0; c<ne; c++ ) for( p=0; p<tot; p++ )
<br> FC->e[c].Q[p] = xx[c*(noe*num)+p];<br></div></blockquote><div><br></div></div><div>You haven't told me about "noe" or "num". Do you mean for this to read xx[c*tot+p]?</div></div></blockquote>
<div><br>noe is the number of equations (equal to number of flow states) and num is the number of shape functions in element (according to the order of accuracy) <br></div><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div class="gmail_quote"><div class="im"><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div><br> ierr = VecRestoreArray( x, &xx ); CHKERRQ( ierr );<br><br> interiorFlux( FC->flw, FC->e );
<br> faceFlux ( FC->flw, FC->f, FC->e );</div></blockquote></div></div><br><div>The first of these should set FC->e (if you are adding into it, then you need to zero it first) and the second should add into it.</div>
</blockquote><div><br>I found that my mistake is that I forgot to set upwind effect in my defined field context. The problem was not stable, it is not related to SNES setup.<br><br>However, I can't find out the difference between "basic" and "basicnonorms" in line search method.<br>
<br>In addition, Is it possible to set the number of line-search corrections or is just decided by the solver?<br><br>Thanks a lot,<br>BehZad<br> <br></div></div><br>