]> https://gitweb.dealii.org/ - dealii.git/commitdiff
add ScaLAPACKMatrix.copy_transposed()
authorBenjamin Brands <benjamin.brands@fau.de>
Tue, 6 Feb 2018 15:41:17 +0000 (16:41 +0100)
committerBenjamin Brands <benjamin.brands@fau.de>
Wed, 14 Feb 2018 20:10:15 +0000 (21:10 +0100)
include/deal.II/lac/scalapack.h
source/lac/scalapack.cc

index 9f701bfbc759ea0a49529b72c528f161f3468a1f..934b2fea7bca7c7845d7559aee7edf0d45786c24 100644 (file)
@@ -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:
    *
index e815bc191184ff9ae0a7826f756ec0c060ca0b16..be2cd3391ff16a7c303315d51b2cfa0cdd49c0ac 100644 (file)
@@ -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);
+}
 
 
 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.