<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Sep 25, 2014 at 1:50 PM, Alletto, John M <span dir="ltr"><<a href="mailto:john.m.alletto@lmco.com" target="_blank">john.m.alletto@lmco.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal">All,<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">In DMDACreate3D I declare a stencil width of 2 for a 13 point star stencil<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Later when I fill the stencil, I am assuming the following order… can you tell me if I am correct in the assumption.</p></div></div></blockquote><div><br></div><div>Are you talking about using MatSetValuesStencil()? If so, there is no inherent ordering. You</div><div>plug in the i,j,k,c values.</div><div><br></div><div> Thanks,</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div lang="EN-US" link="blue" vlink="purple"><div>
<p class="MsoNormal"><span style="background:yellow">My assumption is the each V index maps to the corresponding I,j,k indicies of W</span><u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">W is the Stencil weighting at that position<u></u><u></u></p>
<p class="MsoNormal">v[0] = W (I,j,k-2)<u></u><u></u></p>
<p class="MsoNormal">v[1] =W(I,j-2,k)<u></u><u></u></p>
<p class="MsoNormal">v[2] =W( i-1,j,k)<u></u><u></u></p>
<p class="MsoNormal">v[3] =W( I,j,k-1)<u></u><u></u></p>
<p class="MsoNormal">v[4] =W( I,j-1,k)<u></u><u></u></p>
<p class="MsoNormal">v[5] =W( i-1,j,k)<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">v[7] =W( i+1, j, k)<u></u><u></u></p>
<p class="MsoNormal">v[8] =W( I,j+1,k)<u></u><u></u></p>
<p class="MsoNormal">v[9] =W( I,j,k+1)<u></u><u></u></p>
<p class="MsoNormal">v[10] =W (i+2,j,k)<u></u><u></u></p>
<p class="MsoNormal">v[11] =W(I,j+2,k)<u></u><u></u></p>
<p class="MsoNormal">v[12] =W(I,j,k+2)<u></u><u></u></p>
<p class="MsoNormal">v[6] = sum( v[1..5])+ sum(v[6..12])<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Thanks<span class="HOEnZb"><font color="#888888"><u></u><u></u></font></span></p><span class="HOEnZb"><font color="#888888">
<p class="MsoNormal">John <u></u><u></u></p>
</font></span></div>
</div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener
</div></div>