<div dir="ltr"><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Jan 15, 2014 at 6:24 PM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><div id=":lo" style="overflow:hidden">Hi Praveen,<br>
<br>
To compute continuous weak approximations to du/dy (say), you<br>
would do the following:<br>
<br>
      call gradm1(ux,uy,uz,u)<br>
<br>
      n = nx1*ny1*nz1*nelv<br>
      call col2  (uy,bm1,n)<br>
      call dssum (uy,nx1,ny1,nz1)<br>
      call col2  (uy,binvm1,n)<br>
<br>
assuming all arrays appropriately declared.</div></blockquote></div><br>Thank you. I will try this.</div><div class="gmail_extra"><br></div><div class="gmail_extra">I realized that I am saving the solution in 4 bytes which could lead to loss of precision. I tried to save the averaged solution in 8 bytes with following code in usrchk</div>
<div class="gmail_extra"><br></div><div class="gmail_extra"><p style="margin:0px;font-size:11px;font-family:Menlo">      param(63) = 8.0</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">      call avg_all()</p>
<p style="margin:0px;font-size:11px;font-family:Menlo">      param(63) = 0.0</p><p style="margin:0px;font-size:11px;font-family:Menlo"><br></p>But the "avg*" files still are not saved in 8 bytes. If I set param(63)=8 in rea file then it works.</div>
<div class="gmail_extra"><br></div><div class="gmail_extra">Best</div><div class="gmail_extra">praveen</div></div>