Hi Roy, I spent some time looking through the ProjectVector::operator() code, and have some questions:
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. (?) — 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? — 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? Thanks, -Manav ------------------------------------------------------------------------------ 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
