> If I do a component renumbering of the variable without any cuthill it > seems that my K is not of the form > K=blkdiag(K_1,K_2)
It's hard to tell without seeing your code why it isn't of that form. How did you check? And do you renumber your DoFs by component (y_1, y_2, p) or by block (y, p)? W. -- ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
