deallog.threshold_double(1.e-10);
const double Adata[] =
- { 1, 2, 3, 4, 5, 6 };
+ { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
- FullMatrix<double> A(2, 3);
+ FullMatrix<double> A(3, 3);
A.fill(Adata);
- Vector<double> V(3);
- V(0) = 1;
- V(1) = 2;
- V(2) = 3;
+ Vector<double> V1(3);
+ Vector<double> V2(3);
- checkVmult<double>(A, V, false);
- checkVmult<double>(A, V, true);
+ V1(0) = V2(0) = 1;
+ V1(1) = V2(1) = 2;
+ V1(2) = V2(2) = 3;
+
+ checkVmult<double>(A, V1, false);
+ checkVmult<double>(A, V2, true);
}
DEAL::vmult
-DEAL::1.0000 28.0000
+DEAL::1.0000 28.0000 43.0000
DEAL::vmult
-DEAL::1.0000 32.0000
+DEAL::1.0000 32.0000 50.0000
deallog.threshold_double(1.e-10);
const double Adata[] =
- { 1, 2, 3, 4, 5, 6 };
+ { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
- FullMatrix<double> A(2, 3);
+ FullMatrix<double> A(3, 3);
A.fill(Adata);
- Vector<double> V(2);
- V(0) = 1;
- V(1) = 2;
+ Vector<double> V1(3);
+ Vector<double> V2(3);
- checkTvmult<double>(A, V, false);
- checkTvmult<double>(A, V, true);
+ V1(0) = V2(0) = 1;
+ V1(1) = V2(1) = 2;
+ V1(2) = V2(2) = 3;
+
+ checkTvmult<double>(A, V1, false);
+ checkTvmult<double>(A, V2, true);
}
DEAL::Tvmult
-DEAL::1.0000 10.0000 12.0000
+DEAL::1.0000 34.0000 39.0000
DEAL::Tvmult
-DEAL::1.0000 12.0000 15.0000
+DEAL::1.0000 36.0000 42.0000
deallog.threshold_double(1.e-10);
const double Adata[] =
- { 1, 2, 3, 4, 5, 6 };
+ { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
- FullMatrix<double> A(2, 3);
+ FullMatrix<double> A(3, 3);
A.fill(Adata);
- Vector<double> V(3);
- V(0) = 1;
- V(1) = 2;
- V(2) = 3;
+ Vector<double> V1(3);
+ Vector<double> V2(3);
- checkVmult_Add<double>(A, V, false);
- checkVmult_Add<double>(A, V, true);
+ V1(0) = V2(0) = 1;
+ V1(1) = V2(1) = 2;
+ V1(2) = V2(2) = 3;
+
+ checkVmult_Add<double>(A, V1, false);
+ checkVmult_Add<double>(A, V2, true);
}
DEAL::vmult_add
-DEAL::1.0000 29.0000
+DEAL::1.0000 29.0000 44.0000
DEAL::vmult_add
-DEAL::1.0000 33.0000
+DEAL::1.0000 33.0000 51.0000
deallog.threshold_double(1.e-10);
const double Adata[] =
- { 1, 2, 3, 4, 5, 6 };
+ { 1, 2, 3, 4, 5, 6, 7, 8, 9 };
- FullMatrix<double> A(2, 3);
+ FullMatrix<double> A(3, 3);
A.fill(Adata);
- Vector<double> V(2);
- V(0) = 1;
- V(1) = 2;
+ Vector<double> V1(3);
+ Vector<double> V2(3);
- checkTvmult_Add<double>(A, V, false);
- checkTvmult_Add<double>(A, V, true);
+ V1(0) = V2(0) = 1;
+ V1(1) = V2(1) = 2;
+ V1(2) = V2(2) = 3;
+
+ checkTvmult_Add<double>(A, V1, false);
+ checkTvmult_Add<double>(A, V2, true);
}
DEAL::Tvmult_add
-DEAL::1.0000 11.0000 13.0000
+DEAL::1.0000 35.0000 40.0000
DEAL::Tvmult_add
-DEAL::1.0000 13.0000 16.0000
+DEAL::1.0000 37.0000 43.0000