From: Boyce G. <gri...@ci...> - 2011-10-05 23:31:14
|
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 |