On 10/18/11 3:02 PM, Boyce Griffith wrote: > > > On 10/6/11 3:59 AM, Lorenzo Alessio Botti wrote: >> 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. > > The points that I'm passing to reinit are reference coordinates on the > side element. > > To avoid calling reinit(elem,side,tolerance,pts), I'm now trying: > > AutoPtr<Elem> side_elem = elem->build_side(side); > fe_boundary->reinit(side_elem.get(), &reference_coords); > > but I get an error: > > ERROR: Unsupported 2D element type!: 0
Grrr, I forgot to switch dimension of the FE object... -- 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-d2d-oct _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
