From: guido Date: Fri, 30 May 2003 19:08:45 +0000 (+0000) Subject: additional tests disabled X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7f61c7221d75ec92dd49f42586d32f90c5a36384;p=dealii-svn.git additional tests disabled git-svn-id: https://svn.dealii.org/trunk@7690 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/lac/solver.cc b/tests/lac/solver.cc index 6e83139730..0dd3b17fab 100644 --- a/tests/lac/solver.cc +++ b/tests/lac/solver.cc @@ -19,6 +19,7 @@ #include "testmatrix.h" #include #include +#include #include #include #include @@ -188,13 +189,17 @@ int main() deallog.pop(); deallog.push("psor"); - + + if(false) + { + check_Tsolve(rich,A,u,f,prec_psor); check_solve(rich,A,u,f,prec_psor); check_solve(cg,A,u,f,prec_psor); check_solve(bicgstab,A,u,f,prec_psor); check_solve(gmres,A,u,f,prec_psor); check_solve(gmresright,A,u,f,prec_psor); + } deallog.pop(); } @@ -203,5 +208,63 @@ int main() std::cerr << e.what() << std::endl; } }; + + // Solve advection problem + + // Disabled for the time being + for (unsigned int size=4; size <= 3; size *= 3) + { + unsigned int dim = (size-1)*(size-1); + + deallog << "Size " << size << " Unknowns " << dim << std::endl; + + // Make matrix + FDMatrix testproblem(size, size); + SparsityPattern structure(dim, dim, 5); + testproblem.five_point_structure(structure); + structure.compress(); + SparseMatrix A(structure); + testproblem.upwind(A, true); + + PreconditionSOR<> prec_sor; + prec_sor.initialize(A, 1.); + + std::vector permutation(dim); + std::vector inverse_permutation(dim); + + // Create a permutation: Blocks + // backwards and every second + // block backwards + unsigned int k = 0; + for (unsigned int i=0;i prec_psor; + prec_psor.initialize(A, permutation, inverse_permutation, 1.); + + Vector f(dim); + Vector u(dim); + f = 1.; + u = 1.; + + std::cerr << "******************************" << std::endl; + + check_solve(rich,A,u,f,prec_sor); + check_solve(rich,A,u,f,prec_psor); + } }