From 40cec9b1d6215e667188c245d523ede6684ce61d Mon Sep 17 00:00:00 2001 From: Luca Heltai Date: Sun, 12 Apr 2009 08:58:18 +0000 Subject: [PATCH] Added an image and cleared some details. git-svn-id: https://svn.dealii.org/trunk@18597 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-34/doc/results.dox | 8 ++++++-- deal.II/examples/step-34/step-34.cc | 25 ++++++++++++------------ 2 files changed, 19 insertions(+), 14 deletions(-) diff --git a/deal.II/examples/step-34/doc/results.dox b/deal.II/examples/step-34/doc/results.dox index e6849a3bb0..57e6ab1539 100644 --- a/deal.II/examples/step-34/doc/results.dox +++ b/deal.II/examples/step-34/doc/results.dox @@ -145,7 +145,11 @@ like @image html step-34_2d.png -while in three dimensions we only show the potential on the surface, -together with a countur plot: +while in three dimensions we show first the potential on the surface, +together with a contur plot, @image html step-34_3d.png + +and then the external contur plot of the potential, with opacity set to 25%: + +@image html step-34_3d-2.png diff --git a/deal.II/examples/step-34/step-34.cc b/deal.II/examples/step-34/step-34.cc index b60f78c57a..f23e21a8fa 100644 --- a/deal.II/examples/step-34/step-34.cc +++ b/deal.II/examples/step-34/step-34.cc @@ -1,5 +1,5 @@ //---------------------------- step-34.cc --------------------------- -// $Id: testsuite.html 13373 2006-07-13 13:12:08Z manigrasso $ +// $Id: testsuite.html 13373 2006-07-13 13:12:08Z heltai $ // Version: $Name$ // // Copyright (C) 2009 by the deal.II authors @@ -446,7 +446,7 @@ void BEMProblem::read_parameters (const std::string filename) { // "Run 3d simulation" flag to false. // // This is another example of how to use parameter files in - // dimensio independent programming. + // dimension independent programming. run_in_this_dimension = prm.get_bool("Run " + Utilities::int_to_string(dim) + "d simulation"); @@ -779,11 +779,12 @@ void BEMProblem::assemble_system() { } // One quick way to compute the diagonal matrix of the solid // angles, is to use the Neumann matrix itself. It is enough to - // multiply the matrix with the vector of ones, to get the - // diagonal matrix of the alpha solid angles. + // multiply the matrix with a vector of elements all equal to -1., + // to get the diagonal matrix of the alpha angles, or solid + // angles. Vector ones(dh.n_dofs()); - for(unsigned int i=0; i::output_results(unsigned int cycle) { std::ofstream file(filename.c_str()); dataout.write_vtk(file); - - convergence_table.set_precision("L2(phi)", 3); - convergence_table.set_precision("Linfty(alpha)", 3); - - convergence_table.set_scientific("L2(phi)", true); - convergence_table.set_scientific("Linfty(alpha)", true); if(cycle == n_cycles-1) { + convergence_table.set_precision("L2(phi)", 3); + convergence_table.set_precision("Linfty(alpha)", 3); + + convergence_table.set_scientific("L2(phi)", true); + convergence_table.set_scientific("Linfty(alpha)", true); + convergence_table .evaluate_convergence_rates("L2(phi)", ConvergenceTable::reduction_rate_log2); convergence_table -- 2.39.5