Thanks. I was out of town the end of last week, but it looks like you've
probably found the issue without help...
I'll take a look, but the fix sounds appropriate. BTW, the arrays of
magic numbers map 1D shape functions for tensor-product evaluations, and
there is no guarantee they are right.
[mailto:libmesh-devel-admin@...] On Behalf Of Roy Stogner
Sent: Thursday, October 06, 2005 2:03 PM
Subject: Re: [Libmesh-devel] Problem with 3D cubic Hierarchics?
On Wed, 5 Oct 2005, Roy Stogner wrote:
> On Wed, 5 Oct 2005, Roy Stogner wrote:
>> The HIERARCHIC/HEX27/THIRD specialization of shape() depends on three
>> long arrays of magic numbers. Is the original author out there to
>> help me bughunt?
> This may be a false alarm - after investigating the arrays of magic
> numbers I'm becoming convinced that the bug is in my code after all.
> No need to start poring through fe_hierarchic_shape_etc yet.
Okay, time to pore through fe_hierarchic_shape_3D. On a simple one-cube
mesh, shape functions index 2 and 3 were evaluating identically, and so
giving identical rows in a mass matrix. It appears that the edge
orientation tests which should have only applied to edge DoFs were being
applied to vertex DoFs as well, so I've checked in a 2 line fix that should
correct that. It seems to pass my simple tests now, but I'd appreciate
others taking a look.
This SF.Net email is sponsored by:
Power Architecture Resource Center: Free content, downloads, discussions,
and more. http://solutions.newsforge.com/ibmarch.tmpl
Libmesh-devel mailing list