From 69801d25c6958e047b614a53da062dd5e8173502 Mon Sep 17 00:00:00 2001
From: Andreas Dedner <dedner@dune-project.org>
Date: Sat, 13 Mar 2010 12:37:14 +0000
Subject: [PATCH] fix FS745
[[Imported from SVN: r5928]]
---
dune/common/fmatrix.hh | 22 +++++++++++-----------
1 file changed, 11 insertions(+), 11 deletions(-)
diff --git a/dune/common/fmatrix.hh b/dune/common/fmatrix.hh
index 14f8606b2..2498a15b6 100644
--- a/dune/common/fmatrix.hh
+++ b/dune/common/fmatrix.hh
@@ -1096,77 +1096,77 @@ namespace Dune
template<class X, class Y>
void mv (const X& x, Y& y) const
{
- y.p = a[0] * x.p;
+ y[0] = a[0] * x[0];
}
//! y = A^T x
template<class X, class Y>
void mtv ( const X &x, Y &y ) const
{
- y.p = a[ 0 ] * x.p;
+ y[0] = a[ 0 ] * x[0];
}
//! y += A x
template<class X, class Y>
void umv (const X& x, Y& y) const
{
- y.p += a[0] * x.p;
+ y[0] += a[0] * x[0];
}
//! y += A^T x
template<class X, class Y>
void umtv (const X& x, Y& y) const
{
- y.p += a[0] * x.p;
+ y[0] += a[0] * x[0];
}
//! y += A^H x
template<class X, class Y>
void umhv (const X& x, Y& y) const
{
- y.p += conjugateComplex(a[0]) * x.p;
+ y[0] += conjugateComplex(a[0]) * x[0];
}
//! y -= A x
template<class X, class Y>
void mmv (const X& x, Y& y) const
{
- y.p -= a[0] * x.p;
+ y[0] -= a[0] * x[0];
}
//! y -= A^T x
template<class X, class Y>
void mmtv (const X& x, Y& y) const
{
- y.p -= a[0] * x.p;
+ y[0] -= a[0] * x[0];
}
//! y -= A^H x
template<class X, class Y>
void mmhv (const X& x, Y& y) const
{
- y.p -= conjugateComplex(a[0]) * x.p;
+ y[0] -= conjugateComplex(a[0]) * x[0];
}
//! y += alpha A x
template<class X, class Y>
void usmv (const K& alpha, const X& x, Y& y) const
{
- y.p += alpha * a[0] * x.p;
+ y[0] += alpha * a[0] * x[0];
}
//! y += alpha A^T x
template<class X, class Y>
void usmtv (const K& alpha, const X& x, Y& y) const
{
- y.p += alpha * a[0] * x.p;
+ y[0] += alpha * a[0] * x[0];
}
//! y += alpha A^H x
template<class X, class Y>
void usmhv (const K& alpha, const X& x, Y& y) const
{
- y.p += alpha * conjugateComplex(a[0]) * x.p;
+ y[0] += alpha * conjugateComplex(a[0]) * x[0];
}
//===== norms
--
GitLab