[deal.II] Getting number of hanging support points
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 points. How can I know the total number of these hanging points for the whole mesh? 2. How can I know the number of hanging support points per element? Since I need to know how many dofs per element are included in the system matrix, I need this number as well. I hope my questions are clear and looking forward for a suggestion on solving it. Best regards -- *Deepak K. Gupta* -- 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.
Re: [deal.II] Bug in sparse_matrix.templates.h
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 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.
[deal.II] Re: Possibly wrong format of msh file
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/0ByiwmbgGTAfqazRERTVEUVJKR1k/view?usp=sharing > Your mesh has too many cells to see anything. You should try to make the mesh as coarse as possible and then look if all the cells look good. I don't remember exactly what this error means but the name of the function that throws, I am going to assume that the ordering of the vertices is not consistent. Can you load the two meshes separately? If not, looking at the one which doesn't work could give you an indication on what is happening. Best, Bruno -- 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.
Re: [deal.II] Re: Possibly wrong format of msh file
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 in refining the coarse meshes, > which are not structured? I mean mainly, what the constraint, > that in dealii triangulations the refinement step can be only 1 > between neighbouring meshes. How to establish that for the > case of unstructured meshes. Is there some tutorial step for that? There is no difference if the mesh is structured or not. Deal doesn't know how the mesh looks like. Best, Bruno -- 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.
Re: [deal.II] Re: Possibly wrong format of msh file
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:00 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 in refining the coarse meshes, >> which are not structured? I mean mainly, what the constraint, >> that in dealii triangulations the refinement step can be only 1 >> between neighbouring meshes. How to establish that for the >> case of unstructured meshes. Is there some tutorial step for that? > There is no difference if the mesh is structured or not. Deal doesn't > know how the mesh looks like. > > Best, > > Bruno -- 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.
Re: [deal.II] Re: Possibly wrong format of msh file
> > 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 some helpful features, which diminish this burden? 2016-07-15 17:37 GMT+02:00 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 in refining the coarse meshes, > > which are not structured? I mean mainly, what the constraint, > > that in dealii triangulations the refinement step can be only 1 > > between neighbouring meshes. How to establish that for the > > case of unstructured meshes. Is there some tutorial step for that? > There is no difference if the mesh is structured or not. Deal doesn't > know how the mesh looks like. > > Best, > > Bruno > > -- > 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. > -- 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.
Re: [deal.II] Re: Possibly wrong format of msh file
> > 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 : > 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 > some helpful features, which diminish this burden? > > 2016-07-15 17:37 GMT+02:00 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 in refining the coarse meshes, >> > which are not structured? I mean mainly, what the constraint, >> > that in dealii triangulations the refinement step can be only 1 >> > between neighbouring meshes. How to establish that for the >> > case of unstructured meshes. Is there some tutorial step for that? >> There is no difference if the mesh is structured or not. Deal doesn't >> know how the mesh looks like. >> >> Best, >> >> Bruno >> >> -- >> 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. >> > > -- 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.
Re: [deal.II] Re: Possibly wrong format of msh file
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 only need to use the refine flags for the cells that you want to refine. When you execute the refinement, deal.II will automatically refine extra cells to satisfy the constraint. See the general documentation of Triangulation dealii.org/8.4.1/doxygen/deal.II/classTriangulation.html Bruno -- 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.
Re: [deal.II] Possibly wrong format of msh file
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==0 || n_negative_cells==cells.size() For the record, what this exception means is that there were some cells that had a negative volume (but not *all* cells, in which case deal.II would automatically have inverted every cell). This typically indicates that the generated mesh has a few cells that are degenerate or inverted. The error message is, admittedly, not very useful. I'll write a patch. Best W. -- Wolfgang Bangerth email:bange...@math.tamu.edu www: http://www.math.tamu.edu/~bangerth/ -- 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.
Re: [deal.II] Bug in sparse_matrix.templates.h
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 email:bange...@math.tamu.edu www: http://www.math.tamu.edu/~bangerth/ -- 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.
Re: [deal.II] axis-symmetric problem
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 tutorial, but it is easily possible. deal.II does not care what the first and second and third coordinate axes represent. If you want to call the first 'r' and the second 'z', then that is just fine with deal.II. What matters is how you formulate the bilinear form and right hand side, as well as the operators in them. For example, if you have a problem in which you only have a dependence on r and z, but not on the angle phi in a cylindrical coordinate system, then the Laplace equation will read in weak form \int (dv/dr du/dr + dv/dz du/dz) * 2*pi r dr dz Note that 2*pi*r factor in the integral. As long as you make sure all of the terms you integrate have the appropriate scaling, you will be fine. Best W. -- Wolfgang Bangerth email:bange...@math.tamu.edu www: http://www.math.tamu.edu/~bangerth/ -- 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.
[deal.II] Dirichlet constraint to a specific node
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 dealii. Thank you, Anup. -- 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.
[deal.II] Re: geometry
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 = triangulation.begin_active(),
endc = triangulation.end();
std::set::active_cell_iterator> remove;
// Triangulation<2>::active_cell_iterator cell;
for (cell=triangulation.begin_active(); cell!=triangulation.end(); ++cell)
{
Point<2> v=cell->center();
// If it is in the bottom part of the geometry, remove it
if (v(1) <= 0)
remove.insert(cell);
}
// remove the marked cells
GridGenerator::create_triangulation_with_removed_cells(triangulation,
removal, triangulation);
Unfortunately, It does not create the right geometry. I do really
appreciate your kindness if you let me know which parts of this code go
wrong. Looking forward to hearing from you.
Best,
Benhour
On Wednesday, July 13, 2016 at 6:15:04 AM UTC-5, Jean-Paul Pelteret wrote:
>
> Dear Benhour,
>
> You could use one of the GridGenerator options to create a half circle,
> and then use GridGenerator::create_triangulation_with_removed_cells
>
> to
> remove the excess cells that yo udon't want.
>
> Regards,
> J-P
>
> On Wednesday, July 13, 2016 at 3:40:02 AM UTC+2, benhour@gmail.com
> wrote:
>>
>> Dear All,
>> I want to create a quarter of a circle in deal ii. In this case, I could
>> create half of a circle with half_hyper_ball. I have tried
>> GridTools::delete_unused_vertices or GridTools::transformation to remove or
>> shift the same vertices to eliminate half of the circle. It would be very
>> kind of you if you let me know how to define this geometry in Dealii. It
>> should be noted that I also used quarter_hyper_shell and put inner radius
>> to zero but it did not work for this case.
>>
>> Best,
>> Benhour
>>
>
--
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.
