## [Libmesh-users] Obtaining solution value at non-quadrature point

 [Libmesh-users] Obtaining solution value at non-quadrature point From: David Knezevic - 2006-05-17 11:36:20 ```Hi all, My libMesh is a little rusty, and I need to work out how to get a solution value at a non-quadrature point. I have looped over the elements of the mesh to locate the element that contains (x1,x2) and then I use an inverse map to get the point in the reference element corresponding to (x1,x2), i.e.: Point x_point(x1,x2,0); Point ref_point = inverse_map(elem,x_point); But then I can't recall how one can retrieve the basis function values at ref_point rather than the quadrature points. A pointer would be appreciated :) Regards, David ```

 [Libmesh-users] Obtaining solution value at non-quadrature point From: David Knezevic - 2006-05-17 11:36:20 ```Hi all, My libMesh is a little rusty, and I need to work out how to get a solution value at a non-quadrature point. I have looped over the elements of the mesh to locate the element that contains (x1,x2) and then I use an inverse map to get the point in the reference element corresponding to (x1,x2), i.e.: Point x_point(x1,x2,0); Point ref_point = inverse_map(elem,x_point); But then I can't recall how one can retrieve the basis function values at ref_point rather than the quadrature points. A pointer would be appreciated :) Regards, David ```
 AW: [Libmesh-users] Obtaining solution value at non-quadrature point From: Steffen Petersen - 2006-05-17 11:48:31 ```> Hi all, > > My libMesh is a little rusty, and I need to work out how to get a > solution value at a non-quadrature point. I have looped over the > elements of the mesh to locate the element that contains (x1,x2) and > then I use an inverse map to get the point in the reference element > corresponding to (x1,x2), i.e.: > > Point x_point(x1,x2,0); > Point ref_point = inverse_map(elem,x_point); > > But then I can't recall how one can retrieve the basis > function values > at ref_point rather than the quadrature points. A pointer would be > appreciated :) Just use the funtion FE::reinit(const Elem *elem, const std::vector< Point > *const pts=NULL) with your mapped non-quadrature point stored in some vector and pass the pointer instead of the default NULL, i.e. call something like fe->reinit(elem, &my_point). Steffen > Regards, > David > > > ------------------------------------------------------- > Using Tomcat but need to do more? Need to support web > services, security? > Get stuff done quickly with pre-integrated technology to make > your job easier > Download IBM WebSphere Application Server v.1.0.1 based on > Apache Geronimo > http://sel.as-us.falkag.net/sel?cmd=lnk&kid=120709&bid=263057&; dat=121642 > _______________________________________________ > Libmesh-users mailing list > Libmesh-users@... > https://lists.sourceforge.net/lists/listinfo/libmesh-users > ```