On Wed, Jul 20, 2016 at 12:57 PM, Xujun Zhao <xzha...@gmail.com> wrote:
> I also tried. You are right. It looks like this element and point don't > have the zero det problem. I also tried other mesh and points. They crashed > down and led to the same error messages: > > Assertion `my_det != static_cast<T>(0.)' failed. > my_det = 0 > static_cast<T>(0.) = 0 > > But when I extract the nodal information and construct a one-element mesh, > the function contains_point() works well again. I can't understand why this > happen. > Me either, contains_point() should not depend on whether you are in a 1 element mesh or a Mesh with many elements. Btw, I use mesh.read() to read a mesh file created by CUBIT, then let it > all_second_order(). Will this cause problems? > No, that should work fine unless there has been a regression recently in all_second_order(). -- John ------------------------------------------------------------------------------ What NetFlow Analyzer can do for you? Monitors network bandwidth and traffic patterns at an interface-level. Reveals which users, apps, and protocols are consuming the most bandwidth. Provides multi-vendor support for NetFlow, J-Flow, sFlow and other flows. Make informed decisions using capacity planning reports.http://sdm.link/zohodev2dev _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users