<div class="gmail_quote">On Wed, Nov 23, 2011 at 01:06, behzad baghapour <span dir="ltr"><<a href="mailto:behzad.baghapour@gmail.com">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 id=":3jh"> 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>You haven't told me about "noe" or "num". Do you mean for this to read xx[c*tot+p]?</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":3jh"><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><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>