// an IndexSet where every valid index is part of the set. Note that this
// program can only be run sequentially and will throw an exception if used
// in parallel.
- IndexSet eigenfunction_partitioning = complete_index_set(dof_handler.n_dofs ());
+ IndexSet eigenfunction_index_set = dof_handler.locally_owned_dofs ();
eigenfunctions
.resize (parameters.get_integer ("Number of eigenvalues/eigenfunctions"));
for (unsigned int i=0; i<eigenfunctions.size (); ++i)
- eigenfunctions[i].reinit (eigenfunction_partitioning, MPI_COMM_WORLD);
+ eigenfunctions[i].reinit (eigenfunction_index_set, MPI_COMM_WORLD);
eigenvalues.resize (eigenfunctions.size ());
}
system_matrix.reinit (dsp);
complete_system_matrix.reinit (dsp);
- IndexSet solution_partitioning = complete_index_set(dof_handler.n_dofs());
- solution.reinit (solution_partitioning, MPI_COMM_WORLD);
- system_rhs.reinit (solution_partitioning, MPI_COMM_WORLD);
- complete_system_rhs.reinit (solution_partitioning, MPI_COMM_WORLD);
- contact_force.reinit (solution_partitioning, MPI_COMM_WORLD);
+ IndexSet solution_index_set = dof_handler.locally_owned_dofs();
+ solution.reinit (solution_index_set, MPI_COMM_WORLD);
+ system_rhs.reinit (solution_index_set, MPI_COMM_WORLD);
+ complete_system_rhs.reinit (solution_index_set, MPI_COMM_WORLD);
+ contact_force.reinit (solution_index_set, MPI_COMM_WORLD);
// The only other thing to do here is to compute the factors in the $B$
// matrix which is used to scale the residual. As discussed in the
TrilinosWrappers::SparseMatrix mass_matrix;
mass_matrix.reinit (dsp);
assemble_mass_matrix_diagonal (mass_matrix);
- diagonal_of_mass_matrix.reinit (solution_partitioning);
+ diagonal_of_mass_matrix.reinit (solution_index_set);
for (unsigned int j=0; j<solution.size (); j++)
diagonal_of_mass_matrix (j) = mass_matrix.diag_element (j);
}
* entries in Input_Maps. For this non-distributed vector, the %parallel
* partitioning is not used, just the global size of the partitioner.
*/
- explicit BlockVector (const std::vector<Epetra_Map> &partitioner);
+ explicit BlockVector (const std::vector<Epetra_Map> &partitioner DEAL_II_DEPRECATED);
/**
* Constructor. Generate a block vector with as many blocks as there are
* partitioning is not used, just the global size of the partitioner.
*/
explicit BlockVector (const std::vector<IndexSet> &partitioner,
- const MPI_Comm &communicator = MPI_COMM_WORLD);
+ const MPI_Comm &communicator = MPI_COMM_WORLD) DEAL_II_DEPRECATED;
/**
* Copy-Constructor. Set all the properties of the non-%parallel vector to
* those of the given %parallel vector and import the elements.
*/
- BlockVector (const MPI::BlockVector &V);
+ BlockVector (const MPI::BlockVector &V) DEAL_II_DEPRECATED;
/**
* Copy-Constructor. Set all the properties of the vector to those of the
* given input vector and copy the elements.
*/
- BlockVector (const BlockVector &V);
+ BlockVector (const BlockVector &V) DEAL_II_DEPRECATED;
/**
* Creates a block vector consisting of <tt>num_blocks</tt> components,
* but there is no content in the individual components and the user has
* to fill appropriate data using a reinit of the blocks.
*/
- explicit BlockVector (const size_type num_blocks);
+ explicit BlockVector (const size_type num_blocks) DEAL_II_DEPRECATED;
/**
* Constructor. Set the number of blocks to <tt>n.size()</tt> and
*
* References BlockVector.reinit().
*/
- explicit BlockVector (const std::vector<size_type> &N);
+ explicit BlockVector (const std::vector<size_type> &N) DEAL_II_DEPRECATED;
/**
* Constructor. Set the number of blocks to <tt>n.size()</tt>. Initialize
template <typename InputIterator>
BlockVector (const std::vector<size_type> &n,
const InputIterator first,
- const InputIterator end);
+ const InputIterator end) DEAL_II_DEPRECATED;
/**
* Destructor. Clears memory
<< "local_size = global_size is a necessary condition, but"
<< arg1 << " != " << arg2 << " was given!");
- } DEAL_II_DEPRECATED;
+ };
* function <tt>reinit()</tt> will have to give the vector the correct
* size.
*/
- Vector ();
+ Vector () DEAL_II_DEPRECATED;
/**
* This constructor takes as input the number of elements in the vector.
*/
- explicit Vector (const size_type n);
+ explicit Vector (const size_type n) DEAL_II_DEPRECATED;
/**
* This constructor takes as input the number of elements in the vector.
* ignored, the only thing that matters is the size of the index space
* described by this argument.
*/
- explicit Vector (const Epetra_Map &partitioning);
+ explicit Vector (Const Epetra_Map &partitioning) DEAL_II_DEPRECATED;
/**
* This constructor takes as input the number of elements in the vector.
* size of the index space described by this argument.
*/
explicit Vector (const IndexSet &partitioning,
- const MPI_Comm &communicator = MPI_COMM_WORLD);
+ const MPI_Comm &communicator = MPI_COMM_WORLD) DEAL_II_DEPRECATED;
/**
* This constructor takes a (possibly parallel) Trilinos Vector and
* generates a localized version of the whole content on each processor.
*/
- explicit Vector (const VectorBase &V);
+ explicit Vector (const VectorBase &V) DEAL_II_DEPRECATED;
/**
* Copy-constructor from deal.II vectors. Sets the dimension to that of
* the given vector, and copies all elements.
*/
template <typename Number>
- explicit Vector (const dealii::Vector<Number> &v);
+ explicit Vector (const dealii::Vector<Number> &v) DEAL_II_DEPRECATED;
/**
* Reinit function that resizes the vector to the size specified by
* thus an empty function.
*/
void update_ghost_values () const;
- } DEAL_II_DEPRECATED;
+ };