From: Roy Stogner <roystgnr@ic...> - 2009-10-12 03:44:27
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.
Get latest updates about Open Source Projects, Conferences and News.