data.degree = 3;
prec.initialize(m, data);
+ deallog << "Exact inverse: ";
+ for (unsigned int i=0; i<size; ++i)
+ deallog << in(i)/m(i,i) << " ";
+ deallog << std::endl;
+
deallog << "Check vmult orig: ";
prec.vmult(out, in);
for (unsigned int i=0; i<size; ++i)
{
std::ofstream logfile("output");
deallog << std::fixed;
- deallog << std::setprecision(0);
+ deallog << std::setprecision(2);
deallog.attach(logfile);
deallog.depth_console(0);
deallog.threshold_double(1.e-10);
-DEAL::Check vmult orig: 1. 0. 0. 0. 0. 0. 0. 0. 0. 0.
-DEAL::Check Tvmult orig: 1. 0. 0. 0. 0. 0. 0. 0. 0. 0.
-DEAL::Check vmult diag: 1. 0. 0. 0. 0. 0. 0. 0. 0. 0.
-DEAL::Check Tvmult diag: 1. 0. 0. 0. 0. 0. 0. 0. 0. 0.
+DEAL::Exact inverse: 0.84 0.20 0.26 0.20 0.18 0.03 0.05 0.10 0.03 0.06
+DEAL::Check vmult orig: 1.09 0.26 0.34 0.26 0.24 0.04 0.06 0.12 0.04 0.07
+DEAL::Check Tvmult orig: 1.09 0.26 0.34 0.26 0.24 0.04 0.06 0.12 0.04 0.07
+DEAL::Check vmult diag: 0.82 0.26 0.33 0.21 0.15 0.02 0.03 0.09 0.03 0.07
+DEAL::Check Tvmult diag: 0.82 0.26 0.33 0.21 0.15 0.02 0.03 0.09 0.03 0.07
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