<div dir="ltr"><div class="gmail_default" style="font-size:small"><div class="gmail_default" style="font-size:small">Hello,</div><div class="gmail_default" style="font-size:small">To know more about the subroutines used for dealiasing, I perform a simple exercise.</div><div class="gmail_default" style="font-size:small">I take a simple sine function and project it onto the fine mesh and back to the coarse using intp_rstd.</div><div class="gmail_default" style="font-size:small">When I compare the values of the processed array with the original array, it is scaled by a factor of 1.25204, ie, I get a function like 1.25204*sin(2*x)<br></div><div class="gmail_default" style="font-size:small">Could you please tell me which factor am I missing? Which jacobian has to be multiplied?</div><div class="gmail_default" style="font-size:small">Attached is the usr file and below are the operations I perform in short.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default">      sinax(i,1,1,1)      = sin(2.*x)<br></div><div class="gmail_default"><div class="gmail_default">      </div><div class="gmail_default">     call intp_rstd(uf1,u1,nx1,nxd,if3d,<wbr>0) ! 0 --> forward</div><div class="gmail_default">      nxyzd = lxd*lyd*lzd</div><div class="gmail_default">         do i=1,nxyzd</div><div class="gmail_default">            tr1   = uf1(i)</div><div class="gmail_default">           <span> </span><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:small;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">wf(i)</span><span> </span>= tr1*rx(i,1,e)+tr1*rx(i,2,e)+<wbr>tr1*rx(i,3,e)</div><div class="gmail_default">         enddo<br></div><div class="gmail_default"><br></div><div class="gmail_default">      call intp_rstd(w1,wf,nx1,nxd,if3d,<wbr>1) ! 1 --> back to coarse</div><div class="gmail_default"><br></div><div class="gmail_default">      nxyz = lx1*ly1*lz1           !              </div><div class="gmail_default">         do i=1,nxyz               ! to physical space.</div><div class="gmail_default">            bi = 1./((bm1(i,1,1,e)))   !</div><div class="gmail_default">            w1(i) = bi*w1(i)</div><div class="gmail_default"><div class="gmail_default">            jacm=jacmi(i,e) ! I dont know whether to use jacmi(i,e) or not</div><div>         enddo</div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Thank you in advance.</span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline"><br></span></div><div><span style="color:rgb(34,34,34);font-family:arial,sans-serif;font-size:12.8px;font-style:normal;font-variant-ligatures:normal;font-variant-caps:normal;font-weight:400;letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Sandeep</span></div></div></div><br clear="all"></div><div><div class="gmail_signature"><div dir="ltr"><div dir="ltr"><br></div><div><div style="margin:2px 0px 0px;font-size:12.8px"></div></div></div></div></div>
</div>