Jiaqi,

to create a point source, you can use 
VectorTools::create_point_source_vector 
(http://dealii.org/developer/doxygen/deal.II/namespaceVectorTools.html#ad03b858b1a3b59003a76f6224e67efc7)

Best,

Bruno

On Monday, August 29, 2016 at 10:39:27 AM UTC-4, Jiaqi Zhang wrote:
>
> Dear All,
>
> I was trying to use DG to solve a Hamilton-Jacobi equation, and the 
> problem has a point source condition
>
> for example, the center of a unit square pointCenter = 0.
>
> First I thought of using
> ZeroSourcePoint.add_line(cell->vertex_dof_index(v,0));
>
> but later found that DG elements aren't defined at the vertices of the 
> mesh in FAQ.
>
> So how do I impose this point condition in DEAL II?
>
> Thanks in advance!
>
> Best,
> Jiaqi
>

-- 
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.

Reply via email to