<div dir="ltr">Bill it could be that the vtkPyramid functions were originally written (a long time ago) as degenerate hex shape functions which might explain the apex parametric coordinates. I believe Mathieu found and implemented an alternative formulation and also added the quadratic version; I'm not sure where to find that reference.<div>W</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 20, 2015 at 2: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>
As I continue to work through a comprehensive unit test for vtk cells,<br>
I have found an issue(s) with vtkPyramid and vtkQuadraticPyramid.<br>
<br>
For example, the parametric coordinates for the Pyramid (and Quadratic<br>
Pyramid) place the apex at 0, 0, 1. I believe it should be centered on<br>
the base (which would be .5, .5, 1).<br>
<br>
Changing that causes other issues with the shape functions.<br>
<br>
I am seeking a copy of the reference quoted in the QuadraticPyramid<br>
code so that I can verify the implementation.<br>
// The shape functions and derivatives could be implemented thanks to<br>
// the report Pyramid Solid Elements Linear and Quadratic Iso-P Models<br>
// From Center For Aerospace Structures<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature">William J. Schroeder, PhD<br>Kitware, Inc.<br>28 Corporate Drive<br>Clifton Park, NY 12065<br><a href="mailto:will.schroeder@kitware.com" target="_blank">will.schroeder@kitware.com</a><br><a href="http://www.kitware.com" target="_blank">http://www.kitware.com</a><br>(518) 881-4902</div>
</div>