Christian,
> Following up, here is a test case for deal.II's BlockVector which fails.
Excellent. I could have sworn that Vector::operator= does not resize its left
hand side operand and that consequently BlockVector::operator= should not do
so either. Alas, my recollection is false: The former does, and so should the
latter. The fix for this is now r17380 in the subversion archive. Can you
double check whether that also fixes your SolutionTransfer problem?
> * Where is the definition of BlockVector's operator= ?
It's inline in lac/include/lac/block_vector.h.
> * If it exists, shouldn't it do something like collect_sizes() ?
It does now :-)
Thanks for pointing out where the problem was! Cheers
W.
-------------------------------------------------------------------------
Wolfgang Bangerth email: [EMAIL PROTECTED]
www: http://www.math.tamu.edu/~bangerth/
_______________________________________________
