As an alternative to Martin's suggestion, you could also just add all contributions, and then call constraints.condense (system_matrix, solution, system_rhs); at the very end, rather than eliminating constrained degrees of freedom already during assembly.
W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
