Dear Denis,
Thanks for your replies.
I do not have some data attached to ghost cells. But I have a variable that
changes with time and its change depends on the data of neighbor cells.
For example,
const double value1 = 1, value2 =-1;
for (cell = tr.begin_active(); cell != tr.end(); ++cell)
if (cell -> is_locally_owned()){
bool flag = false;
std::vector<std::shared_ptr<MyData> > qpd = data_storage.template
get_data<MyData>(cell);
for (unsigned int q = 0; q < n_q_points; ++q)
if(qpd->value[q]> 0) {flag = true; break;}
if(!flag)
// finding the qpd values of neighbor cells and check all the values
//if one of them is positive, then flag = true;
//Setting the gpd values of present cell.
for (unsigned int q = 0; q < n_q_points; ++q)
{
if(flag) qpd->value[q] = value1;
else qpd->value[q] = value2;
}
….
}
Thanks!
--
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.
