> thanks. I thought about this as a kind of fall back option but was
> hoping not to have M_y as a part of
> a BlockSparse but to be honest I am sure it will work perfectly.

Sure, you could also just call the function with no coefficient and pass it a 
block sparse matrix. The block you want is then
  mass_matrix.block(0,0)

W.

-------------------------------------------------------------------------
Wolfgang Bangerth                email:            [email protected]
                                 www: http://www.math.tamu.edu/~bangerth/

_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii

Reply via email to