<div dir="ltr">You rock Bill.</div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <span dir="ltr"><<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Folks,<br>
<br>
I reported recently of issues with the shape functions and derivatives<br>
for the pentagonal prism. Since I could not find any references for<br>
the current implementation I used the formulation in this paper:<br>
<a href="http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf" rel="noreferrer" target="_blank">http://dilbert.engr.ucdavis.edu/~suku/nem/papers/polyelas.pdf</a><br>
Since this paper did not have the shape function derivatives, I used<br>
SymPy ( <a href="http://www.sympy.org/en/index.html" rel="noreferrer" target="_blank">http://www.sympy.org/en/index.html</a> ) to compute the<br>
derivatives.<br>
<br>
I'm happy to report that the new formulation does a great job and<br>
passes the new unit test I created for cells.<br>
<br>
I have not submitted a patch yet since I want the code to document the<br>
process. I wrote c++ code and SymPy code to generate the statements<br>
for the new formulation.<br>
<br>
I wanted to give a heads up in case someone else is looking into this.<br>
<br>
Bill<br>
_______________________________________________<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Search the list archives at: <a href="http://markmail.org/search/?q=vtk-developers" rel="noreferrer" target="_blank">http://markmail.org/search/?q=vtk-developers</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/vtk-developers" rel="noreferrer" target="_blank">http://public.kitware.com/mailman/listinfo/vtk-developers</a><br>
<br>
</blockquote></div><br></div>