> i want to go for the matrix-free approach that Timo suggested, but how > about preconditioning? can one build a preconditioner using one of the > standard deal.ii classes and the new Matrix wrapper that Timo gave below?
You can construct a preconditioner in exactly the same way as for the matrix where you implement a vmult matrix that simply calls a function that does the transpose preconditioning. E.g. it could call SparseMatrix::TSOR; of course it could also call the SSOR function which might even work for a transpose matrix. W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
