Re: [deal.II] Unexpected data output with cell data vector

> > 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 > ignored)

Re: [deal.II] Unexpected data output with cell data vector

On 7/7/21 12:08 AM, vachanpo...@gmail.com wrote: Vector temp_vec(gh_vec); // for data output for(auto &cell: 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_v

[deal.II] Unexpected data output with cell data vector

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()