//
// ---------------------------------------------------------------------
-// Test that it is possible to use a operator*() of LinearOperator object for
-// Trilinos matrices and vectors
+// Test that it is possible to use a PackagedOperation created by
+// operator*() of a LinearOperator object for Trilinos matrices and vectors
#include "../tests.h"
typedef TrilinosWrappers::MPI::Vector vector_t;
typedef TrilinosWrappers::SparseMatrix matrix_t;
- matrix_t a;
+ matrix_t a(5U, 5U, 3U);
+ a.compress (VectorOperation::add);
auto op_a = linear_operator<vector_t>(a);
vector_t u,res;
+ op_a.reinit_domain_vector(u, false);
res = op_a * u;
- // ^^ this was not working, whereas
- // op_a.vmult(res,u) did.
+ // ^^ this was not working, whereas op_a.vmult(res,u) did.
deallog << "OK" << std::endl;