Re: [deal.II] How to use CellId class

2019-10-19 Thread Zhidong Brian Zhang
Hi, Dr. W. Bangerth, I have tried the your suggestion that is to use FE_DGQ(0) to build another finite element field. It works even though there is till some small problems. However, I prefer to double check the implementation with you: I built two fem field objects. FESystem fe; // for

Re: [deal.II] How to use CellId class

2019-10-19 Thread Zhidong Brian Zhang
Hi, Dr. W. Bangerth, Thank you so much for your timely help! I will try it and come back if I have further questions. Best regards, Zhidong Brian Zhang On Saturday, October 19, 2019 at 10:56:47 AM UTC-4, Wolfgang Bangerth wrote: > > > Zhidong, > > > In short, my question is how to get the

Re: [deal.II] How to use CellId class

2019-10-19 Thread Wolfgang Bangerth
Zhidong, > In short, my question is how to get the cell_id properly by using CellId > class ? > > I want to assign each cell with a different value (its pesudo density in > topology optimization). All the values are stored in a cell-based vector x, > which is initialized as follows, > >

Re: [deal.II] Loading a .msh mesh of big size into parallel::distributed::triangulation using p4est

2019-10-19 Thread Phạm Ngọc Kiên
Hi praveen, If you run with 1 process(rank = 0) with ./step-40, it's fine. I want to run with MPI on multiple processors with mpirun -np 4 ./step-40. In this case the trouble comes. Thank you very much. Best Kien Vào Th 7, 19 thg 10, 2019 vào lúc 13:45 Praveen C đã viết: > Just an observation