On Tue, 5 May 2015, David Knezevic wrote: > I think the problem is that calling "elem->print_info()" then calls > "elem->volume()", which reinits an FE object, which triggers > negative jacobian error again, and then you keep looping forever.
Ha! We should have seen that coming. The best fix is probably not to revert the print_info(), though, but to use a static bool to make sure we only hit it once. "failing=true; print_info(); libmesh_error();" if it's false; zero out stuff and return garbage (on the theory that *some* info is better than none) if it's true. --- Roy ------------------------------------------------------------------------------ One dashboard for servers and applications across Physical-Virtual-Cloud Widest out-of-the-box monitoring support with 50+ applications Performance metrics, stats and reports that give you Actionable Insights Deep dive visibility with transaction tracing using APM Insight. http://ad.doubleclick.net/ddm/clk/290420510;117567292;y _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
