Documentation: Advertise <deal.II/lac/linear_operator_tools.h>

Let's add some pointers to linear_operator_tools.h for the reader.
Maybe this helps to clarify the situation.

diff --git a/doc/doxygen/headers/laoperators.h b/doc/doxygen/headers/laoperators.h
index 0e09e31a0b4408fca06a1ac535514db416077aa5..d88a81437b754f29939d89b13229edf4ced27117 100644 (file)
--- a/doc/doxygen/headers/laoperators.h
+++ b/doc/doxygen/headers/laoperators.h
@@ -44,6 +44,8 @@
  * LinearOperator <code>op</code> that performs above computation when
  * applied on a vector, one can write:
  * @code
+ * #include <deal.II/lac/linear_operator_tools.h>
+ *
  * dealii::SparseMatrix<double> A, B, C;
  * double k;
  * // Setup and assembly...
@@ -116,6 +118,8 @@
  *
  * As an example consider the addition of multiple vectors:
  * @code
+ *   #include <deal.II/lac/linear_operator_tools.h>
+ *
  *   dealii::Vector<double> a, b, c, d;
  *   // ..
  *   dealii::Vector<double> result = a + b - c + d;
@@ -178,9 +182,16 @@
  *   Vector<double> residual =  b - op_a * x; // computes the residual at this point
  * @endcode
  *
- * The step-20 tutorial program has a detailed usage example of the
+ * @note The step-20 tutorial program has a detailed usage example of the
  * PackagedOperation class.
  *
+ * @note Many use cases of LinearOperator lead to intermediate expressions
+ * requiring a PackagedOperation. In order to include all necessary header
+ * files in one go consider using
+ * @code
+ * #include <deal.II/lac/linear_operator_tools.h>
+ * @endcode
+ *
  *
  * @ingroup LAC
  * @ingroup MATRICES

diff --git a/include/deal.II/lac/linear_operator.h b/include/deal.II/lac/linear_operator.h
index 726ac02564fb3306e6df555fcd55fb19cae4cf18..8e45ec668437670604f703cc2dfa4dd95ccaa9a6 100644 (file)
--- a/include/deal.II/lac/linear_operator.h
+++ b/include/deal.II/lac/linear_operator.h
@@ -113,6 +113,8 @@ identity_operator(const LinearOperator<Range, Domain, Payload> &);
  * one can write:
  *
  * @code
+ * #include <deal.II/lac/linear_operator_tools.h>
+ *
  * dealii::SparseMatrix<double> A, B, C;
  * const double k = ...;
  *
@@ -154,6 +156,15 @@ identity_operator(const LinearOperator<Range, Domain, Payload> &);
  * an interface that produces the result of this composite operation that is
  * compatible with Trilinos operator used by Trilinos solvers.
  *
+ * @note Many use cases of LinearOperator lead to intermediate expressions
+ * requiring a PackagedOperation. In order to include all necessary header
+ * files in one go consider using
+ * @code
+ * #include <deal.II/lac/linear_operator_tools.h>
+ * @endcode
+ *
+ * In order to use the full LinearOperator and PackagedOperation
+ *
  * @note To ensure that the correct payload is provided, wrapper functions
  * for linear operators have been provided within the respective
  * TrilinosWrappers (and, in the future, PETScWrappers) namespaces.