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