From: Guido Kanschat <dr.guido.kanschat@gmail.com>
Date: Wed, 27 Jun 2001 17:10:10 +0000 (+0000)
Subject: shifted matrices in a separate file
X-Git-Tag: v8.0.0~19029
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=8b2ea0b381caae3e978c28368136267710f27d5a;p=dealii.git

shifted matrices in a separate file


git-svn-id: https://svn.dealii.org/trunk@4766 0785d39b-7218-0410-832d-ea1e28bc413d
---

diff --git a/deal.II/lac/include/lac/eigen.h b/deal.II/lac/include/lac/eigen.h
index c2886da3ba..fb43f77a54 100644
--- a/deal.II/lac/include/lac/eigen.h
+++ b/deal.II/lac/include/lac/eigen.h
@@ -14,65 +14,13 @@
 #define __deal2__eigen_h
 
 
+#include <lac/shifted_matrix.h>
 #include <lac/solver.h>
 #include <lac/solver_control.h>
 #include <lac/solver_cg.h>
 #include <lac/vector_memory.h>
 #include <lac/precondition.h>
 
-/**
- * Matrix with shifted diagonal values.
- *
- * Given a matrix $A$, this class implements a matrix-vector product with
- * $A+\sigma I$, where sigma is a provided shift parameter.
- *
- * @author Guido Kanschat, 2000
- */
-template<class MATRIX>
-class ShiftedMatrix
-{
-  public:
-				     /**
-				      * Constructor.
-				      * Provide the base matrix and a shift parameter.
-				      */
-    ShiftedMatrix (const MATRIX& A, const double sigma);
-
-				     /**
-				      * Set the shift parameter.
-				      */
-    void shift (const double sigma);
-
-				     /**
-				      * Access to the shift parameter.
-				      */
-    double shift () const;
-
-				     /**
-				      * Matrix-vector-product.
-				      */
-    template <class VECTOR>
-    void vmult (VECTOR& dst, const VECTOR& src) const;
-
-				     /**
-				      * Residual.
-				      */
-    template <class VECTOR>
-    double residual (VECTOR& dst, const VECTOR& src, const VECTOR& rhs) const;
-    
-  private:
-				     /**
-				      * Storage for base matrix.
-				      */
-    const MATRIX& A;
-
-				     /**
-				      * Shift parameter.
-				      */
-    double sigma;
-};
-    
-
 /**
  * Power method (von Mises).
  *
@@ -230,58 +178,6 @@ class EigenInverse : private Solver<VECTOR>
     AdditionalData additional_data;
 };
 
-//----------------------------------------------------------------------//
-
-template <class MATRIX>
-inline
-ShiftedMatrix<MATRIX>::ShiftedMatrix (const MATRIX& A, const double sigma)
-		:
-		A(A), sigma(sigma)
-{}
-
-
-
-template <class MATRIX>
-inline void
-ShiftedMatrix<MATRIX>::shift (const double s)
-{
-  sigma = s;
-}
-
-
-template <class MATRIX>
-inline double
-ShiftedMatrix<MATRIX>::shift () const
-{
-  return sigma;
-}
-
-
-
-template <class MATRIX>
-template <class VECTOR>
-inline void
-ShiftedMatrix<MATRIX>::vmult (VECTOR& dst, const VECTOR& src) const
-{
-  A.vmult(dst, src);
-  dst.add(sigma, src);
-}
-
-
-template <class MATRIX>
-template <class VECTOR>
-inline double
-ShiftedMatrix<MATRIX>::residual (VECTOR& dst,
-				 const VECTOR& src,
-				 const VECTOR& rhs) const
-{
-  A.vmult(dst, src);
-  dst.add(sigma, src);
-  dst.sadd(-1.,1.,rhs);
-  return dst.l2_norm ();
-}
-
-
 //----------------------------------------------------------------------
 
 
diff --git a/deal.II/lac/include/lac/shifted_matrix.h b/deal.II/lac/include/lac/shifted_matrix.h
new file mode 100644
index 0000000000..57849269c9
--- /dev/null
+++ b/deal.II/lac/include/lac/shifted_matrix.h
@@ -0,0 +1,242 @@
+//----------------------------  eigen.h  ---------------------------
+//    $Id$
+//    Version: $Name$
+//
+//    Copyright (C) 2001 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  eigen.h  ---------------------------
+#ifndef __deal2__shifted_matrix_h
+#define __deal2__shifted_matrix_h
+
+
+#include <base/smartpointer.h>
+
+/**
+ * Matrix with shifted diagonal values.
+ *
+ * Given a matrix p{A}, this class implements a matrix-vector product with
+ * @p{A+\sigma I}, where sigma is a provided shift parameter.
+ *
+ * @author Guido Kanschat, 2000, 2001
+ */
+template<class MATRIX>
+class ShiftedMatrix
+{
+  public:
+				     /**
+				      * Constructor.
+				      * Provide the base matrix and a shift parameter.
+				      */
+    ShiftedMatrix (const MATRIX& A, const double sigma);
+
+				     /**
+				      * Set the shift parameter.
+				      */
+    void shift (const double sigma);
+
+				     /**
+				      * Access to the shift parameter.
+				      */
+    double shift () const;
+
+				     /**
+				      * Matrix-vector-product.
+				      */
+    template <class VECTOR>
+    void vmult (VECTOR& dst, const VECTOR& src) const;
+
+				     /**
+				      * Residual.
+				      */
+    template <class VECTOR>
+    double residual (VECTOR& dst, const VECTOR& src, const VECTOR& rhs) const;
+    
+  private:
+				     /**
+				      * Storage for base matrix.
+				      */
+    SmartPointer<const MATRIX> A;
+
+				     /**
+				      * Shift parameter.
+				      */
+    double sigma;
+};
+    
+
+
+/**
+ * Matrix with shifted diagonal values with respect to a certain scalar product.
+ *
+ * Given a matrix @p{A}, this class implements a matrix-vector product
+ * with @p{A+\sigma M}, where sigma is a provided shift parameter and
+ * @p{M} is the matrix representing the identity
+ *
+ * @author Guido Kanschat, 2001
+ */
+template<class MATRIX, class MASSMATRIX>
+class ShiftedMatrixGeneralized
+{
+  public:
+				     /**
+				      * Constructor.
+				      * Provide the base matrix and a shift parameter.
+				      */
+    ShiftedMatrixGeneralized (const MATRIX& A,
+			      const MASSMATRIX& M,
+			      const double sigma);
+
+				     /**
+				      * Set the shift parameter.
+				      */
+    void shift (const double sigma);
+
+				     /**
+				      * Access to the shift parameter.
+				      */
+    double shift () const;
+
+				     /**
+				      * Matrix-vector-product.
+				      */
+    template <class VECTOR>
+    void vmult (VECTOR& dst, const VECTOR& src) const;
+
+				     /**
+				      * Residual.
+				      */
+    template <class VECTOR>
+    double residual (VECTOR& dst, const VECTOR& src, const VECTOR& rhs) const;
+    
+  private:
+				     /**
+				      * Storage for base matrix.
+				      */
+    SmartPointer<const MATRIX> A;
+				     /**
+				      * Storage for mass matrix.
+				      */
+    SmartPointer<const MASSMATRIX> M;
+
+				     /**
+				      * Shift parameter.
+				      */
+    double sigma;
+};
+    
+
+
+//----------------------------------------------------------------------//
+
+template <class MATRIX>
+inline
+ShiftedMatrix<MATRIX>::ShiftedMatrix (const MATRIX& A, const double sigma)
+		:
+		A(&A), sigma(sigma)
+{}
+
+
+
+template <class MATRIX>
+inline void
+ShiftedMatrix<MATRIX>::shift (const double s)
+{
+  sigma = s;
+}
+
+
+template <class MATRIX>
+inline double
+ShiftedMatrix<MATRIX>::shift () const
+{
+  return sigma;
+}
+
+
+
+template <class MATRIX>
+template <class VECTOR>
+inline void
+ShiftedMatrix<MATRIX>::vmult (VECTOR& dst, const VECTOR& src) const
+{
+  A.vmult(dst, src);
+  dst.add(sigma, src);
+}
+
+
+template <class MATRIX>
+template <class VECTOR>
+inline double
+ShiftedMatrix<MATRIX>::residual (VECTOR& dst,
+				 const VECTOR& src,
+				 const VECTOR& rhs) const
+{
+  A.vmult(dst, src);
+  dst.add(sigma, src);
+  dst.sadd(-1.,1.,rhs);
+  return dst.l2_norm ();
+}
+
+
+//----------------------------------------------------------------------//
+
+template <class MATRIX, class MASSMATRIX>
+inline
+ShiftedMatrixGeneralized<MATRIX, MASSMATRIX>
+::ShiftedMatrixGeneralized (const MATRIX& A,
+			    const MASSMATRIX& M,
+			    const double sigma)
+		:
+		A(&A), M(&M), sigma(sigma)
+{}
+
+
+
+template <class MATRIX, class MASSMATRIX>
+inline void
+ShiftedMatrixGeneralized<MATRIX, MASSMATRIX>::shift (const double s)
+{
+  sigma = s;
+}
+
+
+template <class MATRIX, class MASSMATRIX>
+inline double
+ShiftedMatrixGeneralized<MATRIX, MASSMATRIX>::shift () const
+{
+  return sigma;
+}
+
+
+
+template <class MATRIX, class MASSMATRIX>
+template <class VECTOR>
+inline void
+ShiftedMatrixGeneralized<MATRIX, MASSMATRIX>::vmult (VECTOR& dst,
+						     const VECTOR& src) const
+{
+  A.vmult(dst, src);
+  dst.add(sigma, src);
+}
+
+
+template <class MATRIX, class MASSMATRIX>
+template <class VECTOR>
+inline double
+ShiftedMatrixGeneralized<MATRIX, MASSMATRIX>::residual (VECTOR& dst,
+							const VECTOR& src,
+							const VECTOR& rhs) const
+{
+  A.vmult(dst, src);
+  dst.add(sigma, src);
+  dst.sadd(-1.,1.,rhs);
+  return dst.l2_norm ();
+}
+
+
+#endif