Update method used to nullify Epetra_MultiVectors

diff --git a/source/lac/trilinos_sparse_matrix.cc b/source/lac/trilinos_sparse_matrix.cc
index 7b598a4dc5486ed43feff460591227ea1b887111..c8eaa7f66bf402efef5c239d5e39f7ee0105a05d 100644 (file)
--- 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;