From 1a17cc4024d072d1812d982aa1c3cbe9ccff8e5a Mon Sep 17 00:00:00 2001
From: wolf <wolf@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 2 Jun 2003 22:10:33 +0000
Subject: [PATCH] FullMatrix::operator /=

git-svn-id: https://svn.dealii.org/trunk@7725 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/lac/include/lac/full_matrix.h           |  6 ++++++
 deal.II/lac/include/lac/full_matrix.templates.h | 17 +++++++++++++++++
 2 files changed, 23 insertions(+)

diff --git a/deal.II/lac/include/lac/full_matrix.h b/deal.II/lac/include/lac/full_matrix.h
index 72846570b5..9efe89c773 100644
--- a/deal.II/lac/include/lac/full_matrix.h
+++ b/deal.II/lac/include/lac/full_matrix.h
@@ -239,6 +239,12 @@ class FullMatrix : public Table<2,number>
 				      * fixed factor.
 				      */
     FullMatrix & operator *= (const double factor);
+
+				     /**
+				      * Scale the entire matrix by the
+				      * inverse of the given factor.
+				      */
+    FullMatrix & operator /= (const double factor);
     
 				     /**
 				      * Weighted addition. The matrix
diff --git a/deal.II/lac/include/lac/full_matrix.templates.h b/deal.II/lac/include/lac/full_matrix.templates.h
index 95a6e8c3a4..7d15b5c348 100644
--- a/deal.II/lac/include/lac/full_matrix.templates.h
+++ b/deal.II/lac/include/lac/full_matrix.templates.h
@@ -106,6 +106,23 @@ FullMatrix<number>::operator *= (const double factor)
 
 
 
+template <typename number>
+FullMatrix<number> &
+FullMatrix<number>::operator /= (const double factor)
+{
+  number       *p = &this->el(0,0);
+  const number *e = &this->el(0,0) + n()*m();
+
+  const number factor_inv = 1./factor;
+  
+  while (p != e)
+    *p++ *= factor_inv;
+
+  return *this;
+}
+
+
+
 template <typename number>
 template <typename number2>
 void
-- 
2.39.5