Skip to content
Snippets Groups Projects
Commit da8bed73 authored by Carsten Gräser's avatar Carsten Gräser
Browse files

[amg][bugfix] Fix dirichlet processing in FastAMG

The comment and surroundig code indicate that FastAMG should
solve during pre() for dirichlet entries. This is determined
by checking if a row does contain a diagonal entry while
all other entries are zero.

The check for this was buggy because hasDiagonal was always false.
Now it's set to true if a diagonal entry is found and nonzero.
parent c946134d
No related branches found
No related tags found
No related merge requests found
......@@ -471,7 +471,7 @@ namespace Dune
for(ColIter col=row->begin(); col!=row->end(); ++col) {
if(row.index()==col.index()) {
diag = col;
hasDiagonal = false;
hasDiagonal = (*col != zero);
}else{
if(*col!=zero)
isDirichlet = false;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment