Rajat, It seems that I confused lines and did not express in the code what I wrote above: You have to add that additional index also when initializing the sparsity pattern, i.e. std::cout << " constraints printed " <<std::endl; DynamicSparsityPattern dsp (locally_relevant_dofs); should become std::cout << " constraints printed " <<std::endl; locally_relevant_dofs.add_index(the_dof); DynamicSparsityPattern dsp (locally_relevant_dofs);
Best, Daniel -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.