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

Bugfix.

We need to set the overlap to zero on the coarse levels to avoid
degradation of convergence for coarse levels.

[[Imported from SVN: r1026]]
parent 48dd1d31
No related branches found
No related tags found
No related merge requests found
......@@ -136,6 +136,9 @@ namespace Dune
ParallelInformation<T>& comm)
{
Transfer<V,BlockVector<B>,SequentialInformation>::restrict (aggregates, coarse, fine, SequentialInformation());
// We need this here to avoid it in the smoothers on the coarse level.
// There (in the preconditioner d is const.
comm.project(coarse);
}
template<class V, class B, class T1, class T2>
......@@ -153,8 +156,9 @@ namespace Dune
OwnerOverlapCopyCommunication<T1,T2>& comm)
{
Transfer<V,BlockVector<B>,SequentialInformation>::restrict (aggregates, coarse, fine, SequentialInformation());
comm.copyOwnerToAll(coarse,coarse);
//comm.project(coarse);
// We need this here to avoid it in the smoothers on the coarse level.
// There (in the preconditioner d is const.
comm.project(coarse);
}
#endif
/** @} */
......
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