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

2016-07-16 Thread Jiaqi Zhang
Hey all, I am trying to solve a level set equation with DG. When I use degree=1, everything is ok, but when I change it to degree=0, it cannot output the results, and says: An error occurred in line <957> of file in function void dealii::DataOut_DoFData::add_data_vector(const VectorType&

[deal.II] Re: geometry

2016-07-16 Thread benhour . amirian66
Daniel, In that code, I had a half circle in the right part of the coordinate axis. By writing v(1) <=0, I wanted to choose the part with negative y components from the geometry, however, The attached file has been obtained. I really do not know how I should modify the source you sent to me beca

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

2016-07-16 Thread Daniel Arndt
Anup, if you know which node you want to constrain, you just have to find the corresponding degree of freedom and add a constraint to the ConstraintMatrix. For finding this DoF, DoFTools::map_dofs_to_support_points[1] might be useful. The alternative is to loop through all the cells and ask eac

[deal.II] Re: geometry

2016-07-16 Thread Daniel Arndt
Benhour, It would be helpful if you tell us what the problem with the newly created problem is. Looking at the implementation of GridGenerator::half_hyper_ball [1], you can observe that there are cells which have vertices both with positive and negative y-components. This will likely lead to und