> 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
