Re: [deal.II] A bug of visualization with VISIT

2017-02-24 Thread Jiaqi ZHANG
Hey Prof. Bangerth, Thanks for your help again, I understand now :) Best, Jiaqi 在 2017年2月22日星期三 UTC-5下午11:58:06,Wolfgang Bangerth写道: > > On 02/22/2017 06:41 AM, Jiaqi ZHANG wrote: > > Soluiton2 is fine. For Solution1, I checked the result with MATLAB, the > > contour plot is right. > > I also

[deal.II] 3D mesh with CAD boundary

2017-02-24 Thread Ehsan Nedaaee
Dear Dealii Users I am going to numerically solve a system of PDEs in a 3D domain with complex boundaries. The domain is plotted using CAD, and I would like to import such domain into dealii, then using deal procedures to make 3D mesh and refine it. I think I should first generate a mesh on t

Re: [deal.II] Re: Access specific element within a distributed triangulation

2017-02-24 Thread Wolfgang Bangerth
Seyed, Generally, in FEM you number elements like in the below picture: No, that is patently wrong. First, you are only considering rectangular domains subdivided

[deal.II] Re: Access specific element within a distributed triangulation

2017-02-24 Thread Jean-Paul Pelteret
Hi Seyed, Nope, like I said, there's no such thing as a cell number in deal.II (at least in the sense that a cell knows its own "id" number). The cell->id() function returns a CellId

[deal.II] Re: Access specific element within a distributed triangulation

2017-02-24 Thread 'Seyed Ali Mohseni' via deal.II User Group
Thanks for your suggestions. Generally, in FEM you number elements like in the below picture: I understand of course, that the cells (elements) in deal.II are para

[deal.II] New: Physics module and solid mechanics code gallery example

2017-02-24 Thread Jean-Paul Pelteret
Dear all, It is my pleasure to announce the inception of a Physics module in the deal.II project. This module, and its associated namespaces, is dedicated providing some core functionality for computations commonly required

Re: [deal.II] Reading user_indices from a vtk file.

2017-02-24 Thread Wolfgang Bangerth
On 02/24/2017 04:15 AM, 987093...@qq.com wrote: In the file,I have set up user_indices for each grid cell, but when I use GridIn to generate a mesh, there is no user indices that I have saved on each grid cell. What should i do? I don't actually see anything in the vtk file that would desc

Re: [deal.II] Re: Access specific element within a distributed triangulation

2017-02-24 Thread Timo Heister
To add to this, there might also be no reason for cell "4" (whatever that means) to be available on the current processor. The mesh in a parallel computation is "distributed" after all. On Fri, Feb 24, 2017 at 8:01 AM, Jean-Paul Pelteret wrote: > Dear Seyed, > > There are no cell "numbers" in dea

[deal.II] Re: Access specific element within a distributed triangulation

2017-02-24 Thread Jean-Paul Pelteret
Dear Seyed, There are no cell "numbers" in deal.II, only CellIDs to give some unique identifier to each cell. So its unclear to me as to exactly what you're trying to achieve here. Are you wanting the fourth cell iterator for the

[deal.II] Reading user_indices from a vtk file.

2017-02-24 Thread 987093441
hey guys, In the file,I have set up user_indices for each grid cell, but when I use GridIn to generate a mesh, there is no user indices that I have saved on each grid cell. What should i do? Here is my program. Triangulation<2>triangulation; GridIn<2>grid_in; grid_i

[deal.II] Access specific element within a distributed triangulation

2017-02-24 Thread 'Seyed Ali Mohseni' via deal.II User Group
Hi, Is it possible to access a specific element within a parallel distributed triangulation by means of the following procedure for instance: typename DoFHandler::active_cell_iterator cell = dof_handler.begin_active(), endc = dof_handler.end(); for (; cell != endc; ++cell) { if (cell->is_