From f5964adb4fdfe64d29168511fd65212358d711f1 Mon Sep 17 00:00:00 2001 From: Wolfgang Bangerth Date: Mon, 19 Oct 2020 17:19:08 -0600 Subject: [PATCH] Update documentation in FEFieldFunction. --- include/deal.II/numerics/fe_field_function.h | 17 ++++++++++------- 1 file changed, 10 insertions(+), 7 deletions(-) diff --git a/include/deal.II/numerics/fe_field_function.h b/include/deal.II/numerics/fe_field_function.h index 677c5de61a..6950786715 100644 --- a/include/deal.II/numerics/fe_field_function.h +++ b/include/deal.II/numerics/fe_field_function.h @@ -87,22 +87,25 @@ namespace Functions * DoFHandler dh1 (tria_1); * Vector solution_1; * - * // Do the same with the second + * // On this first domain, set up the various data structures, + * // assemble matrices, solve the linear system, and get a Nobel + * // prize for the work we have done here: + * [...] + * + * // Then create a DoFHandler and solution vector for the second domain: * DoFHandler dh2 (tria_2); * Vector solution_2; * - * // Setup the system, assemble matrices, solve problems and get the - * // nobel prize on the first domain... - * - * // Now project it to the second domain + * // Finally, project the solution on the first domain onto the + * // second domain, assuming that this does not require querying + * // values from outside the first domain: * FEFieldFunction fe_function_1 (dh_1, solution_1); * VectorTools::project (dh_2, constraints_2, quad, * fe_function_1, solution_2); * - * // Or interpolate it... + * // Alternatively, we could have also interpolated it: * Vector solution_3; * VectorTools::interpolate (dh_2, fe_function_1, solution_3); - * * @endcode * * The snippet of code above will work assuming that the second -- 2.39.5