Skip to content
Snippets Groups Projects
Commit a94841ea authored by Markus Blatt's avatar Markus Blatt Committed by Carsten Gräser
Browse files

[release] Do not assume a symmetric sparsity pattern in SymmetricDependency.


Previously, we assumed that if a_{ij} is stored in the sparse matrix
then a_{ji} must be stored also and used mat[i][j] to access. If the entry
was not stored then an exception somewhere in basearray.hh was thrown.
Unfortunately knowing the cause for this exception is insider knowledge,
seldomly leaked ;). In addition there is now reason to throw an expection here.

With this commit we the find method on the row to search for the entry a_{ji}.
If it is not present we simply treat it as being zero.

(cherry picked from commit 05c0a63a)
Signed-off-by: default avatarCarsten Gräser <graeser@dune-project.org>
parent 6e80cede
No related branches found
No related tags found
1 merge request!112Add data() accessors to BlockVector
Loading
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment