From: heltai Date: Sun, 12 Apr 2009 08:58:18 +0000 (+0000) Subject: Added an image and cleared some details. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=12e06c39116b996608778872fdf94c92a56db309;p=dealii-svn.git Added an image and cleared some details. git-svn-id: https://svn.dealii.org/trunk@18597 0785d39b-7218-0410-832d-ea1e28bc413d --- 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/doc/step-34_3d-2.png b/deal.II/examples/step-34/doc/step-34_3d-2.png new file mode 100644 index 0000000000..cbbadc6461 Binary files /dev/null and b/deal.II/examples/step-34/doc/step-34_3d-2.png differ 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