<div dir="ltr"><div><div><div><div><div>ah ok <br></div>i thought your composition was all positive in the nodal points but negative in the integration points after interpolation.<br></div>in that case what it suggested should work nicely<br>

</div>also for steep thermal gradient for instance this is a nice trick<br><br></div>cheers<br></div>Thomas<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Nov 4, 2013 at 6:37 PM, Timo Heister <span dir="ltr"><<a href="mailto:heister@clemson.edu" target="_blank">heister@clemson.edu</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">> this is the result of higher order basis functions no?<br>
<br>
</div>This also happens with linear elements. Without doing anything about<br>
it, you can not expect the composition to stay in [0,1].<br>
<span class="HOEnZb"><font color="#888888"><br>
<br>
--<br>
Timo Heister<br>
<a href="https://urldefense.proofpoint.com/v1/url?u=http://www.math.clemson.edu/~heister/&k=dsHy%2FVymniCD0osh6tze%2Fw%3D%3D%0A&r=WLGUy1n%2F1Vzsqown7KFPcAEzR5DIdbP771jYw4Spy%2Bk%3D%0A&m=Drl%2BpNA39ejCSc2vhWUs%2FSZ74xtbMIIkMV6IHSFmQsk%3D%0A&s=e3eb62ab86b2009902a8015b8bc9703b26b1176554f70d198b895c75601a2071" target="_blank">https://urldefense.proofpoint.com/v1/url?u=http://www.math.clemson.edu/~heister/&k=dsHy%2FVymniCD0osh6tze%2Fw%3D%3D%0A&r=WLGUy1n%2F1Vzsqown7KFPcAEzR5DIdbP771jYw4Spy%2Bk%3D%0A&m=Drl%2BpNA39ejCSc2vhWUs%2FSZ74xtbMIIkMV6IHSFmQsk%3D%0A&s=e3eb62ab86b2009902a8015b8bc9703b26b1176554f70d198b895c75601a2071</a><br>


</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Aspect-devel mailing list<br>
<a href="mailto:Aspect-devel@geodynamics.org">Aspect-devel@geodynamics.org</a><br>
<a href="http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel" target="_blank">http://geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel</a><br>
</div></div></blockquote></div><br></div>