Skip to content
Snippets Groups Projects
Commit 28948f0c authored by Uli Sack's avatar Uli Sack Committed by usack
Browse files

mirror changes in dune/fufem/staticmatrixtools.hh

[[Imported from SVN: r5757]]
parent 0a8d38cf
No related branches found
No related tags found
No related merge requests found
...@@ -95,13 +95,13 @@ class StaticMatrix ...@@ -95,13 +95,13 @@ class StaticMatrix
// } // }
template <typename FieldType, int n, class MatrixB> template <typename FieldType, int n, class MatrixB>
static void axpy(Dune::ScaledIdentityMatrix<FieldType,n>& x, const Dune::ScaledIdentityMatrix<FieldType,n>::field_type a, const MatrixB& y) static void axpy(Dune::ScaledIdentityMatrix<FieldType,n>& x, const typename Dune::ScaledIdentityMatrix<FieldType,n>::field_type a, const MatrixB& y)
{ {
x.scalar() += a * y[0][0]; x.scalar() += a * y[0][0];
} }
template <typename FieldType, int n, class MatrixB> template <typename FieldType, int n, class MatrixB>
static void axpy(Dune::DiagonalMatrix<FieldType,n>& x, const Dune::DiagonalMatrix<FieldType,n>::field_type a, const MatrixB& y) static void axpy(Dune::DiagonalMatrix<FieldType,n>& x, const typename Dune::DiagonalMatrix<FieldType,n>::field_type a, const MatrixB& y)
{ {
for(int i=0; i<n; ++i) for(int i=0; i<n; ++i)
x.diagonal()[i] += a * y[i][i]; x.diagonal()[i] += a * y[i][i];
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment