Is FE::reinit() in fe_boundary.C intended to work when qrule is NULL? Line 203 of fe_boundary.C calls qrule->get_points() whether or not qrule has been set to anything:
// Find where the integration points are located on the // full element. std::vector<Point> qp; //this->inverse_map (elem, xyz, qp, tolerance); this->side_map(elem, side.get(), s, qrule->get_points(), qp); Can this revert to using inverse_map if there are user-supplied points? Or is this not the right solution? Thanks, -- Boyce ------------------------------------------------------------------------------ All the data continuously generated in your IT infrastructure contains a definitive record of customers, application performance, security threats, fraudulent activity and more. Splunk takes this data and makes sense of it. Business sense. IT sense. Common sense. http://p.sf.net/sfu/splunk-d2dcopy1 _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
