From 97c9c7ac51135f96dac04a6c3da4300ebff37501 Mon Sep 17 00:00:00 2001
From: kanschat <kanschat@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Thu, 6 May 2010 22:01:02 +0000
Subject: [PATCH] adding some results

git-svn-id: https://svn.dealii.org/trunk@21081 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/examples/step-39/Makefile             |   2 +-
 deal.II/examples/step-39/doc/results.dox      |  87 ++++++++++++++++++
 .../step-39/doc/step-39-convergence.png       | Bin 0 -> 4223 bytes
 deal.II/examples/step-39/step-39.cc           |   4 +-
 4 files changed, 90 insertions(+), 3 deletions(-)
 create mode 100644 deal.II/examples/step-39/doc/step-39-convergence.png

diff --git a/deal.II/examples/step-39/Makefile b/deal.II/examples/step-39/Makefile
index 4a33ab61d8..01ea1089e2 100644
--- a/deal.II/examples/step-39/Makefile
+++ b/deal.II/examples/step-39/Makefile
@@ -98,7 +98,7 @@ run: $(target)
 	@echo ============================ Running $<
 	@./$(target)$(EXEEXT)
 	@echo ============================ Postprocessing to step-39.dat
-	@$(PERL) postprocess.pl step-39.log > step-39.dat
+	@$(PERL) postprocess.pl deallog > output.dat
 
 # As a last rule to the `make' program, we define what to do when
 # cleaning up a directory. This usually involves deleting object files
diff --git a/deal.II/examples/step-39/doc/results.dox b/deal.II/examples/step-39/doc/results.dox
index a79c9a63f4..236939ef23 100644
--- a/deal.II/examples/step-39/doc/results.dox
+++ b/deal.II/examples/step-39/doc/results.dox
@@ -1,3 +1,90 @@
 <h1>Results</h1>
 
+<h2>Logfile output</h2>
+First, the program produces the usual logfile here stored in <tt>deallog</tt>. It reads (with ommission of intermediate steps)
+<pre>
+DEAL::Element: FE_DGQ<2>(3)
+DEAL::Step 0
+DEAL::Triangulation 4 cells, 1 levels
+DEAL::DoFHandler 64 dofs, level dofs 64
+DEAL::Assemble matrix
+DEAL::Assemble multilevel matrix
+DEAL::Assemble right hand side
+DEAL::Solve
+DEAL:cg::Starting value 27.1275
+DEAL:cg::Convergence step 1 value 1.97998e-14
+DEAL::Error    0.161536
+DEAL::Estimate 1.35839
+DEAL::Writing solution to <sol-00.gnuplot>...
+DEAL::
+DEAL::Step 1
+DEAL::Triangulation 10 cells, 2 levels
+DEAL::DoFHandler 160 dofs, level dofs 64 128
+DEAL::Assemble matrix
+DEAL::Assemble multilevel matrix
+DEAL::Assemble right hand side
+DEAL::Solve
+DEAL:cg::Starting value 35.5356
+DEAL:cg::Convergence step 10 value 1.24790e-13
+DEAL::Error    0.164997
+DEAL::Estimate 1.08528
+DEAL::Writing solution to <sol-01.gnuplot>...
+DEAL::
+DEAL::Step 2
+DEAL::Triangulation 16 cells, 2 levels
+DEAL::DoFHandler 256 dofs, level dofs 64 256
+DEAL::Assemble matrix
+DEAL::Assemble multilevel matrix
+DEAL::Assemble right hand side
+DEAL::Solve
+DEAL:cg::Starting value 37.0552
+DEAL:cg::Convergence step 10 value 4.74024e-14
+DEAL::Error    0.113596
+DEAL::Estimate 0.990460
+DEAL::Writing solution to <sol-02.gnuplot>...
 
+...
+
+DEAL::Step 10
+DEAL::Triangulation 124 cells, 9 levels
+DEAL::DoFHandler 1984 dofs, level dofs 64 256 512 512 256 256 256 256 256
+DEAL::Assemble matrix
+DEAL::Assemble multilevel matrix
+DEAL::Assemble right hand side
+DEAL::Solve
+DEAL:cg::Starting value 38.5798
+DEAL:cg::Convergence step 12 value 7.56337e-14
+DEAL::Error    0.0101279
+DEAL::Estimate 0.0957571
+DEAL::Writing solution to <sol-10.gnuplot>...
+DEAL::
+DEAL::Step 11
+DEAL::Triangulation 163 cells, 10 levels
+DEAL::DoFHandler 2608 dofs, level dofs 64 256 768 576 512 256 256 256 256 256
+DEAL::Assemble matrix
+DEAL::Assemble multilevel matrix
+DEAL::Assemble right hand side
+DEAL::Solve
+DEAL:cg::Starting value 44.1721
+DEAL:cg::Convergence step 12 value 9.49720e-14
+DEAL::Error    0.00716966
+DEAL::Estimate 0.0681646
+DEAL::Writing solution to <sol-11.gnuplot>...
+</pre>
+
+This log for instance shows that the number of conjugate gradient
+iteration steps is constant at approximately 12.
+
+<h2>Postprocessing of the logfile</h2>
+
+<img align="right" alt="Convergence graphs" src="step-39-convergence.png">
+Using the perl script <tt>postprocess.pl</tt>, we extract relevant
+data into <tt>output.dat</tt>, which can be used to plot graphs with
+<tt>gnuplot</tt>. The graph on the right for instance was produced with
+<pre>
+set style data linespoints
+set logscale
+set xrange [50:3000]
+plot "output.dat" using 2:3 title "error", "" using 2:4 title "estimate", \
+     "" using 2:(3000*$2**-1.5) title "3rd order"
+</pre>
diff --git a/deal.II/examples/step-39/doc/step-39-convergence.png b/deal.II/examples/step-39/doc/step-39-convergence.png
new file mode 100644
index 0000000000000000000000000000000000000000..3f85cdf0788236ab6e4c83ed9f8a9859203acb2c
GIT binary patch
literal 4223
zcma)9dpJ~G``$D5Ff+m&oQc6OI?zNpmCVRti0E7@ZBh;u6UkZIkcbc^6`J_c0U5nH
zPi9Z)98pPUlPFE)EJEhn-uL^x|NO4+yWVxJz1Oqu=eeKfUf0@dU)M_aUhAqv)FuJ|
zpybAJS`PqN7yvL5yc~*PRTjTQCkwniHag2>G8EzS`EvsT0I&=JWWc}x!a`pFjfTh&
z8UVBSuncAaGFS#65P%>61S1d(zyM1IGa(p;02bO{GzDstU<iUy?_4?=00Ex_2n&HB
zJ_4aUgpWE<CPYJ{AR5|OEE)?9r$KWxw6SO~4INNP{B{XKqX95l3xE*`%12-Tfgu^H
z1g(b1WT=gL`KTfU6^brF2B07xLHIHVnw!9Yj1RyP2_Nk;8Y&0@5I&8@r^#S6F^x7?
zjRnKdTvixOj-V8ljE{Cm#^=h=tbb;d$&igE!UzCRA-g%TH%1r!c%PotW24IIaPpSP
z9vlAz3tJA)SWOpeMS2}j^}k+t)V06__`l}ye5J!|{ESk6Y)aI0GWKFvF`iWBi!rYi
zTY7&FfJgmt?}U}T8k7q&voD0QzI!2It0ulEH+xNX<b7D{<`64#m~0}LAeP;(@KJ}Y
zB*^Id{$k0vxLmlas&yj7Xf2!-lQ87-p^6%JN^`^8fCnmf@ub`4O}mU)cq7)$AzSIp
zOwDYsJl09>sby%~5&NB46=8K>6X@bldoQ8AOk>HckG)p@jXG?&YNm&%&r=B-K+W{{
z_7YEe<9lD{sS4#nMbm5fD9@f+CLWuzAdBahRtjkvu>x{mH<M|~@mbMzqS~*Gla>kJ
zBBR`AY9;r41MZjiM85fv4Gb@`4IY|`a6fWj^vG;vokJ@i^OfDMt0E_NnnVw;fAqBL
zC8M;;^`)rd>-=s}jcD3oXke6g;mfh-Uk3X!lGHH+i<6ob@kfFSL=IDmP|Ve&uMDI`
zTEFFFX#yfTmMS6{CK^qr6wJJ78x=Qa@xCrety-P&)GBJ*hpF~|#csH!B}vfF0jmSJ
zsvf*armYHx7roSu?8%E}#+jsmpH6DJG7A;x;vg%TEL0z-zw-B<m1s=${>~c{jX9;4
zMVZIQoc8MOq~u@Id5^<B)SphKhJAj{r1x){ts%WUFA*7i>wY$(*Y#VIFf}&qE}3?H
ztVoYfh{~N%8%hQ3-{6OScML{dw&_-gH((#JlAGCviBfh{KE^5nIXvacC2>QMq}<82
z#$Dk9t7?U26y1}hfE(B8I+iXT)?B4_#)W5eFeX@;qMLpP#p&j9bn!tOr5p`G_U!}k
zV&3i{9U4Jd90Aj9fNo4KG$B;r@QOzDX@tXOjxhaIwI&MIm7ukX!mXA5NCuleg=_Cs
zARK0FM!iKrs5h_&$KmZv*F-7wy(oC%3JNlM<WVq97X@jDVY;U87!WhQo;xfTK%Q5b
z47f3Gn4=2o?Fczwsnmo;$-IinW|qv!jy62e8)zWrl>9Z=&C3CezFMFmrO(lj<DuvZ
zH0_j8u-}eF8LJuv+}fK4<ml{fd3+9NuVaYB@s5QR%~MPUqt^`ew|?^TD*y%w&n~R^
z1Wrys3D$x#O-bbNWx`YhCjuv2Tl@R{tN}Z4qqTU;+gUwXE`L$uG~CscD!M@=*!fpu
zeKoXar?B+?)bUmPO&Rq5l+seD@90TzCpT8l`;s&95E2Y4@)V6S`1Uf4`SR?OVC7vR
z_v*R~deEX!SaF-;r_W1C1Rtp(a_NV9-30=jz2ZHte&`l+U?r8q+Z?OR{Th4otO;8u
zuymX6`v(!Q^`P3&RNwIHWU%$}!>L0S=WM;<x6^PHNp@HzAoKD=yGZkMw(Z1%KX~>U
z;xS_I3eT72W#=WY?;0=katl2cdEsArt{G5I1>8^@g$r9mBU`z^A}dP#pxg~rD<^ot
zRZ#TPAYE`GPr;zPGTH&Nw^*Hid$h61J<;-PVA*TE3mtu%M$SJx5}JZ9Crh$QE`Qld
zEGJ2>KMPDt7pP+NO!e183l_8CW7j?hq%o<kVy%U)HTh(VHF+)TxMk8s73O7rUzZ#O
zEQ=a$OYI%^2Iu0LZyO#fT6pBSzzuD_e&fwpix=Z3fYGAwgZpxmlm~5yZHJ}A@`t68
zfMcuYFNshs^co<R@0y<#DfEE*Qhdd~2Dn!~?ZH+U{4(qjA66y^vVwTeSp*7XRt*kL
z!*}Q$cszg-KOd)(OS|mXmUD~$1R#WKe_l*knGQUqto0tbXHP>$ov7av$aF-VBhbdw
zvU#WR0yRp7Bm9HQ+K~k4^uYDHjFR*ID|Vz(J6cadTG;(-?ZnM4Wdv`S%iDRvgUZJj
zZ4xMerYQ;Hm7uAMxp<$eTsdgUk#?@x@TcW&D4tly)p$(95O@1W>SK-@iS1XaUFJ_o
zdD{%b(qk5`eR#H9ac^27|M^>KTSK{`Uaar)PqP<zqB48zi-SL<Jgs7{oEB*|(8-yw
zw$B$szsdak4O=vxsD}mS0Fe)N?&KPOpwhuBo<AFXKDIp~@mF<x3vu3#BeY|ud0}AC
z>Q(HIZ^MoPl^lj<0GZBOLAN3E%_srd_g+a{jXdrY==_EYfYuWUEDu;YQ@Dr&o14=6
z4{gl>igLW#zX)iU#*4MYMn^a+jbR!$x(VN|*3tE}K{wBU71-<RFrlfgTK{nazJ2~h
zjqH`#yj#W`_(Pf~v29Bd*q=T5&DuCjsETPcl1<rqo%h3zx%7=kBogT?ky8|S3*Dt+
zCGKKJxS3$mzZ(@dNT;ec^~GoQPBUT}5_V_y-ls*G9W5U(pkeOO3M2uI6bhueo}itL
z?zFBy;nBPOde}^-Adh<^DNnpHz7>N=%n!!q6IBbY=<p9>=7X=|{X3`5@j%V&7cJM^
zc-2^8VbqA84b8eb<Rsr8{yDy*ERk}mVb5=HcO}fQPfVDUcV=v^q3SY`;N!PQ?p_sA
z6Pu0jqLVuDQUmeDHn~K_Tci(ny}12ASc$8qz4)V?F#SD0fH!i-2b?6u8mQ`>BJWfr
z3>k5g!wK#x`^XaacSnZ?!zc}c{INnTW;X4OB!Y8s*e6<TZ4R)3XP$z6N_O~(!BOdm
z8AfNnL)p+OqWX4u0@qb|)O+N%j|G7uN!pV#uurvFkl$j#fwg><VpY69-u7ulDJv(A
zw#UKx2Pt%9huGzsGj;zf8b-bTzQ*2qhkJcEf&gi>l^MHm1UBbl64<4E3dy+1UpKo+
zd#L+6px|Y767`PO)wpO0<aU$rV!E^KDb6Ol2K(9lOi+y@kr5ppA}~b8uhmV3nXhZ8
zUDy#yHJhn5>fDgLT>NCT>)L%{!JSy^bz`0+FHfc7cM5bQUgP;M`37sY@Tvpts`Ah_
zJ<Q&?V&`Skze|#dorIUgsWCIlt7?0|`fq28J1tCqnWQnw<jYIlzT1c6b01z|cg#x5
z<Qv?Z)sp4=y20Tx8*!@WJdj~(ePzTX&3X@djC?HKsAkAgc_5)Z|KR|C(r{}I4Rh{~
zAa7%Fg?vnvjmlzWw_fU#t@5_S!>dJX2El6P02tl2f?cVqMqK%hLgh2wL|ERiBra(p
zQ23UQs+#kP)<wTQt4(Kx-BhTf2db*qmd%kl5!^N5E`Daqmp@o1SW|W8h;+v4(rnpS
z8x3POwo3bDkawWVNAu;@o^W%z^%bj9=7R5XRzfpj@PUj=)Lv5?^GotIKXGU9Onj02
zIJI|~hpxtnFVaJE7E}9a=K@CmSgpYMV}o5QNMP!st}%2I(^RuKH|IukZKLh$!)9w-
zVDfwE1{*V@(a*c~ftqLTQSbKT5kbv@d&T@iyuvhP&)t@bU+XhUO7ryd$$Bphk^X=S
zYGi%X81naeM1m9y860?L+KL<QZ~}hCa!>e{cbwABw{MM*(03Z_(}Au&AG9%to%Xo+
zg{&xK!S;kntX=W(Z)rVy_GlYF*l_SxhwzF=-gyJ0d57wfTitXv_It*F&a>s59H2>+
zfzCFbdOSmcs>aN3#7Zqh;HxE_;6CfS;490T(P!gTx+rmh$<~QH@YRYhXMPzYvkyhl
zqb-YE@m%1XtI;I~HoUq$PEqg-Je#W3;=FIGz?5C)C*Kp{oeH>KoYFh=3;q6Chq~pK
z4o>5_w2-ExsF*qAl&{r9ejW{z<#N!naBg({1GxhdxU6e9J@l_)wEq3PL3H7}pP6yL
znCauW-j8sh&3$s1EcX~o;ds6ZXk$L!%3a9D>K2H#Cr53ux@QOW79_A{yyFjNl|uA>
zWz!i$q$9X3$e5YyVAE;ash`<<Ct8Wl3i5q1qJVLnD9PW2CwW;BD6|iYDRkE2eze@y
z4)jnDJJ+$WE!LgRsJUh-8|!r3hX<xDy95$Q+c#bhaf4r#aW}2Ow&)GJX<=KGe+9F#
zEzKF<KjJa|-M{y&2}d*jw0}c=t4yU?JVe*>z@&N}h33A@D+gFtb`VXq-m|U+Oe6k4
zt8NTFO~aI|KX0}aQ{tP9RtOa4i=~v&(z9E(Si~mWoZ9iDh!t>VgGR%Xtk0X@T(0eM
zwZ!z~^Vh+ai)d}0O8KBo6*fM8?AA7u+e<fjDyicvpN@v$;uBklY3HR^#d*^(nAX(w
z0vcAySD2hl28K(cw+JSeJ&{($)Tl979+w_NG0WaIozUe;%pKU*b5+PV^u$$yX?5>V
zHI5@$u6R<o1*()|=aDmUJh|dMXiH2LzAs{9Q+PkBH(^tbSI)h2$6hVOIJ(8|KV1s-
z|JsH@gG{W%OlNo~74SR$SAWiAwghXU*MtJY|JQq*md!h8sC^H;2mQAmy-+z7eg0R-
z-_(NCB!e<PDD&^y|J*nkujx<M|3_8+F85CiLmWrKh)x)?80Zo<rgV&!b}_c>3W@Jd
z97r4b5!L!FIVxThnx0@0mpC@nTc`$&D5iTQ9trs%3*TL1ePKF6p_?(&AG33Q8_woS
zB@7k}CxG?i^GCyrP{M5Q#7iyHJhDXSgZbi|@(sXXn$Vrc2mg6xE7p;|*mZXMK|aP^
z4!wagGrM~4gpjaW>$7!LpsThq%9Q*6F%u}d$8zV;hORsH(F?V=QU01-e>@le=fu`Z
yT=NRFX2UgjK~mQfxj?e}C4E~x)_<<HNWm=TdE=5IV)e3$Impd<t<!Z!i1%+sheBTf

literal 0
HcmV?d00001

diff --git a/deal.II/examples/step-39/step-39.cc b/deal.II/examples/step-39/step-39.cc
index 04dfa1bc48..2e6186cea3 100644
--- a/deal.II/examples/step-39/step-39.cc
+++ b/deal.II/examples/step-39/step-39.cc
@@ -711,9 +711,9 @@ Step39<dim>::run(unsigned int n_steps)
 
 int main()
 {
-  std::ofstream logfile("step-39.log");
+  std::ofstream logfile("deallog");
   deallog.attach(logfile);
   FE_DGQ<2> fe1(3);
   Step39<2> test1(fe1);
-  test1.run(20);
+  test1.run(12);
 }
-- 
2.39.5