[cig-commits] commit: Update benchmarks in manual

[Mercurial]

changeset:   473:10fa4d5529a8
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Aug 22 15:03:39 2012 -0700
files:       doc/gale.lyx doc/images/Sphere_Error.eps doc/images/Sphere_v.png doc/images/divergence_2D_error.eps doc/images/divergence_3D_error.eps doc/images/divergence_v.png doc/images/divergence_v_3D.png doc/images/sinusoid.png doc/images/sinusoid_error.eps doc/images/sinusoid_scaled_error.eps doc/images/thermal_diffusion_end.png doc/images/thermal_diffusion_error.eps doc/images/thermal_diffusion_start.png doc/images/thermal_eulear_high.png doc/images/thermal_eulear_low.png doc/images/thermal_eulear_medium.png doc/images/thermal_eulear_start.png doc/images/thermal_euler_high.png doc/images/thermal_euler_low.png doc/images/thermal_euler_medium.png doc/images/thermal_euler_start.png doc/images/thermal_lagrangian_end.png doc/images/thermal_lagrangian_start.png
description:
Update benchmarks in manual


diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/gale.lyx
--- a/doc/gale.lyx	Wed Aug 22 15:03:10 2012 -0700
+++ b/doc/gale.lyx	Wed Aug 22 15:03:39 2012 -0700
@@ -2290,14 +2290,6 @@ In practice, both of these quantities ma
  on the surface, such as landslides.
  A maximum strain rate, in the mean time, will assist in ensuring that the
  code is convergent.
- Both of these parameters are used for the Geomod benchmark (Sections 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
-).
 \end_layout
 
 \begin_layout Subsection
@@ -2311,25 +2303,11 @@ If you have a problem with strong viscos
  viscosities (e.g., Appendix 
 \begin_inset CommandInset ref
 LatexCommand ref
-reference "sec:Circular-Inclusion"
-
-\end_inset
-
- and 
-\begin_inset CommandInset ref
-LatexCommand ref
 reference "sec:Falling-Sphere"
 
 \end_inset
 
-), or when materials yield (e.g., Appendix 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
-).
+), or when materials yield.
 \end_layout
 
 \begin_layout Standard
@@ -3604,13 +3582,6 @@ For some problems where you are comparin
 \begin_inset CommandInset ref
 LatexCommand ref
 reference "sec:Falling-Sphere"
-
-\end_inset
-
-, 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Circular-Inclusion"
 
 \end_inset
 
@@ -17698,7 +17669,7 @@ status open
 \begin_layout Plain Layout
 \align center
 \begin_inset Tabular
-<lyxtabular version="3" rows="8" columns="2">
+<lyxtabular version="3" rows="7" columns="2">
 <features tabularvalignment="middle">
 <column alignment="center" valignment="top" width="0">
 <column alignment="center" valignment="top" width="0">
@@ -17738,13 +17709,6 @@ Stokes solver and interpolate between pa
 \begin_layout Plain Layout
 \begin_inset CommandInset ref
 LatexCommand ref
-reference "sec:Circular-Inclusion"
-
-\end_inset
-
-,
-\begin_inset CommandInset ref
-LatexCommand ref
 reference "sec:Relaxation-of-Topography"
 
 \end_inset
@@ -17756,13 +17720,6 @@ reference "sec:Divergence"
 
 \end_inset
 
- , 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
 
 \end_layout
 
@@ -17822,14 +17779,7 @@ reference "sec:Relaxation-of-Topography"
 
 \end_inset
 
-,  
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
- 
+
 \end_layout
 
 \end_inset
@@ -17862,13 +17812,6 @@ reference "sec:Relaxation-of-Topography"
 
 \end_inset
 
-, 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
 
 \end_layout
 
@@ -17892,39 +17835,6 @@ Free surface
 \begin_inset CommandInset ref
 LatexCommand ref
 reference "sec:Relaxation-of-Topography"
-
-\end_inset
-
-, 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-</cell>
-</row>
-<row>
-<cell alignment="center" valignment="top" topline="true" leftline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-Drucker Prager rheology in 2D
-\end_layout
-
-\end_inset
-</cell>
-<cell alignment="center" valignment="top" topline="true" leftline="true" rightline="true" usebox="none">
-\begin_inset Text
-
-\begin_layout Plain Layout
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "sec:Geomod-2004"
 
 \end_inset
 
@@ -18039,6 +17949,13 @@ Altogether, these benchmarks give us a h
  code.
 \end_layout
 
+\begin_layout Standard
+\begin_inset Newpage pagebreak
+\end_inset
+
+
+\end_layout
+
 \begin_layout Section
 \begin_inset CommandInset label
 LatexCommand label
@@ -18173,7 +18090,7 @@ Setting
 \end_inset
 
 , 
-\begin_inset Formula $\delta\rho=0.01$
+\begin_inset Formula $\delta\rho=1$
 \end_inset
 
 , and 
@@ -18186,7 +18103,7 @@ Setting
 \begin_layout Standard
 \begin_inset Formula 
 \[
-u=0.00222.
+u=0.222.
 \]
 
 \end_inset
@@ -18282,7 +18199,7 @@ For our case with
 \begin_layout Standard
 \begin_inset Formula 
 \[
-u=1.122747319\cdot10^{-3}.
+u=0.1122747319.
 \]
 
 \end_inset
@@ -18327,7 +18244,7 @@ For our case, this gives a speed of
 \begin_layout Standard
 \begin_inset Formula 
 \[
-u=1.12293603939\cdot10^{-3},
+u=0.112293603939,
 \]
 
 \end_inset
@@ -18460,343 +18377,6 @@ Error in computed velocity vs.
 
 \end_layout
 
-\begin_layout Standard
-\begin_inset Newpage pagebreak
-\end_inset
-
-
-\end_layout
-
-\begin_layout Section
-\begin_inset CommandInset label
-LatexCommand label
-name "sec:Circular-Inclusion"
-
-\end_inset
-
-Circular Inclusion
-\end_layout
-
-\begin_layout Standard
-Schmid and Podladchikov 
-\begin_inset CommandInset citation
-LatexCommand cite
-key "Clast"
-
-\end_inset
-
- derived a simple analytic solution for the pressure and velocity fields
- for a circular inclusion under simple shear as in Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:inclusion-setup"
-
-\end_inset
-
-.
-\end_layout
-
-\begin_layout Standard
-\noindent
-\align center
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/inclusion_setup.png
-	lyxscale 50
-	width 50col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:inclusion-setup"
-
-\end_inset
-
-Schematic for the circular inclusion benchmark
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-Because of the symmetry of the problem, we only have to solve over the top
- right quarter of the domain.
- For the velocity boundary conditions, the analytic solution is a bit complicate
-d.
- So we used the simple relation
-\end_layout
-
-\begin_layout Standard
-\begin_inset Formula 
-\begin{eqnarray*}
-v_{x} & = & -\dot{\epsilon}y,\\
-v_{y} & = & \dot{\epsilon}x,
-\end{eqnarray*}
-
-\end_inset
-
- for the boundaries, where 
-\begin_inset Formula $\dot{\epsilon}=1$
-\end_inset
-
- is the magnitude of the shear and 
-\begin_inset Formula $x$
-\end_inset
-
- and 
-\begin_inset Formula $y$
-\end_inset
-
- are the coordinates.
- This induces an error of order 
-\begin_inset Formula $r_{i}^{2}/r^{2}$
-\end_inset
-
-, where 
-\begin_inset Formula $r_{i}=1$
-\end_inset
-
- is the radius of the inclusion, and 
-\begin_inset Formula $r$
-\end_inset
-
- is the radius.
- We have the boundaries at 80 times the radius of the inclusion, giving
- an error of about 0.01%.
- This error dominates near the boundary, but near the inclusion other errors
- dominate.
- 
-\end_layout
-
-\begin_layout Standard
-A characteristic of the analytic solution is that the pressure is zero inside
- the inclusion, while outside it follows the relation
-\end_layout
-
-\begin_layout Standard
-\begin_inset Formula 
-\[
-p_{m}=4\dot{\epsilon}\frac{\mu_{m}\left(\mu_{i}-\mu_{m}\right)}{\mu_{i}+\mu_{m}}\frac{r_{i}^{2}}{r^{2}}\cos\left(2\theta\right),
-\]
-
-\end_inset
-
-where 
-\begin_inset Formula $\mu_{i}=2$
-\end_inset
-
- is the viscosity of the inclusion and 
-\begin_inset Formula $\mu_{m}=1$
-\end_inset
-
- is the viscosity of the background media.
- Many numerical codes that solve Stokes flow (Eq.
- 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "eq:simple momentum conservation"
-
-\end_inset
-
- and 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "eq:continuity"
-
-\end_inset
-
-), including Gale, assume that pressure, velocity, and viscosity are continuous
- within an element.
- The pressure discontinuity at the surface of the inclusion violates that
- assumption, so the error tends to concentrate near the surface of the inclusion.
-\end_layout
-
-\begin_layout Standard
-Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Pressure-inclusion"
-
-\end_inset
-
- plots the error in the pressure along the line 
-\begin_inset Formula $y=x/2$
-\end_inset
-
- for different resolutions.
- The convergence is not uniform.
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/inclusion_pressure.png
-	width 100col%
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Pressure-inclusion"
-
-\end_inset
-
-Pressure along the line 
-\begin_inset Formula $y=x/2$
-\end_inset
-
- for resolutions of 128
-\begin_inset Formula $\times$
-\end_inset
-
-128 (red), 256
-\begin_inset Formula $\times$
-\end_inset
-
-256 (green), and 512
-\begin_inset Formula $\times$
-\end_inset
-
-512 (blue).
- The inclusion has a radius 
-\begin_inset Formula $r_{i}=1.$
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-Outside the inclusion, the error is better behaved.
- Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Pressure-Error"
-
-\end_inset
-
- plots the error in the pressure along the line 
-\begin_inset Formula $y=x/2$
-\end_inset
-
- outside the inclusion for different resolutions.
- This kind of result, where the solution is bad close to the surface, but
- good otherwise, is typical for numerical solutions of this problem 
-\begin_inset CommandInset citation
-LatexCommand cite
-key "FD Stokes"
-
-\end_inset
-
-.
- 
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/inclusion_error.png
-	lyxscale 75
-	width 100col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Pressure-Error"
-
-\end_inset
-
-Error in the pressure outside the inclusion along the line 
-\begin_inset Formula $y=x/2$
-\end_inset
-
- for resolutions of 128
-\begin_inset Formula $\times$
-\end_inset
-
-128 (blue), 256
-\begin_inset Formula $\times$
-\end_inset
-
-256 (red), and 512
-\begin_inset Formula $\times$
-\end_inset
-
-512 (black).
- The inclusion has a radius 
-\begin_inset Formula $r_{i}=1.$
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
 \begin_layout Section
 \begin_inset CommandInset label
 LatexCommand label
@@ -18876,11 +18456,11 @@ reference "fig:Strain-topo"
 
 \end_inset
 
- shows the results of a high-resolution (128
+ shows the results of a high-resolution (256
 \begin_inset Formula $\times$
 \end_inset
 
-256) run.
+512) run.
  Note that we use symmetry to only simulate half of the wavelength.
  
 \end_layout
@@ -18936,16 +18516,17 @@ reference "fig:topo-error"
 \end_inset
 
 .
- Scaling the error with resolution gives Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:scaled-topo-error"
-
-\end_inset
-
-.
- The error decreases linearly with increasing resolution, giving us confidence
- in our ability to accurately track topography.
+ The error behaves a bit erratically because of the damping term applied
+ to the free surface (Section 
+\begin_inset CommandInset ref
+LatexCommand ref
+reference "sub:BouyancyForceTerm"
+
+\end_inset
+
+).
+ Even so, the error decreases linearly with increasing resolution, giving
+ us confidence in our ability to accurately track topography.
 \end_layout
 
 \begin_layout Standard
@@ -18975,58 +18556,6 @@ name "fig:topo-error"
 \end_inset
 
 Error in the height at the peak
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\noindent
-\align center
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/sinusoid_scaled_error.eps
-
-\end_inset
-
-
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:scaled-topo-error"
-
-\end_inset
-
-As in Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:topo-error"
-
-\end_inset
-
-, but with the error scaled with 
-\begin_inset Formula $h$
-\end_inset
-
-.
- So the medium-resolution error is multiplied by 2 and the high-resolution
- error is multiplied by 4.
 \end_layout
 
 \end_inset
@@ -19061,22 +18590,6 @@ reference "eq:divergence"
 .
  In 2D, a constant divergence is applied to a square domain, and the velocity
  on the corners is set to enforce a spreading from the center of the square.
- Figures 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Divergence_v_sri"
-
-\end_inset
-
- and 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Divergence_v_sri3D"
-
-\end_inset
-
- show the velocity and the error in the strain rate invariant for a numerical
- solution.
  For a constant divergence 
 \begin_inset Formula $d$
 \end_inset
@@ -19129,87 +18642,6 @@ linearTolerance
 \end_layout
 
 \begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/divergence_v.png
-	width 50col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Divergence_v_sri"
-
-\end_inset
-
-Velocity and error in the strain rate invariant solution for the 2D Divergence
- benchmark.
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/divergence_v_3D.png
-	width 50col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Divergence_v_sri3D"
-
-\end_inset
-
-Velocity and error in the strain rate invariant solution for the 3D Divergence
- benchmark.
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
 \begin_inset Newpage newpage
 \end_inset
 
@@ -19281,7 +18713,7 @@ reference "fig:Diffusion_end"
 \end_inset
 
  and 
-\begin_inset Formula $t=0.0367647$
+\begin_inset Formula $t=0.0011489$
 \end_inset
 
  for a run with 16
@@ -19297,7 +18729,7 @@ reference "fig:Diffusion_error"
 \end_inset
 
  plots the error in the maximum temperature at 
-\begin_inset Formula $t=0.0367647$
+\begin_inset Formula $t=0.0011489$
 \end_inset
 
  as a function of the grid spacing 
@@ -19305,7 +18737,7 @@ reference "fig:Diffusion_error"
 \end_inset
 
 .
- The error decreases quadratically as the spacing decreases.
+ The error decreases linearly as the spacing decreases.
 \end_layout
 
 \begin_layout Standard
@@ -19426,7 +18858,7 @@ name "fig:Diffusion_error"
 \end_inset
 
 Error in the maximum temperature at 
-\begin_inset Formula $t=0.0367647$
+\begin_inset Formula $t=0.0011489$
 \end_inset
 
  as a function of resolution.
@@ -19674,7 +19106,7 @@ status open
 \begin_layout Plain Layout
 \align center
 \begin_inset Graphics
-	filename images/thermal_eulear_start.png
+	filename images/thermal_euler_start.png
 	width 50col%
 
 \end_inset
@@ -19715,7 +19147,7 @@ status open
 \begin_layout Plain Layout
 \align center
 \begin_inset Graphics
-	filename images/thermal_eulear_low.png
+	filename images/thermal_euler_low.png
 	width 50col%
 
 \end_inset
@@ -19761,7 +19193,7 @@ status open
 \begin_layout Plain Layout
 \align center
 \begin_inset Graphics
-	filename images/thermal_eulear_medium.png
+	filename images/thermal_euler_medium.png
 	width 50col%
 
 \end_inset
@@ -19807,7 +19239,7 @@ status open
 \begin_layout Plain Layout
 \align center
 \begin_inset Graphics
-	filename images/thermal_eulear_high.png
+	filename images/thermal_euler_high.png
 	width 50col%
 
 \end_inset
@@ -19834,196 +19266,6 @@ Temperature at
 \end_inset
 
 64 elements.
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Newpage pagebreak
-\end_inset
-
-
-\end_layout
-
-\begin_layout Section
-\begin_inset CommandInset label
-LatexCommand label
-name "sec:Geomod-2004"
-
-\end_inset
-
-Geomod 2004 Extension
-\end_layout
-
-\begin_layout Standard
-This benchmark simulates a sandbox being extended as in Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Extension-model-setup"
-
-\end_inset
-
-.
- The right side and half of the bottom are translated to the right.
- This creates a velocity discontinuity at the center which is the initial
- seed for localization.
-\end_layout
-
-\begin_layout Standard
-Like half of the codes in the benchmark, boundary friction was not included.
- Rather, the material is held fixed to the bottom boundary, and the velocity
- discontinuity is smoothed over 0.2 cm.
- In addition, the exact background viscosity is not prescribed by the benchmark.
- We have used 
-\begin_inset Formula $10^{12}Pa\cdot s$
-\end_inset
-
-, the same as used in the I2ELVIS calculations.
- This value is somewhere in the middle of the range of values used in the
- calculations for other codes.
- We used a maximum strain rate of 
-\begin_inset Formula $5\cdot10^{4}$
-\end_inset
-
- and a minimum viscosity of 
-\begin_inset Formula $10^{-4}$
-\end_inset
-
-.
- Figure 
-\begin_inset CommandInset ref
-LatexCommand ref
-reference "fig:Comparison-extension"
-
-\end_inset
-
- shows the Gale results along with the other codes.
- While it is difficult to perform exact comparisons, Gale produces similar
- fault patterns.
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/Extension_setup.png
-	lyxscale 50
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Extension-model-setup"
-
-\end_inset
-
-Extension model setup.
- Reproduced, with permission, from Buiter et al.
- 
-\begin_inset CommandInset citation
-LatexCommand cite
-key "Buiter et al Numerical Sandbox"
-
-\end_inset
-
-.
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Standard
-\begin_inset Float figure
-placement H
-wide false
-sideways false
-status open
-
-\begin_layout Plain Layout
-\align center
-\begin_inset Graphics
-	filename images/Extension_comparision.png
-	lyxscale 50
-	width 100col%
-
-\end_inset
-
-
-\end_layout
-
-\begin_layout Plain Layout
-\begin_inset Caption
-
-\begin_layout Plain Layout
-\begin_inset CommandInset label
-LatexCommand label
-name "fig:Comparison-extension"
-
-\end_inset
-
-Strain rate invariant for the numerical extension models after 5 cm of extension.
- The resolutions of the various models are: I2ELVIS: 400
-\begin_inset Formula $\times$
-\end_inset
-
-75, LAPEX-2D: 301
-\begin_inset Formula $\times$
-\end_inset
-
-71, Microfem: 201
-\begin_inset Formula $\times$
-\end_inset
-
-61, SloMo: 401
-\begin_inset Formula $\times$
-\end_inset
-
-71, Sopale: 401
-\begin_inset Formula $\times$
-\end_inset
-
-71, Gale: 1024
-\begin_inset Formula $\times$
-\end_inset
-
-128.
- Upper images reproduced, with permission, from Buiter et al.
- 
-\begin_inset CommandInset citation
-LatexCommand cite
-key "Buiter et al Numerical Sandbox"
-
-\end_inset
-
-.
 \end_layout
 
 \end_inset
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/Sphere_Error.eps
--- a/doc/images/Sphere_Error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ b/doc/images/Sphere_Error.eps	Wed Aug 22 15:03:39 2012 -0700
@@ -1,7 +1,7 @@
 %!PS-Adobe-2.0 EPSF-2.0
 %%Title: Sphere_Error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Wed Nov 23 17:22:15 2011
+%%Creator: gnuplot 4.6 patchlevel 0
+%%CreationDate: Mon Aug 20 16:08:03 2012
 %%DocumentFonts: (atend)
 %%BoundingBox: 50 50 410 302
 %%EndComments
@@ -20,10 +20,12 @@ gnudict begin
 /Level1 false def
 /Rounded false def
 /ClipToBoundingBox false def
+/SuppressPDFMark false def
 /TransparentPatterns false def
 /gnulinewidth 5.000 def
 /userlinewidth gnulinewidth def
 /Gamma 1.0 def
+/BackgroundColor {-1.000 -1.000 -1.000} def
 %
 /vshift -46 def
 /dl1 {
@@ -38,22 +40,6 @@ gnudict begin
 /vpt_ 31.5 def
 /hpt hpt_ def
 /vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (Sphere_Error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Wed Nov 23 17:22:15 2011)
-  /DOCINFO pdfmark
-end
-} ifelse
 /doclip {
   ClipToBoundingBox {
     newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
@@ -61,7 +47,7 @@ end
   } if
 } def
 %
-% Gnuplot Prolog Version 4.4 (January 2010)
+% Gnuplot Prolog Version 4.4 (August 2010)
 %
 %/SuppressPDFMark true def
 %
@@ -73,6 +59,7 @@ end
 /Z {closepath} bind def
 /C {setrgbcolor} bind def
 /f {rlineto fill} bind def
+/g {setgray} bind def
 /Gshow {show} def   % May be redefined later in the file to support UTF-8
 /vpt2 vpt 2 mul def
 /hpt2 hpt 2 mul def
@@ -94,6 +81,7 @@ end
 	dup 1 lt {pop 1} if 10 mul /udl exch def} def
 /PL {stroke userlinewidth setlinewidth
 	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
 % Default Line colors
 /LCw {1 1 1} def
 /LCb {0 0 0} def
@@ -431,8 +419,26 @@ Level1 {Level1PatternFill} {Level2Patter
 /Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
 dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
 currentdict end definefont pop
+Level1 SuppressPDFMark or 
+{} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (Sphere_Error.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 4.6 patchlevel 0)
+  /Author (boo)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Mon Aug 20 16:08:03 2012)
+  /DOCINFO pdfmark
+end
+} ifelse
 end
 %%EndProlog
+%%Page: 1 1
 gnudict begin
 gsave
 doclip
@@ -441,113 +447,98 @@ 0 setgray
 0 setgray
 newpath
 (Helvetica) findfont 140 scalefont setfont
+BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 7200.00 5040.00 BoxColFill} if
 1.000 UL
 LTb
-938 448 M
+854 448 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
-854 448 M
-( 0.002) Rshow
+770 448 M
+( 0) Rshow
 1.000 UL
 LTb
-938 850 M
+854 939 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
-854 850 M
-( 0.003) Rshow
-1.000 UL
-LTb
-938 1252 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
-( 0.004) Rshow
-1.000 UL
-LTb
-938 1654 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
+770 939 M
 ( 0.005) Rshow
 1.000 UL
 LTb
-938 2056 M
+854 1431 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.006) Rshow
-1.000 UL
-LTb
-938 2458 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
-( 0.007) Rshow
-1.000 UL
-LTb
-938 2861 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
-( 0.008) Rshow
-1.000 UL
-LTb
-938 3263 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
-( 0.009) Rshow
-1.000 UL
-LTb
-938 3665 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
+-6114 0 R
 ( 0.01) Rshow
 1.000 UL
 LTb
-938 4067 M
+854 1922 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.011) Rshow
+-6114 0 R
+( 0.015) Rshow
 1.000 UL
 LTb
-938 4469 M
+854 2414 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.012) Rshow
+-6114 0 R
+( 0.02) Rshow
 1.000 UL
 LTb
-938 4871 M
+854 2905 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.013) Rshow
+-6114 0 R
+( 0.025) Rshow
 1.000 UL
 LTb
-938 448 M
+854 3397 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.03) Rshow
+1.000 UL
+LTb
+854 3888 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.035) Rshow
+1.000 UL
+LTb
+854 4380 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.04) Rshow
+1.000 UL
+LTb
+854 4871 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.045) Rshow
+1.000 UL
+LTb
+854 448 M
 0 63 V
 0 4360 R
 0 -63 V
-938 308 M
+854 308 M
 ( 0) Cshow
 1.000 UL
 LTb
-1947 448 M
+1870 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -555,7 +546,7 @@ 0 -4500 R
 ( 0.1) Cshow
 1.000 UL
 LTb
-2955 448 M
+2885 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -563,7 +554,7 @@ 0 -4500 R
 ( 0.2) Cshow
 1.000 UL
 LTb
-3964 448 M
+3901 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -571,7 +562,7 @@ 0 -4500 R
 ( 0.3) Cshow
 1.000 UL
 LTb
-4972 448 M
+4916 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -579,7 +570,7 @@ 0 -4500 R
 ( 0.4) Cshow
 1.000 UL
 LTb
-5981 448 M
+5932 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -587,7 +578,7 @@ 0 -4500 R
 ( 0.5) Cshow
 1.000 UL
 LTb
-6989 448 M
+6947 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -597,20 +588,20 @@ LTb
 LTb
 1.000 UL
 LTb
-938 4871 N
-938 448 L
-6051 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--6051 0 V
+-6093 0 V
 Z stroke
 LCb setrgbcolor
-196 2659 M
+112 2659 M
 currentpoint gsave translate -270 rotate 0 0 M
 (Velocity Error) Cshow
 grestore
 LTb
 LCb setrgbcolor
-3963 98 M
+3900 98 M
 (h) Cshow
 LTb
 1.000 UP
@@ -620,35 +611,35 @@ 1.000 UP
 1.000 UP
 4.000 UL
 LT0
-5981 4163 M
-0 630 V
--31 -630 R
+5932 4037 M
+0 500 V
+-31 -500 R
 62 0 V
--62 630 R
+-62 500 R
 62 0 V
-3459 2003 M
-0 567 V
--31 -567 R
+3393 2261 M
+0 131 V
+-31 -131 R
 62 0 V
--62 567 R
+-62 131 R
 62 0 V
-2199 807 M
-0 505 V
-2168 807 M
+2123 742 M
+0 128 V
+2092 742 M
 62 0 V
--62 505 R
+-62 128 R
 62 0 V
-5981 4478 Pls
-3459 2287 Pls
-2199 1059 Pls
+5932 4287 Pls
+3393 2326 Pls
+2123 806 Pls
 % End plot #1
 1.000 UL
 LTb
-938 4871 N
-938 448 L
-6051 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--6051 0 V
+-6093 0 V
 Z stroke
 1.000 UP
 1.000 UL
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/Sphere_v.png
Binary file doc/images/Sphere_v.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/divergence_2D_error.eps
--- a/doc/images/divergence_2D_error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,768 +0,0 @@
-%!PS-Adobe-2.0 EPSF-2.0
-%%Title: divergence_2D_error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Thu Jul 22 15:35:46 2010
-%%DocumentFonts: (atend)
-%%BoundingBox: 50 50 410 302
-%%EndComments
-%%BeginProlog
-/gnudict 256 dict def
-gnudict begin
-%
-% The following true/false flags may be edited by hand if desired.
-% The unit line width and grayscale image gamma correction may also be changed.
-%
-/Color true def
-/Blacktext false def
-/Solid false def
-/Dashlength 1 def
-/Landscape false def
-/Level1 false def
-/Rounded false def
-/ClipToBoundingBox false def
-/TransparentPatterns false def
-/gnulinewidth 5.000 def
-/userlinewidth gnulinewidth def
-/Gamma 1.0 def
-%
-/vshift -46 def
-/dl1 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
-} def
-/dl2 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul add } if
-} def
-/hpt_ 31.5 def
-/vpt_ 31.5 def
-/hpt hpt_ def
-/vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (divergence_2D_error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Thu Jul 22 15:35:46 2010)
-  /DOCINFO pdfmark
-end
-} ifelse
-/doclip {
-  ClipToBoundingBox {
-    newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
-    clip
-  } if
-} def
-%
-% Gnuplot Prolog Version 4.4 (January 2010)
-%
-%/SuppressPDFMark true def
-%
-/M {moveto} bind def
-/L {lineto} bind def
-/R {rmoveto} bind def
-/V {rlineto} bind def
-/N {newpath moveto} bind def
-/Z {closepath} bind def
-/C {setrgbcolor} bind def
-/f {rlineto fill} bind def
-/Gshow {show} def   % May be redefined later in the file to support UTF-8
-/vpt2 vpt 2 mul def
-/hpt2 hpt 2 mul def
-/Lshow {currentpoint stroke M 0 vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
-  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
-/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
- {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
-/BL {stroke userlinewidth 2 mul setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/AL {stroke userlinewidth 2 div setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/UL {dup gnulinewidth mul /userlinewidth exch def
-	dup 1 lt {pop 1} if 10 mul /udl exch def} def
-/PL {stroke userlinewidth setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-% Default Line colors
-/LCw {1 1 1} def
-/LCb {0 0 0} def
-/LCa {0 0 0} def
-/LC0 {1 0 0} def
-/LC1 {0 1 0} def
-/LC2 {0 0 1} def
-/LC3 {1 0 1} def
-/LC4 {0 1 1} def
-/LC5 {1 1 0} def
-/LC6 {0 0 0} def
-/LC7 {1 0.3 0} def
-/LC8 {0.5 0.5 0.5} def
-% Default Line Types
-/LTw {PL [] 1 setgray} def
-/LTb {BL [] LCb DL} def
-/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
-/LT0 {PL [] LC0 DL} def
-/LT1 {PL [4 dl1 2 dl2] LC1 DL} def
-/LT2 {PL [2 dl1 3 dl2] LC2 DL} def
-/LT3 {PL [1 dl1 1.5 dl2] LC3 DL} def
-/LT4 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
-/LT5 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC5 DL} def
-/LT6 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC6 DL} def
-/LT7 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC7 DL} def
-/LT8 {PL [2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 4 dl2] LC8 DL} def
-/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
-/Dia {stroke [] 0 setdash 2 copy vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke
-  Pnt} def
-/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
-  currentpoint stroke M
-  hpt neg vpt neg R hpt2 0 V stroke
- } def
-/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke
-  Pnt} def
-/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
-  hpt2 vpt2 neg V currentpoint stroke M
-  hpt2 neg 0 R hpt2 vpt2 V stroke} def
-/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke
-  Pnt} def
-/Star {2 copy Pls Crs} def
-/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath fill} def
-/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath fill} def
-/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke
-  Pnt} def
-/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath fill} def
-/DiaF {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath fill} def
-/Pent {stroke [] 0 setdash 2 copy gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore Pnt} def
-/PentF {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath fill grestore} def
-/Circle {stroke [] 0 setdash 2 copy
-  hpt 0 360 arc stroke Pnt} def
-/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
-/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
-/C1 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C2 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C3 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C4 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C5 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc
-	2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc} bind def
-/C6 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C7 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C8 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C9 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 450 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
-	2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C11 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C12 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C13 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C14 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 360 arc closepath fill
-	vpt 0 360 arc} bind def
-/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
-	neg 0 rlineto closepath} bind def
-/Square {dup Rec} bind def
-/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
-/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
-/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
-/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
-/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
-	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
-/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
-	Bsquare} bind def
-/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
-	Bsquare} bind def
-/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
-/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
-/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
-/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
-/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
-/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
-/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
-/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
-/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
-/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
-/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
-/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
-/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
-/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
-/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
-/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
-/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
-/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
-/DiaE {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke} def
-/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke} def
-/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke} def
-/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke} def
-/PentE {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore} def
-/CircE {stroke [] 0 setdash 
-  hpt 0 360 arc stroke} def
-/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
-/DiaW {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V Opaque stroke} def
-/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V Opaque stroke} def
-/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V Opaque stroke} def
-/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V Opaque stroke} def
-/PentW {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  Opaque stroke grestore} def
-/CircW {stroke [] 0 setdash 
-  hpt 0 360 arc Opaque stroke} def
-/BoxFill {gsave Rec 1 setgray fill grestore} def
-/Density {
-  /Fillden exch def
-  currentrgbcolor
-  /ColB exch def /ColG exch def /ColR exch def
-  /ColR ColR Fillden mul Fillden sub 1 add def
-  /ColG ColG Fillden mul Fillden sub 1 add def
-  /ColB ColB Fillden mul Fillden sub 1 add def
-  ColR ColG ColB setrgbcolor} def
-/BoxColFill {gsave Rec PolyFill} def
-/PolyFill {gsave Density fill grestore grestore} def
-/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
-%
-% PostScript Level 1 Pattern Fill routine for rectangles
-% Usage: x y w h s a XX PatternFill
-%	x,y = lower left corner of box to be filled
-%	w,h = width and height of box
-%	  a = angle in degrees between lines and x-axis
-%	 XX = 0/1 for no/yes cross-hatch
-%
-/PatternFill {gsave /PFa [ 9 2 roll ] def
-  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
-  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
-  gsave 1 setgray fill grestore clip
-  currentlinewidth 0.5 mul setlinewidth
-  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
-  0 0 M PFa 5 get rotate PFs -2 div dup translate
-  0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 M 0 PFs V} for
-  0 PFa 6 get ne {
-	0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
- } if
-  stroke grestore} def
-%
-/languagelevel where
- {pop languagelevel} {1} ifelse
- 2 lt
-	{/InterpretLevel1 true def}
-	{/InterpretLevel1 Level1 def}
- ifelse
-%
-% PostScript level 2 pattern fill definitions
-%
-/Level2PatternFill {
-/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
-	bind def
-/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
->> matrix makepattern
-/Pat1 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
-	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
->> matrix makepattern
-/Pat2 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
-	8 8 L 8 0 L 0 0 L fill}
->> matrix makepattern
-/Pat3 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
-	0 12 M 12 0 L stroke}
->> matrix makepattern
-/Pat4 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
-	0 -4 M 12 8 L stroke}
->> matrix makepattern
-/Pat5 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
-	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
->> matrix makepattern
-/Pat6 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
-	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
->> matrix makepattern
-/Pat7 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
-	12 0 M -4 8 L 12 4 M 0 10 L stroke}
->> matrix makepattern
-/Pat8 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
-	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
->> matrix makepattern
-/Pat9 exch def
-/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
-/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
-/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
-/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
-/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
-/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
-/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
-} def
-%
-%
-%End of PostScript Level 2 code
-%
-/PatternBgnd {
-  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
-} def
-%
-% Substitute for Level 2 pattern fill codes with
-% grayscale if Level 2 support is not selected.
-%
-/Level1PatternFill {
-/Pattern1 {0.250 Density} bind def
-/Pattern2 {0.500 Density} bind def
-/Pattern3 {0.750 Density} bind def
-/Pattern4 {0.125 Density} bind def
-/Pattern5 {0.375 Density} bind def
-/Pattern6 {0.625 Density} bind def
-/Pattern7 {0.875 Density} bind def
-} def
-%
-% Now test for support of Level 2 code
-%
-Level1 {Level1PatternFill} {Level2PatternFill} ifelse
-%
-/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
-dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
-currentdict end definefont pop
-end
-%%EndProlog
-gnudict begin
-gsave
-doclip
-50 50 translate
-0.050 0.050 scale
-0 setgray
-newpath
-(Helvetica) findfont 140 scalefont setfont
-1.000 UL
-LTb
-1022 448 M
-63 0 V
-5904 0 R
--63 0 V
-938 448 M
-( 1e-09) Rshow
-1.000 UL
-LTb
-1022 714 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1066 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1247 M
-31 0 V
-5936 0 R
--31 0 V
--5936 86 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 1e-08) Rshow
-1.000 UL
-LTb
-1022 1599 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1951 M
-31 0 V
-5936 0 R
--31 0 V
-1022 2131 M
-31 0 V
-5936 0 R
--31 0 V
--5936 86 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 1e-07) Rshow
-1.000 UL
-LTb
-1022 2483 M
-31 0 V
-5936 0 R
--31 0 V
-1022 2836 M
-31 0 V
-5936 0 R
--31 0 V
-1022 3016 M
-31 0 V
-5936 0 R
--31 0 V
--5936 86 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 1e-06) Rshow
-1.000 UL
-LTb
-1022 3368 M
-31 0 V
-5936 0 R
--31 0 V
-1022 3720 M
-31 0 V
-5936 0 R
--31 0 V
-1022 3901 M
-31 0 V
-5936 0 R
--31 0 V
--5936 85 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 1e-05) Rshow
-1.000 UL
-LTb
-1022 4253 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4605 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4785 M
-31 0 V
-5936 0 R
--31 0 V
--5936 86 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0001) Rshow
-1.000 UL
-LTb
-1022 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 1e-09) Cshow
-1.000 UL
-LTb
-1471 448 M
-0 31 V
-0 4392 R
-0 -31 V
-1734 448 M
-0 31 V
-0 4392 R
-0 -31 V
-1920 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2065 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2183 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2283 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2369 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2445 448 M
-0 31 V
-0 4392 R
-0 -31 V
-2514 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 1e-08) Cshow
-1.000 UL
-LTb
-2963 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3225 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3412 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3556 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3675 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3774 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3861 448 M
-0 31 V
-0 4392 R
-0 -31 V
-3937 448 M
-0 31 V
-0 4392 R
-0 -31 V
-4006 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 1e-07) Cshow
-1.000 UL
-LTb
-4455 448 M
-0 31 V
-0 4392 R
-0 -31 V
-4717 448 M
-0 31 V
-0 4392 R
-0 -31 V
-4904 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5048 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5166 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5266 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5353 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5429 448 M
-0 31 V
-0 4392 R
-0 -31 V
-5497 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 1e-06) Cshow
-1.000 UL
-LTb
-5946 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6209 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6395 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6540 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6658 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6758 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6844 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6921 448 M
-0 31 V
-0 4392 R
-0 -31 V
-6989 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 1e-05) Cshow
-1.000 UL
-LTb
-1.000 UL
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-LCb setrgbcolor
-196 2659 M
-currentpoint gsave translate -270 rotate 0 0 M
-(Maximum StrainRate Error) Cshow
-grestore
-LTb
-LCb setrgbcolor
-4005 98 M
-(Solver Tolerance) Cshow
-LTb
-1.000 UP
-1.000 UL
-LTb
-% Begin plot #1
-1.000 UL
-LT0
-6989 4120 M
-5497 3380 L
-4006 2603 L
-2514 1327 L
-1022 543 L
-% End plot #1
-stroke
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-1.000 UP
-1.000 UL
-LTb
-stroke
-grestore
-end
-showpage
-%%Trailer
-%%DocumentFonts: Helvetica
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/divergence_3D_error.eps
--- a/doc/images/divergence_3D_error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,672 +0,0 @@
-%!PS-Adobe-2.0 EPSF-2.0
-%%Title: divergence_3D_error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Thu Jul 22 15:41:45 2010
-%%DocumentFonts: (atend)
-%%BoundingBox: 50 50 410 302
-%%EndComments
-%%BeginProlog
-/gnudict 256 dict def
-gnudict begin
-%
-% The following true/false flags may be edited by hand if desired.
-% The unit line width and grayscale image gamma correction may also be changed.
-%
-/Color true def
-/Blacktext false def
-/Solid false def
-/Dashlength 1 def
-/Landscape false def
-/Level1 false def
-/Rounded false def
-/ClipToBoundingBox false def
-/TransparentPatterns false def
-/gnulinewidth 5.000 def
-/userlinewidth gnulinewidth def
-/Gamma 1.0 def
-%
-/vshift -46 def
-/dl1 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
-} def
-/dl2 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul add } if
-} def
-/hpt_ 31.5 def
-/vpt_ 31.5 def
-/hpt hpt_ def
-/vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (divergence_3D_error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Thu Jul 22 15:41:45 2010)
-  /DOCINFO pdfmark
-end
-} ifelse
-/doclip {
-  ClipToBoundingBox {
-    newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
-    clip
-  } if
-} def
-%
-% Gnuplot Prolog Version 4.4 (January 2010)
-%
-%/SuppressPDFMark true def
-%
-/M {moveto} bind def
-/L {lineto} bind def
-/R {rmoveto} bind def
-/V {rlineto} bind def
-/N {newpath moveto} bind def
-/Z {closepath} bind def
-/C {setrgbcolor} bind def
-/f {rlineto fill} bind def
-/Gshow {show} def   % May be redefined later in the file to support UTF-8
-/vpt2 vpt 2 mul def
-/hpt2 hpt 2 mul def
-/Lshow {currentpoint stroke M 0 vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
-  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
-/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
- {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
-/BL {stroke userlinewidth 2 mul setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/AL {stroke userlinewidth 2 div setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/UL {dup gnulinewidth mul /userlinewidth exch def
-	dup 1 lt {pop 1} if 10 mul /udl exch def} def
-/PL {stroke userlinewidth setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-% Default Line colors
-/LCw {1 1 1} def
-/LCb {0 0 0} def
-/LCa {0 0 0} def
-/LC0 {1 0 0} def
-/LC1 {0 1 0} def
-/LC2 {0 0 1} def
-/LC3 {1 0 1} def
-/LC4 {0 1 1} def
-/LC5 {1 1 0} def
-/LC6 {0 0 0} def
-/LC7 {1 0.3 0} def
-/LC8 {0.5 0.5 0.5} def
-% Default Line Types
-/LTw {PL [] 1 setgray} def
-/LTb {BL [] LCb DL} def
-/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
-/LT0 {PL [] LC0 DL} def
-/LT1 {PL [4 dl1 2 dl2] LC1 DL} def
-/LT2 {PL [2 dl1 3 dl2] LC2 DL} def
-/LT3 {PL [1 dl1 1.5 dl2] LC3 DL} def
-/LT4 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
-/LT5 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC5 DL} def
-/LT6 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC6 DL} def
-/LT7 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC7 DL} def
-/LT8 {PL [2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 4 dl2] LC8 DL} def
-/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
-/Dia {stroke [] 0 setdash 2 copy vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke
-  Pnt} def
-/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
-  currentpoint stroke M
-  hpt neg vpt neg R hpt2 0 V stroke
- } def
-/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke
-  Pnt} def
-/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
-  hpt2 vpt2 neg V currentpoint stroke M
-  hpt2 neg 0 R hpt2 vpt2 V stroke} def
-/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke
-  Pnt} def
-/Star {2 copy Pls Crs} def
-/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath fill} def
-/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath fill} def
-/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke
-  Pnt} def
-/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath fill} def
-/DiaF {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath fill} def
-/Pent {stroke [] 0 setdash 2 copy gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore Pnt} def
-/PentF {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath fill grestore} def
-/Circle {stroke [] 0 setdash 2 copy
-  hpt 0 360 arc stroke Pnt} def
-/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
-/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
-/C1 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C2 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C3 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C4 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C5 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc
-	2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc} bind def
-/C6 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C7 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C8 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C9 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 450 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
-	2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C11 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C12 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C13 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C14 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 360 arc closepath fill
-	vpt 0 360 arc} bind def
-/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
-	neg 0 rlineto closepath} bind def
-/Square {dup Rec} bind def
-/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
-/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
-/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
-/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
-/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
-	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
-/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
-	Bsquare} bind def
-/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
-	Bsquare} bind def
-/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
-/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
-/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
-/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
-/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
-/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
-/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
-/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
-/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
-/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
-/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
-/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
-/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
-/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
-/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
-/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
-/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
-/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
-/DiaE {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke} def
-/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke} def
-/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke} def
-/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke} def
-/PentE {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore} def
-/CircE {stroke [] 0 setdash 
-  hpt 0 360 arc stroke} def
-/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
-/DiaW {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V Opaque stroke} def
-/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V Opaque stroke} def
-/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V Opaque stroke} def
-/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V Opaque stroke} def
-/PentW {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  Opaque stroke grestore} def
-/CircW {stroke [] 0 setdash 
-  hpt 0 360 arc Opaque stroke} def
-/BoxFill {gsave Rec 1 setgray fill grestore} def
-/Density {
-  /Fillden exch def
-  currentrgbcolor
-  /ColB exch def /ColG exch def /ColR exch def
-  /ColR ColR Fillden mul Fillden sub 1 add def
-  /ColG ColG Fillden mul Fillden sub 1 add def
-  /ColB ColB Fillden mul Fillden sub 1 add def
-  ColR ColG ColB setrgbcolor} def
-/BoxColFill {gsave Rec PolyFill} def
-/PolyFill {gsave Density fill grestore grestore} def
-/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
-%
-% PostScript Level 1 Pattern Fill routine for rectangles
-% Usage: x y w h s a XX PatternFill
-%	x,y = lower left corner of box to be filled
-%	w,h = width and height of box
-%	  a = angle in degrees between lines and x-axis
-%	 XX = 0/1 for no/yes cross-hatch
-%
-/PatternFill {gsave /PFa [ 9 2 roll ] def
-  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
-  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
-  gsave 1 setgray fill grestore clip
-  currentlinewidth 0.5 mul setlinewidth
-  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
-  0 0 M PFa 5 get rotate PFs -2 div dup translate
-  0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 M 0 PFs V} for
-  0 PFa 6 get ne {
-	0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
- } if
-  stroke grestore} def
-%
-/languagelevel where
- {pop languagelevel} {1} ifelse
- 2 lt
-	{/InterpretLevel1 true def}
-	{/InterpretLevel1 Level1 def}
- ifelse
-%
-% PostScript level 2 pattern fill definitions
-%
-/Level2PatternFill {
-/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
-	bind def
-/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
->> matrix makepattern
-/Pat1 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
-	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
->> matrix makepattern
-/Pat2 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
-	8 8 L 8 0 L 0 0 L fill}
->> matrix makepattern
-/Pat3 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
-	0 12 M 12 0 L stroke}
->> matrix makepattern
-/Pat4 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
-	0 -4 M 12 8 L stroke}
->> matrix makepattern
-/Pat5 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
-	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
->> matrix makepattern
-/Pat6 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
-	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
->> matrix makepattern
-/Pat7 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
-	12 0 M -4 8 L 12 4 M 0 10 L stroke}
->> matrix makepattern
-/Pat8 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
-	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
->> matrix makepattern
-/Pat9 exch def
-/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
-/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
-/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
-/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
-/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
-/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
-/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
-} def
-%
-%
-%End of PostScript Level 2 code
-%
-/PatternBgnd {
-  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
-} def
-%
-% Substitute for Level 2 pattern fill codes with
-% grayscale if Level 2 support is not selected.
-%
-/Level1PatternFill {
-/Pattern1 {0.250 Density} bind def
-/Pattern2 {0.500 Density} bind def
-/Pattern3 {0.750 Density} bind def
-/Pattern4 {0.125 Density} bind def
-/Pattern5 {0.375 Density} bind def
-/Pattern6 {0.625 Density} bind def
-/Pattern7 {0.875 Density} bind def
-} def
-%
-% Now test for support of Level 2 code
-%
-Level1 {Level1PatternFill} {Level2PatternFill} ifelse
-%
-/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
-dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
-currentdict end definefont pop
-end
-%%EndProlog
-gnudict begin
-gsave
-doclip
-50 50 translate
-0.050 0.050 scale
-0 setgray
-newpath
-(Helvetica) findfont 140 scalefont setfont
-1.000 UL
-LTb
-1022 448 M
-63 0 V
-5904 0 R
--63 0 V
-938 448 M
-( 0.0001) Rshow
-1.000 UL
-LTb
-1022 892 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1151 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1336 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1479 M
-31 0 V
-5936 0 R
--31 0 V
-1022 1595 M
-31 0 V
-5936 0 R
--31 0 V
--5936 99 R
-31 0 V
-5936 0 R
--31 0 V
--5936 85 R
-31 0 V
-5936 0 R
--31 0 V
--5936 76 R
-31 0 V
-5936 0 R
--31 0 V
--5936 67 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.001) Rshow
-1.000 UL
-LTb
-1022 2366 M
-31 0 V
-5936 0 R
--31 0 V
-1022 2626 M
-31 0 V
-5936 0 R
--31 0 V
-1022 2810 M
-31 0 V
-5936 0 R
--31 0 V
-1022 2953 M
-31 0 V
-5936 0 R
--31 0 V
-1022 3070 M
-31 0 V
-5936 0 R
--31 0 V
--5936 98 R
-31 0 V
-5936 0 R
--31 0 V
--5936 86 R
-31 0 V
-5936 0 R
--31 0 V
--5936 75 R
-31 0 V
-5936 0 R
--31 0 V
--5936 68 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.01) Rshow
-1.000 UL
-LTb
-1022 3840 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4100 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4284 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4427 M
-31 0 V
-5936 0 R
--31 0 V
-1022 4544 M
-31 0 V
-5936 0 R
--31 0 V
--5936 99 R
-31 0 V
-5936 0 R
--31 0 V
--5936 85 R
-31 0 V
-5936 0 R
--31 0 V
--5936 76 R
-31 0 V
-5936 0 R
--31 0 V
--5936 67 R
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.1) Rshow
-1.000 UL
-LTb
-1022 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0) Cshow
-1.000 UL
-LTb
-2215 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 50) Cshow
-1.000 UL
-LTb
-3409 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 100) Cshow
-1.000 UL
-LTb
-4602 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 150) Cshow
-1.000 UL
-LTb
-5796 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 200) Cshow
-1.000 UL
-LTb
-6989 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 250) Cshow
-1.000 UL
-LTb
-1.000 UL
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-LCb setrgbcolor
-196 2659 M
-currentpoint gsave translate -270 rotate 0 0 M
-(Maximum StrainRate Error) Cshow
-grestore
-LTb
-LCb setrgbcolor
-4005 98 M
-(Particles per Cell) Cshow
-LTb
-1.000 UP
-1.000 UL
-LTb
-% Begin plot #1
-1.000 UP
-1.000 UL
-LT0
-1738 3652 M
-716 -761 V
-3886 2347 L
-6750 1935 L
-1738 3652 Pls
-2454 2891 Pls
-3886 2347 Pls
-6750 1935 Pls
-% End plot #1
-1.000 UL
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-1.000 UP
-1.000 UL
-LTb
-stroke
-grestore
-end
-showpage
-%%Trailer
-%%DocumentFonts: Helvetica
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/divergence_v.png
Binary file doc/images/divergence_v.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/divergence_v_3D.png
Binary file doc/images/divergence_v_3D.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/sinusoid.png
Binary file doc/images/sinusoid.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/sinusoid_error.eps
--- a/doc/images/sinusoid_error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ b/doc/images/sinusoid_error.eps	Wed Aug 22 15:03:39 2012 -0700
@@ -1,7 +1,7 @@
 %!PS-Adobe-2.0 EPSF-2.0
 %%Title: sinusoid_error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Thu Nov 24 11:12:32 2011
+%%Creator: gnuplot 4.6 patchlevel 0
+%%CreationDate: Wed Aug 22 12:23:36 2012
 %%DocumentFonts: (atend)
 %%BoundingBox: 50 50 410 302
 %%EndComments
@@ -20,10 +20,12 @@ gnudict begin
 /Level1 false def
 /Rounded false def
 /ClipToBoundingBox false def
+/SuppressPDFMark false def
 /TransparentPatterns false def
 /gnulinewidth 5.000 def
 /userlinewidth gnulinewidth def
 /Gamma 1.0 def
+/BackgroundColor {-1.000 -1.000 -1.000} def
 %
 /vshift -46 def
 /dl1 {
@@ -38,22 +40,6 @@ gnudict begin
 /vpt_ 31.5 def
 /hpt hpt_ def
 /vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (sinusoid_error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Thu Nov 24 11:12:32 2011)
-  /DOCINFO pdfmark
-end
-} ifelse
 /doclip {
   ClipToBoundingBox {
     newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
@@ -61,7 +47,7 @@ end
   } if
 } def
 %
-% Gnuplot Prolog Version 4.4 (January 2010)
+% Gnuplot Prolog Version 4.4 (August 2010)
 %
 %/SuppressPDFMark true def
 %
@@ -73,6 +59,7 @@ end
 /Z {closepath} bind def
 /C {setrgbcolor} bind def
 /f {rlineto fill} bind def
+/g {setgray} bind def
 /Gshow {show} def   % May be redefined later in the file to support UTF-8
 /vpt2 vpt 2 mul def
 /hpt2 hpt 2 mul def
@@ -94,6 +81,7 @@ end
 	dup 1 lt {pop 1} if 10 mul /udl exch def} def
 /PL {stroke userlinewidth setlinewidth
 	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
 % Default Line colors
 /LCw {1 1 1} def
 /LCb {0 0 0} def
@@ -431,8 +419,26 @@ Level1 {Level1PatternFill} {Level2Patter
 /Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
 dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
 currentdict end definefont pop
+Level1 SuppressPDFMark or 
+{} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (sinusoid_error.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 4.6 patchlevel 0)
+  /Author (boo)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Wed Aug 22 12:23:36 2012)
+  /DOCINFO pdfmark
+end
+} ifelse
 end
 %%EndProlog
+%%Page: 1 1
 gnudict begin
 gsave
 doclip
@@ -441,113 +447,90 @@ 0 setgray
 0 setgray
 newpath
 (Helvetica) findfont 140 scalefont setfont
+BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 7200.00 5040.00 BoxColFill} if
 1.000 UL
 LTb
-1022 448 M
+854 448 M
 63 0 V
-5904 0 R
+6030 0 R
 -63 0 V
-938 448 M
+770 448 M
 ( 0) Rshow
 1.000 UL
 LTb
-1022 1080 M
+854 1001 M
 63 0 V
-5904 0 R
+6030 0 R
 -63 0 V
--5988 0 R
-( 0.0002) Rshow
-1.000 UL
-LTb
-1022 1712 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0004) Rshow
-1.000 UL
-LTb
-1022 2344 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0006) Rshow
-1.000 UL
-LTb
-1022 2975 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0008) Rshow
-1.000 UL
-LTb
-1022 3607 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
+-6114 0 R
 ( 0.001) Rshow
 1.000 UL
 LTb
-1022 4239 M
+854 1554 M
 63 0 V
-5904 0 R
+6030 0 R
 -63 0 V
--5988 0 R
-( 0.0012) Rshow
+-6114 0 R
+( 0.002) Rshow
 1.000 UL
 LTb
-1022 4871 M
+854 2107 M
 63 0 V
-5904 0 R
+6030 0 R
 -63 0 V
--5988 0 R
-( 0.0014) Rshow
+-6114 0 R
+( 0.003) Rshow
 1.000 UL
 LTb
-1022 448 M
+854 2660 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.004) Rshow
+1.000 UL
+LTb
+854 3212 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.005) Rshow
+1.000 UL
+LTb
+854 3765 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.006) Rshow
+1.000 UL
+LTb
+854 4318 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.007) Rshow
+1.000 UL
+LTb
+854 4871 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.008) Rshow
+1.000 UL
+LTb
+854 448 M
 0 63 V
 0 4360 R
 0 -63 V
-0 -4500 R
+854 308 M
 ( 0) Cshow
 1.000 UL
 LTb
-1724 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 2) Cshow
-1.000 UL
-LTb
-2426 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 4) Cshow
-1.000 UL
-LTb
-3128 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 6) Cshow
-1.000 UL
-LTb
-3830 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 8) Cshow
-1.000 UL
-LTb
-4532 448 M
+1463 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -555,46 +538,94 @@ 0 -4500 R
 ( 10) Cshow
 1.000 UL
 LTb
-5234 448 M
+2073 448 M
 0 63 V
 0 4360 R
 0 -63 V
 0 -4500 R
-( 12) Cshow
+( 20) Cshow
 1.000 UL
 LTb
-5936 448 M
+2682 448 M
 0 63 V
 0 4360 R
 0 -63 V
 0 -4500 R
-( 14) Cshow
+( 30) Cshow
 1.000 UL
 LTb
-6638 448 M
+3291 448 M
 0 63 V
 0 4360 R
 0 -63 V
 0 -4500 R
-( 16) Cshow
+( 40) Cshow
+1.000 UL
+LTb
+3901 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 50) Cshow
+1.000 UL
+LTb
+4510 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 60) Cshow
+1.000 UL
+LTb
+5119 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 70) Cshow
+1.000 UL
+LTb
+5728 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 80) Cshow
+1.000 UL
+LTb
+6338 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 90) Cshow
+1.000 UL
+LTb
+6947 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 100) Cshow
 1.000 UL
 LTb
 1.000 UL
 LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--5967 0 V
+-6093 0 V
 Z stroke
 LCb setrgbcolor
-196 2659 M
+112 2659 M
 currentpoint gsave translate -270 rotate 0 0 M
 (Error) Cshow
 grestore
 LTb
 LCb setrgbcolor
-4005 98 M
+3900 98 M
 (Time) Cshow
 LTb
 1.000 UP
@@ -604,59 +635,57 @@ 4.000 UL
 4.000 UL
 LT0
 LCb setrgbcolor
-1694 4738 M
-(32x64) Rshow
+1526 4738 M
+(64x128) Rshow
 LT0
-1778 4738 M
+1610 4738 M
 399 0 V
-1022 448 M
-2729 2536 L
-5522 4729 L
-6989 4521 L
+3110 4354 M
+6525 3463 L
+422 -73 V
 % End plot #1
 % Begin plot #2
 stroke
 LT1
 LCb setrgbcolor
-1694 4598 M
-(64x128) Rshow
+1526 4598 M
+(128x256) Rshow
 LT1
-1778 4598 M
+1610 4598 M
 399 0 V
-1022 448 M
-853 540 V
-1060 569 V
-1397 570 V
-2048 620 V
-609 109 V
+1038 561 M
+242 107 V
+354 102 V
+662 81 V
+6947 615 L
 % End plot #2
 % Begin plot #3
 stroke
 LT2
 LCb setrgbcolor
-1694 4458 M
-(128x256) Rshow
+1526 4458 M
+(256x512) Rshow
 LT2
-1778 4458 M
+1610 4458 M
 399 0 V
-1022 448 M
-427 151 V
-472 118 V
-530 148 V
-603 146 V
-699 145 V
-831 144 V
-1025 175 V
-1335 146 V
+936 473 M
+92 31 V
+105 27 V
+121 26 V
+144 27 V
+178 26 V
+231 25 V
+332 19 V
+585 16 V
 % End plot #3
 stroke
 1.000 UL
 LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--5967 0 V
+-6093 0 V
 Z stroke
 1.000 UP
 1.000 UL
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/sinusoid_scaled_error.eps
--- a/doc/images/sinusoid_scaled_error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,677 +0,0 @@
-%!PS-Adobe-2.0 EPSF-2.0
-%%Title: sinusoid_scaled_error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Thu Nov 24 11:12:32 2011
-%%DocumentFonts: (atend)
-%%BoundingBox: 50 50 410 302
-%%EndComments
-%%BeginProlog
-/gnudict 256 dict def
-gnudict begin
-%
-% The following true/false flags may be edited by hand if desired.
-% The unit line width and grayscale image gamma correction may also be changed.
-%
-/Color true def
-/Blacktext false def
-/Solid false def
-/Dashlength 1 def
-/Landscape false def
-/Level1 false def
-/Rounded false def
-/ClipToBoundingBox false def
-/TransparentPatterns false def
-/gnulinewidth 5.000 def
-/userlinewidth gnulinewidth def
-/Gamma 1.0 def
-%
-/vshift -46 def
-/dl1 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul sub dup 0 le { pop 0.01 } if } if
-} def
-/dl2 {
-  10.0 Dashlength mul mul
-  Rounded { currentlinewidth 0.75 mul add } if
-} def
-/hpt_ 31.5 def
-/vpt_ 31.5 def
-/hpt hpt_ def
-/vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (sinusoid_scaled_error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Thu Nov 24 11:12:32 2011)
-  /DOCINFO pdfmark
-end
-} ifelse
-/doclip {
-  ClipToBoundingBox {
-    newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
-    clip
-  } if
-} def
-%
-% Gnuplot Prolog Version 4.4 (January 2010)
-%
-%/SuppressPDFMark true def
-%
-/M {moveto} bind def
-/L {lineto} bind def
-/R {rmoveto} bind def
-/V {rlineto} bind def
-/N {newpath moveto} bind def
-/Z {closepath} bind def
-/C {setrgbcolor} bind def
-/f {rlineto fill} bind def
-/Gshow {show} def   % May be redefined later in the file to support UTF-8
-/vpt2 vpt 2 mul def
-/hpt2 hpt 2 mul def
-/Lshow {currentpoint stroke M 0 vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Rshow {currentpoint stroke M dup stringwidth pop neg vshift R
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/Cshow {currentpoint stroke M dup stringwidth pop -2 div vshift R 
-	Blacktext {gsave 0 setgray show grestore} {show} ifelse} def
-/UP {dup vpt_ mul /vpt exch def hpt_ mul /hpt exch def
-  /hpt2 hpt 2 mul def /vpt2 vpt 2 mul def} def
-/DL {Color {setrgbcolor Solid {pop []} if 0 setdash}
- {pop pop pop 0 setgray Solid {pop []} if 0 setdash} ifelse} def
-/BL {stroke userlinewidth 2 mul setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/AL {stroke userlinewidth 2 div setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-/UL {dup gnulinewidth mul /userlinewidth exch def
-	dup 1 lt {pop 1} if 10 mul /udl exch def} def
-/PL {stroke userlinewidth setlinewidth
-	Rounded {1 setlinejoin 1 setlinecap} if} def
-% Default Line colors
-/LCw {1 1 1} def
-/LCb {0 0 0} def
-/LCa {0 0 0} def
-/LC0 {1 0 0} def
-/LC1 {0 1 0} def
-/LC2 {0 0 1} def
-/LC3 {1 0 1} def
-/LC4 {0 1 1} def
-/LC5 {1 1 0} def
-/LC6 {0 0 0} def
-/LC7 {1 0.3 0} def
-/LC8 {0.5 0.5 0.5} def
-% Default Line Types
-/LTw {PL [] 1 setgray} def
-/LTb {BL [] LCb DL} def
-/LTa {AL [1 udl mul 2 udl mul] 0 setdash LCa setrgbcolor} def
-/LT0 {PL [] LC0 DL} def
-/LT1 {PL [4 dl1 2 dl2] LC1 DL} def
-/LT2 {PL [2 dl1 3 dl2] LC2 DL} def
-/LT3 {PL [1 dl1 1.5 dl2] LC3 DL} def
-/LT4 {PL [6 dl1 2 dl2 1 dl1 2 dl2] LC4 DL} def
-/LT5 {PL [3 dl1 3 dl2 1 dl1 3 dl2] LC5 DL} def
-/LT6 {PL [2 dl1 2 dl2 2 dl1 6 dl2] LC6 DL} def
-/LT7 {PL [1 dl1 2 dl2 6 dl1 2 dl2 1 dl1 2 dl2] LC7 DL} def
-/LT8 {PL [2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 2 dl2 2 dl1 4 dl2] LC8 DL} def
-/Pnt {stroke [] 0 setdash gsave 1 setlinecap M 0 0 V stroke grestore} def
-/Dia {stroke [] 0 setdash 2 copy vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke
-  Pnt} def
-/Pls {stroke [] 0 setdash vpt sub M 0 vpt2 V
-  currentpoint stroke M
-  hpt neg vpt neg R hpt2 0 V stroke
- } def
-/Box {stroke [] 0 setdash 2 copy exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke
-  Pnt} def
-/Crs {stroke [] 0 setdash exch hpt sub exch vpt add M
-  hpt2 vpt2 neg V currentpoint stroke M
-  hpt2 neg 0 R hpt2 vpt2 V stroke} def
-/TriU {stroke [] 0 setdash 2 copy vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke
-  Pnt} def
-/Star {2 copy Pls Crs} def
-/BoxF {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath fill} def
-/TriUF {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath fill} def
-/TriD {stroke [] 0 setdash 2 copy vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke
-  Pnt} def
-/TriDF {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath fill} def
-/DiaF {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath fill} def
-/Pent {stroke [] 0 setdash 2 copy gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore Pnt} def
-/PentF {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath fill grestore} def
-/Circle {stroke [] 0 setdash 2 copy
-  hpt 0 360 arc stroke Pnt} def
-/CircleF {stroke [] 0 setdash hpt 0 360 arc fill} def
-/C0 {BL [] 0 setdash 2 copy moveto vpt 90 450 arc} bind def
-/C1 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C2 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C3 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C4 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C5 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc
-	2 copy moveto
-	2 copy vpt 180 270 arc closepath fill
-	vpt 0 360 arc} bind def
-/C6 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C7 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 270 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C8 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C9 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 270 450 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C10 {BL [] 0 setdash 2 copy 2 copy moveto vpt 270 360 arc closepath fill
-	2 copy moveto
-	2 copy vpt 90 180 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C11 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 180 arc closepath fill
-	2 copy moveto
-	2 copy vpt 270 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C12 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C13 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 0 90 arc closepath fill
-	2 copy moveto
-	2 copy vpt 180 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/C14 {BL [] 0 setdash 2 copy moveto
-	2 copy vpt 90 360 arc closepath fill
-	vpt 0 360 arc} bind def
-/C15 {BL [] 0 setdash 2 copy vpt 0 360 arc closepath fill
-	vpt 0 360 arc closepath} bind def
-/Rec {newpath 4 2 roll moveto 1 index 0 rlineto 0 exch rlineto
-	neg 0 rlineto closepath} bind def
-/Square {dup Rec} bind def
-/Bsquare {vpt sub exch vpt sub exch vpt2 Square} bind def
-/S0 {BL [] 0 setdash 2 copy moveto 0 vpt rlineto BL Bsquare} bind def
-/S1 {BL [] 0 setdash 2 copy vpt Square fill Bsquare} bind def
-/S2 {BL [] 0 setdash 2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S3 {BL [] 0 setdash 2 copy exch vpt sub exch vpt2 vpt Rec fill Bsquare} bind def
-/S4 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S5 {BL [] 0 setdash 2 copy 2 copy vpt Square fill
-	exch vpt sub exch vpt sub vpt Square fill Bsquare} bind def
-/S6 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S7 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt vpt2 Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S8 {BL [] 0 setdash 2 copy vpt sub vpt Square fill Bsquare} bind def
-/S9 {BL [] 0 setdash 2 copy vpt sub vpt vpt2 Rec fill Bsquare} bind def
-/S10 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt Square fill
-	Bsquare} bind def
-/S11 {BL [] 0 setdash 2 copy vpt sub vpt Square fill 2 copy exch vpt sub exch vpt2 vpt Rec fill
-	Bsquare} bind def
-/S12 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill Bsquare} bind def
-/S13 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy vpt Square fill Bsquare} bind def
-/S14 {BL [] 0 setdash 2 copy exch vpt sub exch vpt sub vpt2 vpt Rec fill
-	2 copy exch vpt sub exch vpt Square fill Bsquare} bind def
-/S15 {BL [] 0 setdash 2 copy Bsquare fill Bsquare} bind def
-/D0 {gsave translate 45 rotate 0 0 S0 stroke grestore} bind def
-/D1 {gsave translate 45 rotate 0 0 S1 stroke grestore} bind def
-/D2 {gsave translate 45 rotate 0 0 S2 stroke grestore} bind def
-/D3 {gsave translate 45 rotate 0 0 S3 stroke grestore} bind def
-/D4 {gsave translate 45 rotate 0 0 S4 stroke grestore} bind def
-/D5 {gsave translate 45 rotate 0 0 S5 stroke grestore} bind def
-/D6 {gsave translate 45 rotate 0 0 S6 stroke grestore} bind def
-/D7 {gsave translate 45 rotate 0 0 S7 stroke grestore} bind def
-/D8 {gsave translate 45 rotate 0 0 S8 stroke grestore} bind def
-/D9 {gsave translate 45 rotate 0 0 S9 stroke grestore} bind def
-/D10 {gsave translate 45 rotate 0 0 S10 stroke grestore} bind def
-/D11 {gsave translate 45 rotate 0 0 S11 stroke grestore} bind def
-/D12 {gsave translate 45 rotate 0 0 S12 stroke grestore} bind def
-/D13 {gsave translate 45 rotate 0 0 S13 stroke grestore} bind def
-/D14 {gsave translate 45 rotate 0 0 S14 stroke grestore} bind def
-/D15 {gsave translate 45 rotate 0 0 S15 stroke grestore} bind def
-/DiaE {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V closepath stroke} def
-/BoxE {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V closepath stroke} def
-/TriUE {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V closepath stroke} def
-/TriDE {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V closepath stroke} def
-/PentE {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  closepath stroke grestore} def
-/CircE {stroke [] 0 setdash 
-  hpt 0 360 arc stroke} def
-/Opaque {gsave closepath 1 setgray fill grestore 0 setgray closepath} def
-/DiaW {stroke [] 0 setdash vpt add M
-  hpt neg vpt neg V hpt vpt neg V
-  hpt vpt V hpt neg vpt V Opaque stroke} def
-/BoxW {stroke [] 0 setdash exch hpt sub exch vpt add M
-  0 vpt2 neg V hpt2 0 V 0 vpt2 V
-  hpt2 neg 0 V Opaque stroke} def
-/TriUW {stroke [] 0 setdash vpt 1.12 mul add M
-  hpt neg vpt -1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt 1.62 mul V Opaque stroke} def
-/TriDW {stroke [] 0 setdash vpt 1.12 mul sub M
-  hpt neg vpt 1.62 mul V
-  hpt 2 mul 0 V
-  hpt neg vpt -1.62 mul V Opaque stroke} def
-/PentW {stroke [] 0 setdash gsave
-  translate 0 hpt M 4 {72 rotate 0 hpt L} repeat
-  Opaque stroke grestore} def
-/CircW {stroke [] 0 setdash 
-  hpt 0 360 arc Opaque stroke} def
-/BoxFill {gsave Rec 1 setgray fill grestore} def
-/Density {
-  /Fillden exch def
-  currentrgbcolor
-  /ColB exch def /ColG exch def /ColR exch def
-  /ColR ColR Fillden mul Fillden sub 1 add def
-  /ColG ColG Fillden mul Fillden sub 1 add def
-  /ColB ColB Fillden mul Fillden sub 1 add def
-  ColR ColG ColB setrgbcolor} def
-/BoxColFill {gsave Rec PolyFill} def
-/PolyFill {gsave Density fill grestore grestore} def
-/h {rlineto rlineto rlineto gsave closepath fill grestore} bind def
-%
-% PostScript Level 1 Pattern Fill routine for rectangles
-% Usage: x y w h s a XX PatternFill
-%	x,y = lower left corner of box to be filled
-%	w,h = width and height of box
-%	  a = angle in degrees between lines and x-axis
-%	 XX = 0/1 for no/yes cross-hatch
-%
-/PatternFill {gsave /PFa [ 9 2 roll ] def
-  PFa 0 get PFa 2 get 2 div add PFa 1 get PFa 3 get 2 div add translate
-  PFa 2 get -2 div PFa 3 get -2 div PFa 2 get PFa 3 get Rec
-  gsave 1 setgray fill grestore clip
-  currentlinewidth 0.5 mul setlinewidth
-  /PFs PFa 2 get dup mul PFa 3 get dup mul add sqrt def
-  0 0 M PFa 5 get rotate PFs -2 div dup translate
-  0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 M 0 PFs V} for
-  0 PFa 6 get ne {
-	0 1 PFs PFa 4 get div 1 add floor cvi
-	{PFa 4 get mul 0 2 1 roll M PFs 0 V} for
- } if
-  stroke grestore} def
-%
-/languagelevel where
- {pop languagelevel} {1} ifelse
- 2 lt
-	{/InterpretLevel1 true def}
-	{/InterpretLevel1 Level1 def}
- ifelse
-%
-% PostScript level 2 pattern fill definitions
-%
-/Level2PatternFill {
-/Tile8x8 {/PaintType 2 /PatternType 1 /TilingType 1 /BBox [0 0 8 8] /XStep 8 /YStep 8}
-	bind def
-/KeepColor {currentrgbcolor [/Pattern /DeviceRGB] setcolorspace} bind def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke} 
->> matrix makepattern
-/Pat1 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 8 8 L 0 8 M 8 0 L stroke
-	0 4 M 4 8 L 8 4 L 4 0 L 0 4 L stroke}
->> matrix makepattern
-/Pat2 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 0 M 0 8 L
-	8 8 L 8 0 L 0 0 L fill}
->> matrix makepattern
-/Pat3 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 8 M 8 -4 L
-	0 12 M 12 0 L stroke}
->> matrix makepattern
-/Pat4 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -4 0 M 8 12 L
-	0 -4 M 12 8 L stroke}
->> matrix makepattern
-/Pat5 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 8 M 4 -4 L
-	0 12 M 8 -4 L 4 12 M 10 0 L stroke}
->> matrix makepattern
-/Pat6 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop -2 0 M 4 12 L
-	0 -4 M 8 12 L 4 -4 M 10 8 L stroke}
->> matrix makepattern
-/Pat7 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 8 -2 M -4 4 L
-	12 0 M -4 8 L 12 4 M 0 10 L stroke}
->> matrix makepattern
-/Pat8 exch def
-<< Tile8x8
- /PaintProc {0.5 setlinewidth pop 0 -2 M 12 4 L
-	-4 0 M 12 8 L -4 4 M 8 10 L stroke}
->> matrix makepattern
-/Pat9 exch def
-/Pattern1 {PatternBgnd KeepColor Pat1 setpattern} bind def
-/Pattern2 {PatternBgnd KeepColor Pat2 setpattern} bind def
-/Pattern3 {PatternBgnd KeepColor Pat3 setpattern} bind def
-/Pattern4 {PatternBgnd KeepColor Landscape {Pat5} {Pat4} ifelse setpattern} bind def
-/Pattern5 {PatternBgnd KeepColor Landscape {Pat4} {Pat5} ifelse setpattern} bind def
-/Pattern6 {PatternBgnd KeepColor Landscape {Pat9} {Pat6} ifelse setpattern} bind def
-/Pattern7 {PatternBgnd KeepColor Landscape {Pat8} {Pat7} ifelse setpattern} bind def
-} def
-%
-%
-%End of PostScript Level 2 code
-%
-/PatternBgnd {
-  TransparentPatterns {} {gsave 1 setgray fill grestore} ifelse
-} def
-%
-% Substitute for Level 2 pattern fill codes with
-% grayscale if Level 2 support is not selected.
-%
-/Level1PatternFill {
-/Pattern1 {0.250 Density} bind def
-/Pattern2 {0.500 Density} bind def
-/Pattern3 {0.750 Density} bind def
-/Pattern4 {0.125 Density} bind def
-/Pattern5 {0.375 Density} bind def
-/Pattern6 {0.625 Density} bind def
-/Pattern7 {0.875 Density} bind def
-} def
-%
-% Now test for support of Level 2 code
-%
-Level1 {Level1PatternFill} {Level2PatternFill} ifelse
-%
-/Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
-dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
-currentdict end definefont pop
-end
-%%EndProlog
-gnudict begin
-gsave
-doclip
-50 50 translate
-0.050 0.050 scale
-0 setgray
-newpath
-(Helvetica) findfont 140 scalefont setfont
-1.000 UL
-LTb
-1022 448 M
-63 0 V
-5904 0 R
--63 0 V
-938 448 M
-( 0) Rshow
-1.000 UL
-LTb
-1022 1001 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0002) Rshow
-1.000 UL
-LTb
-1022 1554 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0004) Rshow
-1.000 UL
-LTb
-1022 2107 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0006) Rshow
-1.000 UL
-LTb
-1022 2660 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0008) Rshow
-1.000 UL
-LTb
-1022 3212 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.001) Rshow
-1.000 UL
-LTb
-1022 3765 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0012) Rshow
-1.000 UL
-LTb
-1022 4318 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0014) Rshow
-1.000 UL
-LTb
-1022 4871 M
-63 0 V
-5904 0 R
--63 0 V
--5988 0 R
-( 0.0016) Rshow
-1.000 UL
-LTb
-1022 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0) Cshow
-1.000 UL
-LTb
-1724 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 2) Cshow
-1.000 UL
-LTb
-2426 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 4) Cshow
-1.000 UL
-LTb
-3128 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 6) Cshow
-1.000 UL
-LTb
-3830 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 8) Cshow
-1.000 UL
-LTb
-4532 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 10) Cshow
-1.000 UL
-LTb
-5234 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 12) Cshow
-1.000 UL
-LTb
-5936 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 14) Cshow
-1.000 UL
-LTb
-6638 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 16) Cshow
-1.000 UL
-LTb
-1.000 UL
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-LCb setrgbcolor
-196 2659 M
-currentpoint gsave translate -270 rotate 0 0 M
-(Scaled Error) Cshow
-grestore
-LTb
-LCb setrgbcolor
-4005 98 M
-(Time) Cshow
-LTb
-1.000 UP
-1.000 UL
-LTb
-% Begin plot #1
-4.000 UL
-LT0
-LCb setrgbcolor
-1694 4738 M
-(32x64) Rshow
-LT0
-1778 4738 M
-399 0 V
-1022 448 M
-2729 2275 L
-5522 4194 L
-6989 4012 L
-% End plot #1
-% Begin plot #2
-stroke
-LT1
-LCb setrgbcolor
-1694 4598 M
-(64x128) Rshow
-LT1
-1778 4598 M
-399 0 V
-1022 448 M
-853 945 V
-1060 995 V
-1397 999 V
-6380 4471 L
-609 192 V
-% End plot #2
-% Begin plot #3
-stroke
-LT2
-LCb setrgbcolor
-1694 4458 M
-(128x256) Rshow
-LT2
-1778 4458 M
-399 0 V
-1022 448 M
-427 530 V
-472 413 V
-530 517 V
-603 512 V
-699 506 V
-831 503 V
-1025 614 V
-1335 511 V
-% End plot #3
-stroke
-1.000 UL
-LTb
-1022 4871 N
-0 -4423 V
-5967 0 V
-0 4423 V
--5967 0 V
-Z stroke
-1.000 UP
-1.000 UL
-LTb
-stroke
-grestore
-end
-showpage
-%%Trailer
-%%DocumentFonts: Helvetica
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_diffusion_end.png
Binary file doc/images/thermal_diffusion_end.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_diffusion_error.eps
--- a/doc/images/thermal_diffusion_error.eps	Wed Aug 22 15:03:10 2012 -0700
+++ b/doc/images/thermal_diffusion_error.eps	Wed Aug 22 15:03:39 2012 -0700
@@ -1,7 +1,7 @@
 %!PS-Adobe-2.0 EPSF-2.0
 %%Title: thermal_diffusion_error.eps
-%%Creator: gnuplot 4.4 patchlevel 0
-%%CreationDate: Fri Nov 25 11:22:24 2011
+%%Creator: gnuplot 4.6 patchlevel 0
+%%CreationDate: Wed Aug 22 13:27:21 2012
 %%DocumentFonts: (atend)
 %%BoundingBox: 50 50 410 302
 %%EndComments
@@ -20,10 +20,12 @@ gnudict begin
 /Level1 false def
 /Rounded false def
 /ClipToBoundingBox false def
+/SuppressPDFMark false def
 /TransparentPatterns false def
 /gnulinewidth 5.000 def
 /userlinewidth gnulinewidth def
 /Gamma 1.0 def
+/BackgroundColor {-1.000 -1.000 -1.000} def
 %
 /vshift -46 def
 /dl1 {
@@ -38,22 +40,6 @@ gnudict begin
 /vpt_ 31.5 def
 /hpt hpt_ def
 /vpt vpt_ def
-Level1 {} {
-/SDict 10 dict def
-systemdict /pdfmark known not {
-  userdict /pdfmark systemdict /cleartomark get put
-} if
-SDict begin [
-  /Title (thermal_diffusion_error.eps)
-  /Subject (gnuplot plot)
-  /Creator (gnuplot 4.4 patchlevel 0)
-  /Author (boo)
-%  /Producer (gnuplot)
-%  /Keywords ()
-  /CreationDate (Fri Nov 25 11:22:24 2011)
-  /DOCINFO pdfmark
-end
-} ifelse
 /doclip {
   ClipToBoundingBox {
     newpath 50 50 moveto 410 50 lineto 410 302 lineto 50 302 lineto closepath
@@ -61,7 +47,7 @@ end
   } if
 } def
 %
-% Gnuplot Prolog Version 4.4 (January 2010)
+% Gnuplot Prolog Version 4.4 (August 2010)
 %
 %/SuppressPDFMark true def
 %
@@ -73,6 +59,7 @@ end
 /Z {closepath} bind def
 /C {setrgbcolor} bind def
 /f {rlineto fill} bind def
+/g {setgray} bind def
 /Gshow {show} def   % May be redefined later in the file to support UTF-8
 /vpt2 vpt 2 mul def
 /hpt2 hpt 2 mul def
@@ -94,6 +81,7 @@ end
 	dup 1 lt {pop 1} if 10 mul /udl exch def} def
 /PL {stroke userlinewidth setlinewidth
 	Rounded {1 setlinejoin 1 setlinecap} if} def
+3.8 setmiterlimit
 % Default Line colors
 /LCw {1 1 1} def
 /LCb {0 0 0} def
@@ -431,8 +419,26 @@ Level1 {Level1PatternFill} {Level2Patter
 /Symbol-Oblique /Symbol findfont [1 0 .167 1 0 0] makefont
 dup length dict begin {1 index /FID eq {pop pop} {def} ifelse} forall
 currentdict end definefont pop
+Level1 SuppressPDFMark or 
+{} {
+/SDict 10 dict def
+systemdict /pdfmark known not {
+  userdict /pdfmark systemdict /cleartomark get put
+} if
+SDict begin [
+  /Title (thermal_diffusion_error.eps)
+  /Subject (gnuplot plot)
+  /Creator (gnuplot 4.6 patchlevel 0)
+  /Author (boo)
+%  /Producer (gnuplot)
+%  /Keywords ()
+  /CreationDate (Wed Aug 22 13:27:21 2012)
+  /DOCINFO pdfmark
+end
+} ifelse
 end
 %%EndProlog
+%%Page: 1 1
 gnudict begin
 gsave
 doclip
@@ -441,65 +447,114 @@ 0 setgray
 0 setgray
 newpath
 (Helvetica) findfont 140 scalefont setfont
+BackgroundColor 0 lt 3 1 roll 0 lt exch 0 lt or or not {BackgroundColor C 1.000 0 0 7200.00 5040.00 BoxColFill} if
 1.000 UL
 LTb
-938 448 M
+854 448 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
-854 448 M
-( 0.005) Rshow
+770 448 M
+( 0.002) Rshow
 1.000 UL
 LTb
-938 1333 M
+854 1080 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
+-6114 0 R
+( 0.004) Rshow
+1.000 UL
+LTb
+854 1712 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.006) Rshow
+1.000 UL
+LTb
+854 2344 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
+( 0.008) Rshow
+1.000 UL
+LTb
+854 2975 M
+63 0 V
+6030 0 R
+-63 0 V
+-6114 0 R
 ( 0.01) Rshow
 1.000 UL
 LTb
-938 2217 M
+854 3607 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.015) Rshow
+-6114 0 R
+( 0.012) Rshow
 1.000 UL
 LTb
-938 3102 M
+854 4239 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.02) Rshow
+-6114 0 R
+( 0.014) Rshow
 1.000 UL
 LTb
-938 3986 M
+854 4871 M
 63 0 V
-5988 0 R
+6030 0 R
 -63 0 V
--6072 0 R
-( 0.025) Rshow
+-6114 0 R
+( 0.016) Rshow
 1.000 UL
 LTb
-938 4871 M
-63 0 V
-5988 0 R
--63 0 V
--6072 0 R
-( 0.03) Rshow
-1.000 UL
-LTb
-938 448 M
+854 448 M
 0 63 V
 0 4360 R
 0 -63 V
-938 308 M
+854 308 M
+( 0.03) Cshow
+1.000 UL
+LTb
+1463 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 0.04) Cshow
+1.000 UL
+LTb
+2073 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 0.05) Cshow
+1.000 UL
+LTb
+2682 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
 ( 0.06) Cshow
 1.000 UL
 LTb
-1543 448 M
+3291 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 0.07) Cshow
+1.000 UL
+LTb
+3901 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -507,7 +562,15 @@ 0 -4500 R
 ( 0.08) Cshow
 1.000 UL
 LTb
-2148 448 M
+4510 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 0.09) Cshow
+1.000 UL
+LTb
+5119 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -515,7 +578,15 @@ 0 -4500 R
 ( 0.1) Cshow
 1.000 UL
 LTb
-2753 448 M
+5728 448 M
+0 63 V
+0 4360 R
+0 -63 V
+0 -4500 R
+( 0.11) Cshow
+1.000 UL
+LTb
+6338 448 M
 0 63 V
 0 4360 R
 0 -63 V
@@ -523,78 +594,30 @@ 0 -4500 R
 ( 0.12) Cshow
 1.000 UL
 LTb
-3358 448 M
+6947 448 M
 0 63 V
 0 4360 R
 0 -63 V
 0 -4500 R
-( 0.14) Cshow
-1.000 UL
-LTb
-3963 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.16) Cshow
-1.000 UL
-LTb
-4569 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.18) Cshow
-1.000 UL
-LTb
-5174 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.2) Cshow
-1.000 UL
-LTb
-5779 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.22) Cshow
-1.000 UL
-LTb
-6384 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.24) Cshow
-1.000 UL
-LTb
-6989 448 M
-0 63 V
-0 4360 R
-0 -63 V
-0 -4500 R
-( 0.26) Cshow
+( 0.13) Cshow
 1.000 UL
 LTb
 1.000 UL
 LTb
-938 4871 N
-938 448 L
-6051 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--6051 0 V
+-6093 0 V
 Z stroke
 LCb setrgbcolor
-196 2659 M
+112 2659 M
 currentpoint gsave translate -270 rotate 0 0 M
 (Error) Cshow
 grestore
 LTb
 LCb setrgbcolor
-3963 98 M
+3900 98 M
 (h) Cshow
 LTb
 1.000 UP
@@ -605,20 +628,20 @@ 4.000 UL
 4.000 UL
 LT0
 /Helvetica findfont 140 scalefont setfont
-6686 4372 M
-2905 1622 L
-1014 680 L
-6686 4372 Pls
-2905 1622 Pls
-1014 680 Pls
+6642 4659 M
+2834 2425 L
+930 989 L
+6642 4659 Pls
+2834 2425 Pls
+930 989 Pls
 % End plot #1
 1.000 UL
 LTb
-938 4871 N
-938 448 L
-6051 0 V
+854 4871 N
+854 448 L
+6093 0 V
 0 4423 V
--6051 0 V
+-6093 0 V
 Z stroke
 1.000 UP
 1.000 UL
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_diffusion_start.png
Binary file doc/images/thermal_diffusion_start.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_eulear_high.png
Binary file doc/images/thermal_eulear_high.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_eulear_low.png
Binary file doc/images/thermal_eulear_low.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_eulear_medium.png
Binary file doc/images/thermal_eulear_medium.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_eulear_start.png
Binary file doc/images/thermal_eulear_start.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_euler_high.png
Binary file doc/images/thermal_euler_high.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_euler_low.png
Binary file doc/images/thermal_euler_low.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_euler_medium.png
Binary file doc/images/thermal_euler_medium.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_euler_start.png
Binary file doc/images/thermal_euler_start.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_lagrangian_end.png
Binary file doc/images/thermal_lagrangian_end.png has changed
diff -r 38aa5c76ef3f -r 10fa4d5529a8 doc/images/thermal_lagrangian_start.png
Binary file doc/images/thermal_lagrangian_start.png has changed