Dear Wolfgang, 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.
Thanks for your help. Martin On Thu, Jun 3, 2010 at 2:25 PM, Wolfgang Bangerth <[email protected]> wrote: > >> I have a FESystem (Q_y,Q_u) and I am assembling a matrix [K_y >> N_yu;N_uy M_u] but I also need a mass matrix M_y. >> For this I thought that I could simply call something like >> MatrixCreator::create_mass_matrix(dof_handler,...,mass_matrix) >> by specifying the dof_handler only for the y-components. Is there >> anyway that this is possible by only accessing the y_components of the >> dof_handler? > > You can build the matrix > [M_y 0; 0 0] > by calling MatrixCreator::create_mass_matrix on the joint DoFHandler object > and passing a coefficient function object as last argument to > create_mass_matrix. This coefficient needs to be vector valued with two > components > [1 ; 0] > > Does this help? > > W. > > > ------------------------------------------------------------------------- > Wolfgang Bangerth email: [email protected] > www: http://www.math.tamu.edu/~bangerth/ > > _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
