<div dir="ltr">Brilliant work Bill,<div> That is a really interesting paper. I'm glad you found a solution.</div><div>Your work here just underscores the fundamental importance of good testing. I'm looking forward to seeing it in VTK.</div><div><br></div><div>I had never heard of SymPy, I had always used Maple but it is too expensive now that I have retired. I'll give it a SymPy a go.</div><div><br></div><div>Nothing wrong with punch cards ... except when you drop a whole box. I used them to do Tempromandibular joint (TMJ) analysis using FORTRAN, way back in the 70's. Then I graduated to using a Teletype and APL!</div><div><br></div><div>Regards</div><div> Andrew<br><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Jun 16, 2015 at 6:26 PM, <span dir="ltr"><<a href="mailto:vtk-developers-request@vtk.org" target="_blank">vtk-developers-request@vtk.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send vtk-developers mailing list submissions to<br>
<a href="mailto:vtk-developers@vtk.org">vtk-developers@vtk.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<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>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:vtk-developers-request@vtk.org">vtk-developers-request@vtk.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:vtk-developers-owner@vtk.org">vtk-developers-owner@vtk.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of vtk-developers digest..."<br>
<br>Today's Topics:<br>
<br>
1. Pentagonal prism bug resolved! (Bill Lorensen)<br>
2. Re: Pentagonal prism bug resolved! (Berk Geveci)<br>
3. Re: Pentagonal prism bug resolved! (Bill Lorensen)<br>
4. Re: Pentagonal prism bug resolved! (Berk Geveci)<br>
5. Re: VTK for android example: missing vtkRenderingOpenGL2<br>
(Lonni Besan?on)<br>
6. android_native_app_glue.c doesn't exist (Lonni Besan?on)<br>
<br><br>---------- Forwarded message ----------<br>From: Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>To: VTK Developers <<a href="mailto:vtk-developers@vtk.org">vtk-developers@vtk.org</a>><br>Cc: <br>Date: Mon, 15 Jun 2015 13:31:58 -0400<br>Subject: [vtk-developers] Pentagonal prism bug resolved!<br>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>
<br><br>---------- Forwarded message ----------<br>From: Berk Geveci <<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>><br>To: Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>Cc: VTK Developers <<a href="mailto:vtk-developers@vtk.org">vtk-developers@vtk.org</a>><br>Date: Mon, 15 Jun 2015 14:17:00 -0400<br>Subject: Re: [vtk-developers] Pentagonal prism bug resolved!<br><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>
<br><br>---------- Forwarded message ----------<br>From: Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>To: Berk Geveci <<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>><br>Cc: VTK Developers <<a href="mailto:vtk-developers@vtk.org">vtk-developers@vtk.org</a>><br>Date: Mon, 15 Jun 2015 14:27:34 -0400<br>Subject: Re: [vtk-developers] Pentagonal prism bug resolved!<br>Brings back memories of my finite element days, 1968-1978, at<br>
Watervliet Arsenal. At that time we worked with NASTRAN and we added a<br>
few elements for crack tip analysis....<br>
<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic<br>
Isoparametric Element as a Singular Element for Crack Problems,<br>
International Journal for Numerical Methods in Engineering, vol. 12,<br>
no. 11, pp. 1727-1742, 1978.<br>
<br>
Hussain, M. A. and W. E. Lorensen, Isoparametric Elements as Singular<br>
Elements for Crack Problems, AR-LCB-TR-77041, Benet Weapons<br>
Laboratory, Watervliet, NY, November 1977.<br>
<br>
Hussain, M. A.; Pu, S. L.; Lorensen, W. E., Singular Plastic Element:<br>
NASTRAN Implementation and Application, Sixth NASTRAN Users' Colloq.,<br>
pp. 257-274, October, 1977, (NASA CP-2018).<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, Collapsed 12-Node<br>
Triangular Elements as Crack Tip Elements for Elastic Fracture,<br>
ARLCB-TR-77047, Watervliet Arsenal, Watervliet, NY, December 1977.<br>
<br>
Lorensen, W. E. and M. A. Hussain, "The Stability of Isoparametric<br>
Elements as a Singular Element for Crack Problems," in Developments in<br>
Mechanics, ed. T. C. T. Ting, 1977.<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic<br>
Isoparametric Element as a Singular Element for Crack Problems,<br>
AR-LCB-TR77047, Watervliet Arsenal, 1977.<br>
<br>
Hussain, M. A., W. E. Lorensen, and G. Pflegl, The Quarter-Point<br>
Quadratic Isoparametric Element as A Singular Element for Crack<br>
Problems, TM-X-3428, p. 419, NASA, 1976.<br>
<br>
<br>
On Mon, Jun 15, 2015 at 2:17 PM, Berk Geveci <<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>> wrote:<br>
> You rock Bill.<br>
><br>
> On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>
> wrote:<br>
>><br>
>> 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<br>
>> <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>
><br>
<br>
<br>
<br>
--<br>
Unpaid intern in BillsBasement at noware dot com<br>
<br>
<br><br>---------- Forwarded message ----------<br>From: Berk Geveci <<a href="mailto:berk.geveci@kitware.com">berk.geveci@kitware.com</a>><br>To: Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com">bill.lorensen@gmail.com</a>><br>Cc: VTK Developers <<a href="mailto:vtk-developers@vtk.org">vtk-developers@vtk.org</a>><br>Date: Mon, 15 Jun 2015 15:03:09 -0400<br>Subject: Re: [vtk-developers] Pentagonal prism bug resolved!<br><div dir="ltr">Finite element analysis with punch cards. That's hard core!<div><br></div><div>-berk</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Jun 15, 2015 at 2:27 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">Brings back memories of my finite element days, 1968-1978, at<br>
Watervliet Arsenal. At that time we worked with NASTRAN and we added a<br>
few elements for crack tip analysis....<br>
<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic<br>
Isoparametric Element as a Singular Element for Crack Problems,<br>
International Journal for Numerical Methods in Engineering, vol. 12,<br>
no. 11, pp. 1727-1742, 1978.<br>
<br>
Hussain, M. A. and W. E. Lorensen, Isoparametric Elements as Singular<br>
Elements for Crack Problems, AR-LCB-TR-77041, Benet Weapons<br>
Laboratory, Watervliet, NY, November 1977.<br>
<br>
Hussain, M. A.; Pu, S. L.; Lorensen, W. E., Singular Plastic Element:<br>
NASTRAN Implementation and Application, Sixth NASTRAN Users' Colloq.,<br>
pp. 257-274, October, 1977, (NASA CP-2018).<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, Collapsed 12-Node<br>
Triangular Elements as Crack Tip Elements for Elastic Fracture,<br>
ARLCB-TR-77047, Watervliet Arsenal, Watervliet, NY, December 1977.<br>
<br>
Lorensen, W. E. and M. A. Hussain, "The Stability of Isoparametric<br>
Elements as a Singular Element for Crack Problems," in Developments in<br>
Mechanics, ed. T. C. T. Ting, 1977.<br>
<br>
Pu, S. L., M. A. Hussain, and W. E. Lorensen, The Collapsed Cubic<br>
Isoparametric Element as a Singular Element for Crack Problems,<br>
AR-LCB-TR77047, Watervliet Arsenal, 1977.<br>
<br>
Hussain, M. A., W. E. Lorensen, and G. Pflegl, The Quarter-Point<br>
Quadratic Isoparametric Element as A Singular Element for Crack<br>
Problems, TM-X-3428, p. 419, NASA, 1976.<br>
<div><div><br>
<br>
On Mon, Jun 15, 2015 at 2:17 PM, Berk Geveci <<a href="mailto:berk.geveci@kitware.com" target="_blank">berk.geveci@kitware.com</a>> wrote:<br>
> You rock Bill.<br>
><br>
> On Mon, Jun 15, 2015 at 1:31 PM, Bill Lorensen <<a href="mailto:bill.lorensen@gmail.com" target="_blank">bill.lorensen@gmail.com</a>><br>
> wrote:<br>
>><br>
>> 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<br>
>> <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>
><br>
<br>
<br>
<br>
</div></div><span><font color="#888888">--<br>
Unpaid intern in BillsBasement at noware dot com<br>
</font></span></blockquote></div><br></div>
<br><br></blockquote></div><div><br></div>-- <br><div class="gmail_signature">___________________________________________<br>Andrew J. P. Maclean<br><br>___________________________________________</div>
</div></div></div>