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
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
