GrowingVectorMemory<Vector<double> > mem;
ProductSparseMatrix<number, number> product(m1, m2, mem);
+
product.vmult(w,v);
-
for (unsigned int i=0;i<w.size();++i)
deallog << ' ' << w(i);
deallog << std::endl;
-
+
+ product.vmult_add(w,v);
+ for (unsigned int i=0;i<w.size();++i)
+ deallog << ' ' << w(i);
+ deallog << std::endl;
+
product.Tvmult(v,w);
+ for (unsigned int i=0;i<v.size();++i)
+ deallog << ' ' << v(i);
+ deallog << std::endl;
+ product.Tvmult_add(v,w);
for (unsigned int i=0;i<v.size();++i)
deallog << ' ' << v(i);
deallog << std::endl;
DEAL::Sizes 4 2
DEAL:: 1. 2. 3. 4.
DEAL:: -14. -8.
-DEAL:: -86. -44. 86. 44.
-DEAL::GrowingVectorMemory:Overall allocated vectors: 2
+DEAL:: -28. -16.
+DEAL:: -172. -88. 172. 88.
+DEAL:: -344. -176. 344. 176.
+DEAL::GrowingVectorMemory:Overall allocated vectors: 4
DEAL::GrowingVectorMemory:Maximum allocated vectors: 1
In the beginning the Universe was created. This has made a lot of
people very angry and has been widely regarded as a bad move.
Douglas Adams