I found VectorTools::create_point_source_vector

This should be what I need. I guess I can scale this vector by any scalar /a/ and obtain /f(p) = a * \delta (p - p_s);
/Thanks/.
/
On 15.04.2012 14:50, Alexander Grayver wrote:
Hello,

What is the best way to represent delta source in deal.II?
More precisely I want my RHS to be
/f(p) = a * \delta (p - p_s);/
/f(p) != 0/ iif /|p - p_s| = 0/;

Doing this properly would mean that I have to define element within this function happens not to be zero and then distribute contribution among its q points?
--
Regards,
Alexander


--
Regards,
Alexander

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to