add ScaLAPACKMatrix.copy_transposed()

diff --git a/include/deal.II/lac/scalapack.h b/include/deal.II/lac/scalapack.h
index 9f701bfbc759ea0a49529b72c528f161f3468a1f..934b2fea7bca7c7845d7559aee7edf0d45786c24 100644 (file)
--- a/include/deal.II/lac/scalapack.h
+++ b/include/deal.II/lac/scalapack.h
@@ -191,6 +191,14 @@ public:
                const std::pair<unsigned int,unsigned int> &submatrix_size) const;
 
   /**
+   * Transposing assignment: <i>A = B<sup>T</sup></i>
+   *
+   * The matrices <tt>A</tt> and <tt>B</tt> must have the same process grid
+   *
+   * Following alignment conditions have to be fulfilled: MB_A = NB_B and NB_A = MB_B
+   */
+  void copy_transposed(const ScaLAPACKMatrix<NumberType> &B);
+
   /**
    * Matrix-addition:
    *

diff --git a/source/lac/scalapack.cc b/source/lac/scalapack.cc
index e815bc191184ff9ae0a7826f756ec0c060ca0b16..be2cd3391ff16a7c303315d51b2cfa0cdd49c0ac 100644 (file)
--- a/source/lac/scalapack.cc
+++ b/source/lac/scalapack.cc
@@ -441,6 +441,10 @@ ScaLAPACKMatrix<NumberType>::copy_to (ScaLAPACKMatrix<NumberType> &dest) const
 
 
 template <typename NumberType>
+void ScaLAPACKMatrix<NumberType>::copy_transposed(const ScaLAPACKMatrix<NumberType> &B)
+{
+  add(B,0,1,true);
+}