On 4/15/11 4:06 PM, Roy Stogner wrote: > > On Fri, 15 Apr 2011, Boyce Griffith wrote: > >> The reason I ask is that I have a variable, say x(X), restricted to a >> particular side of a particular element for which I would like to >> identify any positions at which x(X0) = x0. If any such positions >> exist, I then need to evaluate some other field variables at those >> locations, which I can easily accomplish if I know the reference >> coordinates corresponding to position X0, but not if I only know the >> value(s) of X0. > > Just use inverse_map()? That way you can also handle off-node > positions.
Thanks, I had managed to totally overlook inverse_map()... It seems like it would at least provide an easy way to get the element coordinates. In fact, it seems like inverse_map() comes close to doing everything that I want to do. I'm trying to do is to find any intersection(s) between a line and a side (in 3D) or edge (in 2D). Conveniently, the lines in question are always in the coordinate directions, so I think that maybe I can simply project the sides in the direction of the line and then call inverse_map(). The only complication is that the coordinate system being used for the mesh is a Lagrangian (material) coordinate system, and the line is in the current (physical) coordinate system, but I suppose that I can call Elem::build_side() and then reset the nodal positions of the generated side element to correspond to the current positions. Thanks, -- Boyce ------------------------------------------------------------------------------ Benefiting from Server Virtualization: Beyond Initial Workload Consolidation -- Increasing the use of server virtualization is a top priority.Virtualization can reduce costs, simplify management, and improve application availability and disaster protection. Learn more about boosting the value of server virtualization. http://p.sf.net/sfu/vmware-sfdev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
