applications/linearFEM1d.doxy

Go to the documentation of this file.

string  author  "(empty)"
string  comment "(empty)"
string  title   "(empty)"
int     gauss_p 1
int     level   3
end
// output starts here
double  exEnergie       0.0246385
array double energieFehler {
        1       0.100759
        2       0.0318543
}
array double h {
        1       0.5
        2       0.25
}
array double maxNodalFehler {
        1       0.0644745
        2       0.0464661
} @endcode

    Error diagrams:
    <img src="linearFEM1d-conv.png" width="640" height="480" alt="error diagram">
    The diagram above was created using the following commands: @code
set logscale x
set logscale y
set grid
set grid mxtics
set grid mytics
set key top left
set xlabel 'h'
set ylabel 'error'
plot 'linearFEM1d-1.gnuplot' using 3:2 title 'relative energy error, gauss_p = 1', \
  '' using 3:4 title 'max nodal error, gauss_p = 1', \
  'linearFEM1d-2.gnuplot' using 3:2 title 'relative energy error, gauss_p = 2', \
  '' using 3:4 title 'max nodal error, gauss_p = 2' @endcode

    Exact and numerical solution:
    <img src="linearFEM1d-solution.png" width="640" height="480"
  alt="exact and numerical solution">
    The diagram above was created using the following commands: @code
set xlabel 'x'
set ylabel 'u(x)'
set key bottom left
plot 'solution1d-1.gnuplot', 'solution1d-2.gnuplot', 'solution1d-3.gnuplot', -3.0*cosh(x)/cosh(1.0)+x*x+2.0 @endcode

    @section complete Complete Source Code
    @author Philipp Frauenfelder, 2004
*/