>
> You want to use cell->active_cell_index() as the index into the vector. The
> vector should have
>triangulation.n_active_cells()
> as its size. This corresponds to the *local* number of active cells,
> including
> ghost and artificial cells (for which vector entries are then just
>
On 7/7/21 12:08 AM, vachanpo...@gmail.com wrote:
Vector temp_vec(gh_vec); // for data output
for(auto : dof_handler.active_cell_iterators()){
if(!(cell->is_locally_owned())) continue;
temp_vec[cell->index()] = gh_vec[cell->global_active_cell_index()];
}
data_out.add_data_vector(temp_vec,
Dear all,
I am having some confusion regarding the index used in DataOut for a cell
vector.
I have created a distributed vector for storing cell data using
p::d::Triangulation::global_active_cell_index_partitioner()