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

Fixes overwritten right hand side for FastAMG with one level.

Unfortunately, SuperLU always overwrites the right hand side with
the solution. Therefore we now copy it if there is only one level
in the matrix hierarchy of AMG.
parent 0f850d8b
No related branches found
No related tags found
No related merge requests found
......@@ -415,7 +415,12 @@ namespace Dune
pinfo->copyOwnerToAll(*lhs, *lhs);
v=*lhs;
}else{*/
mgc(v, d);
if(matrices_->maxlevels()==1){
// The coarse solver might modify the d!
Range b(d);
mgc(v, b);
}else
mgc(v, d);
if(postSteps_==0||matrices_->maxlevels()==1)
pinfo->copyOwnerToAll(v, v);
//}
......
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