[deal.II] Getting number of hanging support points

2016-07-15 Thread Deepak Gupta
Dear All, I am currently dealing with hanging support points in my problem and I have a couple of questions. These are: 1. For my case, different elements of the mesh have different polynomial orders of the basis. Due to support points at the element edges, I believe there will be hanging support

Re: [deal.II] Bug in sparse_matrix.templates.h

2016-07-15 Thread Michał Wichrowski
Thanks Wolfgang. I think it would be helpful if there was something about it in documentation. Michał -- 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 s

[deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Bruno Turcksin
Marek, On Friday, July 15, 2016 at 9:27:02 AM UTC-4, Marek Čapek wrote: > P.S.: I have tried finer meshes (refined in gmsh) and merged meshes of > the geometries in the files in the attachment. > Here is the picture of the whole chamber > https://drive.google.com/file/d/0ByiwmbgGTAfqazRERTVEUVJK

Re: [deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Bruno Turcksin
Marek 2016-07-15 11:27 GMT-04:00 Vinetou Incucuna : > i have redesigned my *.geo files and got > much coarser meshes, which work in dealii GridIn Great. > I dont know much about adaptive refinement in deal.ii. > (Till now I used only the uniform refinement featur). > Are there some obstacles i

Re: [deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Bruno Turcksin
Marek, One more thing. Because the domain is curved, you will need to attach a manifold to the Triangulation (http://dealii.org/8.4.1/doxygen/deal.II/step_53.html) So that when you refine your mesh, the new points on the boundary are put at the right place. Best, Bruno 2016-07-15 11:37 GMT-04:0

Re: [deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Vinetou Incucuna
> > There is no difference if the mesh is structured or not. Deal doesn't > know how the mesh looks like. Yes, but I have thought, that for a programmer it is easier to refine meshes (using the refine flags), as he or she has better control over the above mentioned constraint. Or are in dealii som

Re: [deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Vinetou Incucuna
> > One more thing. Because the domain is curved, you will need to attach > a manifold to the Triangulation Thank You Bruno for the warning, I have never used it (I only assigned boundary objects to mesh). This feature is assumably more general 2016-07-15 17:45 GMT+02:00 Vinetou Incucuna : > Th

Re: [deal.II] Re: Possibly wrong format of msh file

2016-07-15 Thread Bruno Turcksin
2016-07-15 11:45 GMT-04:00 Vinetou Incucuna : > Yes, but I have thought, that for a programmer it is easier to > refine meshes (using the refine flags), as he or she has better > control over the above mentioned constraint. Or are in dealii > some helpful features, which diminish this burden? You o

Re: [deal.II] Possibly wrong format of msh file

2016-07-15 Thread Wolfgang Bangerth
On 07/15/2016 08:27 AM, Marek Čapek wrote: An error occurred in line <1713> of file in function static void dealii::GridReordering::invert_all_cells_of_negative_grid(const std::vector >&, std::vector >&) [with int dim = 3; int spacedim = 3] The violated condition was: n_negative_cells

Re: [deal.II] Bug in sparse_matrix.templates.h

2016-07-15 Thread Wolfgang Bangerth
On 07/15/2016 07:56 AM, Michał Wichrowski wrote: I think it would be helpful if there was something about it in documentation. Good point -- see https://github.com/dealii/dealii/pull/2836 Cheers W> -- Wolfgang Bangerth

Re: [deal.II] axis-symmetric problem

2016-07-15 Thread Wolfgang Bangerth
Is it possible to define a problem in spherical or cylindrical coordinate system instead of Cartesian coordinate system?It would be very kind of you if you refer me to a tutorial in which there is a problem that has been solved in spherical or cylindrical coordinate system. There is no such tu

[deal.II] Dirichlet constraint to a specific node

2016-07-15 Thread Anup Basak
Hello, I have a vector valued (say, displacement) problem and I want to apply Dirichlet boundary condition corresponding to one component of the displacement vector to a single node (whose position is know) on a boundary. I shall be thankful if someone can tell me how can I implement it in de

[deal.II] Re: geometry

2016-07-15 Thread benhour . amirian66
Dear J-P, Thanks for your response. I used the code for creating a quarter of a circle that comes as follow: Triangulation<2> triangulation; const Point<2> center; const double radius = 1.; GridGenerator::half_hyper_ball(triangulation, center, radius); Triangulation<2>::active_cell_iterator cell