Skip to content
Snippets Groups Projects
Commit 4fd6b18a authored by Robert Klöfkorn's avatar Robert Klöfkorn
Browse files

Jesus, check that the code compiles before commiting.

[[Imported from SVN: r6206]]
parent 92ea0264
Branches
Tags
No related merge requests found
......@@ -1079,16 +1079,16 @@ namespace Dune
//! calculates ret = matrix * x
template <typename MAT, typename V1, typename V2>
static inline void multAssign(const DenseMatrix<MAT> &matrix, const DenseVector<V2> & x, DenseVector<V2> & ret)
static inline void multAssign(const DenseMatrix<MAT> &matrix, const DenseVector<V1> & x, DenseVector<V2> & ret)
{
assert(x.size() == matrix.rows());
assert(ret.size() == matrix.cols());
typedef typename DenseMatrix<MAT>::size_type size_type;
for(size_type i=0; i<matrix.size(); ++i)
for(size_type i=0; i<matrix.rows(); ++i)
{
ret[i] = 0.0;
for(size_type j=0; j<matrix[i].size(); ++j)
for(size_type j=0; j<matrix.cols(); ++j)
{
ret[i] += matrix[i][j]*x[j];
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment