Skip to content
Snippets Groups Projects
Commit d628dfb1 authored by Christian Engwer's avatar Christian Engwer
Browse files

[amg,bugfix] fix template arguments for MatrixHierarchy::coarsenVector

if the AMG had a non-standard allocator, the template parameters didn't
match, as the hierarchy was assumend to use the default allocator.
parent 1e30e207
No related branches found
No related tags found
1 merge request!87[multirhstest] Check with AlignedNumber.
......@@ -378,8 +378,8 @@ namespace Dune
* @brief Coarsen the vector hierarchy according to the matrix hierarchy.
* @param hierarchy The vector hierarchy to coarsen.
*/
template<class V, class TA>
void coarsenVector(Hierarchy<BlockVector<V,TA> >& hierarchy) const;
template<class V, class BA, class TA>
void coarsenVector(Hierarchy<BlockVector<V,BA>, TA>& hierarchy) const;
/**
* @brief Coarsen the smoother hierarchy according to the matrix hierarchy.
......@@ -1106,8 +1106,8 @@ namespace Dune
}
template<class M, class IS, class A>
template<class V, class TA>
void MatrixHierarchy<M,IS,A>::coarsenVector(Hierarchy<BlockVector<V,TA> >& hierarchy) const
template<class V, class BA, class TA>
void MatrixHierarchy<M,IS,A>::coarsenVector(Hierarchy<BlockVector<V,BA>, TA>& hierarchy) const
{
assert(hierarchy.levels()==1);
typedef typename ParallelMatrixHierarchy::ConstIterator Iterator;
......
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