On 9/23/2010 7:16 PM, fname lname wrote:
Thanks Timo,I've tried that already, it gets me the same matrix as if I do: sparsity_pattern.block(i,j).reinit(dof_handler.n_dofs(),dof_handler.n_dofs(),0);
Can't you call sparsity_pattern.block(i,j).compress() right there? G _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
