[deal.II] is there a way to refine mesh only in one direction

2017-03-29 Thread Jaekwang Kim
Hi all, I am running deal.ii for study non-newtonian fluids. I have a mesh which is attached. I am using AMR, but, when I use just conventional function, "cell->refine_flag()" for refinement, the flagged cell divide into 4. but as my problem only need high resolution in r-direction, I wa

[deal.II] Re: is there a way to refine mesh only in one direction

2017-03-29 Thread Bruno Turcksin
Jaekwang, you can do anisotropic refinement (http://dealii.org/developer/doxygen/deal.II/classCellAccessor.html#afb6cc537720a5b6381c237abe0887de2) but only in serial. So depending if your code is parallel or not this may be what you are looking for. Best, Bruno On Wednesday, March 29, 2017

Re: [deal.II] Re: is there a way to refine mesh only in one direction

2017-03-29 Thread Timo Heister
To add to Bruno's answer: If you don't want to use anisotropic refinement you have two other options: 1. Just accept that you get additional refinement in areas/directions you don't particularly care about. If you have a realistic computation, you likely need to refine many more times and thus the

[deal.II] Fe_values->shape_grad() wrt to reference mesh

2017-03-29 Thread RAJAT ARORA
Hello all, I am using deal.ii to solve a 3D solid mechanics problem. I am using p:d:triangulation and Petsc in the code. I am updating my mesh after every time-increment. So, cell->vertex(v) contains the current coordinates of the mesh. I know the coordinates of the cell (and whole mesh) at t

[deal.II] Re: is there a way to refine mesh only in one direction

2017-03-29 Thread Jaekwang Kim
thanks, I would need it and I am still running serially I am sorry but I am not that good at c++, so I was bit confused how to use this function. I have just set flag cell like *cell->set_refine_flag ();* Reading the document pages, you linked, I can understand that I should put some argum

Re: [deal.II] Re: is there a way to refine mesh only in one direction

2017-03-29 Thread Wolfgang Bangerth
On 03/29/2017 06:37 PM, Jaekwang Kim wrote: I am sorry but I am not that good at c++, so I was bit confused how to use this function. I have just set flag cell like *cell->set_refine_flag ();* Reading the document pages, you linked, I can understand that I should put some argument inside par