[deal.II] interpolating boundary values from a previously found FE solution

2018-03-29 Thread Mustafa Aggul
Hello all, We usually interpolate boundary values from an explicitly defined function. I have a good understanding of how to deal with this kind of cases with a constrain matrix. However, in my current project I need to interpolate boundary values from a previously found FE solution. In other

Re: [deal.II] functions to measure quality of mesh?

2018-03-29 Thread Wolfgang Bangerth
On 03/29/2018 07:23 PM, Weixiong Zheng wrote: Just a quick question: are there functions to measure mesh quality? I didn't see anything in FAQ and only saw one old post without an answer. There isn't anything immediately that produces mesh quality measures, but all of the pieces are there.

[deal.II] functions to measure quality of mesh?

2018-03-29 Thread Weixiong Zheng
Hi Just a quick question: are there functions to measure mesh quality? I didn't see anything in FAQ and only saw one old post without an answer. Thanks in advance! Weixiong -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Weixiong Zheng
Both idea work fine now. Thank you so much! 在 2018年3月29日星期四 UTC-7上午11:31:23,Jean-Paul Pelteret写道: > > Hi Weixiong, > > So, have you looked at the FAQ as the error message suggests? There’s also > actually a discussion on this point in the very first tutorial. Here’s the > relevant snippet for

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Weixiong Zheng
Thanks, the former doesn't work but the latter way works fine. 在 2018年3月29日星期四 UTC-7上午11:31:23,Jean-Paul Pelteret写道: > > Hi Weixiong, > > So, have you looked at the FAQ as the error message suggests? There’s also > actually a discussion on this point in the very first tutorial. Here’s the >

Re: [deal.II] Load force is different with different number of processes

2018-03-29 Thread Wolfgang Bangerth
On 03/28/2018 03:48 PM, Qing Yin wrote: I would like to post that snippet of code. To be honest, I verified the solution in a naive way, that is, I compare the two different number of processes solutions in paraview. The comparison(including data range, plot curve over time at a point) shows

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Jean-Paul Pelteret
Hi Weixiong, So, have you looked at the FAQ as the error message suggests? There’s also actually a discussion on this point in the very first tutorial. Here’s the relevant snippet for you: > At this point, all objects created in this function will be destroyed in > reverse order.

Re: [deal.II] Any easy way to output cellwise average solution

2018-03-29 Thread Weixiong Zheng
Thanks, this looks much prettier! 在 2018年3月27日星期二 UTC-7下午5:04:18,Jean-Paul Pelteret写道: > > Dear Weixiong, > > There’s an interpolation function that should do this for you: > > https://dealii.org/8.5.1/doxygen/deal.II/namespaceVectorTools.html#a5e3af70a47cedfaf361cf5c621e94e3d > The associated

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Weixiong Zheng
Thanks, I changed the number to 2 and straight boundary comes out. But I got errors that I could not understand saying An error occurred in line <128> of file in function void dealii::Subscriptor::check_no_subscribers() const The violated condition was: counter == 0

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Weixiong Zheng
Hi Thanks so much for the fast reply. I actually tried changing to: const Point<2> center (0,0); const SphericalManifold<2> boundary_description(center); //const StraightBoundary<2> boundary_description; //tria.set_all_manifold_ids_on_boundary(0); tria.set_manifold (2,

Re: [deal.II] set manifold for merged boundaries

2018-03-29 Thread Weixiong Zheng
Hi It does give me straight boundary after changing to: const Point<2> center (0,0); const SphericalManifold<2> boundary_description(center); //const StraightBoundary<2> boundary_description; //tria.set_all_manifold_ids_on_boundary(0); tria.set_manifold (2, boundary_description);

Re: [deal.II] Re: Any easy way to output cellwise average solution

2018-03-29 Thread Weixiong Zheng
Nice catch! I forgot to mention it is bilinear. But thanks for all the explanation. 在 2018年3月28日星期三 UTC-7下午12:42:27,Wolfgang Bangerth写道: > > > > | > > // The vector I want to output is the solution for vector-valued > problem. > > // I want to output the cell-average value for the 0th