Gabriel, But I think I did a mistake somewhere in the function. Using the definiion > of the function from your pull request > I call > > *DoFTools::make_sparsity_pattern > (dof_handler_velocity,dof_handler_pressure,dsp, > constraints_vel,false);* > > and afterwards I call > > *SparsityTools::distribute_sparsity_pattern (dsp, > dof_handler_velocity.n_locally_owned_dofs_per_processor(),mpi_communicator,locally_relevant_dofs_vel);* > > When I run the program on more than 2 processors I get the following error > during the *distribute_sparsity... *function: > * The violated condition was: local_lines.is_element(line_n)* > > * Additional information: The index set given to this constraints > object indicates constraints for degree of freedom 16 should not be stored > by this object, but a constraint is being added.* > > Do you have an idea how I can fix this? >
You should first try create a unit test for this function, i.e. a small program that tests the function. Just create a small (non-trivial) mesh where you can compute the sparsity pattern by hand and compare with the output of the fiunction. I suspect that you need a second AffineConstraints object, such that you have one for the rows and one for the columns. 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 [email protected]. For more options, visit https://groups.google.com/d/optout.
