← Back to branch summary

~ubuntu-branches/ubuntu/trusty/rheolef/trusty-proposed

« back to all changes in this revision

Viewing changes to doc/usrman/crack-compare-uniform-adapt-P2.plot

  • Committer: Bazaar Package Importer
  • Author(s): Christophe Prud'homme
  • Date: 2010-06-12 09:08:59 UTC
  • Revision ID: james.westby@ubuntu.com-20100612090859-8gpm2gc7j3ab43et
Tags: upstream-5.89
Import upstream version 5.89

Show diffs side-by-side

added added

removed removed

Lines of Context:
doc/usrman/crack-compare-uniform-adapt-P2.plot
 
1
# crack P1 / P2
 
2
# comparaison
 
3
#   uniforme de pas h 
 
4
#   adapte avec hmin = coef^2 et phi=u pour P1
 
5
#   adapte avec hmin = coef^4 et phi=grad u pour P2
 
6
#
 
7
# representation des erreurs L2, Linf et H1
 
8
# en fct de n_elt
 
9
#
 
10
# format adapt.dat:
 
11
# coef niter n_node n_elt hmin hmax err_l2 err_infty err_h1
 
12
#  1     2     3     4     5    6    7       8        9
 
13
 
 
14
# format uniform.dat:
 
15
# i n_node n_elt h err_l2 err_infty err_h1
 
16
# 1   2      3   4   5       6        7
 
17
 
 
18
set logscale
 
19
set size square
 
20
 
 
21
set xlabel "n_dof"
 
22
set title "erreur H1"
 
23
plot \
 
24
        "crack-adapt-P1.dat" u 3:9 w lp, \
 
25
        "crack-uniform-P1.dat" u 2:7 w lp, \
 
26
        "crack-adapt-P2.dat" u (($3)+($4)):9 w lp, \
 
27
        "crack-uniform-P2.dat" u (($2)+($3)):7 w lp, \
 
28
        0.1*x**(-0.25)
 
29
pause -1 "<retour>"
 
30
 
 
31
set title "erreur L2"
 
32
plot \
 
33
        "crack-adapt-P1.dat" u 3:7 w lp, \
 
34
        "crack-uniform-P1.dat" u 2:5 w lp, \
 
35
        "crack-adapt-P2.dat" u (($3)+($4)):7 w lp, \
 
36
        "crack-uniform-P2.dat" u (($2)+($3)):5 w lp, \
 
37
        0.1*x**(-3./4)
 
38
pause -1 "<retour>"
 
39
 
 
40
set title "erreur Linf"
 
41
plot \
 
42
        "crack-adapt-P1.dat" u 3:8 w lp, \
 
43
        "crack-uniform-P1.dat" u 2:6 w lp, \
 
44
        "crack-adapt-P2.dat" u (($3)+($4)):8 w lp, \
 
45
        "crack-uniform-P2.dat" u (($2)+($3)):6 w lp, \
 
46
        0.1*x**(-0.25)
 
47
pause -1 "<retour>"
 
48
 
 
49