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

2018-03-30 Thread Jean-Paul Pelteret
Hi Weixiong, Great, I’m glad that you managed to fix your problem! Best, Jean-Paul > On 29 Mar 2018, at 23:47, Weixiong Zheng wrote: > > 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

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 yo

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 > rele

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. Unfortunately,

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 Additional

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, boundary_descript

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); f

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

2018-03-28 Thread Jean-Paul Pelteret
Hi Weixiong, The issue that you’re seeing here is actually related to a discussion that we’re currently having regarding the behaviour of boundaries and manifolds (and will hopefully be fixed soon in the developer version). You’re hitting the issue discussed here https://github.com/dealii/deali

[deal.II] set manifold for merged boundaries

2018-03-28 Thread Weixiong Zheng
Hi, I just started to try manifold and tria merging functionalities. My question comes from the problem that I have to merge two curved triangulations. When refining, how do I make sure the merged boundaries are still curved? I currently set manifold before merging triangulation but it seems n