Ben... did you get a reply in there that didn't go to the list? I
don't see a reply from you.
On 7/12/06, Roy Stogner <roystgnr@...> wrote:
> On Wed, 12 Jul 2006, Kirk, Benjamin (JSC-EG) wrote:
> > Have a look at the embedding matrices in src/geom/face_quad4.C, for
> > example. These define how children's new nodes are created from old
> > nodes.
Thanks for pointing me there... I think that's the "right" way to do
it... I have cludged up something that works for Quad4's with a level
one rule... but it will totally destroy itself for anything else
(including triangles). I have written it as a seperate mesh_tools
method... if anyone is interested let me know.... not that cludgy
narrow minded code is all that interesting....
> Good catch - I completely forgot about that.
> > Of course, some day we might be able to
> > solve the smoothing linear system the same way we solve all others, then
> > the dof_map constraints could be applied directly.
This would, of course, be the end goal... but it's still a long way off.
> Only for isogeometric Lagrange elements. I've updated more than
> enough old code that used to only work with Lagrange elements; could
> we please try to make new code more generalized from the start?
I have tried to write my code as generally as I can... but even when I
think it is general it almost always breaks with anything other than
plain quads or tris... the reason being that I just don't have enough
experience with all the different FE types yet... at the end of my
time here I hope to leave a good clean base that doesn't make any
assumptions. We'll see...
One of my big problems is that smoothing is inherently tied to the
physical mesh (meaning vertexes) reconciling that with nodes and dofs
has been tough for me... but I'm coming around.