75 cell_iter != rCellPopulation.
End();
78 auto p_cell_edge_model =
static_cast<CellSrnModel*
>(cell_iter->GetSrnModel());
81 std::vector<double> notch_vec;
82 std::vector<double> delta_vec;
83 for (
unsigned i = 0 ; i < p_cell_edge_model->GetNumEdgeSrn(); i++)
85 boost::shared_ptr<DeltaNotchEdgeSrnModel> p_model
86 = boost::static_pointer_cast<DeltaNotchEdgeSrnModel>(p_cell_edge_model->GetEdgeSrn(i));
87 double this_delta = p_model->GetDelta();
88 double this_notch = p_model->GetNotch();
89 delta_vec.push_back(this_delta);
90 notch_vec.push_back(this_notch);
93 cell_iter->GetCellEdgeData()->SetItem(
"edge notch", notch_vec);
94 cell_iter->GetCellEdgeData()->SetItem(
"edge delta", delta_vec);
98 cell_iter->GetCellData()->SetItem(
"interior delta", 0);
99 cell_iter->GetCellData()->SetItem(
"interior notch", 0);
104 cell_iter != rCellPopulation.
End();
107 auto p_cell_edge_model =
static_cast<CellSrnModel*
>(cell_iter->GetSrnModel());
108 const unsigned num_cell_edges = p_cell_edge_model->
GetNumEdgeSrn();
109 std::vector<double> neigh_mean_delta(num_cell_edges);
111 for (
unsigned i=0; i<num_cell_edges; ++i)
115 double mean_delta = 0;
116 for (
auto neighbourIndex: elemNeighbours)
119 std::vector<double> neighbour_delta_vec = neighbourCell->GetCellEdgeData()->GetItem(
"edge delta");
120 mean_delta += neighbour_delta_vec[neighbourIndex.second];
122 if (elemNeighbours.size() > 0)
124 mean_delta = mean_delta/elemNeighbours.size();
126 neigh_mean_delta[i] = mean_delta;
128 cell_iter->GetCellEdgeData()->SetItem(
"neighbour delta", neigh_mean_delta);