// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii8SolverCGINS_16TrilinosWrappers6VectorEEE
DEAL:cg::Starting value 31.0000
-DEAL:cg::Convergence step 128 value 0.000807536
-DEAL::Solver stopped after 128 iterations
+DEAL:cg::Convergence step 43 value 0.000818932
+DEAL::Solver stopped after 43 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii14SolverBicgstabINS_16TrilinosWrappers6VectorEEE
DEAL:Bicgstab::Starting value 31.0000
-DEAL:Bicgstab::Convergence step 157 value 0.000791865
-DEAL::Solver stopped after 157 iterations
+DEAL:Bicgstab::Convergence step 29 value 0.000873504
+DEAL::Solver stopped after 29 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii11SolverGMRESINS_16TrilinosWrappers6VectorEEE
DEAL:GMRES::Starting value 31.0000
-DEAL:GMRES::Convergence step 1216 value 0.000999125
-DEAL::Solver stopped after 1216 iterations
+DEAL:GMRES::Convergence step 75 value 0.000809583
+DEAL::Solver stopped after 75 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii12SolverMinResINS_16TrilinosWrappers6VectorEEE
DEAL:minres::Starting value 31.0000
-DEAL:minres::Convergence step 33 value 0.000792239
-DEAL::Solver stopped after 33 iterations
+DEAL:minres::Convergence step 5 value 0.000393678
+DEAL::Solver stopped after 5 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii10SolverQMRSINS_16TrilinosWrappers6VectorEEE
DEAL:QMRS::Starting value 31.0000
-DEAL:QMRS::Convergence step 135 value 0.000563443
-DEAL::Solver stopped after 135 iterations
+DEAL:QMRS::Convergence step 46 value 0.000906349
+DEAL::Solver stopped after 46 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
DEAL::Solver type: N6dealii16SolverRichardsonINS_16TrilinosWrappers6VectorEEE
DEAL:Richardson::Starting value 31.0000
-DEAL:Richardson::Failure step 304 value nan
-DEAL::Iterative method reported convergence failure in step 304 with residual nan
-DEAL::Solver stopped after 304 iterations
+DEAL:Richardson::Failure step 374 value nan
+DEAL::Iterative method reported convergence failure in step 374 with residual nan
+DEAL::Solver stopped after 374 iterations
{
- SolverControl control(200, 1.e-3);
+ SolverControl control(2000, 1.e-3);
const unsigned int size = 32;
unsigned int dim = (size-1)*(size-1);
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
-DEAL::Solver type: N6dealii13TrilinosWrappers8SolverCGE
-DEAL::Starting value 7.750
-DEAL::Convergence step 41 value 0.0006730
-DEAL::Solver stopped after 41 iterations
+DEAL::Solver type: N6dealii16TrilinosWrappers8SolverCGE
+DEAL::Convergence step 43 value 0.000818932
+DEAL::Solver stopped after 43 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
-DEAL::Solver type: N6dealii13TrilinosWrappers11SolverGMRESE
-DEAL::Starting value 7.750
-DEAL::Convergence step 48 value 0.0009841
-DEAL::Solver stopped after 48 iterations
+DEAL::Solver type: N6dealii16TrilinosWrappers11SolverGMRESE
+DEAL::Convergence step 71 value 0.000916839
+DEAL::Solver stopped after 71 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);
DEAL::Size 32 Unknowns 961
-DEAL::Solver type: N6dealii13TrilinosWrappers9SolverCGSE
-DEAL::Starting value 7.750
-DEAL::Convergence step 40 value 0.0007519
-DEAL::Solver stopped after 40 iterations
+DEAL::Solver type: N6dealii16TrilinosWrappers9SolverCGSE
+DEAL::Convergence step 41 value 0.000549329
+DEAL::Solver stopped after 41 iterations
// Make matrix
FDMatrix testproblem(size, size);
- TrilinosWrappers::SparseMatrix A(dim, dim, 5U);
+ CompressedSimpleSparsityPattern csp (dim, dim);
+ testproblem.five_point_structure(csp);
+ TrilinosWrappers::SparseMatrix A;
+ A.reinit(csp);
testproblem.five_point(A);
TrilinosWrappers::Vector f(dim);