From: Timo Heister Date: Thu, 30 May 2013 20:06:25 +0000 (+0000) Subject: add PETSc matrix copy_from X-Git-Tag: v8.0.0~120^2~54 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=4513a2cea9dac91329061547a16c2b89da4c680e;p=dealii.git add PETSc matrix copy_from git-svn-id: https://svn.dealii.org/branches/branch_unify_linear_algebra@29670 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/include/deal.II/lac/petsc_parallel_sparse_matrix.h b/deal.II/include/deal.II/lac/petsc_parallel_sparse_matrix.h index 7eb8b9d087..93c04f4371 100644 --- a/deal.II/include/deal.II/lac/petsc_parallel_sparse_matrix.h +++ b/deal.II/include/deal.II/lac/petsc_parallel_sparse_matrix.h @@ -294,6 +294,13 @@ namespace PETScWrappers */ SparseMatrix &operator = (const value_type d); + + /** + * Make a copy of the PETSc matrix @p other. It is assumed that both matrices have + * the same SparsityPattern. + */ + void copy_from(const SparseMatrix & other); + /** * Throw away the present matrix and * generate one that has the same diff --git a/deal.II/source/lac/petsc_parallel_sparse_matrix.cc b/deal.II/source/lac/petsc_parallel_sparse_matrix.cc index 405bc0c2bf..ecc5dac544 100644 --- a/deal.II/source/lac/petsc_parallel_sparse_matrix.cc +++ b/deal.II/source/lac/petsc_parallel_sparse_matrix.cc @@ -99,7 +99,17 @@ namespace PETScWrappers return *this; } + void + SparseMatrix::copy_from (const SparseMatrix &other) + { + if (&other == this) + return; + + this->communicator = other.communicator; + int ierr = MatCopy(other.matrix, matrix, SAME_NONZERO_PATTERN); + AssertThrow (ierr == 0, ExcPETScError(ierr)); + } void SparseMatrix::reinit (const MPI_Comm &communicator,