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
If I print out the element type for the side proxy element, it seems
like I always get 0, regardless of the parent element type. Are proxy
elements not supposed to know their type?
-- 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
