]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
adjust to new structure
authorkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 11 Oct 2007 16:24:08 +0000 (16:24 +0000)
committerkanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Thu, 11 Oct 2007 16:24:08 +0000 (16:24 +0000)
git-svn-id: https://svn.dealii.org/trunk@15290 0785d39b-7218-0410-832d-ea1e28bc413d

tests/deal.II/filtered_matrix.cc

index 68f2c57331b9c6c2b5a1526fdc335172078509b4..19e245e22d307a4845ec22da9e5a903f141f2012 100644 (file)
@@ -50,20 +50,22 @@ solve_filtered (std::map<unsigned int,double> &bv,
                Vector<double>                &u,
                Vector<double>                &f)
 {
-  FilteredMatrix<SparseMatrix<double> > A1 (A);
+  FilteredMatrix<Vector<double> > A1 (A);
   A1.add_constraints (bv);
   
   SolverControl control (1000, 1.e-10, false, false);
   PrimitiveVectorMemory<Vector<double> > mem;
   SolverCG<Vector<double> > solver (control, mem);
-  PreconditionJacobi<FilteredMatrix<SparseMatrix<double> > > prec;
-  prec.initialize (A1, 1.2);
+  PreconditionJacobi<SparseMatrix<double> > prec;
+  FilteredMatrix<Vector<double> > fprec;
+  prec.initialize (A, 1.2);
+  fprec.initialize(prec);
   
   Vector<double> f1 (f.size());
   f1 = f;
   A1.apply_constraints (f1, true);
   
-  solver.solve (A1, u, f1, prec);
+  solver.solve (A1, u, f1, fprec);
 
   for (std::map<unsigned int,double>::const_iterator i=bv.begin();
        i!=bv.end(); ++i)

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


Typeset in Trocchi and Trocchi Bold Sans Serif.