Greeting Timo,
> I don't understand what you mean. Why can you not use it? I am using
> the deal solvers with PETSc block matrices for sure. There _is_ a
> PETScWrappers::MPI::BlockSparseMatrix::vmult (BlockVector &dst, const
> BlockVector &src) function. Or is this a misunderstanding from my
> side?
You understood me correctly. Maybe it is my err. I remebe that the
SolverGMRES asked for a
PETScWrappers::PreconditionerBase::vmult (BlockVector &dst, const
BlockVector &src). I'll go back to it soon, I can see how he should work
for sure; too many days in front of the green screen methinks.
> It was an eye-opener for me when I found that out. In those moments I
> am impressed by the power of abstraction. Somebody wrote a simple,
> templated GMRES algorithm without thinking about parallelization and
> then you can use it to solve distributed block systems with thousands
> of CPUs without a single change. :-)
Similar thoughts went through my head; moments before my eyes popped out
and I passed out from the cold shock of it. That is indeed a cunning
little gem hidden away there in deal.II; and it leads to more....
Cheers,
Toby
-----
Toby D. Young
Assistant Professor
Institute of Fundamental Technological Problems
Polish Academy of Sciences
ul Adolfa Pawinskiego 5b
02-106 Warsaw
Poland
www: http://www.ippt.gov.pl/~tyoung
skype: stenografia
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
