Hello
I want to use the FE<dim,XYZ> space which is a discontinuous space of
polynomials.
Suppose I am at element "elem" and let
nbr_elem = neighbour of elem
I want to evaluate the shape functions of nbr_elem at the quadrature points
of elem. Is this possible and how can I do it. I suppose I would have to do
something like
const std::vector<Point>& q_point = fe->get_xyz();
for ( ; el != end_el; ++el)
{
const Elem* elem = *el;
fe->reinit(elem);
// get iterator for neighbour: nbr_elem
fe_nbr->reinit(nbr_elem);
I want to evaluate fe_nbr at q_point but do not know how to do this.
Thanks
praveen
------------------------------------------------------------------------------
Introducing Performance Central, a new site from SourceForge and
AppDynamics. Performance Central is your source for news, insights,
analysis and resources for efficient Application Performance Management.
Visit us today!
http://pubads.g.doubleclick.net/gampad/clk?id=48897511&iu=/4140/ostg.clktrk
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
