On Fri, 6 Dec 2013, Manav Bhatia wrote:
K U = F solution based projection is performed for h-refined
elements. For p-refined, however, the dof based assignment is
performed.
— Is there a reason for why the same K U = F based projection is
not performed for p-refined elements? I understand that this would
be more expensive than what is currently being done for p-refined,
but it will likely also take care of the non-hierarchich basis. (?)
IIRC it simply would have made the code more complicated, and I was
more worried about bugs from overly-complicated code than about
efficiency or generality.
— Is the dof based assignment for p-refined elements the source of
the error for hp-refinement that you had referred to in your initial
email?
Not the "source" exactly, but it's the location of the error, yes.
— I noticed that a similar operation are performed in
FEBase::coarsened_dof_values() and
FEBase::compute_proj_constraints(). So, these too would need
attention for hp- studies?
I'm not sure if those are wrong, but they'd definitely need
attention, yes.
---
Roy
------------------------------------------------------------------------------
Sponsored by Intel(R) XDK
Develop, test and display web and hybrid apps with a single code base.
Download it for free now!
http://pubads.g.doubleclick.net/gampad/clk?id=111408631&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
