->cell_dof_indices_cache[accessor.present_index * n_dofs];
// distribute cell vector
- global_destination.add(n_dofs, dofs,
- &(*const_cast<dealii::Vector<number>*>(&local_source))(0));
+ global_destination.add(n_dofs, dofs,local_source.begin());
}
template <typename Number>
inline
void
-BlockVectorBase<VectorType> ::add (const std::vector<unsigned int> &indices,
- const std::vector<Number> &values)
+BlockVectorBase<VectorType>::add (const std::vector<unsigned int> &indices,
+ const std::vector<Number> &values)
{
Assert (indices.size() == values.size(),
ExcDimensionMismatch(indices.size(), values.size()));
template <typename Number>
inline
void
-BlockVectorBase<VectorType> ::add (const std::vector<unsigned int> &indices,
- const Vector<Number> &values)
+BlockVectorBase<VectorType>::add (const std::vector<unsigned int> &indices,
+ const Vector<Number> &values)
{
Assert (indices.size() == values.size(),
ExcDimensionMismatch(indices.size(), values.size()));
template <typename Number>
inline
void
-BlockVectorBase<VectorType> ::add (const unsigned int n_indices,
- const unsigned int *indices,
- const Number *values)
+BlockVectorBase<VectorType>::add (const unsigned int n_indices,
+ const unsigned int *indices,
+ const Number *values)
{
for (unsigned int i=0; i<n_indices; ++i)
(*this)(indices[i]) += values[i];
{
Assert (indices.size() == values.size(),
ExcMessage ("Function called with arguments of different sizes"));
- do_set_add_operation(indices.size(), &indices[0],
- &(*const_cast<dealii::Vector<PetscScalar>*>(&values))(0),
- true);
+ do_set_add_operation(indices.size(), &indices[0], values.begin(), true);
}