<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Wed, Feb 26, 2014 at 12:29 PM, Jones,Martin Alexander <span dir="ltr"><<a href="mailto:MAJones2@mdanderson.org" target="_blank">MAJones2@mdanderson.org</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 lang="EN-US" link="blue" vlink="purple">
<div>
<p class="MsoNormal">How do we input an auxillary field so that we can use the a[] array in the f0_u function in ex12.c ?<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Void f0_u(const PetscScalar u[], const PetscScalar gradU[], const PetscScalar a[], const PetscScalar gradA[],const PetscReal x[], PetscScalar f0[])<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">I would like a=constant within a bounded volume and a=0 outside the bounded volume.<u></u><u></u></p>
<p class="MsoNormal"><u></u> <u></u></p>
<p class="MsoNormal">Any help on this is much appreciated. Thanks.</p></div></div></blockquote><div><br></div><div>This mode sounds like -variable_coefficient field, which creates an a[] in some FE space. Then you use</div>
<div>-mat_petscspace_order k for order k (in the example I test 0 and 1). The function SetupMaterial() makes<br></div><div>the auxiliary field values, and uses the pointwise function nu_2d() to do it. You can substitute your own</div>
<div>pointwise function there, or take over at a higher level and directly put your values in the Vec.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div><br></div><div><br></div><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 lang="EN-US" link="blue" vlink="purple"><div><p class="MsoNormal"><span style="color:rgb(136,136,136)"> </span></p><span class=""><font color="#888888"><p class="MsoNormal"><u></u></p>
<p class="MsoNormal">Martin Jones<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>