Skip to content
Snippets Groups Projects
Commit 969fa8a4 authored by Markus Blatt's avatar Markus Blatt
Browse files

Merge branch 'feature/fix-fastamg-dirichlet-processing' into 'master'

[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.
This merge commits sets it to true if a diagonal entry is found and nonzero.

See merge request !32
parents ca9342c0 da8bed73
No related branches found
No related tags found
1 merge request!32[amg][bugfix] Fix dirichlet processing in FastAMG
......@@ -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