/**
* @relates LinearOperator
*
- * This function implement a forward substitution argument to invert a lower
+ * This function implements a forward substitution argument to invert a lower
* block triangular matrix.
- * It takes as argement an array of array of LinearOperators @p block_matrix
+ * As argument, it takes an array of array of LinearOperators @p block_matrix
* representing a lower block triangular matrix and an array of LinearOperators
* @p inverse_diagonal representing inverses of digonal blocks of @p block_matrix.
*
- * Let us assume we have a linear system where each coefficient of the system is a
- * matrix:
+ * Let us assume we have a linear system with the following block structure:
+ *
* A00 x0 + ... = y0
* A01 x0 + A11 x1 + ... = y1
* ... ...
* and therefore:
* xn = Ann^-1 ( yn - A0n x0 - ... - A(n-1)n x(n-1) )
*
- * Notice that we are not using the whole matrix: just the lower triangular block
- * matrix obtained from @p block_matrix it is used.
+ * @note Notice that we are not using the whole matrix:
+ * just the lower triangular block matrix obtained
+ * from @p block_matrix it is used.
*
- * Caveat: Tvmult and Tvmult_add have not been implemented, yet. This may lead to mistakes.
* @ingroup LAOperators
*/
template < unsigned int n,
/**
* @relates LinearOperator
*
- * This function implement a back substitution argument to invert an upper
+ * This function implements a back substitution argument to invert an upper
* block triangular matrix.
- * It takes as argement an array of array of LinearOperators @p block_matrix
+ * As argument, it takes an array of array of LinearOperators @p block_matrix
* representing an upper block triangular matrix and an array of LinearOperators
* @p inverse_diagonal representing inverses of digonal blocks of @p block_matrix.
*
- * Let us assume we have a linear system where each coefficient of the system is a
- * matrix:
+ * Let us assume we have a linear system with the following block structure:
*
* A00 x0 + A01 x1 + ... + A0n xn = yn
* A11 x1 + ... = y1
* and therefore:
* x0 = A00^-1 ( y0 - A0n xn - ... - A01 x1 )
*
- * Notice that we are not using the whole matrix: just the upper triangular block
- * matrix obtained from @p block_matrix it is used.
+ * @note Notice that we are not using the whole matrix:
+ * just the upper triangular block matrix obtained
+ * from @p block_matrix it is used.
*
- * Caveat: Tvmult and Tvmult_add have not been implemented, yet. This may lead to mistakes.
* @ingroup LAOperators
*/
template < unsigned int n,