Skip to content
Snippets Groups Projects
Commit 15f4dfc1 authored by Jorrit Fahlke's avatar Jorrit Fahlke
Browse files

Dont use a matrix allocator to allocate vectors

[[Imported from SVN: r1192]]
parent 27e0d930
No related branches found
No related tags found
No related merge requests found
......@@ -55,7 +55,11 @@ namespace Dune
*/
template<typename T, typename A, int n, int m>
class SuperLU<BCRSMatrix<FieldMatrix<T,n,m>,A > >
: public InverseOperator<BlockVector<FieldVector<T,m>,A>,BlockVector<FieldVector<T,n>,A> >
: public InverseOperator<
BlockVector<FieldVector<T,m>,
typename A::template rebind<FieldVector<T,m> >::other>,
BlockVector<FieldVector<T,n>,
typename A::template rebind<FieldVector<T,n> >::other> >
{
public:
/* @brief The matrix type. */
......@@ -63,9 +67,13 @@ namespace Dune
/* @brief The corresponding SuperLU Matrix type.*/
typedef Dune::SuperLUMatrix<Matrix> SuperLUMatrix;
/** @brief The type of the domain of the solver. */
typedef Dune::BlockVector<FieldVector<T,m>,A> domain_type;
typedef Dune::BlockVector<
FieldVector<T,m>,
typename A::template rebind<FieldVector<T,m> >::other> domain_type;
/** @brief The type of the range of the solver. */
typedef Dune::BlockVector<FieldVector<T,n>,A> range_type;
typedef Dune::BlockVector<
FieldVector<T,n>,
typename A::template rebind<FieldVector<T,n> >::other> range_type;
/**
* @brief Constructs the SuperLU solver.
*
......
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