]> https://gitweb.dealii.org/ - dealii.git/commit
Fix bug in distribute_local_to_global 7726/head
authorSebastian Stark <stark.sebastian@gmx.de>
Wed, 13 Feb 2019 11:33:09 +0000 (13:33 +0200)
committerSebastian Stark <stark.sebastian@gmx.de>
Thu, 14 Feb 2019 13:09:13 +0000 (15:09 +0200)
commit5cca395fd110d6d30112aca0f157133c92cc9dfd
tree90a1e50b529a5a39f7f36f804d66df7310844eab
parent4f7b5902333cbcfdb031d33406268720c79bee47
Fix bug in distribute_local_to_global

If a local matrix with all diagonal elements equal to zero is distributed to a global matrix, the l1 norm of the local matrix divided by the size of the local matrix is now added to those diagonal elements of the global matrix which correspond to a constrained dof. In case the entire local matrix is zero, 1 is added. Previously zero was added for both cases, possibly resulting in singular global matrices. Additionally, a test for the patched version of distribute_local_to_global has been added.
Fixes #7658
doc/news/changes/minor/20190213SebastianStark [new file with mode: 0644]
include/deal.II/lac/affine_constraints.templates.h
tests/lac/constraint_matrix_distribute_local_global.cc [new file with mode: 0644]
tests/lac/constraint_matrix_distribute_local_global.output [new file with mode: 0644]

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.