IndexSet
VectorBase::locally_owned_elements() const
{
+ Assert(owned_elements.size()==size(),
+ ExcMessage("The locally owned elements have not been properly initialized!"
+ " This happens for example if this object has been initialized"
+ " with exactly one overlapping IndexSet."));
return owned_elements;
}
if (v.nonlocal_vector.get() != 0)
nonlocal_vector.reset(new Epetra_MultiVector(v.nonlocal_vector->Map(), 1));
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
last_action = Zero;
vector.reset (new Epetra_FEVector(*v.vector));
has_ghosts = v.has_ghosts;
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
vector.reset (new Epetra_FEVector(v.vector->Map()));
has_ghosts = v.has_ghosts;
last_action = Zero;
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
else if (omit_zeroing_entries == false)
{
my_global_elements(v.block(block).vector_partitioner());
for (size_type i=0; i<v.block(block).local_size(); ++i)
global_ids[added_elements++] = glob_elements[i] + block_offset;
- owned_elements.add_indices(v.block(block).locally_owned_elements(),
+ owned_elements.add_indices(v.block(block).owned_elements,
block_offset);
block_offset += v.block(block).size();
}
vector.reset (new Epetra_FEVector(*v.vector));
last_action = Zero;
has_ghosts = v.has_ghosts;
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
if (v.nonlocal_vector.get() != 0)
v.vector->Map().IndexBase(),
v.vector->Comm());
vector.reset (new Epetra_FEVector(map));
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
else if (omit_zeroing_entries)
{
v.vector->Map().IndexBase(),
v.vector->Comm());
vector.reset (new Epetra_FEVector(map));
- owned_elements = v.locally_owned_elements();
+ owned_elements = v.owned_elements;
}
const int ierr = vector->Update(1.0, *v.vector, 0.0);