On Fri, Apr 15, 2011 at 2:02 PM, Boyce Griffith <griff...@cims.nyu.edu> wrote:
> Hi, Folks --
>
> Is there a function that gives the reference coordinates for each node
> of a (Lagrange) element?
>
> 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.
>
> Doing this in general seems like it could be a real mess, but for lower
> order Lagrange elements it does not seem too bad.  However, I think that
> I need to know the reference coordinates of the nodes, and I have been
> unable to find an API that seems to give access to these values.  (It is
> my impression that functions such as Elem::point() give the current
> coordinates of the node, not the reference coordinates.)
>
> Also, I realize this probably seems like kind of a bizarre thing to be
> computing, but does there happen to be any existing code in libMesh that
> could be adapted to do this kind of thing?  :-)
>
> Thanks!

FE::inverse_map()

-- 
John

------------------------------------------------------------------------------
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
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to