On Sun, 11 Oct 2009, Anthony Han wrote: > Can anybody tell me how to implement a point source within the > domain using libMesh? Is there a function that can realize it?
Use a PointLocator to find which element your source is located in (this answers your other email's question), then use FE::inverse_map() to get the point's corresponding location on the master element, and finally use FE::reinit(Elem*, std::vector<Point>*) to calculate shape functions/derivatives at that point. Then all that's left is the arithmetic corresponding to your weak formulation, same as in the example assembly routines. --- Roy ------------------------------------------------------------------------------ Come build with us! The BlackBerry(R) Developer Conference in SF, CA is the only developer event you need to attend this year. Jumpstart your developing skills, take BlackBerry mobile applications to market and stay ahead of the curve. Join us from November 9 - 12, 2009. Register now! http://p.sf.net/sfu/devconference _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
