From 815cac0f4103c1275df145a16489ee7c0641bc95 Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Thu, 12 Dec 2013 21:52:20 +0000
Subject: [PATCH] add axpy and operator*=
[[Imported from SVN: r9586]]
---
dune/gfe/symmetricmatrix.hh | 12 +++++++++++-
1 file changed, 11 insertions(+), 1 deletion(-)
diff --git a/dune/gfe/symmetricmatrix.hh b/dune/gfe/symmetricmatrix.hh
index 05c39c6e..3ce94c90 100644
--- a/dune/gfe/symmetricmatrix.hh
+++ b/dune/gfe/symmetricmatrix.hh
@@ -37,7 +37,13 @@ public:
return *this;
}
- /** \brief Matrix style random read/write access to components
+ SymmetricMatrix<T,N>& operator*=(const T& s)
+ {
+ data_ *= s;
+ return *this;
+ }
+
+ /** \brief Matrix style random read/write access to components
* \param i line index
* \param j column index
* \note You need to know what you are doing: You can only access the lower
@@ -71,6 +77,10 @@ public:
return result;
}
+ void axpy(const T& a, const SymmetricMatrix<T,N>& other)
+ {
+ data_.axpy(a,other.data_);
+ }
private:
Dune::FieldVector<T,N*(N+1)/2> data_;
--
GitLab