Hey Daniel,
thanks for the pull, with this instantiation everything compiles fine and I can use it.
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);
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.
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?
(When I run the code on 1 or 2 processors, this error doesnt come up.
But I get divergence in the Gmres mehod,
which might imply that the boundaryconditions are not
applied in the correct way yet, but I first want to focus on the upper error.)
Thanks and best regards
Gabriel
Gesendet: Mittwoch, 03. April 2019 um 03:03 Uhr
Von: "Daniel Arndt" <daniel.ar...@iwr.uni-heidelberg.de>
An: "deal.II User Group" <dealii@googlegroups.com>
Betreff: Re: [deal.II] Re: Applying boundary values in parll::distr:triang setting for two dof_handler Sparsematrux
Von: "Daniel Arndt" <daniel.ar...@iwr.uni-heidelberg.de>
An: "deal.II User Group" <dealii@googlegroups.com>
Betreff: Re: [deal.II] Re: Applying boundary values in parll::distr:triang setting for two dof_handler Sparsematrux
Gabriel,
Now I cloned a more recent version, and added the Constraintsmatrix and changed the instantiations.
That looks much better!
But somehow it still doesn't work. I still get the error "undefined reference..."
Admittedly, getting the instantiations right in this case might not be straightforward.
I created a pull request to your branch at https://github.com/gabpeters/dealii/pull/1 to fix some issues such that the branch at least compiles for me.
If you are satisfied with everything, please create a pull request to the master branch.
I hope I dont bother you too much with all these questions.
Not at all! We are happily helping people to contribute.
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.
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.
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.
