reference_rhs);
deallog << " Reference matrix nonzeros: " << reference_matrix.n_nonzero_elements()
- << ", actually: " << reference_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << reference_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
}
test_all_constraints.condense (test_matrix, test_rhs);
deallog << " Test matrix 1 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
test_equality();
test_rhs);
}
deallog << " Test matrix 2 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
test_equality();
}
logfile << std::setprecision (2);
deallog.attach(logfile);
deallog.depth_console(0);
- deallog.threshold_double(1.e-12);
+ deallog.threshold_double(1.e-10);
LaplaceProblem<2> laplace_problem;
laplace_problem.run ();
DEAL:: Number of active cells: 192
DEAL:: Number of degrees of freedom: 864
DEAL:: Number of constraints : 0
-DEAL:: Reference matrix nonzeros: 12672, actually: 8974
-DEAL:: Test matrix 1 nonzeros: 12672, actually: 8974
+DEAL:: Reference matrix nonzeros: 12672, actually: 8480
+DEAL:: Test matrix 1 nonzeros: 12672, actually: 8480
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 12672, actually: 8974
+DEAL:: Test matrix 2 nonzeros: 12672, actually: 8480
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
DEAL:cg::Starting value 20.
DEAL:: Number of active cells: 309
DEAL:: Number of degrees of freedom: 1575
DEAL:: Number of constraints : 187
-DEAL:: Reference matrix nonzeros: 25901, actually: 17333
-DEAL:: Test matrix 1 nonzeros: 25901, actually: 17333
+DEAL:: Reference matrix nonzeros: 25901, actually: 16531
+DEAL:: Test matrix 1 nonzeros: 25901, actually: 16531
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 25901, actually: 17333
+DEAL:: Test matrix 2 nonzeros: 25901, actually: 16531
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
DEAL:cg::Starting value 22.
DEAL:: Number of active cells: 450
DEAL:: Number of degrees of freedom: 2692
DEAL:: Number of constraints : 475
-DEAL:: Reference matrix nonzeros: 51910, actually: 33908
-DEAL:: Test matrix 1 nonzeros: 51910, actually: 33908
+DEAL:: Reference matrix nonzeros: 51910, actually: 32730
+DEAL:: Test matrix 1 nonzeros: 51910, actually: 32730
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 51910, actually: 33908
+DEAL:: Test matrix 2 nonzeros: 51910, actually: 32730
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
DEAL:cg::Starting value 25.
+++ /dev/null
-
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 192
-DEAL:: Number of degrees of freedom: 864
-DEAL:: Number of constraints : 0
-DEAL:: Reference matrix nonzeros: 12672, actually: 8974
-DEAL:: Test matrix 1 nonzeros: 12672, actually: 8974
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 12672, actually: 8974
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 20.
-DEAL:cg::Convergence step 25 value 2.5e-07
-DEAL::Distribute error: 0
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 309
-DEAL:: Number of degrees of freedom: 1575
-DEAL:: Number of constraints : 187
-DEAL:: Reference matrix nonzeros: 25901, actually: 17333
-DEAL:: Test matrix 1 nonzeros: 25901, actually: 17333
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 25901, actually: 17333
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 22.
-DEAL:cg::Convergence step 42 value 3.2e-07
-DEAL::Distribute error: 0
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 450
-DEAL:: Number of degrees of freedom: 2692
-DEAL:: Number of constraints : 475
-DEAL:: Reference matrix nonzeros: 51910, actually: 33906
-DEAL:: Test matrix 1 nonzeros: 51910, actually: 33906
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 51910, actually: 33906
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 25.
-DEAL:cg::Convergence step 61 value 3.8e-07
-DEAL::Distribute error: 0
+++ /dev/null
-
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 192
-DEAL:: Number of degrees of freedom: 864
-DEAL:: Number of constraints : 0
-DEAL:: Reference matrix nonzeros: 12672, actually: 8974
-DEAL:: Test matrix 1 nonzeros: 12672, actually: 8974
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 12672, actually: 8974
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 20.
-DEAL:cg::Convergence step 25 value 2.5e-07
-DEAL::Distribute error: 0
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 309
-DEAL:: Number of degrees of freedom: 1575
-DEAL:: Number of constraints : 187
-DEAL:: Reference matrix nonzeros: 25901, actually: 17333
-DEAL:: Test matrix 1 nonzeros: 25901, actually: 17333
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 25901, actually: 17333
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 22.
-DEAL:cg::Convergence step 42 value 3.2e-07
-DEAL::Distribute error: 0
-DEAL::
-DEAL::
-DEAL:: Number of active cells: 450
-DEAL:: Number of degrees of freedom: 2692
-DEAL:: Number of constraints : 475
-DEAL:: Reference matrix nonzeros: 51910, actually: 33908
-DEAL:: Test matrix 1 nonzeros: 51910, actually: 33908
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 51910, actually: 33908
-DEAL:: Matrix difference norm: 0
-DEAL:: RHS difference norm: 0
-DEAL:cg::Starting value 25.
-DEAL:cg::Convergence step 61 value 3.8e-07
-DEAL::Distribute error: 0
reference_rhs);
deallog << " Reference matrix nonzeros: " << reference_matrix.n_nonzero_elements()
- << ", actually: " << reference_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << reference_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
}
test_all_constraints.condense (test_matrix, test_rhs);
deallog << " Test matrix 1 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
test_equality();
test_rhs);
}
deallog << " Test matrix 2 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-10)
<< std::endl;
test_equality();
}
logfile << std::setprecision (2);
deallog.attach(logfile);
deallog.depth_console(0);
- deallog.threshold_double(1.e-12);
+ deallog.threshold_double(1.e-10);
{
AdvectionProblem<2> advection_problem;
DEAL:: Number of active cells: 22
DEAL:: Number of degrees of freedom: 228
DEAL:: Number of constraints : 18
-DEAL:: Reference matrix nonzeros: 6700, actually: 1402
-DEAL:: Test matrix 1 nonzeros: 6700, actually: 1402
+DEAL:: Reference matrix nonzeros: 6700, actually: 1254
+DEAL:: Test matrix 1 nonzeros: 6700, actually: 1254
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 6700, actually: 1396
+DEAL:: Test matrix 2 nonzeros: 6700, actually: 1254
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
DEAL::
DEAL:: Number of active cells: 78
DEAL:: Number of degrees of freedom: 1856
DEAL:: Number of constraints : 174
-DEAL:: Reference matrix nonzeros: 189772, actually: 28840
-DEAL:: Test matrix 1 nonzeros: 189772, actually: 28840
+DEAL:: Reference matrix nonzeros: 189772, actually: 24098
+DEAL:: Test matrix 1 nonzeros: 189772, actually: 24194
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 189772, actually: 28832
+DEAL:: Test matrix 2 nonzeros: 189772, actually: 24108
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
reference_rhs);
deallog << " Reference matrix nonzeros: " << reference_matrix.n_nonzero_elements()
- << ", actually: " << reference_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << reference_matrix.n_actually_nonzero_elements (1e-12)
<< std::endl;
}
test_all_constraints.condense (test_matrix, test_rhs);
deallog << " Test matrix 1 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-12)
<< std::endl;
test_equality();
test_rhs);
}
deallog << " Test matrix 2 nonzeros: " << test_matrix.n_nonzero_elements()
- << ", actually: " << test_matrix.n_actually_nonzero_elements ()
+ << ", actually: " << test_matrix.n_actually_nonzero_elements (1e-12)
<< std::endl;
test_equality();
}
DEAL:: Number of active cells: 22
DEAL:: Number of degrees of freedom: 114
DEAL:: Number of constraints : 9
-DEAL:: Reference matrix nonzeros: 1618, actually: 701
-DEAL:: Test matrix 1 nonzeros: 1618, actually: 701
+DEAL:: Reference matrix nonzeros: 1618, actually: 627
+DEAL:: Test matrix 1 nonzeros: 1618, actually: 627
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 1618, actually: 698
+DEAL:: Test matrix 2 nonzeros: 1618, actually: 627
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
DEAL::
DEAL:: Number of active cells: 78
DEAL:: Number of degrees of freedom: 928
DEAL:: Number of constraints : 87
-DEAL:: Reference matrix nonzeros: 46980, actually: 14420
-DEAL:: Test matrix 1 nonzeros: 46980, actually: 14420
+DEAL:: Reference matrix nonzeros: 46980, actually: 12049
+DEAL:: Test matrix 1 nonzeros: 46980, actually: 12097
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0
-DEAL:: Test matrix 2 nonzeros: 46980, actually: 14416
+DEAL:: Test matrix 2 nonzeros: 46980, actually: 12054
DEAL:: Matrix difference norm: 0
DEAL:: RHS difference norm: 0