From: Jean-Paul Pelteret Date: Wed, 1 Feb 2017 15:05:42 +0000 (+0100) Subject: Update method used to nullify Epetra_MultiVectors X-Git-Tag: v8.5.0-rc1~161^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F3890%2Fhead;p=dealii.git Update method used to nullify Epetra_MultiVectors --- diff --git a/source/lac/trilinos_sparse_matrix.cc b/source/lac/trilinos_sparse_matrix.cc index 7b598a4dc5..c8eaa7f66b 100644 --- a/source/lac/trilinos_sparse_matrix.cc +++ b/source/lac/trilinos_sparse_matrix.cc @@ -2805,27 +2805,37 @@ namespace TrilinosWrappers return_op.vmult = [](Range &tril_dst, const Domain &) { - tril_dst.Scale(0); + const int ierr = tril_dst.PutScalar(0.0); + + AssertThrow (ierr == 0, ExcTrilinosError(ierr)); }; return_op.Tvmult = [](Domain &tril_dst, const Range &) { - tril_dst.Scale(0); + const int ierr = tril_dst.PutScalar(0.0); + + AssertThrow (ierr == 0, ExcTrilinosError(ierr)); }; return_op.inv_vmult = [](Domain &tril_dst, const Range &) { AssertThrow(false, ExcMessage("Cannot compute inverse of null operator")); - tril_dst.Scale(0); + + const int ierr = tril_dst.PutScalar(0.0); + + AssertThrow (ierr == 0, ExcTrilinosError(ierr)); }; return_op.inv_Tvmult = [](Range &tril_dst, const Domain &) { AssertThrow(false, ExcMessage("Cannot compute inverse of null operator")); - tril_dst.Scale(0); + + const int ierr = tril_dst.PutScalar(0.0); + + AssertThrow (ierr == 0, ExcTrilinosError(ierr)); }; return return_op;