On Mon, 8 Nov 2010, Roy Stogner wrote: > Indeed it does! I can replicate that (seeing an exception in the > ghosted vector support) just by switching ex14 to cubic hierarchics. > > I'll look into this ASAP.
I rooted out the problem I found, and the bug fix for that is in the svn head now. Saurabh, would you test this with your own code in debug or devel mode and make sure that no further assertions are thrown? I can't reproduce any other errors myself yet, but, as you've had the misfortune to find out, it's been a little while since anyone's given the higher order HIERARCHICs a full shaking out. Some more info on the bug, for other developers interested: During the work on ghosted vectors, we did a little optimization on the DofMap send_list, so now it should contain only ghost dof IDs and not all semilocal IDs. When we removed the addition of local IDs, it also inadvertently removed the only place where a certain category of ghost DoF ID was added: IDs on hanging nodes might be inadvertently omitted if the finite element type had vertex DoFs that did not overlap their side/edge DoFs. So the bug affected Clough-Tocher elements, Hierarchics, and 4th-order-or-higher Hermites. The DofMap bug fix is in there now, but I'm going to wait for Saurabh's reply before putting out a libMesh 0.7.0.3. We also construct a different sort of send_list in independent code for System::project_vector, and although I couldn't find any bugs in that code, Saurabh first reported seeing a problem with projections, so maybe there's a failure case that I just haven't noticed or hit yet. This bug would have been a silent data-corrupting error if not for the extra tests in the ghosted vector code. Thanks again to Tim: if you're ever in Austin dinner is on me. --- Roy ------------------------------------------------------------------------------ The Next 800 Companies to Lead America's Growth: New Video Whitepaper David G. Thomson, author of the best-selling book "Blueprint to a Billion" shares his insights and actions to help propel your business during the next growth cycle. Listen Now! http://p.sf.net/sfu/SAP-dev2dev _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
