From b2d2554a5d188b8762c1904fccc44ee21f70cd66 Mon Sep 17 00:00:00 2001
From: Christian Engwer <christi@dune-project.org>
Date: Fri, 14 Dec 2007 00:16:05 +0000
Subject: [PATCH] added missing method mv
[[Imported from SVN: r5029]]
---
common/fmatrix.hh | 19 +++++++++++++++++++
1 file changed, 19 insertions(+)
diff --git a/common/fmatrix.hh b/common/fmatrix.hh
index 23f27f1f2..302310a75 100644
--- a/common/fmatrix.hh
+++ b/common/fmatrix.hh
@@ -235,6 +235,19 @@ namespace Dune {
//===== linear maps
+ //! y = A x
+ template<class X, class Y>
+ void mv (const X& x, Y& y) const
+ {
+#ifdef DUNE_FMatrix_WITH_CHECKING
+ if (x.N()!=M()) DUNE_THROW(FMatrixError,"index out of range");
+ if (y.N()!=N()) DUNE_THROW(FMatrixError,"index out of range");
+#endif
+ for (size_type i=0; i<n; i++)
+ for (size_type j=0; j<m; j++)
+ y[i] = (*this)[i][j] * x[j];
+ }
+
//! y += A x
template<class X, class Y>
void umv (const X& x, Y& y) const
@@ -955,6 +968,12 @@ namespace Dune {
//===== linear maps
+ //! y = A x
+ void mv (const FieldVector<K,1>& x, FieldVector<K,1>& y) const
+ {
+ y.p = a[0] * x.p;
+ }
+
//! y += A x
void umv (const FieldVector<K,1>& x, FieldVector<K,1>& y) const
{
--
GitLab