[deal.II] Re: geometry

2016-07-17 Thread benhour . amirian66
Daniel, Thanks very much for your help. For improving my knowledge and information, It would be very kind of you if you let me know why my code does not produce right geometry and which part is wrong. Sincerely, Benhour On Sunday, July 17, 2016 at 7:10:41 AM UTC-5, Daniel Arndt wrote: > > Benho

[deal.II] Re: Output problem of DG with degree =0

2016-07-17 Thread Jiaqi ZHANG
Hello Jean-Paul, It works, thank you very much for you help. Jiaqi 在 2016年7月17日星期日 UTC-4上午6:03:37,Jean-Paul Pelteret写道: > > Dear Jiaqui, > > The issue here is that when you choose a piece-wise constant ansatz you > end up with exactly the same number of DoFs as there are cells. > deal.II usuall

[deal.II] Fwd: Getting number of hanging support points

2016-07-17 Thread Deepak Gupta
Dear All, Since I didn't receive any response to my question, I thought of raising it up as a reminder if someone can help. Best Deepak On Friday, July 15, 2016 at 2:02:43 PM UTC+2, Deepak Gupta wrote: > > Dear All, > > I am currently dealing with hanging support points in my problem and I >

[deal.II] Re: geometry

2016-07-17 Thread Daniel Arndt
Benhour, you can find a modification of step-3 in which the domain is a quarter circle attached. The four cells from GridGenerator::half_hyper_ball have been replaced by just three. The vertex (0,-radius) has been removed and all the remaining vertices with a negative y-component have been move

[deal.II] Re: Dirichlet constraint to a specific node

2016-07-17 Thread Jean-Paul Pelteret
Dear Anup, In addition to Daniel's approach, if you're using just linear shape functions then vertex_dof_index can also help here. See this post and this list of search results

[deal.II] Re: Output problem of DG with degree =0

2016-07-17 Thread Jean-Paul Pelteret
Dear Jiaqui, The issue here is that when you choose a piece-wise constant ansatz you end up with exactly the same number of DoFs as there are cells. deal.II usually treats the output of scalar cell data (e.g. if you were to output the material ID of each cell) in a special way, so you just need