Hello Timo, Does the replacement code serve the same purpose as the original code? The purpose of this line is to sync the error computed by each process and summarized in the vector estimated_error_per_cell, and we should get the same estimated_error_per_cell on each process.
I wrote this part of the code myself since there is no tutorial on how to write this part... Step-18 only has a guide on how to write this part with PESTc. Best, Yiyang On Sunday, October 29, 2017 at 6:03:12 PM UTC-5, Timo Heister wrote: > > Yiyang, > > I am not quite sure why, but the code > > estimated_error_per_cell = distributed_error_per_cell; > > generates an incorrect result. You can replace it with > > for(auto cell = tria.begin_active(); cell != tria.end(); ++cell) > if(cell->is_locally_owned()) > estimated_error_per_cell(cell->active_cell_index()) = > distributed_error_per_cell(cell->active_cell_index()); > > Can you try if that works for you? > > I will continue to investigate what is going on. > > > On Sat, Oct 28, 2017 at 1:00 PM, Yiyang Zhang <zyyce...@gmail.com > <javascript:>> wrote: > > Hello Prof. Bangerth, > > > > I have attached a test case. > > > > Best, > > Yiyang > > > > On Friday, October 27, 2017 at 5:12:49 PM UTC-5, Wolfgang Bangerth > wrote: > >> > >> On 10/27/2017 02:13 PM, Yiyang Zhang wrote: > >> > > >> > Yes I think I am setting them in a way that is exactly same > throughout > >> > all processes. > >> > Since I can check the n_active_cells() before the refinement, and > also > >> > the refine_flags and coarsen_flags for each process. They are the > same > >> > for each process. But after the refinement, the number of > active_cells > >> > suddenly becomes different. > >> > > >> > I attached the refine_mesh() code. > >> > >> Can you fabricate a complete testcase that we can run that shows this > >> phenomenon? > >> > >> Best > >> W. > >> > >> -- > >> > ------------------------------------------------------------------------ > >> Wolfgang Bangerth email: bang...@colostate.edu > >> www: > https://urldefense.proofpoint.com/v2/url?u=http-3A__www.math.colostate.edu_-7Ebangerth_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=Avq4CtF3shqaTzaJnB2SqcYaMgNDwAkFHlvEyo45f24&s=GT8h1eKbhfid-4BHzYVtaoniDOC2gT0gnXEmZiLsVLY&e= > > > > > -- > > The deal.II project is located at > https://urldefense.proofpoint.com/v2/url?u=http-3A__www.dealii.org_&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=Avq4CtF3shqaTzaJnB2SqcYaMgNDwAkFHlvEyo45f24&s=3zUo9GeZ-rDnZ1tlZXSuEvw4gcM279OP2eh_1qCYcZs&e= > > > For mailing list/forum options, see > > > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_forum_dealii-3Fhl-3Den&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=Avq4CtF3shqaTzaJnB2SqcYaMgNDwAkFHlvEyo45f24&s=4KjaGXKFUbi23g-CuQL1kAAPFr7IXyw1mGiKnk6eROs&e= > > > --- > > You received this message because you are subscribed to the Google > Groups > > "deal.II User Group" group. > > To unsubscribe from this group and stop receiving emails from it, send > an > > email to dealii+un...@googlegroups.com <javascript:>. > > For more options, visit > https://urldefense.proofpoint.com/v2/url?u=https-3A__groups.google.com_d_optout&d=DwIBaQ&c=Ngd-ta5yRYsqeUsEDgxhcqsYYY1Xs5ogLxWPA_2Wlc4&r=4k7iKXbjGC8LfYxVJJXiaYVu6FRWmEjX38S7JmlS9Vw&m=Avq4CtF3shqaTzaJnB2SqcYaMgNDwAkFHlvEyo45f24&s=Vsm-Y6HviYtRxDZP6E1hzM3JP9jGx_a8utwaysqI8a0&e= > > . > > > > -- > Timo Heister > http://www.math.clemson.edu/~heister/ > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.
