From: guido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 29 Jan 2001 11:32:51 +0000 (+0000)
Subject: interpolation matrices in FE_Q, starting FE_DGQ
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=1948152b7e3f04ab73db4f219a55538b185d93ee;p=dealii-svn.git

interpolation matrices in FE_Q, starting FE_DGQ


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

diff --git a/deal.II/lac/include/lac/full_matrix.h b/deal.II/lac/include/lac/full_matrix.h
index 19c9b94fd6..ed0150ecba 100644
--- a/deal.II/lac/include/lac/full_matrix.h
+++ b/deal.II/lac/include/lac/full_matrix.h
@@ -19,6 +19,7 @@
 #include <base/exceptions.h>
 #include <base/subscriptor.h>
 
+#include <vector>
 
 // forward declarations
 template<typename number> class Vector;
@@ -157,6 +158,29 @@ class FullMatrix : public Subscriptor
 	       const unsigned int         j=0);
     
 
+				     /**
+				      * Fill with permutation of another matrix.
+				      *
+				      * The matrix @p{src} is copied
+				      * into the target. The two
+				      * permutation @p{p_r} and
+				      * @p{p_c} operate in a way, such
+				      * that @{result(i,j) =
+				      * src(p_r[i], p_c[j])}.
+				      *
+				      * The vectors may also be a
+				      * selection from a larger set of
+				      * integers, if the matrix
+				      * @p{src} is bigger. It is also
+				      * possible to duplicate rows or
+				      * columns by this method.
+				      */
+    template<typename number2>
+    void fill_permutation (const FullMatrix<number2> &src,
+	       const std::vector<unsigned int>& p_rows,
+	       const std::vector<unsigned int>& p_cols);
+    
+
                                      /**
 				      * Fill matrix with an array of numbers.
 				      * The array is arranged line by line. No
diff --git a/deal.II/lac/include/lac/full_matrix.templates.h b/deal.II/lac/include/lac/full_matrix.templates.h
index 53aea79ee2..58bf9e8eca 100644
--- a/deal.II/lac/include/lac/full_matrix.templates.h
+++ b/deal.II/lac/include/lac/full_matrix.templates.h
@@ -457,6 +457,24 @@ void FullMatrix<number>::fill (const FullMatrix<number2>& src,
 }
 
 
+template <typename number>
+template <typename number2>
+void FullMatrix<number>::fill_permutation (const FullMatrix<number2> &src,
+					   const std::vector<unsigned int>& p_rows,
+					   const std::vector<unsigned int>& p_cols)
+{
+  Assert (p_rows.size() == dim_image,
+	  ExcDimensionMismatch (p_rows.size(), dim_image));
+  Assert (p_cols.size() == dim_range,
+	  ExcDimensionMismatch (p_cols.size(), dim_range));
+
+  for (unsigned int i=0;i<dim_image;++i)
+    for (unsigned int j=0;j<dim_range;++j)
+      el(i,j) = src(p_rows[i], p_cols[j]);
+}
+
+
+
 template <typename number>
 template <typename number2>
 void FullMatrix<number>::fill (const number2* entries)