// ---------------------------------------------------------------------
-// Test vmult and Tvmult of PointerMatrix and TransposeMatrix
+// Test vmult and Tvmult of PointerMatrix
#include "../tests.h"
#include <deal.II/base/logstream.h>
-#include <deal.II/lac/pointer_matrix.h>
-#include <deal.II/lac/transpose_matrix.h>
#include <deal.II/lac/full_matrix.h>
+#include <deal.II/lac/pointer_matrix.h>
#include <deal.II/lac/vector.h>
#include <fstream>
A(i,j) = ++k;
PointerMatrix<FullMatrix<double>, Vector<double> > P(&A, "P");
- TransposeMatrix<FullMatrix<double>, Vector<double> > T(&A, "T");
Vector<double> x(5);
Vector<double> y(5);
diff = y;
diff.add(-1., y2);
deallog << "P.vmult: diff " << diff.l2_norm() << std::endl;
- T.Tvmult(y2,x);
- diff = y;
- diff.add(-1., y2);
- deallog << "T.Tvmult: diff " << diff.l2_norm() << std::endl;
A.Tvmult(y,x);
P.Tvmult(y2,x);
diff = y;
diff.add(-1., y2);
deallog << "P.Tvmult: diff " << diff.l2_norm() << std::endl;
- T.vmult(y2,x);
- diff = y;
- diff.add(-1., y2);
- deallog << "T.vmult: diff " << diff.l2_norm() << std::endl;
-
-
}
}
DEAL::P.vmult: diff 0.00000
-DEAL::T.Tvmult: diff 0.00000
DEAL::P.Tvmult: diff 0.00000
-DEAL::T.vmult: diff 0.00000
DEAL::P.vmult: diff 0.00000
-DEAL::T.Tvmult: diff 0.00000
DEAL::P.Tvmult: diff 0.00000
-DEAL::T.vmult: diff 0.00000
DEAL::P.vmult: diff 0.00000
-DEAL::T.Tvmult: diff 0.00000
DEAL::P.Tvmult: diff 0.00000
-DEAL::T.vmult: diff 0.00000
DEAL::P.vmult: diff 0.00000
-DEAL::T.Tvmult: diff 0.00000
DEAL::P.Tvmult: diff 0.00000
-DEAL::T.vmult: diff 0.00000
DEAL::P.vmult: diff 0.00000
-DEAL::T.Tvmult: diff 0.00000
DEAL::P.Tvmult: diff 0.00000
-DEAL::T.vmult: diff 0.00000