Hi Boyce, > // 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? > >
I introduced the above change in order avoid inverse maps and speed up quadrature evaluation on sides. It should be reverted to inverse map if the quadrature points you are providing are defined in the physical frame (xyz), if the quadrature points are defined on the reference side the code should work. However, if your points are in the physical frame the best option (you don't need to modify the reinit behavior) is probably to find the reference points in your code with inverse map and then use reinit providing a vector of points FEInterface::inverse_map (elem->dim(), fe->get_fe_type(), elem, your_xyz_points, your_reference_frame_points); fe_neighbor_face->reinit(neighbor, &your_reference_frame_points); This assuming that you are using basis functions defined in the reference frame. Lorenzo > ------------------------------------------------------------------------------ > 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 ------------------------------------------------------------------------------ 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
