Bruno, Thanks for your reply. However, I Don't know how to use this function. In the documentation of the VectorTools::create_point_source_vector function it says:
void VectorTools::create_point_source_vector ( const Mapping <http://www.dealii.org/8.1.0/doxygen/deal.II/classMapping.html>< dim, spacedim > & mapping, const DoFHandler <http://www.dealii.org/8.1.0/doxygen/deal.II/classDoFHandler.html>< dim, spacedim > & dof, const Point <http://www.dealii.org/8.1.0/doxygen/deal.II/classPoint.html>< spacedim > & p, Vector <http://www.dealii.org/8.1.0/doxygen/deal.II/classVector.html>< double <http://www.dealii.org/8.1.0/doxygen/deal.II/classdouble.html> > & rhs_vector ) It Creates a right hand side vector for a point source at point p. In other > words, it creates a vector F so that Fi=∫Ωδ(x−p)ϕi(x)dx. Prior content of > the given rhs_vector vector is deleted. > I was wondering if I should use the rhs_vector when assembling the cell_rhs as follows: cell_rhs(i) += (rhs_vector(i) * fe_values.JxW (q)); or do something else. would you be able to give me part of a sample code where this function has been used? Thanks -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
