void
ConstraintMatrix::distribute (TrilinosWrappers::MPI::Vector &vec) const
{
+ Assert (sorted==true, ExcMatrixIsClosed());
+
//TODO: not implemented yet, we need to fix
//LocalRange() first to only include
//"owned" indices. For this we need to keep
void
ConstraintMatrix::distribute (TrilinosWrappers::MPI::BlockVector &vec) const
{
+ Assert (sorted==true, ExcMatrixIsClosed());
+
IndexSet my_indices (vec.size());
for (unsigned int block=0; block<vec.n_blocks(); ++block)
{
void
ConstraintMatrix::distribute (PETScWrappers::MPI::Vector &vec) const
{
+ Assert (sorted==true, ExcMatrixIsClosed());
+
typedef std::vector<ConstraintLine>::const_iterator constraint_iterator;
ConstraintLine index_comparison;
index_comparison.line = vec.local_range().first;
void
ConstraintMatrix::distribute (PETScWrappers::MPI::BlockVector &/*vec*/) const
{
+ Assert (sorted==true, ExcMatrixIsClosed());
AssertThrow (false, ExcNotImplemented());
}