Skip to content
Snippets Groups Projects
Commit fd613e1a authored by Christoph Gersbacher's avatar Christoph Gersbacher
Browse files

Fix segmentation fault in matrix interface check.

[[Imported from SVN: r7104]]
parent 157c684f
No related branches found
No related tags found
No related merge requests found
......@@ -164,7 +164,7 @@ namespace CheckMatrixInterface
static void checkSolve ( const Matrix &matrix )
{
typename Traits::domain_type x;
typename Traits::domain_type x = Traits::domain( matrix );
const typename Traits::range_type b = Traits::range( matrix );
matrix.solve( x, b );
}
......@@ -267,16 +267,16 @@ namespace CheckMatrixInterface
typename Traits::value_type alpha( 1 );
matrix.mv( domain, range );
matrix.mtv( domain, range );
matrix.mtv( range, domain );
matrix.umv( domain, range );
matrix.umtv( domain, range );
matrix.umhv( domain, range );
matrix.umtv( range, domain );
matrix.umhv( range, domain );
matrix.mmv( domain, range );
matrix.mmtv( domain, range );
matrix.mmhv( domain, range );
matrix.mmtv( range, domain );
matrix.mmhv( range, domain );
matrix.usmv( alpha, domain, range );
matrix.usmtv( alpha, domain, range );
matrix.usmhv( alpha, domain, range );
matrix.usmtv( alpha, range, domain );
matrix.usmhv( alpha, range, domain );
}
// check iterator methods
......
......@@ -572,7 +572,7 @@ void test_interface()
typedef CheckMatrixInterface::UseFieldVector< K, rows, cols > Traits;
typedef Dune::FieldMatrix< K, rows, cols > FMatrix;
FMatrix m;
FMatrix m( 1 );
checkMatrixInterface< FMatrix >( m );
checkMatrixInterface< FMatrix, Traits >( m );
}
......
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