From: Martin Kronbichler Date: Fri, 19 Jun 2009 09:01:53 +0000 (+0000) Subject: Update comments for mmult. X-Git-Tag: v8.0.0~7609 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f96d3e9ef721c5ccf0f7d053e8f819a4f06bce38;p=dealii.git Update comments for mmult. git-svn-id: https://svn.dealii.org/trunk@18941 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 c73046ef6c..1b60fed51a 100644 --- a/deal.II/lac/include/lac/full_matrix.h +++ b/deal.II/lac/include/lac/full_matrix.h @@ -678,7 +678,7 @@ class FullMatrix : public Table<2,number> * inverse of the given factor. */ FullMatrix & operator /= (const number factor); - + /** * Simple addition of a scaled * matrix, i.e. *this += @@ -694,7 +694,7 @@ class FullMatrix : public Table<2,number> template void add (const number a, const FullMatrix &A); - + /** * Multiple addition of scaled * matrices, i.e. *this += @@ -713,7 +713,7 @@ class FullMatrix : public Table<2,number> const FullMatrix &A, const number b, const FullMatrix &B); - + /** * Multiple addition of scaled * matrices, i.e. *this += @@ -734,7 +734,7 @@ class FullMatrix : public Table<2,number> const FullMatrix &B, const number c, const FullMatrix &C); - + /** * Add rectangular block. * @@ -759,7 +759,7 @@ class FullMatrix : public Table<2,number> const unsigned int dst_offset_j = 0, const unsigned int src_offset_i = 0, const unsigned int src_offset_j = 0); - + /** * Weighted addition of the * transpose of B to @@ -770,7 +770,7 @@ class FullMatrix : public Table<2,number> template void Tadd (const number s, const FullMatrix &B); - + /** * Add transose of a rectangular block. * @@ -1023,9 +1023,17 @@ class FullMatrix : public Table<2,number> * if (!adding) * C = A*B * - * Assumes that A and B have - * compatible sizes and that C - * already has the right size. + * Assumes that A and + * B have compatible sizes and + * that C already has the + * right size. + * + * This function uses the BLAS function + * Xgemm if the calling matrix has more + * than 15 rows and BLAS was detected + * during configuration. Using BLAS + * usually results in considerable + * performance gains. */ template void mmult (FullMatrix &C, @@ -1051,6 +1059,13 @@ class FullMatrix : public Table<2,number> * B have compatible * sizes and that C * already has the right size. + * + * This function uses the BLAS function + * Xgemm if the calling matrix has more + * than 15 columns and BLAS was + * detected during configuration. Using + * BLAS usually results in considerable + * performance gains. */ template void Tmmult (FullMatrix &C,